使用coxme对象的前向程序(R)

时间:2013-06-27 13:10:27

标签: r survival-analysis mixed-models cox-regression

我正在研究心脏移植数据,我必须在R中建立一个Cox混合模型。因为它是为此目的而设计的,所以我选择使用coxme包,但是我很难选择我的模型。我想使用正向过程,但它不适用于coxme个对象。我在这里问了另一个问题,他的答案允许我访问step()的代码进行修改,因此它适用于coxme模型。然而,这对于我来说非常困难(至少对我而言)因为函数非常复杂(它绝对不仅仅是修改exctractAIC(),还有许多其他函数需要适应)。

所以这是我的问题:

1)是否有更简单的方法在coxme模型上运行前进?

2)我应该使用另一个套餐吗?

1 个答案:

答案 0 :(得分:0)

在尝试了很长时间以适应step()之后,我决定从头开始创建一个简单但精确前进的功能,这更简单。因此,如果有人想到和我一样的问题,我真的推荐他:

1)从开始就不怕编码整个功能! (原始功能的80%对于正向程序无用,不需要重新编码)

2)使用coxme包。

然而,如果你这样做,这里有一件事可以防止混淆:coxme模型中所谓的AIC / BIC实际上是测试模型与null模型相比的改进。这意味着如果您决定使用它,您应该最大化它而不是最小化它(如果您愿意,可以将计算更改为真正的AIC / BIC)。