Hi, You may want to try the double exponential transformation on the numerator and the denominator on this one. The method is described in detail here: http://projecteuclid.org/DPubS?service=UI&version=1.0&verb=Display&handle=euclid.prims/1145474600 If you want to give it a shot outside R there are a couple of C/C++ implementations out there: http://www.codeproject.com/KB/recipes/FastNumericalIntegration.aspx http://www.crbond.com/download/deint.c I have used the second one in a project of mine that is similar to the one you describe (calculation of posterior expectations in Bayesian inference), and I can say that it works reasonably well. As a side note, the DE method is one (actually the default) of the numerical integration methods implemented by Mathematica's NIntegrate function. It would be nice if someone found the time to augment R's integration facilities with this (and possibly other) numerical integration methods. Christos Argyropoulos University of Pittsburgh Medical Center _________________________________________________________________
s. It's easy! aspx&mkt=en-us ______________________________________________ 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.
