R语言中如何保存与加载SBCK包的拟合后R6类对象？

问题根源 #

直接用saveRDS()保存拟合后的dOTC R6对象时，R6的**活动绑定（active bindings）**会被序列化为普通属性，加载后丢失原R6类的结构，变成仅包含活动绑定的普通对象，无法直接调用原类的方法。

解决方案：提取核心拟合参数重新构建实例 #

R6对象的拟合结果实际存储在其内部的静态属性中（而非活动绑定），我们可以提取这些核心参数保存，加载后重新初始化R6实例并赋值参数，恢复拟合状态。

1. 保存拟合后的核心参数

先查看拟合后dOTC对象的关键属性（用str(dotc_fitted)或ls(dotc_fitted, all.names = TRUE)），提取fit()方法生成的静态参数：

# 假设拟合后的对象是 dotc_fitted
# 提取核心拟合参数（根据SBCK::dOTC的实际结构调整）
dotc_save_data <- list(
  params = dotc_fitted$params,
  X0 = dotc_fitted$X0,
  Y0 = dotc_fitted$Y0,
  fitted_flag = TRUE  # 标记对象已完成拟合
)

# 保存到磁盘
saveRDS(dotc_save_data, "dOTC_fitted_model.rds")

2. 加载并恢复R6实例

加载保存的参数后，重新初始化dOTC对象，将参数赋值给实例，并标记为已拟合：

# 加载保存的参数
dotc_load_data <- readRDS("dOTC_fitted_model.rds")

# 初始化新的dOTC实例
dotc_loaded <- SBCK::dOTC$new()

# 赋值拟合参数
dotc_loaded$params <- dotc_load_data$params
dotc_loaded$X0 <- dotc_load_data$X0
dotc_loaded$Y0 <- dotc_load_data$Y0
dotc_loaded$fitted <- dotc_load_data$fitted_flag

# 验证：调用预测方法测试
# predict_result <- dotc_loaded$predict(X_new)

补充说明 #

• 如果不确定哪些参数是必要的，可以查看SBCK包中dOTC类的源代码，重点关注fit()方法中修改的属性；也可以用R6::is.active_binding(dotc_fitted, "属性名")区分活动绑定和静态属性，只保存静态属性。
• 避免直接保存整个R6对象，活动绑定是动态计算的属性，不需要单独保存，加载后R6类会自动重建这些绑定。

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