R中save()与load()操作致数据不一致的问题排查与解决求助

问题背景 #

在R 4.2.1（Windows 10/CentOS 7环境）中，使用save()/load()序列保存加载lme4构建的线性混合效应模型后，模型对象无法通过identical()验证，all.equal()提示pp和resp组件的类定义不一致，导致结果无法精准复现。

可复现代码 #

library(lme4)

# 构建混合模型
model <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
model2 <- model

# 初始状态验证：完全一致
identical(model, model2) # 返回TRUE
all.equal(model, model2) # 返回TRUE

# 保存并重新加载
save(model2, file="./test.RData")
load("./test.RData")

# 加载后验证：不一致
identical(model, model2) # 返回FALSE
all.equal(model, model2) 
# 输出：
# [1] "Attributes: < Component “pp”: Class definitions are not identical >"  
# [2] "Attributes: < Component “resp”: Class definitions are not identical >"

解决方案/绕过方法 #

1. 使用saveRDS()/readRDS()替代save()/load() #

这是最稳妥的方案，saveRDS()针对单个对象进行序列化，能更精准地保存对象的类结构和属性，加载后可完全通过identical()验证：

# 替代save
saveRDS(model2, file="./test.rds")
# 替代load
model2 <- readRDS("./test.rds")

# 验证
identical(model, model2) # 返回TRUE
all.equal(model, model2) # 返回TRUE

2. 手动修复加载后的类定义（临时绕过） #

如果必须使用save()/load()，加载后可以手动同步pp和resp组件的类定义：

load("./test.RData")
# 同步类定义
class(model2$pp) <- class(model$pp)
class(model2$resp) <- class(model$resp)

# 验证
identical(model, model2) # 返回TRUE

注意：此方法仅针对当前问题场景，若模型其他组件存在类似问题可能需要额外处理。

3. 更新lme4包至最新稳定版 #

旧版本lme4可能存在序列化相关的bug，更新到最新版（如1.1-34及以上）可修复部分此类问题：

install.packages("lme4")

新版本R的情况验证 #

在R 4.3.x及以上版本中，配合最新版lme4：

• 使用saveRDS()/readRDS()完全不存在该问题；
• save()/load()的类定义不一致问题也已被修复，加载后的模型可通过identical()验证。

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