急！！！请教高手关于NLMIXED和GLIMMIX程序包的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我正在处理多级数据，需要建立多层三分类LOGIT模型，在用到NLMIXED GLIMMIX两个程序包时，总是出错，我的程序如下：

/*空模型 nlmixed*/
PROC NLMIXED data=my.data2008100404;
PARMS a0=0 V_u0a=1 b0=0 V_u0b=1;
z1=a0 + u0a;
z2=b0 + u0b;
if(x1tran=3) then p=1/(1 + exp(z1) + exp(z2));
else if(x1tran=1) then p=exp(z1)/(1 +exp(z1) + exp(z2));
else if(x1tran=2) then p=exp(z2)/(1 +exp(z1) + exp(z2));
ll=log(p);
Model x1tran~general(ll);
RANDOM u0a,u0b~normal([0,0],[V_u0a,0,V_u0b]) subject=villageid;
ESTIMATE 'ICCa' V_u0a/(V_u0a+3.289868134);
Run;

似乎RANDOM语句有错，请教高人指点，应该如何同时设定几个模型的随机截距方差的分布呢？我急用在我的毕业论文中，如有高手能解答，我将万分感激！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：glimmix nlmixed Mixed 请教高手 GLIM general normal 程序 模型

不会，望高手解答

可不可以把NLMIXED和GLIMMIX程序包copy一份给我？
邮箱：

xiangqfeng@163.com

1# zerohero

RANDOM u0a,u0b~

no comma here.