* in das Unterverzeichnis mit den Daten wechseln * den Datensatz oeffnen use golf.dta , clear * Grafik-Layout auf schwarz-weiss setzen set scheme s1mono * orthogonale Polynome berechnen orthpoly kilstand , generate(kilstandop1 kilstandop2 kilstandop3) degree(3) orthpoly alter , generate(alterop1 alterop2 alterop3) degree(3) * Modelldiagnose bestes Modell nach All subset selection regress preis alterop1 alterop2 kilstandop1 kilstandop2 matrix beta = e(b) * berechne den Alterseffekt und den Effekt des Kilometerstands gen alterseffekt = beta[1,1]*alterop1+beta[1,2]*alterop2 gen kilstandseffekt = beta[1,3]*kilstandop1+ beta[1,4]*kilstandop2 * Berechne vorhergesagte Werte, standardisierte Residuen und studentistierte Residuen predict preishat predict resstand , rstandard predict rstudent , rstudent * obere und untere Grenze der t-Verteilung gen t = invttail(172-6,0.01/344) gen t2 = -invttail(172-6,0.01/344) * Grafiken der studentisierten Residuen sort preishat graph twoway scatter rstudent t t2 preishat , yline(0) /* */ xsize(4.5) msymbol(oh i i) connect(. l l) clpattern(. l l) legend(off) /* */ xtitle("geschätzter Verkaufspreis") ytitle("studentisierte Residuen") /* */ subtitle("studentisierte Residuen vs geschätzter Verkaufspreis") ylab(-5(2)5) sort kilstand graph twoway scatter rstudent t t2 kilstand , yline(0) legend(off) /* */ xsize(4.5) msymbol(oh i i) connect(. l l) clpattern(. l l) /* */ connect(i l l i) clpattern(i l l i) ytitle("studentisierte Residuen") /* */ subtitle("studentisierte Residuen vs Kilometerstand") ylab(-5(2)5) sort alter graph twoway scatter rstudent t t2 alter , yline(0) xlabel(65(10)140) /* */ xsize(4.5) msymbol(oh i i) connect(. l l) clpattern(. l l) /* */ xtitle("Alter in Monaten") ytitle("studentisierte Residuen") /* */ subtitle("studentisierte Residuen vs Alter") legend(off) ylab(-5(2)5) sort tuev graph twoway scatter rstudent t t2 tuev , yline(0) /* */ xsize(4.5) msymbol(oh i i) connect(. l l) clpattern(. l l) /* */ xtitle("Monate bis zum nächsten Tüv") ytitle("studentisierte Residuen") ylab(-5(2)5) /* */ subtitle("studentisierte Residuen vs Monate bis Tüv") xlabel(0(3)24) legend(off) * QQ-Plot der standardisierten Residuen qnorm resstand , xsize(4.5) xtitle("Quantile der Normalverteilung") ytitle("Quantile der Residuen") xlab(-2(1)2)