Я использовал пакет coxme для оценки модели пропорциональной опасности Кокса со случайными эффектами (просто перехватами). Однако у меня возникают проблемы, когда мне нужно исправить дисперсию двух случайных перехватов. С одним перехватом отлично работает следующее:
require(coxme)
y = sample(1:10, 500, replace=TRUE)
cens = sample(0:1, 500, replace=TRUE)
id.1 = sample(1:20, 500, replace=TRUE)
id.2 = sample(1:20, 500, replace=TRUE)
data = cbind(y,cens,id.1,id.2)
f.values = list(id.1=c(0.3))
fit.1 = coxme(Surv(y, cens) ~ (1|id.1), data = data.red, vfixed = f.values)
summary(fit.1)
Однако я просто не могу найти информацию о том, что делать, когда я хочу исправить два случайных эффекта. Я получаю следующую ошибку:
f.values = list(id.1=c(0.3),id.2=c(0.1))
fit.2 = coxme(Surv(y, cens) ~ (1|id.1) + (1|id.2), data = data.red, vfixed = f.values)
Error in vfixed[[i]] : subscript out of bounds
Я предполагаю, что я неправильно использую здесь объект f.values, но не знаю, как это сделать в противном случае.
Любая помощь приветствуется, Алекс