So extrahieren Sie spezifische Restdaten aus einem linearen Modell in R

Question

Wie würde ich Extract Restdaten für ein bestimmtes Baseball-Team im folgenden linearen Modell extrahieren? Zum Beispiel, wie würde ich die Residuen für "CLE" extrahieren?

library(Lahman) 
library(dplyr) 
library(broom) 

# create baseball team data 
data(Teams) 
teams <- Teams 
teams <- teams %>% mutate(win_percentage = (W/(W + L)) * 100) 

# summarize baseball team salary by year 
salaries <- Salaries 
salaries <- salaries %>% 
    group_by(teamID, yearID, lgID) %>% 
    summarise(payroll_M = sum(as.numeric(salary))/10^6) %>% 
    ungroup() 

# add winning percentage to the salary table 
salaries <- teams %>% 
    select(yearID, teamID, win_percentage) %>% 
    right_join(salaries, by = c("yearID", "teamID")) 

# compute linear model of winning vs team salary 
model <- salaries %>% 
    group_by(yearID) %>% 
    do(fit = augment(lm(win_percentage ~ payroll_M, data = .))) 

# extract residuals for Cleveland ?????? 

Answer 1

Sie sind in der Nähe, aber müssen zwei Änderungen an der augment Linie.

1. Sie den resultierenden (augmented) Datenrahmen auf eine Säule fit genannt zu speichern. Versuchen Sie es stattdessen direkt an do (entfernen Sie die fit =).

2. Die Augment-Funktion muss die Spalte teamID als Teil der resultierenden Daten beibehalten, obwohl sie nicht im Modell enthalten ist. Beachten Sie, dass augment ein zweites Argument data für genau diesen Zweck verwendet (siehe help(augment.lm) für mehr).

Somit ist die neue Linie würde wie folgt aussehen:

do(augment(lm(win_percentage ~ payroll_M, data = .), data = .)) 

Der resultierende Datenrahmen eine Zeile pro ursprüngliche Beobachtung haben wird, und wird die teamID zusammen mit den Residuen und angepaßten Werte umfassen (was ermöglicht Sie für filtern).