Displaying 1 result from an estimated 1 matches for "pln2".
Did you mean:
pin2
2004 Mar 02
2
Problem with Integrate
..., mu=2, sigma2=2.8)
Error in integrate(PLN1, -1000, 1000, Count = 73, mu = 2, sigma2 = 2.8):
non-finite function value
In addition: Warning message:
NaNs produced in: dpois(x, lambda, log)
So, the integral gets smaller, and then gives an error.
I then tried entering the formula directly:
PLN2 <- function(LL, Count, mu, sigma2) {
exp(-LL-(log(LL)-mu)^2/(2*sigma2))*LL^(Count-1)/(gamma(Count+1)*sqrt(2*pi*sigma2))
}
> integrate(PLN2, 0, 100, Count=73, mu=2, sigma2=2.8)
0.001287821 with absolute error < 2.6e-10
> integrate(PLN2, 0, 1000, Count=73, mu=2, sigma2=2.8)
0.001289726...