Displaying 1 result from an estimated 1 matches for "figa3".
Did you mean:
figa
2006 Aug 26
1
Capture of iterative integration output
...n.
Any help on y$value capture would be appreciated. I would be ecstatic
for help on implementation of the optim function using finite
differences for the following code.
By the way in case anyone is interested function(x) is the
Poisson-lognormal pdf.
q <- read.table(file="c:/Bayes/FigA3-1.prn",header=TRUE)
attach(q)
V <- c(volume)
c <- c(count)
n <- length(count)
mu <- -8.202900475
sigma <- 1.609437912
for(i in c(1:n))
{
integrand <- function(x)
{(x*V[[i]])^c[[i]]*exp(-x*V[[i]])/factorial(c[[i]])*0.3989422804014327/
(x*sigma)*exp(...