全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
3397 8
2012-02-28
我用SAS计算零膨胀负二项模型(zero inflated negative binomial model):
proc nlmixed data=one;
parms j1=-3 j2=0.5 j3=-3
      k1=-2 k2=0.8 k3=-2.8 k4=-3.2 ;
mort=j1*A1+j2*Rs1+j3*H1;
pmort=1/(1+exp(-mort));
lamda=exp(k1+k2*Rs1+k3*Dm1+k4*Pw);
m=1/alpha ;
p=1/(1+alpha*lamda);
if y=0 then
     ll=log(pmort+(1-pmort)*(p**m));
else ll = log(1-pmort) + log(gamma(m + y)) - log(gamma(y + 1))
       - log(gamma(m)) + m*log(p) + y*log(1-p);
   model y ~ general(ll);out=pp;
run;

出现以下错误,能否帮忙解决下,谢谢了!

NOTE: Execution error for observation 48.
NOTE: PROCEDURE NLMIXED used (Total process time):
      real time           0.03 seconds
      cpu time            0.03 seconds

二维码

扫码加我 拉你入群

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

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

全部回复
2012-2-28 14:44:51
第48条观测是不是异常值?或者该值不满足模型参数计算的条件.
二维码

扫码加我 拉你入群

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

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

2012-2-28 19:42:04
应该不是的。我用proc countreg 这个模块就能做呢。
二维码

扫码加我 拉你入群

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

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

2012-2-29 09:37:46
是在写毕业论文吧?近几年零膨胀模型热起来了啊~哈哈。我是用proc mixed过程去做的
二维码

扫码加我 拉你入群

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

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

2012-2-29 10:22:58
it must be a lot of pains of using proc mixed to do zero-inflated...jingu
二维码

扫码加我 拉你入群

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

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

2012-2-29 10:29:56
bridog 发表于 2012-2-28 19:42
应该不是的。我用proc countreg 这个模块就能做呢。
good question.
jingju
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群