Как извлечь конкретные остаточные данные из линейной модели в R

Question

Как извлечь остаточные данные для конкретной бейсбольной команды в следующей линейной модели? Например, как я могу извлечь остатки для «CLE»?

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

Вы близки, но нужно два изменения в augment линии.

1. Вы экономите результирующее (дополненный) кадр данных столбца называется fit. Вместо этого попробуйте дать его непосредственно do (удалите fit =).

2. Функция увеличения должна содержать столбец teamID как часть полученных данных, даже если это не в модели. Обратите внимание, что augment принимает второй аргумент data для этой цели (см. help(augment.lm) для получения дополнительной информации).

Таким образом, новая строка будет выглядеть следующим образом:

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

Полученный кадр данных будет иметь одну строку для каждого первоначального наблюдения, и будет включать в себя teamID вместе с остатками и подобранными значениями (что позволяет вы должны фильтровать для CLE).