Tilbake
5.4

5.4 Drill: ML-utledning med estimatorvurdering (sjanger G)

Sjangerdrill på storoppgavens første halvdel: utled ML, vis egenskaper, sammenlign kandidater — i alle arkivets varianter.

90 min
11 oppgaver
DrillML-utledning med estimatorvurdering (sjanger G)
Din fremgang i kapitlet
0 / 11 oppgaver
Forkunnskaper: Dette drillkapitlet forutsetter hele Del 5: estimatorvurdering (kap. 5.1 — forventningsretthet, varians, MSE, sammenligningsprotokollen), ML-ritualet (kap. 5.2 — inkludert kombinert likelihood, bibetingelse og tilfellet uten lukket løsning) og eventuelt asymptotisk normalitet/deltametoden (kap. 5.3). Tetthetene hentes fra kap. 2.22.3.

Matematisk må du kunne derivere logaritmer og enkle uttrykk, bruke logaritmeregler (ln(ab)=lna+lnb\ln(ab)=\ln a+\ln b, ln(ab)=blna\ln(a^b)=b\ln a), og regne forventning og varians av veide summer. Ingen ny teknikk introduseres — bare drill på eksamensnivå.

Løsningsoppskriften (åtte trinn)

Enhver sjanger G-oppgave følger samme åtte trinn. De fem første er ML-ritualet fra kap. 5.2; de tre siste er vurderingen fra kap. 5.1. Lær rekkefølgen utenat.

1. Likelihood over ALLE observasjoner. L(θ)=f(xi;θ)L(\theta)=\prod f(x_i;\theta) — og ved to kilder produktet over begge utvalg. Ved bibetingelse: substituer bindingen inn , før derivasjon.

2. Log. l(θ)=lnL(θ)=lnf(xi;θ)l(\theta)=\ln L(\theta)=\sum\ln f(x_i;\theta) — produkt blir sum.

3. Deriver. l(θ)l'(\theta), sett lik 0 (likelihoodligningen).

4. Løs fullstendig. Kom helt frem til θ^=\hat{\theta}=\ldots. Er ligningen transcendental, sett den opp ryddig og si at den løses numerisk.

5. Maksimumskontroll. l(θ^)<0l''(\hat{\theta})<0 eller monotoni — trinnet de fleste hopper over.

6. Forventningsrett? Regn E(θ^)E(\hat{\theta}) og sjekk om den blir θ\theta.

7. Varians. Regn Var(θ^)\text{Var}(\hat{\theta}).

8. Ved flere kandidater: MSE-protokollen. Sjekk forventningsretthet for alle, regn varianser (eller MSE ved skjevhet), og velg med begrunnelse.

Samlet tid ≈ 90 min. Ta gjerne pause etter den gjennomregnede casen — den viser alle åtte trinnene på ett sted.

✏️Gjennomregnet eksamenscase (med sensornotater)

Levetidene x1,,xnx_1,\ldots,x_n (i år) til nn pumper er uavhengige og eksponensialfordelte med ukjent forventning β\beta, f(x;β)=1βex/β\displaystyle f(x;\beta)=\frac{1}{\beta}e^{-x/\beta}.

a) Utled ML-estimatoren β^\hat{\beta} og kontroller at det er et maksimum.

b) Er β^\hat{\beta} forventningsrett? Finn variansen.

c) En kollega foreslår i stedet β~=X1\tilde{\beta}=X_1 (bruk kun den første målingen). Sammenlign de to estimatorene og velg med begrunnelse.

a) Trinn 1–4 (ML).
L(β)=βnexp ⁣(1βxi),l(β)=nlnβ1βxi,l(β)=nβ+1β2xi=0.L(\beta)=\beta^{-n}\exp\!\Big(-\tfrac{1}{\beta}\sum x_i\Big),\quad l(\beta)=-n\ln\beta-\tfrac{1}{\beta}\sum x_i,\quad l'(\beta)=-\tfrac{n}{\beta}+\tfrac{1}{\beta^2}\sum x_i=0.
Løsning: β^=1nxi=Xˉ\displaystyle \hat{\beta}=\frac{1}{n}\sum x_i=\bar{X}.

> Sensornotat (2 p): likelihooden må stå som produkt over alle nn observasjoner, og ligningen løses helt frem til β^=Xˉ\hat{\beta}=\bar{X}.

Trinn 5 (maksimumskontroll). l(β)=nβ22β3xi\displaystyle l''(\beta)=\frac{n}{\beta^2}-\frac{2}{\beta^3}\sum x_i; med xi=nβ^\sum x_i=n\hat{\beta} blir l(β^)=nβ^2<0\displaystyle l''(\hat{\beta})=-\frac{n}{\hat{\beta}^2}<0. Maksimum.

> Sensornotat (1 p): dette trinnet skiller toppbesvarelsen. Én linje holder.

