R 預測擬合的殘差？

Question

我知道當您生成線性 model 時，您可以像這樣從擬合中提取殘差：

# Model
model <- lm(y ~ x, data = coolstuff)

# Residuals
myresids <- model$residuals

我進一步了解您可以使用此 model 來預測基於第二個數據集（例如，驗證數據集）的值，如下所示：

mypreds <- predict(model, newdata = coolvalid)

我迷路的地方就是我可以從預測中找到殘差的地方。 predict不會生成數據框或小標題 - 只是一個命名的數字列表。

我在哪里可以找到預測的殘差？

Answer 1

正如 DanY 在評論中指出的那樣，殘差只是觀察值-預測值。 有關內置數據的簡單示例，請參見下文。

# sample data
set.seed(1)
split_indicies <- sample(nrow(mtcars), nrow(mtcars)/2)
train <- mtcars[split_indicies,]
test <- mtcars[-split_indicies,]

# model
model <- lm(mpg ~ disp, data = train)

# residuals of prediction are actual - predicted
test$mpg - predict(model, data = test)
#> Pontiac Firebird   Hornet 4 Drive       Duster 360        Mazda RX4 
#>        9.8318659       -0.1503095        3.4749151        1.4901610 
#>    Mazda RX4 Wag      AMC Javelin        Merc 280C      Merc 450SLC 
#>       -0.1098390       -0.5448161       -5.2950183       -7.7129664 
#>         Fiat 128      Honda Civic   Ford Pantera L    Toyota Corona 
#>      -15.8775915      -11.7018628        0.5021011      -11.2626474 
#>         Merc 280       Volvo 142E   Toyota Corolla     Ferrari Dino 
#>        4.7049817        1.4746340        3.8075878       -8.5311955

^{由reprex package (v2.0.1) 創建於 2022-02-07}