你可以在 R 的 cox 回归中使用两次集群函数吗？答案

【问题标题】：Can you use the cluster function twice in a cox regression in R?你可以在 R 的 cox 回归中使用两次集群函数吗？
【发布时间】：2015-09-20 22:58:10
【问题描述】：

所以我有一个 cox 查询 - 所以我正在尝试做一个 cox 回归，我想知道我是否可以将数据聚类两次：

cox1<-coxph(formula = Surv(time.to.arrive..trap, Arrived.1.or.0) ~     
            Is.feeder.control.or.stimulus + cluster(id.location) + cluster(New.ID), 
            data = all.data)
cox1

我有一个数据集，其中有 200 只鸡，我在树林中的五个位置一个接一个地设置了成对的喂食器（控制和刺激）。（在任何时候只有一个控制和一个刺激馈线）

我正在尝试确定到达控制或刺激喂食器的时间是否存在差异 - 但是我想考虑到鸟类可以在我的数据中出现多次 - （在控制和刺激时，或者五个地点中的每一个）。两次使用集群功能有意义吗？在上面的公式中，我根据鸟的 ID 和 location.Id 对它进行了聚类。

这似乎是一个显而易见的问题 - 但给我带来了无穷无尽的麻烦！非常感谢您对此事的任何启发！

亲切的问候

【问题讨论】：

标签： r cox-regression

【解决方案1】：

我不知道您关于在coxph 中对两个变量进行聚类的问题的具体答案，但也许您可以在id.location 和New.ID 之间创建一个交互因子，然后对该变量进行聚类（并调用cluster只有一次）

例如，

all.data$newvar <- interaction(all.data$id.location, all.data$New.ID)
cox1<-coxph(formula = Surv(time.to.arrive..trap, Arrived.1.or.0) ~     
        Is.feeder.control.or.stimulus + cluster(newvar), data = all.data)

【讨论】：