On Apr 13, 2013, at 06:47 , Thomas Lumley wrote:
> On Sat, Apr 13, 2013 at 6:31 AM, peter dalgaard <pda...@gmail.com> wrote:
> But is it supposed to be t^{-3/2} or t^{-0.5}?? The formula has the former
> and the code the latter, and the integral is clearly divergent with the
> former.
>
>
> I didn't look at that -- I just assumed that the code would reproduce the
> reported error.
>
So did I, but I looked at the formula instead, and got confused when you
obtained convergence.
For the benefit of other recipients:
Of course it makes quite a difference since t^-alpha is integrable near zero
iff alpha < 1, and it is the dominating term when x==0.
Numerically, even theoretically convergent integrals can confuse the algorithm,
basically because it tries to infer global behavior from a few evaluations of
the integrand.
-pd
--
Peter Dalgaard, Professor,
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: pd....@cbs.dk Priv: pda...@gmail.com
______________________________________________
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.
