我正在使用logit实现一个模型。我如何查找数据的方差来衡量,如果预测是好的,那么测量的数据的方差是否良好输入代码?
m <- glm(a~ f+ EN*RT, family = binomial(link = "logit"),data=df)`
我读过我可以使用置信区间做这样的事情:
p <- predict(m, newdata=df, type="link",se.fit = T)
critval <- 1.96 ## approx 95% CI
upr <- p$fit + (critval * p$se.fit)
lwr <- p$fit - (critval * p$fit)
fit <- p$fit
mean(((upr-fit)/1.96)^2, na.rm=T)
但我得到的价值非常巨大。我将非常感谢你的帮助。