我在R中创建的模型是:
fit <- lm(hired ~ educ + exper + sex, data=data)
我不确定的是如何适应模型来预测感兴趣的概率,其中p = pr(雇用= 1).
任何帮助将不胜感激,谢谢
编辑:glm在我的模型中扮演什么角色呢？(我的答案如下)基于Jason对Greg的回答所做的编辑我不明白它具体做了什么.
我的答案是否分析了被录用的几率？
对于估计的模型glm,您可以使用该predict函数为数据集中的每个观测值提取线性预测变量.然后,您可以简单地使用适当的概率分布函数来获得预测概率.例如,在逻辑回归的情况下,使用plogis.换句话说,如果mod您的模型适合glm:
> plogis(predict(mod))
假设您估算了一个逻辑模型,它将返回数据集中每个观测值的预测概率.如果您需要计算不在数据集中的点的预测概率,请参阅newdata选项predict.请注意,predict还可以在每个点提供标准错误.predict.glm有关更多信息,请查看文档.
编辑:根据Greg的建议,您可以type="response"在电话中使用免费predict获得plogis:
> predict(mod, type="response")