Sam Wong <qxsr <at> yahoo.com> writes:
>
>
> Hi, R Users;
>
> I am trying to maximize a likelihood function which
> contains an integral. The integral contains the
> unknown parameter as well. I am trying to use the
> following code to do the maximization:
>
[snip]
> The error message is:
> Error in log(x) : Non-numeric argument to mathematical
> function
>
The only obvious (?) mistake in the code is that
you refer to integrate(...)$integral -- the help page
for ?integrate says that $value is the estimate of
the integral, so you might be getting a NULL here.
However, the code you gave us is not reproducible/self-contained -- since I
don't
know what z1, z2, n1, rho are, I can't run it and see if that's the only
problem.
cheers
Ben Bolker
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.