Tilbake
4.2

4.2 Drill: eksponensialens kjikvadrat-apparat

Den hyppigste vanskelige gjengangeren i hele arkivet: vis at $2X/\beta \sim \chi^2_2$, løft til $2\sum X_i/\beta \sim \chi^2_{2n}$ — og gjør det på autopilot.

80 min
10 oppgaver
Drilleksponensialens kjikvadrat-apparat
Din fremgang i kapitlet
0 / 10 oppgaver
Forkunnskaper: Dette drillkapitlet forutsetter de tre teknikkene fra kap. 4.1: fordelingsfunksjonsteknikken (CDF-metoden), MGF-produktregelen med entydighetsargumentet, og — som råstoff — eksponensial- og kjikvadratfordelingen fra kap. 2.3 (husk at 2X/βχ222X/\beta \sim \chi^2_2 ble forhåndsvarslet der; her utleder vi det).

Du bør kunne at χ22\chi^2_2 er en eksponensialfordeling med forventning 2 (så Fχ22(y)=1ey/2F_{\chi^2_2}(y) = 1 - e^{-y/2}), og at MGF-en til eksponensial(β)(\beta) er (1βt)1(1 - \beta t)^{-1} mens MGF-en til χν2\chi^2_\nu er (12t)ν/2(1 - 2t)^{-\nu/2}. Ingen ny teori — bare tempo og presisjon på den samme utledningen, på mange varianter.

Løsningsoppskriften (fire trinn)

Uansett hvordan oppgaven er kledd, følger kjikvadrat-apparatet samme fire trinn:

1. Vis 2X/βχ222X/\beta \sim \chi^2_2 med CDF-metoden. For Xeksp(β)X \sim \text{eksp}(\beta), sett Y=2X/βY = 2X/\beta. Da er, for y>0y > 0,
FY(y)=P ⁣(2Xβy)=P ⁣(Xβy2)=1e(βy/2)/β=1ey/2.F_Y(y) = P\!\left(\frac{2X}{\beta} \le y\right) = P\!\left(X \le \frac{\beta y}{2}\right) = 1 - e^{-(\beta y/2)/\beta} = 1 - e^{-y/2}.
Dette er fordelingsfunksjonen til en eksponensial med forventning 2, altså χ22\chi^2_2. Pass på ulikhetsretningeng(x)=2x/βg(x) = 2x/\beta er voksende, så ulikheten snur ikke.

2. Løft til summen med MGF-produktregelen. Hver 2Xi/βχ222X_i/\beta \sim \chi^2_2 har MGF (12t)1(1 - 2t)^{-1}. For uavhengige X1,,XnX_1, \ldots, X_n:
M2Xi/β(t)=i=1n(12t)1=(12t)n,M_{2\sum X_i/\beta}(t) = \prod_{i=1}^{n}(1 - 2t)^{-1} = (1 - 2t)^{-n},
som er MGF-en til χ2n2\chi^2_{2n}. Ved entydighet: 2i=1nXi/βχ2n22\sum_{i=1}^{n} X_i/\beta \sim \chi^2_{2n}. Frihetsgradene er 2n2n, ikke nn — det er faktoren 2 fra hver observasjon.

3. Generaliser via potens-/log-løft. Er dataene Weibull med F(x)=1e(x/θ)kF(x) = 1 - e^{-(x/\theta)^k}, så er (X/θ)keksp(1)(X/\theta)^k \sim \text{eksp}(1) — potens-transformér først, så er du tilbake til eksponensialtilfellet: 2(Xi/θ)kχ2n22\sum (X_i/\theta)^k \sim \chi^2_{2n}. Samme grep for Gumbel og Pareto med et oppgitt hint om transformasjonen.

4. Pek fremover. Dette er nettopp pivotalen som gir et eksakt konfidensintervall for β\beta i kap. 6.2 og en eksakt test i Del 7-kjeden — fordi χ2n2\chi^2_{2n} ikke inneholder den ukjente β\beta. Utledningen du driller her, er selve broen inn i inferensen.

✏️Gjennomregnet eksamenscase (med sensornotater)

