
R code here:
lodpval_uncons = function(lod) { loglike = 2 * log(10^lod) pchisq(loglike, df=2, lower.tail=F) } lodpval_pt = function(lod, dom_ignore) { loglike = 2 * log(10^lod) p2 = pchisq(loglike, df=2, lower.tail=F) p1 = pchisq(loglike, df=1, lower.tail=F) if(! dom_ignore) { 0.1 * p2 + 0.5 * p1 } else { 0.5 * p1 } }
0 comments:
Post a Comment