lm 使用插入符号火车的结果

Question

我正在使用 caret R 包作为一个非常方便的建模包装器。尽管这是一种奇怪的用法，但在使用模型类型 =“lm”和“none”的交叉验证时，我在从模型中提取结果时遇到了一些麻烦。见下面的例子：

library(caret)
## Make data
Xs <- matrix(rnorm(300*20), nrow = 300, ncol = 20)
Yvec <- rnorm(300)
## Make traincontrol, cv of "none"
tcontrol <- trainControl(method = "none")
## Fit lm model using train
fit <- train(x= Xs, y = Yvec, method = "lm", metric = "RMSE", trControl = tcontrol)

fit$results
[1] RMSE      Rsquared  parameter
<0 rows> (or 0-length row.names)

知道为什么 fit$results 是空的吗？对于所有其他模型和 cv 类型，这似乎都有效。例如使用 2 倍 CV：

tcontrol2 <- trainControl(method = "cv", number = 2)
fit2      <- train(x= Xs, y = Yvec, method = "lm", metric = "RMSE", trControl = tcontrol2)
fit2$results
  parameter     RMSE     Rsquared      RMSESD   RsquaredSD
1      none 1.023666 0.0008921194 0.006499062 0.0003463973

我明白这是一个人为的例子，但这个模型和 cv 方法只是我正在测试的一个更大数字的组合（因此插入符号包装器是理想的）。

Answer 1

其实我想我可能有办法解决。当没有重采样时，我相信你仍然可以从 train 的输出中得到拟合值。如果为真，这将意味着单个模型（即没有 cross-validation/resampling）的 RMSE/Rsq 只是：

> caret::RMSE(pred = fit$finalModel$fitted.values, obs = Yvec)
[1] 0.9348365
> caret::R2(pred = fit$finalModel$fitted.values, obs = Yvec)
[1] 0.04692012

而这个 R^2 对应于如果你只是手动 lm 你会得到什么：

> dat <- cbind(Yvec, Xs) %>% as.data.frame()
> lm(Yvec ~., data=dat) %>% summary

Call:
lm(formula = Yvec ~ ., data = dat)

...

...Multiple R-squared:  0.04692...

我认为这已经结案，但如果其他人觉得它有帮助，我会保留它，我仍然很想听听是否有人可以 confirm/refute 这个。

lm 使用插入符号火车的结果

results from lm using caret train

r

lm

r-caret