Re: [R] Random numbers from noncentral t-distribution

From: Thomas Lumley <tlumley_at_u.washington.edu>
Date: Sat 01 Jul 2006 - 03:24:32 EST

On Sat, 1 Jul 2006, Long Qu wrote:

> Thank you very much for your kind reply. It solved the problem of rt( ). :D
>
> But it seems that the qt( ) also have problems:

Yes, there does seem to be a problem near zero. A clearer version is

   curve(qt(x,df=20,ncp=3),from=0,to=0.004)    curve(qt(10^x,df=20,ncp=3),from=-10,to=-2,n=1000)

The fix is less obvious here. I'll file it as a bug.

         -thomas

Thomas Lumley			Assoc. Professor, Biostatistics
tlumley@u.washington.edu	University of Washington, Seattle

______________________________________________
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 Sat Jul 01 03:37:28 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 Sat 01 Jul 2006 - 04:15:58 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.