我正在寻找一种方法,在使用Stata中的blogit命令估算分组数据的logit模型后,绘制二元结果的预测概率。 该图需要显示y轴上的预测概率和x轴上的连续变量(天)。
我已尝试过UC LA的xi3和postgr3命令,但它们似乎不适用于blogit命令。
非常感谢帮助!如果您需要更多信息,请与我们联系。
谢谢!
答案 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