我希望使用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))
我确定我没有正确指定随机效果结构。衷心感谢任何建议。