使用lme4::lmer拟合混合效应模型时遇dgTMatrix维度赋值错误

最近用lme4包的lmer()拟合混合效应模型时踩了个特别隐蔽的坑，报错信息完全摸不着头脑，谷歌搜遍了都没找到完整匹配的结果，最后靠回滚代码对比才揪出问题，特意记录下来帮大家避坑：

遇到的错误信息 #

Error in (function (cl, name, valueClass) : assignment of an object of class "numeric" is not valid for @'Dim' in an object of class "dgTMatrix"; is(value, "integer") is not TRUE

问题根源 #

这个错误的核心原因是混合效应公式中的分组变量被意外转换成了numeric类型（正确的类型应该是factor或integer）。当分组变量是numeric时，lme4在构建稀疏矩阵（dgTMatrix类型）时，会错误地将矩阵的Dim属性设置为numeric类型，而dgTMatrix要求Dim必须是integer类型，因此触发了这个赋值错误。

复现代码 #

错误触发场景（模拟预处理时的误操作） #

library(lme4)
data(sleepstudy)

# 模拟数据处理时的错误：将分组变量Subject转为numeric
sleepstudy$Subject <- as.numeric(sleepstudy$Subject)

# 拟合模型时触发目标错误
fm1 <- lmer(Reaction ~ Days + (0 | Subject), data = sleepstudy)

正确运行场景 #

library(lme4)
data(sleepstudy)

# 确保分组变量为factor（sleepstudy默认已为factor，显式转换更稳妥）
sleepstudy$Subject <- as.factor(sleepstudy$Subject)

# 正常拟合模型
fm1 <- lmer(Reaction ~ Days + (0 | Subject), data = sleepstudy)

排查思路 #

因为错误信息完全没关联到分组变量类型，一开始根本想不到问题出在哪。后来通过逐步回滚代码，对比能正常运行的版本和报错版本的差异，才发现是分组变量被不小心转成了numeric。如果大家遇到类似的dgTMatrix相关错误，优先检查：

• 混合效应中的分组变量类型是否为factor或integer
• 数据预处理过程中有没有误转换变量类型

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