Можно ли использовать кластерную функцию дважды в регрессии Кокса в R?

Итак, у меня есть запрос 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
person E.D.D    schedule 20.09.2015    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Я не знаю конкретного ответа на ваш вопрос о кластеризации двух переменных в 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)
person Whitebeard    schedule 20.09.2015