我正在尝试在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)的值?
答案 0 :(得分:2)
您必须将CVlm
的结果存储在变量中并访问ms
属性:
houseprices.lm <- lm(sale.price ~ area, data=houseprices)
cv <- CVlm(houseprices, houseprices.lm, plotit=TRUE)
attr(cv, "ms")
# [1] 3934