Levetidene X1,,X12X_1, \ldots, X_{12} (i tusen timer) til et kretskort er uavhengige og eksponensialfordelte med forventning β\beta.

a) Vis ved fordelingsfunksjonsteknikken at 2X1/β2X_1/\beta er kjikvadratfordelt med 2 frihetsgrader.

b) Angi med begrunnelse fordelingen til 2i=112Xi/β2\sum_{i=1}^{12} X_i/\beta.

c) Forklar med én setning hvorfor denne størrelsen er en pivotal.

a) CDF-metoden (trinn 1). La Y=2X1/βY = 2X_1/\beta. Funksjonen g(x)=2x/βg(x) = 2x/\beta er voksende, så ulikheten snur ikke. For y>0y > 0:
FY(y)=P ⁣(2X1βy)=P ⁣(X1βy2)=1e(βy/2)/β=1ey/2.F_Y(y) = P\!\left(\frac{2X_1}{\beta} \le y\right) = P\!\left(X_1 \le \frac{\beta y}{2}\right) = 1 - e^{-(\beta y/2)/\beta} = 1 - e^{-y/2}.
Dette er fordelingsfunksjonen til en eksponensial med forventning 2, som er χ22\chi^2_2. Altså 2X1/βχ222X_1/\beta \sim \chi^2_2.

> Sensornotat: poenget sitter i overgangen P(2X1/βy)=P(X1βy/2)P(2X_1/\beta \le y) = P(X_1 \le \beta y/2) — vis at gg er voksende, så leseren ser at ulikheten ikke snur. Legg merke til at β\beta kanselleres i eksponenten; det er derfor størrelsen blir parameterfri.

b) MGF-produktregelen (trinn 2). Hver 2Xi/βχ222X_i/\beta \sim \chi^2_2 med MGF (12t)1(1 - 2t)^{-1}, og variablene er uavhengige, så
M2Xi/β(t)=i=112(12t)1=(12t)12.M_{2\sum X_i/\beta}(t) = \prod_{i=1}^{12}(1 - 2t)^{-1} = (1 - 2t)^{-12}.
Dette er MGF-en til χ242\chi^2_{24}. Ved entydighet (MGF-en bestemmer fordelingen): 2i=112Xi/βχ2422\sum_{i=1}^{12} X_i/\beta \sim \chi^2_{24}.

> Sensornotat: frihetsgraden er 2n=212=242n = 2 \cdot 12 = 24, ikke n=12n = 12. Dette er den vanligste enkeltfeilen. Og entydighetsargumentet må stå eksplisitt — ellers har du bare regnet ut en funksjon.

c) Størrelsen 2Xi/β2\sum X_i/\beta er en pivotal fordi den er bygget av data og β\beta, men fordelingen χ242\chi^2_{24} ikke inneholder β\beta (eller andre ukjente) — nettopp det som lar oss snu den til et konfidensintervall for β\beta i kap. 6.2.

— naturlig pausepunkt — nå driller du selv, på nye varianter.

📝Oppgave 1

(Innstegsoppgave — kjernesteget.) XX er eksponensialfordelt med forventning β\beta. La Y=2X/βY = 2X/\beta.

a) Er g(x)=2x/βg(x) = 2x/\beta voksende eller avtakende? Hva betyr det for ulikheten?

b) Regn ut FY(y)=P(2X/βy)F_Y(y) = P(2X/\beta \le y) og vis at den er 1ey/21 - e^{-y/2}.

c) Hvilken fordeling er dette?

📝Oppgave 2

Det er kjent at 2Xi/βχ222X_i/\beta \sim \chi^2_2 for hver av n=8n = 8 uavhengige eksponensiale.

a) Hvilken fordeling har 2i=18Xi/β2\sum_{i=1}^{8} X_i/\beta?

b) Hva blir frihetsgraden hvis n=20n = 20 i stedet?

📝Oppgave 3

Ventetidene X1,,X6X_1, \ldots, X_6 mellom feil på en fiberkabel er uavhengige og eksponensialfordelte med forventning β\beta.

a) Vis kort at 2X1/βχ222X_1/\beta \sim \chi^2_2.

