存储均方由CVlm产生

时间:2013-07-11 12:34:28

标签: r

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

的输出格式为

弃1 测试集中的观察结果:......

折叠2 测试集中的观察结果:......

总体来说 2023

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

1 个答案:

答案 0 :(得分:2)

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

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