search for: probh

...Gp = 2*sum(obs*log(obs/exp)) dfp = max(as.integer(levels(cov)))-1 probp = pchisq(Gp, dfp, lower.tail=F) #total G Gt = sum(G) dft = sum(dfs) probt = pchisq(Gt, dft, lower.tail=F) #heterogeneity G Gh = Gt-Gp dfh = dft - dfp probh = pchisq(Gh, dfh, lower.tail=F) #F ratio Fratio = (Gp/dfp)/(Gh/dfh) probf = pf(Fratio, dfp, dfh, lower.tail=F) res = list(TotG=c(Gt, dft, probt), PooledG=c(Gp, dfp, probp), HetG=c(Gh, dfh, probh), Fratio=c(Fratio, dfp, dfh, probf)) res } #end hetg...