1

所以我有一个 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 对它进行了聚类。

这似乎是一个显而易见的问题——但给我带来了无穷无尽的麻烦！非常感谢您提前为您提供任何关于此事的信息！！

亲切的问候

4

1 回答 1

0

我不知道您关于在 内聚类两个变量的问题的具体答案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)
于 2015-09-20T23:27:26.750 回答