SAS中文论坛

标题: 请教SAS/STAT 中的PROC GENMOD的问题 [打印本页]

作者: shiyiming 时间: 2003-10-7 00:39
标题: 请教SAS/STAT 中的PROC GENMOD的问题
<img src="{SMILIES_PATH}/icon_sad.gif" alt=":(" title="Sad" /> 最近我用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 时间: 2003-10-7 20:42
不知是否方便提供你的程序和部分数据？

似乎是程序写的有问题。

作者: shiyiming 时间: 2003-10-7 21:10
标题: to Jerry
谢谢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 时间: 2003-10-8 12:11
I guess some of your variables are highly correlated, you may need to delete some variables.

作者: shiyiming 时间: 2003-10-10 17:53
前四个asso是哪些变量的协方差阵？

作者: shiyiming 时间: 2003-10-11 14:55
asso并不是变量的协方差阵，它表示类内观测对子间共享的协变量。

欢迎光临 SAS中文论坛 (http://mysas.net/forum/)