使用DHARMa检验Gamma分布glmmTMB模型时遇family方法报错求助

问题详情 #

我在分析延迟数据时，用Gamma分布构建了glmmTMB模型m1，代码如下：

m1 <- glmmTMB(formula = Lat.Min ~ Fruta,
              family  = Gamma(link = "log"),
              data    = datos)

library(DHARMa)
simulationOutput <- simulateResiduals(fittedModel = m1, refit = T, plot = T) # 注：原代码注释为西班牙语“no converge para”（未收敛）
hist(simulationOutput)
testDispersion(simulationOutput)
simulate

运行DHARMa的残差分析步骤时，触发如下错误：

Error in UseMethod("family") : no applicable method for 'family' applied to an object of class "glmmTMB"

已知DHARMa官方支持glmmTMB模型，但不清楚报错原因，希望得到解决方案。

可能的解决方案 #

1. 检查并更新包版本 #

DHARMa对glmmTMB的支持需要版本匹配，旧版本的包可能存在兼容性问题。执行以下命令更新到最新稳定版：

install.packages(c("glmmTMB", "DHARMa"))

如果使用开发版功能，可从GitHub安装最新迭代版本：

# 若未安装devtools，先执行 install.packages("devtools")
devtools::install_github("glmmTMB/glmmTMB")
devtools::install_github("florianhartig/DHARMa")

更新后重启R会话，重新拟合模型并尝试DHARMa分析。

2. 调整包的加载顺序 #

包的加载顺序可能导致S3方法注册冲突。尝试先加载DHARMa，再加载glmmTMB，之后重新拟合模型：

# 重启R会话后执行
library(DHARMa)
library(glmmTMB)

# 重新拟合模型
m1 <- glmmTMB(formula = Lat.Min ~ Fruta,
              family  = Gamma(link = "log"),
              data    = datos)

# 再次运行DHARMa分析
simulationOutput <- simulateResiduals(fittedModel = m1, refit = T, plot = T)

3. 验证模型拟合有效性 #

报错可能源于模型本身未正确拟合。先检查模型的收敛状态和基本信息：

# 查看模型摘要
summary(m1)

# 检查模型是否收敛
is_converged <- glmmTMB:::isConv(m1)
print(is_converged)

如果模型未收敛（原代码注释也提到“No converge para”），需要先解决模型收敛问题：比如检查数据是否存在异常值、调整模型结构（如添加随机效应）、尝试其他链接函数或分布，确保模型拟合成功后再进行残差分析。

4. 手动注册family方法（临时 workaround） #

如果上述方法都无效，可以手动为glmmTMB类对象注册family方法：

# 定义family方法
family.glmmTMB <- function(object, ...) {
  object$modelInfo$family
}

# 注册S3方法
registerS3method("family", "glmmTMB", family.glmmTMB)

注册完成后，重新运行DHARMa的残差分析命令。

内容的提问来源于stack exchange，提问作者Franco Pacelli