The R Project SVN R

Rev

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

Rev 8431 Rev 10470
Line 43... Line 43...
43 
	return qchisq(p, n1, lower_tail, log_p) / n1;
 43 
	return qchisq(p, n1, lower_tail, log_p) / n1;
44 
 
 44 
 
45 
    if (n1 > 4e5)
 45 
    if (n1 > 4e5)
46 
	return 1/qchisq(p, n2, !lower_tail, log_p) * n2;
 46 
	return 1/qchisq(p, n2, !lower_tail, log_p) * n2;
47 
 
 47 
 
48 
    p = (1. / qbeta(R_DT_CIv(p), n2/2, n1/2, LTRUE, LFALSE) - 1.) * (n2 / n1);
 48 
    p = (1. / qbeta(R_DT_CIv(p), n2/2, n1/2, TRUE, FALSE) - 1.) * (n2 / n1);
49 
    return ML_VALID(p) ? p : ML_NAN;
 49 
    return ML_VALID(p) ? p : ML_NAN;
50 
}
 50 
}