零膨胀模型,sas代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我用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

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：sas代码 proc nlmixed observation Procedure EXECUTION alpha 模型 general

第48条观测是不是异常值？或者该值不满足模型参数计算的条件.

应该不是的。我用proc countreg 这个模块就能做呢。

是在写毕业论文吧？近几年零膨胀模型热起来了啊~哈哈。我是用proc mixed过程去做的

it must be a lot of pains of using proc mixed to do zero-inflated...jingu

bridog 发表于 2012-2-28 19:42
应该不是的。我用proc countreg 这个模块就能做呢。

good question.
jingju

呵呵!谢谢了，问题已经解决了！

bridog 发表于 2012-3-1 10:45
呵呵!谢谢了，问题已经解决了！

请问你最后是怎么解决的 我也出现了类似的问题

零膨胀模型请参考http://en.wikipedia.org/wiki/Zero-inflated_models