SAS proc mixed 过程 运行中日志提示满足收敛准则，但最终的 Hessian 不是正定的

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位好：        在下请教诸位。
       我的研究目的是想分析各国gdp, 老年人口比例等因素对于脑血管病发病风险影响。
       使用的变量：
     1，时间范围：1990-2016年，year;
                          4个国家，country；
      2.连续型变量：人均GDP，GDP_per;
                          老龄人口比例，oldman；
                           人均BMI，BMI；
                           人均吸烟率Smoking
                           脑血管病发病率，incidence。
       共108条记录， 分析前我们对连续型变量进行了标准化正态变换处理。

       使用proc mixed 过程拟合模型，程序如下：
proc mixed data = isdata2 method=reml covtest ;
      class country;
      Model incidence = country year gdp_per   oldman  bmi  smoking /S ddfm=KR;
      random int year/ subject=country G type=UN;
run;

    运行后能够出模型结果，但日志总是提示： 满足收敛准则，但最终的 Hessian 不是正定的。
    请问这是怎么回事？该如何解决？

请诸位赐教，谢谢！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：proc mixed Hessian Mixed HESS ROC

Try adding a PARMS statement to your PROC MIXED code. Either use the OLS option to specify ordinary least squares starting values (rather than the default MIVQUE0 values) or specify your own starting values. The syntax requesting the OLS starting values is:

   parms / ols;

The syntax for specifying your own starting values is a little more complex. See the syntax section of the PROC MIXED chapter in the SAS/STAT User's Guide.

If you still get this message with new starting values, then there are probably linear dependencies in the parameters. This means that the model is too complex for your data to support. Try simplifying the covariance structure until you no longer receive this message.

The results from a model that has a nonpositive definite Hessian are not reliable and should be interpreted with extreme caution.

ArisQ 发表于 2020-4-1 12:40
Try adding a PARMS statement to your PROC MIXED code. Either use the OLS option to specify ordinary  ...

谢谢你的建议。