b) Bruk MGF-produktregelen til å identifisere fordelingen til 2i=16Xi/β2\sum_{i=1}^{6} X_i/\beta, med eksplisitt entydighetsargument.

📝Oppgave 4
XX er eksponensialfordelt med forventning β\beta, og 2X/βχ222X/\beta \sim \chi^2_2. Bruk at χ22\chi^2_2 har fordelingsfunksjon F(y)=1ey/2F(y) = 1 - e^{-y/2}.

a) Finn P(2X/β>3)P(2X/\beta > 3).

b) Finn P(2X/β>5)P(2X/\beta > 5).

c) (Sanity-sjekk.) Uttrykk hendelsen {2X/β>3}\{2X/\beta > 3\} som en hendelse om XX, og bekreft at svaret i a) også kan fås direkte fra P(X>3β/2)P(X > 3\beta/2).

📝Oppgave 5

Vindstyrken XX (i m/s) på en offshore-plattform er Weibull-fordelt med F(x)=1e(x/θ)kF(x) = 1 - e^{-(x/\theta)^k}, med kjent form kk og skala θ\theta. La W=(X/θ)kW = (X/\theta)^k.

a) Vis ved CDF-metoden at Weksp(1)W \sim \text{eksp}(1).

b) Hvilken fordeling har da 2W2W?

📝Oppgave 6
nn vindturbin-lagre har Weibull-fordelte levetider med kjent form kk og skala θ\theta: F(x)=1e(x/θ)kF(x) = 1 - e^{-(x/\theta)^k}. Du har allerede vist at (Xi/θ)keksp(1)(X_i/\theta)^k \sim \text{eksp}(1).

a) Angi fordelingen til 2i=1n(Xi/θ)k2\sum_{i=1}^{n}(X_i/\theta)^k, med begrunnelse.

b) Med n=10n = 10: hva blir frihetsgraden?

📝Oppgave 7

Bruddstyrkene X1,,X15X_1, \ldots, X_{15} (i kN) til offshore-bolter er uavhengige og eksponensialfordelte med forventning β\beta.

a) Utled fra bunnen av at 2X1/βχ222X_1/\beta \sim \chi^2_2 (vis CDF-steget).

b) Angi med fullt MGF-argument fordelingen til V=2i=115Xi/βV = 2\sum_{i=1}^{15} X_i/\beta.

c) Uttrykk Xˉ\bar{X} ved VV, og bruk det til å argumentere for at 2nXˉ/βχ2n22n\bar{X}/\beta \sim \chi^2_{2n}.

📝Oppgave 8

En oppskriftsfordeling: nedbørstoppene XX (over et referansenivå x0=1x_0 = 1) følger en Pareto-fordeling med F(x)=1xαF(x) = 1 - x^{-\alpha} for x>1x > 1, med kjent α>0\alpha > 0. La W=αlnXW = \alpha \ln X.

a) Vis ved CDF-metoden at Weksp(1)W \sim \text{eksp}(1).

b) Angi fordelingen til 2i=1nαlnXi2\sum_{i=1}^{n}\alpha \ln X_i for et uavhengig utvalg.

📝Oppgave 9

En oppskriftsfordeling: ekstreme temperaturer XX følger en Gumbel-fordeling med F(x)=ee(xa)/bF(x) = e^{-e^{-(x-a)/b}} for alle reelle xx, med kjente aa og b>0b > 0. La W=e(Xa)/bW = e^{-(X-a)/b}.

a) Er g(x)=e(xa)/bg(x) = e^{-(x-a)/b} voksende eller avtakende? Vis ved CDF-metoden at Weksp(1)W \sim \text{eksp}(1).

b) Angi fordelingen til 2i=1ne(Xia)/b2\sum_{i=1}^{n} e^{-(X_i-a)/b}.

📝Oppgave 10

(Frihetsgrad-fellen, Des23-stil.) For et utvalg X1,,XnX_1, \ldots, X_n av eksponensiale med forventning β\beta er β^=Xˉ\hat{\beta} = \bar{X} maximum-likelihood-estimatoren (kap. 5.2).

a) Skriv 2nXˉ/β2n\bar{X}/\beta ved hjelp av Xi\sum X_i, og angi fordelingen.

