在blogit命令之后绘制预测概率

时间:2015-07-28 12:32:02

标签: graph stata logistic-regression

我正在寻找一种方法,在使用Stata中的blogit命令估算分组数据的logit模型后,绘制二元结果的预测概率。 该图需要显示y轴上的预测概率和x轴上的连续变量(天)。

我已尝试过UC LA的xi3和postgr3命令,但它们似乎不适用于blogit命令。

非常感谢帮助!如果您需要更多信息,请与我们联系。

谢谢!

2 个答案:

答案 0 :(得分:1)

如果我理解正确,一个简单的解决方案是在pop_var为1的情况下创建一个反事实数据集,然后只需使用predict命令获取概率,如下所示:

webuse xmpl2
blogit deaths pop agecat exposed
replace pop = 1
predict probability

然后,您可以使用各种twoway命令对此进行图形化,无论您想要什么。虽然我的例子没有连续的x轴,但它可能是这样的:

twoway scatter probability agecat, by(exposed)

答案 1 :(得分:1)

你也可以在glm框架中使用这种类型的模型:

webuse xmpl2, clear
list, clean noobs
blogit deaths pop agecat exposed
glm deaths agecat exposed, link(logit) family(binomial pop) nolog
replace pop = 1
margins, at(agecat = (0 1)) 
marginsplot