【求助】关于underestimate from NLMIXED Model 的问题

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大家好，我有一个关于NLMIXED Model 的问题，想请教一下，希望可以得到宝贵意见。 任何建议都将十分感激。这里我将用"Fixedeffects regression methods fro longitudinal data" (by Paul D. Allison) 里的例子来说明我遇见的问题。

PROC
NLMIXED DATA=teenyrs5;
eta=b0 + byr1*(year=1)+byr2*(year=2) + byr3*(year=3) +byr4*(year=4) + bmother*mother+bspouse*spouse +bschool*inschool + bhours*hours + alpha;
p=1/(1+ EXP(-eta));
MODEL pov~BINARY(p);
RANDOM alpha~ NORMAL(0,s2) SUBJECT=id out=abc;
PARMS b0=-.29 byr1=-.06 byr2=.16 byr3=.09 byr4=.09 bmother=.99 bspouse=-1.26 bschool=.24 bhours=-.03 s2=1 ;
predict P out=p;
predict 1/(1+ EXP(-(b0 + byr1*(year=1)+ byr2*(year=2) + byr3*(year=3) +
byr4*(year= 4)+ bmother*mother + bspouse*spouse +
bschool*inschool + bhours*hours))) out=p2;
RUN;

当我用上面的程序执行了PROC NLMIXED 以后，估算的FIXED 和 RANDOM EFFECT 的PHAT 都存在了数据 P2和P里面。 然后我计算出了每一年的 POV% 和实际值比较。发现无论是FIXED还是RANDOM的都比实际值低，而且是平行移动。

year	actual_rt	RT_random	RT_Fixed
1	34.84%	33.97%	31.21%
2	39.88%	39.49%	37.45%
3	37.97%	37.43%	35.34%
4	38.92%	38.51%	36.54%
5	36.84%	36.27%	34.25%

这是我第一次用NLMIXED MODEL，不知道是不是中间有什么步骤用的不对，或者MIXED MODEL不可以这样估值。同时，还想请教一下可否把ALPHA的MEAN也设成RANDOM的变量，而不是常数0.


十分希望并且感激大家可以提供宝贵的意见。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：mixed model estimate nlmixed model Mixed methods