Displaying 1 result from an estimated 1 matches for "prob_a".
Did you mean:
prob_1
2008 Mar 02
2
difference between lrm's "Model L.R." and anova's "Chi-Square"
...e" in the anova results, but not quite; the latter value
is always slightly smaller.
anova() calculates the p-value based on "Chi-Square", but I have
independent evidence that "Model L.R." is the actual -2*log(LR), so should
I be using that?
Why are the values different?
prob_a <- inv.logit(rnorm(1,0,1))
prob_b <- inv.logit(rnorm(1,0,1))
data <- data.frame(
factor=c(rep("a",500),rep("b",500)),
outcome=c(sample(c(1,0),100,replace=T,prob=c(prob_a,1-prob_a)),
sample(c(1,0),100,replace=T,prob=c(prob_b,1-prob_b))))
fit <- lrm(outcome~f...