Displaying 1 result from an estimated 1 matches for "probp".
Did you mean:
probe
2007 May 23
0
Replicated LR goodness-of-fit tests, heterogeneity G, with loglm?
...ength(weightexp)) {
sum = sum +
weightexp[[i]]
}
expproppooled =
sum/sum(weights)
#Now the pooled G
obs =
tapply(obscounts,cov,sum)
exp =
expproppooled*sum(obs)
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...