1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
#include "pdfs.hpp"
float Binomial(float p, int x, int n) {
return pow(p,n) * pow(1-p,x-n) * Choose(x,n);
}
float Poisson(float lambda, int n) {
if (n < FactorialTableLength) {
float lambdaToN = pow(lambda, n);
if (lambdaToN == HUGE_VALF)
return 1;
return (pow(lambda, n) * exp(-lambda)) / FactorialTable[n];
}
else {
return 0;
}
}
float Gamma(float lambda, float t, int n) {
if (n <= FactorialTableLength) {
return (pow(lambda, n) * pow(t,n-1) * exp(-lambda*t)) / FactorialTable[n-1];
}
else {
return 0;
}
}
float Exponential(float lambda, int t) {
return lambda * exp(-lambda*t);
}
|
