[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


Hi

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)
      )
      dx

    )
  dalpha
)

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?

Thanks,

Rainer

-- 
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
https://stat.ethz.ch/mailman/listinfo/r-help
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.