proc nlmixed（warning：The final Hessian matrix is full rank ）

zixinye

用proc nlmixed拟合模型时，日志出现警告：The final Hessian matrix is full rank but has at least one negative eigenvalue. Second-order optimality condition violated.

这怎么解释啊？程序是照着参考书写的。

micro@

Stuck at a saddle point

zixinye

回复第2楼的 micro@：我拟合只有截距项的模型。程序如下

proc nlmixed data=chf.censor;

parms b0=0 sigma=1;

y=b0+u0j;

p=1/(1+exp(-y));

model e1~binary(p);

random u0j~normal(0,sigma) subject=hospital;

estimate 'icc' sigma/(sigma+3.289868134);

run;

你的回答我还是不太明白，能不能说的具体点？非常感谢！