请教SAS/STAT 中的PROC GENMOD的问题

shiyiming

<!-- s:( --><img src="{SMILIES_PATH}/icon_sad.gif" alt=":(" title="Sad" /><!-- s:( --> 最近我用SAS中的PROC GENMOD 过程分析一批类似重复测量的资料，我的兴趣主要是类内的相关结构，由于因变量为二分类资料，所以我用了ALR的方法求odds ratio值作为相关结构的测量。但Log窗口中出现：
WARNING:The generalized Hessian matrix is not positive definite.Iteration will be terminated.
Error:Error in parameter estimate covariance computation.
Error:Error in estimation routine.
或者：
WARNING: The generalized Hassian matrix for the log odds ratio is
singular.

我医学出身，数学基础不好。请教高人，出现这种情况的可能原因是什么？这是否提示了我的资料有什么特殊的情况？

非常感谢任何宝贵的建议！

shiyiming

不知是否方便提供你的程序和部分数据？

似乎是程序写的有问题。

shiyiming

谢谢Jerry对该问题的关注！

程序应该没有问题，因为只是某些关联结构的组合时才出现错误。
程序示意如下：

proc genmod data=analyse desc;
      class center;
      model y= x1 x2  x3 /  dist=bin ;
      repeated  subject=center / maxiter=1000
                                     logor=zfull
                                     zdata=z_matrix
                                     zrow =
          (asso1 asso2 asso3 asso4  )
                                     ypair=(t1 t2);
run;

其中analyse是要分析的数据集，包括因变量y（1，0变量），自变量x1,x2,x3...，分组变量center等等。z_matrix是关联结构矩阵的数据集，包括asso1，asso2，asso3 ，asso4...等等表示类内不同的关联结构，t1 和t2是指示变量。
举例来说，我发现如果关联结构只选择asso1－asso4，程序还是可以正常运行的并得到参数估计的。如果添加asso5或者asso6，就会出现错误。另外，最近发现，选择不同的自变量组合有时也会出现相同的错误。

由于一些非个人的原因，不能提供原始数据，望谅。不过这种非独立结构的资料很普通，也可以在SAS中模拟一些。如果这个问题很常见的话，应该在模拟数据集的分析中也会碰到。

再次感谢！

shiyiming

I guess some of your variables are highly correlated, you may need to delete some  variables.

shiyiming

前四个asso是哪些变量的协方差阵？

shiyiming

asso并不是变量的协方差阵，它表示类内观测对子间共享的协变量。