b) Trinn 6–7. E(β^)=E(Xˉ)=βE(\hat{\beta})=E(\bar{X})=\beta (siden E(Xi)=βE(X_i)=\beta), så forventningsrett. Varians: Var(β^)=Var(Xˉ)=Var(X)n=β2n\displaystyle \text{Var}(\hat{\beta})=\text{Var}(\bar{X})=\frac{\text{Var}(X)}{n}=\frac{\beta^2}{n}.

> Sensornotat (1 p): skriv estimatoren med stor bokstav (β^=Xˉ\hat{\beta}=\bar{X}) — det er estimatoren (stokastisk) som har varians, ikke tallet.

c) Trinn 8 (MSE-protokollen). Begge kandidater er forventningsrette: E(β~)=E(X1)=βE(\tilde{\beta})=E(X_1)=\beta. Da avgjør variansen. Var(β^)=β2/n\text{Var}(\hat{\beta})=\beta^2/n, mens Var(β~)=Var(X1)=β2\text{Var}(\tilde{\beta})=\text{Var}(X_1)=\beta^2. For n>1n>1 er β2/n<β2\beta^2/n<\beta^2, så β^=Xˉ\hat{\beta}=\bar{X} velges — den bruker all informasjonen i utvalget, mens β~\tilde{\beta} kaster n1n-1 målinger.

> Sensornotat (2 p): begge forventningsrette → sammenlign varians; begrunn valget. Hadde én vært skjev, måtte du regnet MSE =Var+bias2=\text{Var}+\text{bias}^2 for begge.

— naturlig pausepunkt — nå driller du selv, på alle variantene.

📝Oppgave 1

(Innstegsoppgave.) Antall x1,,xnx_1,\ldots,x_n hendelser i uavhengige intervaller er Poisson-fordelt med forventning λ\lambda.

a) Skriv opp log-likelihooden l(λ)l(\lambda) (dropp konstantledd).

b) Utled λ^\hat{\lambda} og gjør maksimumskontrollen.

📝Oppgave 2

Bruddstyrken til nn prøver er eksponensialfordelt med forventning β\beta. Utled ML-estimatoren og oppgi E(β^)E(\hat{\beta}) og Var(β^)\text{Var}(\hat{\beta}).

📝Oppgave 3

En kvalitetstest på nn enheter finner xx defekte; andelen defekte er pp (binomisk).

a) Utled ML-estimatoren p^\hat{p}.

b) Vis at p^\hat{p} er forventningsrett og finn Var(p^)\text{Var}(\hat{p}).

📝Oppgave 4

Bruddstyrken til nn fiberprøver er gammafordelt med kjent formparameter α=3\alpha=3 og ukjent skala θ\theta: f(x)=1Γ(3)θ3x2ex/θ\displaystyle f(x)=\frac{1}{\Gamma(3)\theta^3}x^2 e^{-x/\theta}, med E(X)=3θE(X)=3\theta.

a) Utled ML-estimatoren θ^\hat{\theta}.

b) Er den forventningsrett?

📝Oppgave 5

En «oppskriftsfordeling» (Pareto-type) har tetthet f(x;θ)=θx(θ+1)f(x;\theta)=\theta\,x^{-(\theta+1)} for x>1x>1, θ>0\theta>0. Du har uavhengige observasjoner x1,,xnx_1,\ldots,x_n.

a) Skriv opp log-likelihooden l(θ)l(\theta).

b) Utled ML-estimatoren θ^\hat{\theta} og gjør maksimumskontrollen.

📝Oppgave 6

Antall mikrosprekker per prøve er Poisson-fordelt med forventning λ\lambda. Ett laboratorium har talt sprekker i nn prøver (x1,,xnx_1,\ldots,x_n), et annet i mm prøver med dobbel prøvestørrelse, så disse har forventning 2λ2\lambda (y1,,ymy_1,\ldots,y_m). Alt uavhengig.

a) Skriv opp den kombinerte log-likelihooden l(λ)l(\lambda).

b) Utled ML-estimatoren λ^\hat{\lambda} basert på begge datasettene.

c) Vis at λ^\hat{\lambda} er forventningsrett.

📝Oppgave 7

To måleserier av feilmellomrom er eksponensialfordelte. Serie A har forventning μA\mu_A, serie B har forventning μB\mu_B, og det er kjent at μA=3μB\mu_A=3\mu_B (A-utstyret er tregere). Vi har pp observasjoner a1,,apa_1,\ldots,a_p fra A og qq observasjoner b1,,bqb_1,\ldots,b_q fra B.

a) Substituer bibetingelsen og skriv opp log-likelihooden i den frie parameteren.

b) Utled ML-estimatoren og oppgi μ^A\hat{\mu}_A og μ^B\hat{\mu}_B.

📝Oppgave 8

Antall defekter per produksjonsbatch modelleres med null-trunkert Poisson (hver batch har minst én defekt), p(x)=eλλxx!(1eλ)p(x)=\dfrac{e^{-\lambda}\lambda^x}{x!(1-e^{-\lambda})} for x1x\ge 1.

