r - 带有 coxme 对象的转发过程 (R)

Question

我正在研究心脏移植数据，我必须在R. 由于它是为此目的而设计的，因此我选择使用该coxme软件包，但在选择模型时遇到了很多麻烦。我想使用转发过程，但它不适用于coxme对象。我在这里问了另一个问题，他的答案允许我访问step()'s 的代码来修改它，以便它适用于coxme模型。但是，这非常困难（至少对我而言），因为功能相当复杂（绝对不仅仅是修改exctractAIC()，还有许多其他功能需要适应）。

所以这是我的问题：

1) 有没有更简单的方法在模型上运行前锋coxme？

2）我应该使用另一个包吗？

Answer 1

在尝试了很长时间来适应step()之后，我决定从一开始就创建一个简化但精确的函数，这要简单得多。所以如果有人想知道和我一样的问题，我真的推荐他：

1）从一开始就不怕编码整个功能！（原始函数的 80% 对前向过程没有用处，不需重新编码）

2）使用coxme包。

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