从Probit模型 - R命令计算概率？

Question

我正在使用Probit模型，并希望从我的模型计算DF中每个观察的概率。我知道我可以使用公式来计算这个，但是我想知道是否有一种快速输出概率的方法并将它们附加到我的DF。

我正在运行以下模型：

attach(non.part.2)
y <- cbind(E)
x1 <- cbind(tech.ems, med.com, tech.nonemerg)
probit <- glm(y ~ x1, family = binomial (link = "probit"))
summary(probit)

我正在运行几个模型，所以能够让R吐出概率并允许我在我的DF（非.part.2）中命名它们会很好 - 比如p_x1 - 这样我就可以运行以后各种型号的摘要统计。

非常感谢任何帮助！

Answer 1

以下情况应该有效。

non.part2$p_x1 <- predict(probit, yourDataToPredictOn, type = "response")