a) Skriv opp log-likelihooden l(λ)l(\lambda) for nn batcher.

b) Deriver og reduser likelihoodligningen til en ligning i λ\lambda og xˉ\bar{x}.

c) Kan den løses eksplisitt? Hva gjør du?

📝Oppgave 9
To estimatorer for forventningen μ\mu til n=5n=5 uavhengige, likt fordelte målinger (varians σ2\sigma^2) vurderes:
μ^1=Xˉ,μ^2=14X1+14X2+14X3+14X4.\hat{\mu}_1=\bar{X},\qquad \hat{\mu}_2=\tfrac{1}{4}X_1+\tfrac{1}{4}X_2+\tfrac{1}{4}X_3+\tfrac{1}{4}X_4.

a) Vis at begge er forventningsrette.

b) Regn variansene og velg med begrunnelse.

📝Oppgave 10
(Regresjons-forberedelse.) I en modell uten konstantledd er Yi=βxi+εiY_i=\beta x_i+\varepsilon_i, der xix_i er kjente tall og εi\varepsilon_i er uavhengige med forventning 0 og varians σ2\sigma^2. To forventningsrette estimatorer for β\beta vurderes:
β^=xiYixi2(minste kvadrater),β~=1nYixi.\hat{\beta}=\frac{\sum x_iY_i}{\sum x_i^2}\qquad(\text{minste kvadrater}),\qquad \tilde{\beta}=\frac{1}{n}\sum\frac{Y_i}{x_i}.

a) Vis at begge er forventningsrette.

b) Finn variansen til hver.

c) For x=(1,2,3,4,5)x=(1,2,3,4,5) og σ2=1\sigma^2=1: regn ut begge variansene og velg.

📝Oppgave 11

Levetider x1,,xnx_1,\ldots,x_n er eksponensialfordelte med forventning β\beta. Vi vil estimere medianlevetiden m=βln2m=\beta\ln 2.

a) Bruk invariansprinsippet til å finne ML-estimatoren m^\hat{m}.

b) Bruk deltametoden til å finne en tilnærmet varians for m^\hat{m}, gitt at Var(β^)=β2/n\text{Var}(\hat{\beta})=\beta^2/n.

Oppskrifts- og formelkort

De fem kortene du bør kunne på autopilot før eksamen — selve algoritmen for sjanger G.

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

Oppskriftskort: de åtte trinnene

Sjanger G, i rekkefølge: (1) likelihood over ALLE observasjoner (begge kilder; bibetingelse inn nå); (2) log; (3) deriver; (4) løs fullstendig (transcendental → numerisk); (5) maksimumskontroll l<0l''<0; (6) E(θ^)E(\hat{\theta}) — forventningsrett?; (7) Var(θ^)\text{Var}(\hat{\theta}); (8) ved flere kandidater: MSE-protokollen og VELG med begrunnelse.

Formelkort: standard-ML-resultater

Eksponensial: β^=Xˉ\hat{\beta}=\bar{X} (forventningsrett, Var=β2/n\text{Var}=\beta^2/n). Poisson: λ^=Xˉ\hat{\lambda}=\bar{X}. Andel: p^=X/n\hat{p}=X/n (Var=p(1p)/n\text{Var}=p(1-p)/n). Normal: μ^=Xˉ\hat{\mu}=\bar{X}, σ^2=1n(XiXˉ)2\displaystyle \hat{\sigma}^2=\frac1n\sum(X_i-\bar X)^2 (skjev). Gamma kjent form α\alpha: θ^=Xˉ/α\hat{\theta}=\bar{X}/\alpha. Oppskriftsfordeling: bruk ritualet på oppgitt ff.

Teknikkort: to kilder og bibetingelse
To kilder: L=f(xi)g(yj)L=\prod f(x_i)\cdot\prod g(y_j); deriver den samlede ll. Ulik forventning → vektet nevner (n+2mn+2m når kilde 2 har 2λ2\lambda). Bibetingelse: substituer bindingen (μA=kμB\mu_A=k\mu_B) inn FØR derivasjon, så én fri parameter; sett de andre tilbake med invarians.
Oppskriftskort: MSE-protokollen (trinn 8)

Ved flere kandidater: (1) sjekk forventningsretthet (regn E(θ^)E(\hat{\theta})); (2) regn varianser; (3) alle forventningsrette → minst varians; ellers MSE =Var+bias2=\text{Var}+\text{bias}^2 for alle; (4) VELG og begrunn. Forventningsrett slår ikke automatisk skjev.

Feilkort: transcendental likelihoodligning

Gir l(θ)=0l'(\theta)=0 en ligning som blander θ\theta og eθe^{-\theta} (null-trunkert Poisson: xˉ=λ/(1eλ)\bar{x}=\lambda/(1-e^{-\lambda})), finnes ingen lukket løsning. Sett opp ligningen ryddig og oppgi at θ^\hat{\theta} løses numerisk. Ikke lat som det finnes en formel.

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.