我正在处理多级数据，需要建立多层三分类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语句有错，请教高人指点，应该如何同时设定几个模型的随机截距方差的分布呢？我急用在我的毕业论文中，如有高手能解答，我将万分感激！