Die Funktion rlm (MASS) erlaubt sowohl die M- als auch die MM-Schätzung für eine robuste Regression. Ich möchte den Smoother der MM-Robust-Regression in ggplot2 darstellen, denke jedoch, dass bei Auswahl von method = "rlm" in stat_smooth die automatisch ausgewählte Schätzmethode der M-Typ ist.

Gibt es eine Möglichkeit, die MM-Typ-Schätzungstechnik für die rlm-Funktion über ggplot2 auszuwählen?

Hier ist mein Code:

df <-  data.frame("x"=c(119,118,144,127,78.8,98.4,108,50,74,30.4,
50,72,99,155,113,144,102,131,105,127,120,85,153,40.6,133),
"y"=c(1.56,2.17,0.81,1.07,1.12,2.03,0.90,1.48,0.64,
0.91,0.85,0.41,0.55,2.18,1.49,1.56,0.82,0.93,0.84,1.84,
0.78,1.15,3.85,3.30,0.94))      

library(ggplot2)
library(MASS)      

ggplot(df,aes(x=x,y=y))+geom_point()+
stat_smooth(method="rlm",fullrange=TRUE)+xlim(0,160)

Ich habe die Ergebnisse mit der rlm-Zusammenfassung selbst überprüft und bin mir ziemlich sicher, dass ggplot2 die (Standard-?) M-Schätzung verwendet.

Wie kann ich die MM-Schätzung aus der rlm-Funktion verwenden?

rlm(formula, ...,method = "MM")

Vielen Dank im Voraus!