library(nlme)

rd<-read.csv(“data.csv”)

fm.nlm<-nlme(total~a*db^b,

data=rd,

fixed = list(a+b ~ 1),

random=?

start=c(a=0.01,b=2)

)

数据是region、species和sample的三重拟合，我想设定的随机效应有四种情况：

1         只region作为随机效应

2         只species作为随机效应

3         对region和 species作为随机效应

4         对region和 species作为随机效应并认为之间有交互

不知对着四种情况怎么做随机效应（希望能讲解清楚详细一点，刚学混合模型）数据在附件中