从 R 中的模型摘要中导出输出

Question

有没有办法从 R 中的摘要过程中导出输出（类似于 SAS 中的 ods 输出）？

我有这个时间序列：

summary(ets(ts(c(100, 110, 120, 200, 300, 100, 200, 100, 120, 300), start = c(2009, 1), frequency = 12)))

当您运行此程序时，您会获得各种重要信息，例如 MAPE 和 RMSE。我想有选择地将其中一些指标分配给变量，以便在以后的其他流程中使用。

MAPE.ets <- [some bit of code that takes the MAPE from the output and pops it right here]

R 中有没有办法做到这一点？

编辑：经过一些回答后，我意识到 stlf ets 的行为与 ets 不同，我也需要一个适用于 stlf ets 的解决方案。

summary(stlf(ts(c(100, 110, 120, 200, 300, 100, 200, 100, 120, 300), start = c(2009, 1), frequency = 4)))

Answer 1

只需将摘要分配给一个新对象。这将创建一个矩阵，您可以在其中访问您提到的指标：

summ <- summary(ets(ts(c(100, 110, 120, 200, 300, 100, 200, 100, 120, 300), start = c(2009, 1), frequency = 12)))

summ[,"MAPE"]
[1] 44.51514

编辑：

预测包中的预期方式似乎是使用accuracy 方法（ets 类的summary 方法实际上调用accuracy）；这对ets 和stlf 都有效：

res2 <- stlf(ts(c(100, 110, 120, 200, 300, 100, 200, 100, 120, 300), start = c(2009, 1), frequency = 4))
acc <- accuracy(res2)
acc[,"MAPE"]

Answer 2

摘要为您提供了一个向量。如果你想要 colnames 和 rownames，你可以像这样得到 MAPE。

out <-summary(ets(ts(c(100, 110, 120, 200, 300, 100, 200, 100, 120, 300), 
              start = c(2009, 1), frequency = 12)))
out[,5,drop=FALSE]

             MAPE
Training set 44.51514

我只想要号码，删除drop=FALSE

out[5]
[1] 44.51514

EDIT 使用stfl，最好使用accuracy 函数而不是summary。要获得 MAPE，您可以：

res <- stlf(ts(c(100, 110, 120, 200, 300, 100, 200, 100, 120, 300), 
            start = c(2009, 1), frequency = 4))
accuracy(res)[5]
[1] 44.5994