b) En student skriver at 2nXˉ/β^χ2n22n\bar{X}/\hat{\beta} \sim \chi^2_{2n}. Forklar hvorfor det er meningsløst, og hva frihetsgraden i pivotalen faktisk avhenger av.

Oppskrifts- og formelkort

De fem kortene du bør kunne på autopilot — selve kjikvadrat-apparatet.

Kortene er flashcard-/repetisjonsstoff — de gjentar oppskriften. Hopp trygt over ved førstegangslesing; tidsanslaget gjelder drillen.

Oppskriftskort: kjikvadrat-apparatet i fire trinn
(1) Vis 2X/βχ222X/\beta \sim \chi^2_2 med CDF-metoden (F(y)=P(Xβy/2)=1ey/2F(y) = P(X \le \beta y/2) = 1 - e^{-y/2}, β\beta kanselleres). (2) Løft til summen med MGF-produktregelen: 2Xi/βχ2n22\sum X_i/\beta \sim \chi^2_{2n} (entydighet). (3) For Weibull/Gumbel/Pareto: potens-/log-transformér til eksponensial(1) først. (4) Husk at dette er pivotalen for eksakt KI (kap. 6.2) og eksakt test (Del 7).
Formelkort: 2X/βχ222X/\beta \sim \chi^2_2 per observasjon

For Xeksp(β)X \sim \text{eksp}(\beta): 2X/βχ222X/\beta \sim \chi^2_2. Utledning med CDF-metoden: P(2X/βy)=P(Xβy/2)=1ey/2P(2X/\beta \le y) = P(X \le \beta y/2) = 1 - e^{-y/2}, som er fordelingsfunksjonen til χ22\chi^2_2 (en eksponensial med forventning 2). Nøkkelen er at β\beta kanselleres i eksponenten — det er derfor størrelsen blir parameterfri.

Formelkort: 2Xi/βχ2n22\sum X_i/\beta \sim \chi^2_{2n} — frihetsgrad 2n2n
For uavhengige X1,,Xneksp(β)X_1, \ldots, X_n \sim \text{eksp}(\beta):

2i=1nXiβ=2nXˉβχ2n2.\frac{2\sum_{i=1}^{n} X_i}{\beta} = \frac{2n\bar{X}}{\beta} \sim \chi^2_{2n}.

Fra MGF-produktregelen: hver 2Xi/βχ222X_i/\beta \sim \chi^2_2 med MGF (12t)1(1-2t)^{-1}, og produktet (12t)n(1-2t)^{-n} er χ2n2\chi^2_{2n}. Frihetsgraden er 2n2n, aldri nn eller n1n-1, i eksponensialtilfellet.

Formelkort: potens-/log-løftet

Apparatet krever eksponensiale data. Transformér først: Weibull F(x)=1e(x/θ)kF(x) = 1 - e^{-(x/\theta)^k} gir (X/θ)keksp(1)(X/\theta)^k \sim \text{eksp}(1); Pareto F(x)=1xαF(x) = 1 - x^{-\alpha} gir αlnXeksp(1)\alpha \ln X \sim \text{eksp}(1); Gumbel F(x)=ee(xa)/bF(x) = e^{-e^{-(x-a)/b}} gir e(Xa)/beksp(1)e^{-(X-a)/b} \sim \text{eksp}(1). Deretter er 2(transformert)χ2n22\sum(\text{transformert}) \sim \chi^2_{2n}.

Feilkort: pivotalen krever ukjent β\beta

En pivotal er bygget av data og den ukjente parameteren, men har en fordeling uten ukjente (χ2n2\chi^2_{2n}). Derfor må β\beta (ikke estimatet β^\hat{\beta}) stå i nevneren: 2nXˉ/β2n\bar{X}/\beta. Setter du β^=Xˉ\hat{\beta} = \bar{X} inn, blir det 2n2n — en konstant, ubrukelig. Frihetsgraden 2n2n avhenger av antall observasjoner, ikke av om β\beta er kjent.

Symbol- og formelliste

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 Norges teknisk-naturvitenskapelige universitet. Dette er ikke offisielt studiemateriell. Les mer.