求教：sas nlmixed 命令拟合多水平模型报错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位好：

    我想用sas   nlmixed 命令拟合多水平模型，反映AGE（中心化定量变量）, GENDER，ECONOMY1， ECONOMY2（此3个为2分类变量）对HBP（2分类变量）的影响大小，其中ECONOMY1 ECONOMY2为二水平协变量，AGE, GENDER 为一水平协变量。

sas程序如下：

proc nlmixed data=anlyzem;
parms B0=0 B1=0 B2=0 B3=0 B4=0 V_u0=0 C_u03=0 V_u3=0 C_u04=0 C_u34=0 V_u4=0;
Z=B0+B1*AGE+B2*GENDER+B3*ECONOMY1+B4*ECONOMY2+u0j+u3j*ECONOMY1+u4j*ECONOMY2;
If (hbp=1) then P=1/(1+exp(-Z));
Else P=1-(1/(1+EXP(-Z)));
LL=LOG(P);
MODEL HBP~GENERAL(LL);
RANDOM u0j u3j u4j~normal([0,0,0],[V_u0,C_u03,0,C_u04,0,0]) SUBJECT=HOSPITAL_N;
RUN;

运行后sas日志提示为非正定矩阵。


上述程序去掉parms语句后，个参数默认为1，运行后sas日志提示： Quadrature accuracy of 0.000100 could not be achieved with 31 points.  The achieved accuracy was 1.000000.


我用mlwin分析数据后得到个参数值，以此为基础设定parms语句中各参数值，结果也是提示：Quadrature accuracy of 0.000100 could not be achieved with 31 points.  The achieved accuracy was 1.000000.

请问这是什么意思，我该怎么做。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：nlmixed 多水平模型 Mixed 多水平 mix 模型 SAS 命令 拟合 nlmixed

一个变量,一个变量地加.
proc nlmixed data=anlyzem;
parms B0=1 B1=1;
Z=B0+B1*AGE;
.......
如果不行,从B0开始,每次加一个变量.