Well if you define
F(x,y) = I_{x > g(y)} f(x,y),
where I is an indicator function giving a 0 if x < g(y) and a 1 of x > g(y) then
\int_a^b \int_{g(y)}^Inf f(x,y) dx dy = \int_a^b \int_{-Inf}^Inf F(x,y) dx dy
There are several things you can do to make the integrals easier to be
evaluated. The first thing I would is to make a substitution in x to map the
interval (-Inf, Inf) to e.g. (0, 1). Then you only need to integrate over (0,1)
in the inner integral. I think you can read about that in Abramowitz, Milton
and Stegun, Irene A. (Eds.). Handbook of Mathematical Functions with Formulas,
Graphs, and Mathematical Tables ( or is it Gradshteyn, I.S. and Ryzhik, I.M.
Tables of Integrals, Series, and Products??). Perhaps you or other can provide
other references.
Med venlig hilsen
Frede Aakmann Tøgersen
> -----Oprindelig meddelelse-----
> Fra: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] På vegne af Chris Rhoads
> Sendt: 23. oktober 2007 22:58
> Til: r-help@r-project.org
> Emne: [R] 2-D numerical integration over odd region
>
> Hello all,
>
> I'm hoping to find a way to evaluate the following sort of
> integral in R.
>
> \int_a^b \int_{g(y)}^Inf f(x,y) dx dy.
>
> The integral has no closed form and so must be evaluated
> numerically. The "adapt" package provides for
> multidimensional integration but does not appear to allow the
> limits of integration to be a function. I need to evaluate a
> number of integrals of this sort and I need the evaluations
> to be fairly precise (something like 10^-6 or 10^-7 would be
> sufficient) so I'd prefer to avoid time-consuming MCMC methods.
>
> Any ideas?
>
> Thanks,
>
> Chris Rhoads
>
> ______________________________________________
> R-help@r-project.org 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.
>
______________________________________________
R-help@r-project.org 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.
