如何在R中用glmmTMB拟合随机效应模型？

Question

我想用R中的glmmTMB函数拟合一个随机效应模型，我指定模型为模型= glmmTMB(Y ~ (1+x1+x2|group))。估计结果仍然包含存储在obj$env$last.par.best1中的作为beta的通用截距，即使我没有将1作为随机部分之外的通用截距。有人知道如何用glmmTMB来拟合随机效应模型吗？谢谢!

我之所以想要一个随机效应模型，是因为我有一个包含很多组的面板数据，每个组的y和x之间有一个非常不同的关系，所以基本上我只想为每个组建立一个模型。我以为glmmTMB能帮我做到这一点。

Answer 1

我不明白其中的动机，但在公式中使用-1将常量从固定效果设计矩阵中删除的任何常见包中实现这一点是非常简单的。下面是一个使用来自lme4::sleepstudy的数据的示例

library('glmmTMB')
df <- lme4::sleepstudy
names(df) <- tolower(names(df))

m1 <- glmmTMB(formula = reaction ~ -1 + (days|subject), data = df)
summary(m1)
# Family: gaussian  ( identity )
# Formula:          reaction ~ -1 + (days | subject)
# Data: df
# 
# AIC      BIC   logLik deviance df.resid 
# 1840.8   1853.6   -916.4   1832.8      176 
# 
# Random effects:
#   
#   Conditional model:
#   Groups   Name        Variance Std.Dev. Corr 
# subject  (Intercept) 63770.1  252.53        
# days          142.2   11.93   0.88 
# Residual               654.9   25.59        
# Number of obs: 180, groups:  subject, 18
# 
# Dispersion estimate for gaussian family (sigma^2):  655