The R Project SVN R

Rev

Rev 68947 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 68947 Rev 77685
Line 35... Line 35...
35
 
35
 
36
double rf(double n1, double n2)
36
double rf(double n1, double n2)
37
{
37
{
38
    double v1, v2;
38
    double v1, v2;
39
    if (ISNAN(n1) || ISNAN(n2) || n1 <= 0. || n2 <= 0.)
39
    if (ISNAN(n1) || ISNAN(n2) || n1 <= 0. || n2 <= 0.)
40
	ML_ERR_return_NAN;
40
	ML_WARN_return_NAN;
41
 
41
 
42
    v1 = R_FINITE(n1) ? (rchisq(n1) / n1) : 1;
42
    v1 = R_FINITE(n1) ? (rchisq(n1) / n1) : 1;
43
    v2 = R_FINITE(n2) ? (rchisq(n2) / n2) : 1;
43
    v2 = R_FINITE(n2) ? (rchisq(n2) / n2) : 1;
44
    return v1 / v2;
44
    return v1 / v2;
45
}
45
}