各位好：

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

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