SAS中文论坛

标题: LOGISTIC回归模型不收敛如何解决 [打印本页]

作者: chenli9045 时间: 2014-7-24 14:06
标题: LOGISTIC回归模型不收敛如何解决
求帮忙~~~
我在做LOGISTIC回归，结局变量y是二分类：1=有效，2=无效。
自变量有三个，分别为：x1-组别，取值1、2、3；x2-中心，取值1到11；x1*x2-表示中心和组别的交互效应。
可是结果显示如下，请问该如何解决？

作者: mono 时间: 2014-7-25 08:24
1. Lositic 结局变量最好写成0,1的格式
2. 其实你的自变量不是三个，应该是35个，X1， X2，还有33个 X1*X2，如果你的样本量很少，就很难converge了。
3.看看能不能合并一下你的变量
4. 11中心用dummy variale (指定一个reference 中心）
5. 用 PENALIZED MAXIMUM LIKELIHOOD ESTIMATION
6. 用Baysian.
7. 最后一招，你也可以强制SAS covernge, 比如让 Covergence > 1E-4 (默认值）

作者: chenli9045 时间: 2014-7-25 09:17

mono 发表于 2014-7-25 08:24
1. Lositic 结局变量最好写成0,1的格式
2. 其实你的自变量不是三个，应该是35个，X1， X2，还有33个 X1*X2 ...

谢谢你的回答。我一会试一下你说的解决办法。。。
没错的，自变量是35个，总的样本量是自变量个数的10倍吧。
我找到不收敛的原因了，因为有一个中心某一组结局变量的取值全是1，所以极大似然估计的变量的标准误特别大，不可用，出现了模型不收敛。说是可以看似然比检验的结果，或者是用精确LOGISTIC回归。

作者: mono 时间: 2014-7-25 12:14

chenli9045 发表于 2014-7-25 09:17
谢谢你的回答。我一会试一下你说的解决办法。。。
没错的，自变量是35个，总的样本量是自变量个数的10倍 ...

你这种研究应该使用random effect logistic regression 而不是一般的logistic regression

作者: chenli9045 时间: 2014-7-28 09:09

mono 发表于 2014-7-25 12:14
你这种研究应该使用random effect logistic regression 而不是一般的logistic regression ...

random effect Logistic 没研究过，能否推荐点参考资料之类的？

作者: mono 时间: 2014-7-28 13:50
本帖最后由 mono 于 2014-7-28 13:52 编辑

chenli9045 发表于 2014-7-28 09:09
random effect Logistic 没研究过，能否推荐点参考资料之类的？

因为你是多中心的研究，要考虑各中心的cluster effect, 你这东西如果写成论文，只要提到多中心，肯定会有人建议你用random effect model, 比如GLIMMIX 或GEE。具体材料自己搜搜吧。

作者: chenli9045 时间: 2014-7-28 16:55

mono 发表于 2014-7-28 13:50
因为你是多中心的研究，要考虑各中心的cluster effect, 你这东西如果写成论文，只要提到多中心，肯定会有 ...

原来你说的是这个方法啊，以前研究的还真是少。
总之，谢谢你的回答。

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