The R Project SVN R

Rev

Rev 8431 | Rev 42302 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 8431 Rev 23878
Line 37... Line 37...
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_ERR_return_NAN;
41 
 
 41 
 
42 
    v1 = R_FINITE(n1) ? (rchisq(n1) / n1) : norm_rand();
 42 
    v1 = R_FINITE(n1) ? (rchisq(n1) / n1) : 1;
43 
    v2 = R_FINITE(n2) ? (rchisq(n2) / n2) : norm_rand();
 43 
    v2 = R_FINITE(n2) ? (rchisq(n2) / n2) : 1;
44 
    return v1 / v2;
 44 
    return v1 / v2;
45 
}
 45 
}