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

From: Thomas Lumley <>
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 Lumley			Assoc. Professor, Biostatistics	University of Washington, Seattle

______________________________________________ mailing list PLEASE do read the posting guide! 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 Please read the posting guide before posting to the list.