3.7 Drill: regresjonsoppgaven fra R-utskrift (sjanger H–K)
Den tunge storoppgaven trent ledd for ledd: fra R-utskrift til koeffisientinferens, R²/justert R², konfundering, diagnostikk og prediksjonsintervall — enkel og multippel regresjon.
Kjeden ruller gjennom fire faste sjangre (oppgavetyper):
- Sjanger H — enkel regresjon: koeffisienter, t-test og KI for stigningstallet (kap. 3.1–3.2).
- Sjanger I — multippel regresjon: matriseform, justert , konfundering (kap. 3.3).
- Sjanger J — diagnostikk: de tre residualplottene (kap. 3.4).
- Sjanger K — prediksjon: KI for forventet respons vs. prediksjonsintervall (kap. 3.5).
Formelsamlingen dekker KI-/PI-formlene og -definisjonene — tren valg og oppslag. Det som må regnes aktivt: tolkning i kontekst, riktige frihetsgrader, konfunderingsargumentet og mønster-lesing av plott. Prioritet: må sitte perfekt.
Tidsbruk: ~80 min. Del gjerne i to økter: oppskriften + eksamenscasen + oppgave 1–5 i én, oppgave 6–10 i neste. Anslaget er lesetid; regner du for hånd, gang med ca. 1,5.
Alle nødvendige kvantiler legges — som på eksamen — rett i oppgavetekstene. R-utskriftene er rene teksttabeller: du skal lese dem, aldri skrive kode.
Alt i dette kapitlet er kombinasjoner av disse fire — pluss ord: tolkning i kontekst, konfunderingsargument og mønster-lesing av plott.
Løsningsoppskrift: syv trinn fra utskrift til full besvarelse
Hver regresjonsoppgave — uansett hvilke deloppgaver den stiller — er et utvalg av de samme syv trinnene. Kjør dem i rekkefølge:
Trinn 1 — Definer modellen og antakelsene. Skriv med uavhengige . Navngi variablene i kontekst.
Trinn 2 — Les/utled koeffisientene og tolk i kontekst. Fra utskriften (Estimate) eller via minste kvadrater. Tolk ALLTID med enhet: «per enhet endres forventet med …» — i multippel modell: «… med de øvrige variablene holdt konstant».
Trinn 3 — Inferens for en koeffisient. t-test og KI — med riktig frihetsgrad: (enkel), (multippel). Ved samlet spørsmål: F-testen ( alle helninger ).
Trinn 4 — vs. justert . Tolk som forklart andel. Ved modellvalg brukes ALLTID justert — øker automatisk for hver variabel som legges til.
Trinn 5 — Konfundering. Endrer en koeffisient seg fra enkel til multippel modell? Forklar: den nye variabelen er korrelert med både responsen og den gamle variabelen; i multippel modell måles effekten med den nye holdt konstant.
Trinn 6 — Diagnostikk i mønstre. Residual mot tilpasset (kurvatur → ikke-linearitet), scale-location (trend → heteroskedastisitet), QQ-plott (avvik fra linje → ikke-normalitet). Aldri bare «ser greit ut» — nevn mønsteret du ser eller ikke ser.
Trinn 7 — Prediksjon. KI for forventet respons vs. prediksjonsintervall for ny observasjon — PI har «» under rota og er alltid bredere. Kommenter alltid om ligger i dataområdet (ekstrapolering).
Kjede-regelen: står du fast på ett ledd, bruk oppgitt mellomresultat og gå videre — et blankt deltrinn sperrer aldri resten av kjeden.
| Du trenger | Sted i utskriften | Regel |
|---|---|---|
| Estimate-kolonnen | tolk i kontekst, med enhet | |
| Std. Error-kolonnen | til t-test og KI | |
| t-test av | t value / Pr(>|t|) | forkast når P-verdi |
| test av | må regnes selv | |
| og | Residual standard error-linja | eller |
| forklart andel | Multiple R-squared | tolkning, ikke modellvalg |
| modellvalg | Adjusted R-squared | størst justert vinner |
| samlet test | F-statistic-linja |
I enkel regresjon er F-testen og t-testen for stigningstallet samme test: .
For bilmodeller er drivstofforbruket (liter per 100 km) regressert på vekt (tonn), og deretter på vekt og motorvolum (liter). Vektene i utvalget ligger mellom 0,9 og 2,1 tonn, med og .
Modell 1 (enkel):
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.1000 0.3424 6.134 1.28e-06
vekt 3.2000 0.2357 13.576 7.68e-14
---
Residual standard error: 0.5 on 28 degrees of freedom
Multiple R-squared: 0.8681, Adjusted R-squared: 0.8634
F-statistic: 184.3 on 1 and 28 DF, p-value: 7.68e-14Modell 2 (multippel):
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.7500 0.5200 3.365 0.00231
vekt 2.5300 0.3400 7.441 5.27e-08
volum 0.9500 0.4100 2.317 0.02832
---
Residual standard error: 0.463 on 27 degrees of freedom
Multiple R-squared: 0.891, Adjusted R-squared: 0.8829
F-statistic: 110.4 on 2 and 27 DF, p-value: 1.01e-13a) Tolk vekt-koeffisienten i modell 1 i kontekst.
b) Test i modell 1 om vekt har effekt, og gi et 95 %-KI for stigningstallet. Bruk .
c) Hvilken modell foretrekker du, og hvorfor?
d) Vekt-koeffisienten falt fra 3,20 til 2,53 da volum kom inn. Forklar.
e) Residualplottene for modell 2 viser: ingen kurvatur i residual-mot-tilpasset, flat scale-location, QQ-punkter nær linja. Vurder antakelsene.
f) Bruk modell 1 til et 95 % prediksjonsintervall for forbruket til en NY bil på 1,6 tonn, og sammenlign med KI-et for forventet forbruk ved samme vekt.
b) Trinn 3. mot . Fra utskriften: og Pr(>|t|) : forkast — vekt har klar effekt. KI: . Sensor-margnotat: df er — hentes rett fra «on 28 degrees of freedom»-linja; her er den klassiske feilen.
c) Trinn 4. Modellvalg gjøres med justert : (modell 2) mot (modell 1) — modell 2 foretrekkes. I tillegg er volum-koeffisienten signifikant (). Sensor-margnotat: å sammenligne Multiple R-squared (0,891 mot 0,868) gir ikke uttelling — den øker uansett når en variabel legges til.
d) Trinn 5 — konfundering. Tyngre biler har som regel større motor, så vekt og volum er positivt korrelert, og volum påvirker forbruket. I modell 1 fanger vekt-koeffisienten derfor også noe av volum-effekten. I modell 2 måles vekt-effekten med volum holdt konstant — den rendyrkede vekt-effekten er mindre (2,53).
e) Trinn 6. Ingen kurvatur → lineariteten står; flat scale-location → konstant varians (homoskedastisitet) er forenlig med data; QQ nær linja → normalantakelsen er rimelig. Antakelsene virker oppfylt — konklusjonen begrunnes i mønstrene, ikke med et rent «ser greit ut». Sensor-margnotat: navngi hvilket plott som sjekker hva.
f) Trinn 7. . Med , , , :
KI for forventet forbruk: .
PI for ny bil: .
PI-et er mye bredere: den nye bilens eget støyledd kommer i tillegg («» under rota). ligger trygt i dataområdet [0,9, 2,1], så ingen ekstrapoleringsfare. Sensor-margnotat: KI/PI-forvekslingen er den vanligste feilen i sjanger K — pek eksplisitt på «».
Drill: ti a–f-kjeder på eksamensnivå
Oppgavene under ruller gjennom sjanger H, I, J og K — enkle og multiple modeller om hverandre, alle med R-utskrift. Kjør syvtrinnsoppskriften hver gang. To oppgaver er merket «(krevende)», og helt sist ligger en kald oppgave uten hint, der fasiten bare er en momentliste — som generalprøve.
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.4900 0.52200 6.685 8.06e-07
areal 0.0820 0.00919 8.920 6.30e-09
---
Residual standard error: 0.65 on 23 degrees of freedom
Multiple R-squared: 0.7758, Adjusted R-squared: 0.766
F-statistic: 79.57 on 1 and 23 DF, p-value: 6.30e-09a) Tolk stigningstallet i kontekst.
b) Verifiser t value for areal.
c) Er arealeffekten signifikant på 5 %-nivå?
d) Gi et 95 %-KI for . Bruk .
e) Angi og frihetsgradene, og forklar hvor kan leses ut.
f) Tolk Multiple R-squared.
a) Finn minste kvadraters estimater og .
b) Estimer med riktige frihetsgrader.
c) Test mot på 5 %-nivå. Bruk og at .
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 4.0000 1.61000 2.485 0.0230
tid 1.8000 0.15179 11.859 6.11e-10
---
Residual standard error: 2.4 on 18 degrees of freedom
Multiple R-squared: 0.8865, Adjusted R-squared: 0.8802Bruk .
a) Angi et 95 %-KI for forventet strekkfasthet ved døgn.
b) Angi et 95 % prediksjonsintervall for en NY prøve ved , og forklar hvorfor det er bredere.
c) En kollega vil predikere ved døgn. Regn ut prediksjonsintervallet og forklar hvorfor det ikke bør brukes.
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.2000 0.3500 3.429 0.00211
areal 0.0280 0.0060 4.667 8.84e-05
alder -0.0150 0.0070 -2.143 0.04205
---
Residual standard error: 0.31 on 25 degrees of freedom
Multiple R-squared: 0.712, Adjusted R-squared: 0.689
F-statistic: 30.90 on 2 and 25 DF, p-value: 1.75e-07a) Hvor mange boliger er med, og hvordan ser du det?
b) Tolk alder-koeffisienten i kontekst.
c) Gi et 95 %-KI for alder-koeffisienten. Bruk .
d) Hvilken hypotese tester F-linja, og hva konkluderer du?
Modell A (enkel): utdanning 42.0 (Std. Error 6.0)
Adjusted R-squared: 0.41
Modell B (multippel): utdanning 28.0 (Std. Error 7.0)
erfaring 15.0 (Std. Error 4.0)
Adjusted R-squared: 0.55a) Utdanningskoeffisienten falt fra 42 til 28. Gi konfunderingsforklaringen.
b) Hvilken modell foretrekkes, og med hvilket kriterium?
c) Under hvilken betingelse ville koeffisienten IKKE endret seg fra A til B?
- Residual mot tilpasset: tydelig omvendt U-form (parabel) — positive residualer i midten, negative i begge ender.
- Scale-location: omtrent flat trend.
- QQ-plott: punktene følger linja godt.
a) Hva sjekker hvert av de tre plottene?
b) Hvilken antakelse er brutt her, og hva betyr mønsteret faglig?
c) Foreslå en konkret modellutvidelse, og forklar hvorfor konklusjoner fra den lineære modellen er upålitelige.
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.1000 1.9000 2.684 0.0139
kanal1 0.8400 0.6400 1.312 0.2035
kanal2 0.7100 0.6800 1.044 0.3083
---
Residual standard error: 3.4 on 21 degrees of freedom
Multiple R-squared: 0.706, Adjusted R-squared: 0.678
F-statistic: 25.21 on 2 and 21 DF, p-value: 2.62e-06a) Ingen av kanal-koeffisientene er signifikante, men F-testen er sterkt signifikant. Forklar hvordan dette henger sammen.
b) Betyr t-resultatene at annonsering er uten verdi? Formuler den korrekte, betingede tolkningen.
c) Hva ville du gjort videre med modellen?
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.3100 0.4200 0.738 0.4662
referanse 0.9140 0.0380 24.053 < 2e-16
---
Residual standard error: 1.02 on 30 degrees of freedoma) Forklar hvorfor Pr(>|t|)-kolonnen IKKE besvarer kalibreringsspørsmålet.
b) Test mot på 5 %-nivå. Bruk og .
c) Kryssjekk med et 95 %-KI for .
Modell A (bare temperatur):
Coefficients:
Estimate Std. Error t value Pr(>|t|)
temp 1.4200 0.1800 7.889 4.27e-09
---
Residual standard error: 2.32 on 33 degrees of freedom
Multiple R-squared: 0.658, Adjusted R-squared: 0.6476Modell B (temperatur + nedbør):
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 8.2000 2.1000 3.905 0.00046
temp 1.2100 0.1600 7.562 1.30e-08
nedbor -0.6500 0.2100 -3.095 0.00407
---
Residual standard error: 2.05 on 32 degrees of freedom
Multiple R-squared: 0.741, Adjusted R-squared: 0.7248
F-statistic: 45.78 on 2 and 32 DF, p-value: 4.10e-10a) Tolk temp-koeffisienten i modell B i kontekst.
b) Test på 5 %-nivå om nedbør har effekt gitt temperatur, med alle ledd i konklusjonen.
c) Gi et 95 %-KI for nedbør-koeffisienten. Bruk .
d) Velg modell med begrunnelse.
e) Temp-koeffisienten endret seg fra 1,42 til 1,21. Forklar mekanismen.
f) En festivalarrangør vil bruke modell B ved 35 °C. Kommenter.
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 13.2000 1.00700 13.109 7.19e-09
gjodsel 0.1200 0.02373 5.056 0.000220
---
Residual standard error: 1.3 on 13 degrees of freedom
Multiple R-squared: 0.6629, Adjusted R-squared: 0.6369
F-statistic: 25.56 on 1 and 13 DF, p-value: 0.000220Bruk .
a) Tolk stigningstallet i kontekst.
b) Test på 5 %-nivå om gjødsel har effekt, og verifiser t value.
c) Gi et 95 %-KI for .
d) Angi , frihetsgradene og sammenhengen mellom F og t her.
e) Gi et 95 % prediksjonsintervall for en NY rute med , og kommenter dataområdet.
2. i stedet for justert til modellvalg. øker alltid når du legger til en prediktor — bare justert straffer unødvendige ledd.
3. Konfidens- og prediksjonsintervall forvekslet. KI gjelder forventet respons ; PI gjelder én ny observasjon og har det ekstra -leddet — PI er alltid bredere.
4. Koeffisient tolket som «hele effekten». I multippel regresjon er effekten av når de andre prediktorene holdes faste — utelater du det, tolker du feil (konfundering).
5. Lese R-utskriften feil og «ser greit ut» i diagnostikken. Vit hvilken kolonne som er Estimate/Std. Error/t/P-verdi, og tolk residualplottene i mønstre (vifteform, kurvatur, tunge haler) — aldri bare «ser greit ut».
Begrepsbank
Den faste kjeden i enhver regresjonsoppgave: (1) modell + antakelser; (2) koeffisienter + tolkning i kontekst; (3) t-test/KI med riktig df (evt. F-test); (4) vs. justert ; (5) konfundering; (6) diagnostikk i mønstre; (7) prediksjon med KI/PI-skillet og ekstrapoleringssjekk.
Startpunktet for full uttelling: med uavhengige . Antakelsene (linearitet, konstant varians, normalitet, uavhengighet) er nettopp det diagnostikken i trinn 6 etterprøver.
Sensorkravet for koeffisienter: tallfest MED enhet, og i multippel modell alltid med forbeholdet «med de øvrige variablene holdt konstant». En naken tallverdi uten kontekst gir ikke full uttelling.
Observatoren under . Utskriftens t value og Pr(>|t|) gjelder nullverdien 0; mot enhver annen nullverdi må regnes selv.
Konfidensintervallet bygges rett fra kolonnene: Estimate Std. Error, med df fra Residual standard error-linja. Dualitet: nullverdier utenfor intervallet forkastes av den tosidige testen.
F-statistic-linja tester — at ingen av prediktorene har effekt. I enkel regresjon er , så F- og t-testen er samme test der.
df i enkel regresjon (to parametere estimert), i multippel med prediktorer. Leses direkte av «on … degrees of freedom»; feil df straffes hardt.
At en koeffisient endrer seg når en korrelert variabel legges til: uten den nye variabelen fanger den gamle koeffisienten også noe av dens effekt. Med ukorrelerte (ortogonale) prediktorer skjer ingen endring.
Residual mot tilpasset → linearitet (kurvatur = feil form); scale-location → konstant varians (trend = heteroskedastisitet); QQ-plott → normalitet (avvik fra linja). Tolkes alltid i MØNSTRE, aldri «ser greit ut».
Intervallet for : — usikkerheten i selve linja ved .
Intervallet for en NY observasjon ved : samme formel med . «»-leddet er den nye observasjonens eget støyledd — derfor alltid bredere enn KI-et.
Før enhver prediksjon: ligger i dataområdet? Utenfor er lineariteten udokumentert og intervallene eser ut med — prediksjonen skal frarådes eller rapporteres med kraftig forbehold.
En t-test i multippel modell er betinget på de øvrige variablene: «ikke signifikant» betyr «ingen påvisbar tilleggseffekt», ikke «uten verdi». Ved kollineære prediktorer kan alle t-tester være svake selv om F-testen er sterk.
Eksamens a–f-kjeder er bygget slik at ett blankt deltrinn ikke sperrer resten: bruk det oppgitte mellomresultatet og regn videre. Å stoppe helt ved første hinder er den dyreste strategifeilen i storoppgaven.
Dette kapitlet er skrevet av Anthropics toppmodeller (Claude Opus og Claude Fable) og er foreløpig ikke manuelt gjennomgått — kvalitetskontrollen gjøres av uavhengige KI-agenter, og innmeldte feil rettes fortløpende. Funnet en feil? Meld fra, så retter vi den. Les mer om hvordan innholdet lages.
Skolesaga er en uavhengig læringsressurs og er ikke tilknyttet eller godkjent av Universitetet i Oslo. Dette er ikke offisielt studiemateriell. Les mer.