我阅读了 glmmTMB 包小插曲(https://cran.r-project.org/web/packages/glmmTMB/vignettes/covstruct.html)并有以下问题:
- 在小插图中,他们使用以下方法拟合模型
glmmTMB(y~ar1(times+0|group),data=dat0)
并提到time+0
对应于连接观察向量 y(行)和随机效应向量 u(列)的设计矩阵 Z"
是什么意思+0
?和有什么区别(times|group)
吗?(times+1|group)
(1|group)
关于协方差结构的语法有什么全面的总结吗?
- 如果我想在 R 函数生成结果 y_{ij} 时拟合负二项式模型:
rnbinom(mu=x_{ij}beta+b_i+e_{ij},size=1)
,其中 i 是组索引,j 是个体索引,并且 b_i~N(0,1), e_{ij}~ N(0,1)。以下代码会正确指定模型吗?
dat <- data.frame(y,x,group)
glmmTMB(y~x+(1|group),data=dat,family=nbinom2)
任何建议和帮助表示赞赏。在此先感谢您的帮助!