r - 存储来自 CVlm 的均方

Question

我正在尝试在 R 中创建一个脚本，用于自动评估各种可能的线性模型的预测能力。为了评估模型的预测能力，我使用来自交叉验证的整体均方作为质量指标，为此我使用包 DAAG 中的函数 CVlm。我的问题是如何以自动方式检索由 CVlm 产生的整体均方值（无需观察 CVlm 的文本输出）？

例如来自http://maths-people.anu.edu.au/~johnm/r-book/3edn/scripts/reg1.R的以下代码

houseprices.lm <- lm(sale.price ~ area, data=houseprices) 
CVlm(houseprices, houseprices.lm, plotit=TRUE) 

有一个形式的输出

fold 1 测试集中的观察结果：...

折叠测试集中的 2 个观察值：...

总体女士 2023

如何访问/存储每次运行的 ms (2023) 值？

Answer 1

您必须将结果存储CVlm在变量中并访问ms属性：

houseprices.lm <- lm(sale.price ~ area, data=houseprices) 
cv <- CVlm(houseprices, houseprices.lm, plotit=TRUE) 
attr(cv, "ms")
# [1] 3934