我正在研究心脏移植数据,我必须在R
中建立一个Cox混合模型。因为它是为此目的而设计的,所以我选择使用coxme
包,但是我很难选择我的模型。我想使用正向过程,但它不适用于coxme
个对象。我在这里问了另一个问题,他的答案允许我访问step()
的代码进行修改,因此它适用于coxme
模型。然而,这对于我来说非常困难(至少对我而言)因为函数非常复杂(它绝对不仅仅是修改exctractAIC()
,还有许多其他函数需要适应)。
所以这是我的问题:
1)是否有更简单的方法在coxme
模型上运行前进?
2)我应该使用另一个套餐吗?
答案 0 :(得分:0)
在尝试了很长时间以适应step()
之后,我决定从头开始创建一个简单但精确前进的功能,这更简单。因此,如果有人想到和我一样的问题,我真的推荐他:
1)从开始就不怕编码整个功能! (原始功能的80%对于正向程序无用,不需要重新编码)
2)使用coxme
包。
然而,如果你这样做,这里有一件事可以防止混淆:coxme模型中所谓的AIC / BIC实际上是测试模型与null模型相比的改进。这意味着如果您决定使用它,您应该最大化它而不是最小化它(如果您愿意,可以将计算更改为真正的AIC / BIC)。