Robuste Regression in R

Arndt Regorz, Dipl. Kfm. & MSc. Psychologie, 19.10.2021

Dieses ist eine Begleitseite zum Video-Tutorial über Robuste Regression in R.

(Hinweis: Mit Anklicken des Videos wird ein Angebot des Anbieters YouTube genutzt.)

Code aus dem Video-Tutorial:

Hier ist der gesamte Code aus dem Video-Tutorial zur robusten Regression:

library(robustbase)

data(package="robustbase")

# Datensatz aus robustbase
head(CrohnD)

# Gewöhnliche multiple Regression
fit.reg <- lm(nrAdvE ~ BMI + height, data= CrohnD)
summary(fit.reg)

# Robuste Regression
fit.rob <- lmrob(nrAdvE ~ BMI + height, data= CrohnD)
summary(fit.rob)

# Streudiagramm
plot(CrohnD$BMI, CrohnD$nrAdvE)
plot(CrohnD$height, CrohnD$nrAdvE)

#Ausreißerdiagnose Outlier&Leverage
library(olsrr)
ols_plot_resid_lev(fit.reg)

Weitere Info zum MM-Schätzer:

Der in dem robustbase-Package für die robuste Regression verwendete MM-Schätzer basiert auf folgendem Grundlagenpaper:

Yohai, V. J. (1987). High breakdown-point and high efficiency robust estimates for regression. The Annals of Statistics, 15(2), 642-656. https://www.jstor.org/stable/2241331