Modelpassing visualiseren

Introductie tot regressie in R

Richie Cotton

Data Evangelist at DataCamp

Gewenste eigenschappen van residuen

Residuen zijn normaal verdeeld
Het gemiddelde van de residuen is nul

Brasem en baars opnieuw

Brasem: het “goede” model

mdl_bream <- lm(mass_g ~ length_cm, data = bream)

Het eerder getoonde spreidingsdiagram van brasem‑massa’s versus lengte, met een trendlijn.

Baars: het “slechte” model

mdl_perch <- lm(mass_g ~ length_cm, data = perch)

Het eerder getoonde spreidingsdiagram van baars‑massa’s versus lengte, met een trendlijn.

Residuen vs. voorspelde waarden

Brasem

Een spreidingsdiagram van brasem‑modelresiduen vs. voorspelde waarden, met een LOESS-trendlijn. De trendlijn blijft dicht bij de x‑as.

Baars

Een spreidingsdiagram van baars‑modelresiduen vs. voorspelde waarden, met een LOESS-trendlijn. De trendlijn vormt een V‑vorm.

Q-Q-plot

Brasem

Een Q-Q-plot van gestandaardiseerde brasem‑modelresiduen vs. theoretische kwantielen. De punten volgen de lijn vrijwel, op twee uitschieters na.

Baars

Een Q-Q-plot van gestandaardiseerde baars‑modelresiduen vs. theoretische kwantielen. De punten volgen de gelijkheidslijn niet goed, vooral rechts in de plot.

Schaal-locatie

Een spreidingsdiagram van bream-model wortel-gestandaardiseerde residuen vs. voorspelde waarden, met een LOESS-trendlijn. De trendlijn gaat licht omhoog naarmate de voorspelde waarden toenemen.

Een spreidingsdiagram van perch-model wortel-gestandaardiseerde residuen vs. voorspelde waarden, met een LOESS-trendlijn. De trendlijn gaat meerdere keren op en neer naarmate de voorspelde waarden toenemen.

autoplot()

library(ggplot2)
library(ggfortify)

autoplot(model_object, which = ???)

Waarden voor which

1 residuen vs. voorspelde waarden
2 Q-Q-plot
3 schaal-locatie

autoplot() met het perch‑model

autoplot(
  mdl_perch, 
  which = 1:3, 
  nrow = 3, 
  ncol = 1
)

De drie diagnostische plots voor het perch‑model in één figuur met drie panelen.

'Autoplots, uitrollen!' -Plotimus Prime

Introductie tot regressie in R