为何无法将SSMcustom状态空间对象传入KFAS的SSModel公式？

为什么会报错？ #

R的公式系统（~表达式）处理变量时，默认会尝试从当前环境提取变量并构建model.frame——这个流程只支持向量、因子这类适合做预测变量的简单类型，而SSMcustom()返回的是状态空间组件列表，属于复杂对象，自然过不了model.frame的类型检查。

而直接在SSModel()公式里写SSMcustom(...)时，SSModel的内部逻辑会特殊识别这类KFAS包的状态空间构造函数，跳过普通的公式变量解析流程，直接把函数返回的组件整合到模型里，所以不会触发错误。

解决方法 #

方法1：用I()包裹变量 #

I()是R的“恒等函数”，作用是告诉公式解析器：这个对象不要当成普通变量处理，直接原样传递。修改代码如下：

ssm <- SSMcustom(Z = Z_t, T = T_t, R = R_t, Q = Q_t, a1 = a1_t, P1 = P1_t)
model_fixed <- KFAS::SSModel(ts ~ -1 + I(ssm))

方法2：用substitute()动态构造公式 #

如果需要更灵活的变量插入，可以用substitute()生成包含变量的公式对象：

ssm <- SSMcustom(Z = Z_t, T = T_t, R = R_t, Q = Q_t, a1 = a1_t, P1 = P1_t)
formula <- substitute(ts ~ -1 + ssm, list(ssm = ssm))
model_fixed <- KFAS::SSModel(formula)

验证 #

运行上述修改后的代码，就能得到和model_good一致的状态空间模型。

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