Я пытаюсь распечатать «остатки между событиями» (которые являются условными режимами случайного эффекта) из модели смешанных эффектов в файл .csv. То, что я не смог получить, это то, что длина вектора ранфа не совпадает с факторным уровнем моих данных, я верю, что приведенный ниже код объяснит все, что я не смог.
> fm_rjb <- lmer(log(data[[periods[[8]]]]) ~ 1 + data$M1 + data$M2 + data$M3
+ data$Nflag + data$Rflag + data$Dist1_RJB + data$Dist2_RJB + data$VN
+ data$depth1 + data$depth2 + ( 1 | data$STID) + ( 1 | data$EQID), data,
REML = FALSE)
> length(ranef(fm_rjb)[[2]][,1])
[1] 155
> length(as.numeric(levels(factor(data$EQID))))
[1] 159
Мне нужноэти два должны быть одинаковой длины.