Я надеюсь создать нелинейную модель со смешанным эффектом, используя пакет nlme.У меня есть три случайных эффекта: shoot вложен в tree, который вложен в site.У меня также есть три фиксированных эффекта: a0, b1 и b2
Это то, что я собрал:
fun_1 <- n ~ a0 * length^b1 * exp(b2 * age)
mm1 <- nlme(fun_1,
data = branches_new_counts,
fixed = a0 + b1 + b2 ~ 1,
random = pdDiag(a0 + b1 + b2 ~ 1|dataset/tree/shoot_no),
start = coef(fixed_model_final))
Я уверен, что я указал случайныйнеправильно структурирует структуруЛюбые советы искренне приветствуются.