[R] How to use mle or similar with integrate?

From: Rainer M Krug <rkrug_at_sun.ac.za>
Date: Fri 23 Jun 2006 - 23:29:07 EST


I have the following formula (I hope it is clear - if no, I can try to do better the next time)

h(x, a, b) =
integral(0 to pi/2)
    integral(D/sin(alpha) to Inf)

        f(x, a, b)


and I want to do an mle with it.
I know how to use mle() and I also know about integrate(). My problem is to give the parameter values a and b to the integrate function.

In other words, how can I write

h <- function...

so that I can estimate a and b?



Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation
Biology (UCT)

Department of Conservation Ecology and Entomology
University of Stellenbosch
Matieland 7602
South Africa

R-help@stat.math.ethz.ch mailing list
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Fri Jun 23 23:31:49 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sun 25 Jun 2006 - 06:11:55 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.