本文共 522 字，大约阅读时间需要 1 分钟。

书里代码少了一行，这里修改过了

library(ggplot2)library(tidyverse)mod <- loess(hwy~displ,data=mpg)grid <- data.frame(displ=seq(min(mpg$displ),max(mpg$displ),length=50))grid$hwy <- predict(mod,newdata = grid)std_resid <- resid(mod)/mod$smpg <- mutate(mpg,std_resid)outlier <- filter(mpg,abs(std_resid)>2)ggplot(mpg,aes(displ,hwy))+geom_point()+geom_line(data=grid,colour="blue",size=1.5)+geom_text(data=outlier,aes(label=model))ggplot(mpg,aes(displ,hwy))+geom_point(data = mpg)+geom_line(data=grid)+geom_text(data=outlier,aes(label=model))