glmmTMB拟合报错：初始参数值下负对数似然为NaN求排查

针对你用glmmTMB拟合果蝇发育时间混合模型时遇到的报错，可按以下步骤排查解决：

• 检查数据分布与模型类型适配性
  发育时间为连续型数据，但如果存在零值、极端偏态或不符合正态假设，默认高斯模型易出现数值问题。先查看数据特征：

  summary(tddcVZ$age)
  hist(tddcVZ$age)
  

  若数据右偏明显，尝试对数转换后拟合：

  model_tddVZ <- glmmTMB(log(age) ~ treatment + (1|vial), tddcVZ)
  

  若存在零值，改用截断高斯模型：

  model_tddVZ <- glmmTMB(age ~ treatment + (1|vial), tddcVZ, family = truncated_gaussian())
  

• 验证随机效应分组的样本量
  若单个vial分组仅含1个观测，或分组总数过少（如少于5组），会导致随机效应方差估计不稳定。统计每个分组的观测数：

  table(tddcVZ$vial)
  

  对单样本分组，可考虑合并或移除后重试。

• 手动指定初始参数
  初始参数值不合理可能触发似然值为NaN，可手动设置随机效应方差的初始值（theta为方差的对数形式）：

  model_tddVZ <- glmmTMB(age ~ treatment + (1|vial), tddcVZ,
                         start = list(theta = log(0.1)))
  

• 清理数据中的异常值与缺失值
  NA、Inf或极端值会干扰似然计算，先清理数据：

  # 移除缺失值
  tddcVZ_clean <- na.omit(tddcVZ)
  # 移除age的极端异常值（基于四分位距）
  q1 <- quantile(tddcVZ_clean$age, 0.25)
  q3 <- quantile(tddcVZ_clean$age, 0.75)
  iqr <- q3 - q1
  tddcVZ_clean <- subset(tddcVZ_clean, age >= q1 - 1.5*iqr & age <= q3 + 1.5*iqr)
  # 重新拟合模型
  model_tddVZ <- glmmTMB(age ~ treatment + (1|vial), tddcVZ_clean)
  

• 更新glmmTMB至最新版本
  即使重装过包，也可能未安装最新版，尝试更新：

  install.packages("glmmTMB", dependencies = TRUE)
  

  若CRAN安装失败，安装GitHub开发版：

  install.packages("devtools")
  devtools::install_github("glmmTMB/glmmTMB")
  

内容的提问来源于stack exchange，提问作者Lucas Kreiman