R语言lm线性回归中如何剔除极值?

R语言中做lm线性拟合如何剔除极值？请教各位大神，car包中的influencePlot函数可以显示极值，但不知道如何处理，有大神指导嘛？比如用什么包？或者函数？
回归方程如下：
fit<-lm(zd~xd+rd+mdy+sd+wd,data=na.omit(train_data),model=TRUE)