From tlumley at u.washington.edu  Mon Jul  3 20:04:43 2006
From: tlumley at u.washington.edu (Thomas Lumley)
Date: Mon, 3 Jul 2006 11:04:43 -0700 (PDT)
Subject: [Rd] noncentral F-distributed random numbers (PR#9055)
In-Reply-To: <20060701191347.47E8950A21@slim.kubism.ku.dk>
References: <20060701191347.47E8950A21@slim.kubism.ku.dk>
Message-ID: <Pine.LNX.4.64.0607031103200.21683@homer22.u.washington.edu>

On Sat, 1 Jul 2006, longor at iastate.edu wrote:
>
> The rf() function reads:
>> rf
> function (n, df1, df2, ncp = 0)
> {
>    if (ncp == 0)
>        .Internal(rf(n, df1, df2))
>    else rchisq(n, df1, ncp = ncp)/rchisq(n, df2)
> }
> <environment: namespace:stats>
>
> where I believe both the numerator and the denominator should be divided by
> their corresponding degrees of freedom.
>

Yes. Fixed for 2.4.0.

I didn't add a denominator non-centrality parameter, since we don't have 
code for qf, pf, or df for the doubly non-central distribution.

 	-thomas



[Rd] noncentral F-distributed random numbers (PR#9055)
Thomas Lumley