3

我正在尝试在 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) 值?

4

1 回答 1

2

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

houseprices.lm <- lm(sale.price ~ area, data=houseprices) 
cv <- CVlm(houseprices, houseprices.lm, plotit=TRUE) 
attr(cv, "ms")
# [1] 3934
于 2013-07-11T12:48:02.840 回答