data(wvs)
library(MASS)
lmod <- lm(secpay ~ gini_disp + democrat + log(pop), data=wvs)
e1_m <- rlm(secpay ~ gini_disp + democrat + log(pop),
data=wvs, method="M")$residuals
e1_mm <- rlm(secpay ~ gini_disp + democrat + log(pop),
data=wvs, method="MM")$residuals
e1dat <- data.frame(OLS = lmod$residuals,
M = e1_m,
MM = e1_mm)
rrPlot(OLS ~ M + MM, data=e1dat)
Run the code above in your browser using DataLab