使用nlme(非线性混合效应模型)的适当的嵌套随机效应结构

时间:2019-01-18 01:39:49

标签: r nlme random-effects

我希望使用nlme包创建一个非线性混合效果模型。我有三种随机效果:shoot嵌套在tree中,嵌套在site中。我还具有三个固定效果:a0b1b2

这是我整理的:

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))  

我确定我没有正确指定随机效果结构。衷心感谢任何建议。

0 个答案:

没有答案