小弟在做完回归模型后,想要看看符合特定条件的个案所对应的概率。但是用margin 和 predict 求得的概率不一致;还请知道的坛友帮忙解答一下。
为了便于解释,我用stata自带的数据做演示。
webuse margex
reg y i.sex i.group
margin sex
得到的结果是:
------------------------------------------------------------------------------
| Delta-method
| Margin Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
sex |
0 | 60.56034 .5781782 104.74 0.000 59.42713 61.69355
1 | 78.88236 .5772578 136.65 0.000 77.75096 80.01377
------------------------------------------------------------------------------
可以看到,sex==1 时,概率为78.88
如果我用predict 预测 sex==1的概率的话,
predict p1 if sex==1
tab p1 // p1有三个取值,因为group有三个取值。
Fitted |
values | Freq. Percent Cum.
------------+-----------------------------------
71.64348 | 984 65.51 65.51
79.6811 | 452 30.09 95.61
90.28271 | 66 4.39 100.00
------------+-----------------------------------
Total | 1,502 100.00
然后sum一下
sum p1
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
y1 | 1502 74.88129 4.923423 71.64348 90.28271
——————————————————————————————————————————————
这样,用margin得出的sex==1的概率值就是78.88;
而用predict得出的sex==1时的概率值就是74.88。
那么请问,如果我想知道性别为1的概率,到底该用哪一个呢?
谢谢大家!


雷达卡



京公网安备 11010802022788号







