Tilbake
2.3

2.3 Drill: den egendefinerte og «eksotiske» tettheten (sjanger C)

Standardtreningen i kontinuerlig sannsynlighet: fra oppgitt tetthet til normering, CDF, median og momenter — og beredskapen mot et ukjent fordelingsnavn.

80 min
10 oppgaver
Drillden egendefinerteeksotiske» tettheten (sjanger C)
Din fremgang i kapitlet
0 / 10 oppgaver
Forkunnskaper: Alt fra kap. 2.1 (normeringskonstant, stykkevis FF, median ved F(m)=pF(m) = p, momenter ved integrasjon, eksponensial) og kap. 2.2 (gammaintegralet 0xa1ex/bdx=baΓ(a)\int_0^\infty x^{a-1}e^{-x/b}dx = b^a\Gamma(a), gammafunksjonens regneregler).

Sist du var her — de tre grepene du bruker mest:
- Normering: f=1\int f = 1 bestemmer konstanten.
- Median: løs F(m)=12F(m) = \tfrac12 (potens-, log- eller andregradslikning).
- Momenter: E(X)=xfdxE(X) = \int x f\,dx, ofte via gammaintegralet.

Nye fordelingsnavn (Laplace, Pareto, Rayleigh, Weibull, lognormal) forklares der de dukker opp — du trenger ikke ha sett dem før.

Løsningsoppskriften for sjanger C

Møter du en oppgitt tetthet f(x)f(x) (kjent eller ukjent fordelingsnavn), kjør denne faste algoritmen — den løser praktisk talt hver sjanger C-oppgave:

1. Normer. Er tettheten gitt «proporsjonal med» noe, finn konstanten av f=1\int f = 1. Er ff allerede komplett, hopp over — men sjekk gjerne at den integrerer til 1.

2. Utled F(x)F(x). Integrer tettheten fra bunnen av støtten: F(x)=Axf(t)dtF(x) = \int_{A}^{x} f(t)\,dt. Skriv den stykkevis — 0 under støtten, formel på støtten, 1 over. (For fordelinger definert ved halen, som Pareto/Weibull, er det ofte lettere å regne P(X>x)=1F(x)P(X > x) = 1 - F(x) direkte.)

3. Regn sannsynligheter fra FF: P(a<Xb)=F(b)F(a)P(a < X \le b) = F(b) - F(a), P(X>x)=1F(x)P(X > x) = 1 - F(x), og betingede via definisjonen.

4. Median og kvantiler: løs F(m)=12F(m) = \tfrac12 (og F(xp)=pF(x_p) = p). Dette blir en potenslikning (Pareto), logaritmelikning (Weibull/Rayleigh) eller andregradslikning — regn nøyaktig.

5. Momenter: E(X)=xfdxE(X) = \int x f\,dx, E(X2)=x2fdxE(X^2) = \int x^2 f\,dx, Var=E(X2)[E(X)]2\text{Var} = E(X^2) - [E(X)]^2. Bruk gammaintegralet og Γ\Gamma-reglene fra kap. 2.2 der integranden har form xa1ex^{a-1}e^{-\ldots}.

6. Pek fremover: kvantilfunksjonen F1F^{-1} er råstoffet i inversjonsmetoden for simulering (kap. 7.1), og en passende potens- eller logaritmetransformasjon gjør ofte den eksotiske fordelingen om til en eksponensial/gamma (kap. 4.3) — nøkkelen til estimering senere.

✏️Gjennomregnet eksamenscase: Pareto-fordelt inntekt (med sensornotater)
Årsinntekten XX (i tusen kroner) over en minsteinntekt k=250k = 250 er Pareto-fordelt med tetthet

f(x)=θkθxθ+1,xk,θ=2,5.f(x) = \frac{\theta\,k^\theta}{x^{\theta+1}}, \qquad x \ge k, \quad \theta = 2{,}5.

a) Vis at ff er en gyldig tetthet og finn fordelingsfunksjonen FF.
b) Finn medianinntekten.
c) Finn forventet inntekt E(X)E(X).
d) Vis at 2θln(X/k)2\theta\ln(X/k) er kjikvadratfordelt, og angi frihetsgradene.

a) Sensornotat: start med å bekrefte at tettheten er gyldig — 1 poeng nesten gratis. Integrer fra kk:

F(x)=kxθkθtθ+1dt=θkθ[tθθ]kx=kθ(kθxθ)=1(kx)θ,xk.F(x) = \int_k^x \frac{\theta k^\theta}{t^{\theta+1}}\,dt = \theta k^\theta\left[\frac{t^{-\theta}}{-\theta}\right]_k^x = k^\theta\left(k^{-\theta} - x^{-\theta}\right) = 1 - \left(\frac{k}{x}\right)^{\theta}, \quad x \ge k.

Stykkevis: F(x)=0F(x) = 0 for x<kx < k og F(x)=1(k/x)θF(x) = 1 - (k/x)^\theta for xkx \ge k. At F(k)=0F(k) = 0 og F()=1F(\infty) = 1 bekrefter at tettheten er gyldig. Sensornotat: her holder det å vise F()=1F(\infty) = 1.

b) Sensornotat: median = løs F(m)=12F(m) = \tfrac12, ikke f(m)f(m).

1(km)θ=12    (km)θ=12    km=21/θ    m=k21/θ.1 - \left(\frac{k}{m}\right)^\theta = \tfrac12 \;\Rightarrow\; \left(\frac{k}{m}\right)^\theta = \tfrac12 \;\Rightarrow\; \frac{k}{m} = 2^{-1/\theta} \;\Rightarrow\; m = k\cdot 2^{1/\theta}.

Med k=250k = 250, θ=2,5\theta = 2{,}5: m=25020,42501,320=329,9m = 250\cdot 2^{0{,}4} \approx 250\cdot 1{,}320 = 329{,}9 tusen kroner.

c) Sensornotat: her kreves θ>1\theta > 1 for at integralet konvergerer — nevn betingelsen.

E(X)=kxθkθxθ+1dx=θkθkxθdx=θkθ[x1θ1θ]k.E(X) = \int_k^\infty x\cdot\frac{\theta k^\theta}{x^{\theta+1}}\,dx = \theta k^\theta\int_k^\infty x^{-\theta}\,dx = \theta k^\theta\left[\frac{x^{1-\theta}}{1-\theta}\right]_k^\infty.

For θ>1\theta > 1 går den øvre grensen mot 0, og

E(X)=θkθk1θθ1=θkθ1=2,52501,5416,7 tusen kroner.E(X) = \theta k^\theta\cdot\frac{k^{1-\theta}}{\theta - 1} = \frac{\theta k}{\theta - 1} = \frac{2{,}5\cdot 250}{1{,}5} \approx 416{,}7 \text{ tusen kroner}.

Merk: E(X)>mE(X) > m — fordelingen er høyreskjev (rike-halen drar snittet opp).

d) Sensornotat: dette er A-differensiatoren — bruk CDF-metoden. La Y=2θln(X/k)Y = 2\theta\ln(X/k). Siden XkX \ge k er Y0Y \ge 0. For y0y \ge 0:

P(Yy)=P ⁣(2θlnXky)=P ⁣(Xkey/(2θ))=F ⁣(key/(2θ)).P(Y \le y) = P\!\left(2\theta\ln\tfrac{X}{k} \le y\right) = P\!\left(X \le k\,e^{y/(2\theta)}\right) = F\!\left(k e^{y/(2\theta)}\right).

Sett inn FF: (k/(key/(2θ)))θ=ey/2\left(k/(k e^{y/(2\theta)})\right)^\theta = e^{-y/2}, så

P(Yy)=1ey/2.P(Y \le y) = 1 - e^{-y/2}.

Det er fordelingsfunksjonen til en eksponensial med forventning 2, altså gamma(1,2)=χ22\text{gamma}(1, 2) = \chi^2_2. Y=2θln(X/k)χ22Y = 2\theta\ln(X/k) \sim \chi^2_2 (2 frihetsgrader). Sensornotat: nettopp denne transformasjonen er broen til den eksakte pivotalen for θ\theta i Del 6.

> — naturlig pausepunkt — nå har du sett hele oppskriften på ett case. Under følger drilloppgavene.

📝Oppgave 1

En variabel har tetthet f(x)=cx3f(x) = c\,x^3 for 0x10 \le x \le 1 og 0 ellers.

a) Finn cc.
b) Finn F(x)F(x) på støtten og medianen.

📝Oppgave 2

For tettheten f(x)=4x3f(x) = 4x^3[0,1][0,1] (fra oppgave 1), finn E(X)E(X).

📝Oppgave 3

Laplace-fordelingen (dobbelt-eksponensial) har tetthet f(x)=12bex/bf(x) = \tfrac{1}{2b}e^{-|x|/b} for alle reelle xx, med b=2b = 2.

a) Forklar hvorfor medianen er 0.
b) Finn P(X>2)P(X > 2).

📝Oppgave 4

Bølgehøyden XX (meter) er Rayleigh-fordelt med tetthet f(x)=xσ2ex2/(2σ2)f(x) = \tfrac{x}{\sigma^2}e^{-x^2/(2\sigma^2)}, x0x \ge 0, med σ=2\sigma = 2.

a) Utled F(x)F(x).
b) Finn medianen.
c) Finn P(X>3)P(X > 3).

📝Oppgave 5

Levetiden XX (timer) til en komponent er Weibull-fordelt med f(x)=αβ(xβ)α1e(x/β)αf(x) = \tfrac{\alpha}{\beta}\left(\tfrac{x}{\beta}\right)^{\alpha-1}e^{-(x/\beta)^\alpha}, x0x \ge 0, med formparameter α=2\alpha = 2 og skala β=100\beta = 100.

a) Vis at F(x)=1e(x/β)αF(x) = 1 - e^{-(x/\beta)^\alpha}.
b) Finn medianlevetiden.
c) Finn P(X>120)P(X > 120).

📝Oppgave 6

Alkoholforbruket XX (liter per år) er lognormalfordelt: lnXN(μ,σ2)\ln X \sim N(\mu, \sigma^2) med μ=1\mu = 1 og σ=0,5\sigma = 0{,}5. Da er F(x)=Φ ⁣(lnxμσ)F(x) = \Phi\!\left(\tfrac{\ln x - \mu}{\sigma}\right), der Φ\Phi er standardnormal fordelingsfunksjon.

a) Finn medianen.
b) Finn P(X>5)P(X > 5) (bruk Φ(1,22)0,889\Phi(1{,}22) \approx 0{,}889).

📝Oppgave 7

Formuen XX (millioner) er Pareto-fordelt med f(x)=θkθ/xθ+1f(x) = \theta k^\theta/x^{\theta+1}, xkx \ge k, med k=1k = 1 og θ=3\theta = 3.

a) Finn F(x)F(x) og medianen.
b) Finn E(X)E(X).

📝Oppgave 8

Inntekten XX er Pareto-fordelt med kjent minsteinntekt kk og ukjent formparameter θ\theta: F(x)=1(k/x)θF(x) = 1 - (k/x)^\theta, xkx \ge k.

a) Vis at Y=ln(X/k)Y = \ln(X/k) er eksponensialfordelt, og angi forventningen.
b) Bruk dette til å vise at 2θln(X/k)χ222\theta\ln(X/k) \sim \chi^2_2.

📝Oppgave 9

For en Weibull-variabel XX med F(x)=1e(x/β)αF(x) = 1 - e^{-(x/\beta)^\alpha} (α,β\alpha, \beta kjent), vis at Z=(X/β)αZ = (X/\beta)^\alpha er eksponensialfordelt med forventning 1, og finn deretter medianen til ZZ.

📝Oppgave 10

En stykkevis tetthet er f(x)=cx(2x)f(x) = c\,x(2 - x) for 0x20 \le x \le 2 og 0 ellers.

a) Finn cc.
b) Finn E(X)E(X), og forklar hvorfor svaret er 1 uten å regne (symmetri).

Oppskriftskort — flashcards

Flashcard-/repetisjonsstoff — hopp trygt over ved førstegangslesing. Formel- og oppskriftskort for sjanger C, til rask oppslag under drilling.

Sjanger C — oppskriften

Oppgitt tetthet: (1) normer (f=1\int f = 1), (2) utled stykkevis FF, (3) regn PP-er fra FF, (4) median/kvantil ved F(m)=pF(m) = p, (5) momenter ved integrasjon/gammaintegralet, (6) pek fremover (F1F^{-1}, transformasjon til gamma).

Normeringskortet

Er f=cg(x)f = c\,g(x): c=1/gdxc = 1/\int g\,dx over støtten. Alltid første steg ved «proporsjonal med». Uten riktig cc blir alt etterpå galt.

CDF-kortet
F(x)=Axf(t)dtF(x) = \int_A^x f(t)\,dt, stykkevis (0 under støtten, 1 over). Sjekk F(A)=0F(A) = 0, F(B)=1F(B) = 1. For hale-definerte fordelinger: regn P(X>x)=1F(x)P(X > x) = 1 - F(x) direkte.
Median-kortet

Løs F(m)=12F(m) = \tfrac12 — aldri f(m)f(m). Typiske former: potens ((k/m)θ=12(k/m)^\theta = \tfrac12), logaritme ((m/β)α=ln2(m/\beta)^\alpha = \ln 2), andregrad.

Momentkortet
E(X)=xfdxE(X) = \int x f\,dx, E(X2)=x2fdxE(X^2) = \int x^2 f\,dx, Var=E(X2)[E(X)]2\text{Var} = E(X^2) - [E(X)]^2. Har integranden form xa1ex/bx^{a-1}e^{-x/b}, bruk gammaintegralet baΓ(a)b^a\Gamma(a).
Laplace (dobbelt-eksponensial)
f(x)=12bexμ/bf(x) = \tfrac{1}{2b}e^{-|x-\mu|/b}, alle reelle xx. Symmetrisk om μ\mu, så median =μ= \mu og E(X)=μE(X) = \mu; Var(X)=2b2\text{Var}(X) = 2b^2. To eksponensialhaler speilet.
Pareto (tung hale)
f(x)=θkθ/xθ+1f(x) = \theta k^\theta/x^{\theta+1}, xkx \ge k. F(x)=1(k/x)θF(x) = 1 - (k/x)^\theta. Median =k21/θ= k\cdot 2^{1/\theta}. E(X)=θk/(θ1)E(X) = \theta k/(\theta-1) for θ>1\theta > 1; endelig varians krever θ>2\theta > 2.
Pareto → gamma

For Pareto er ln(X/k)\ln(X/k) eksponensial med rate θ\theta, og 2θln(X/k)χ22=gamma(1,2)2\theta\ln(X/k) \sim \chi^2_2 = \text{gamma}(1,2). Broen til den eksakte pivotalen for θ\theta (Del 6).

Rayleigh
f(x)=xσ2ex2/(2σ2)f(x) = \tfrac{x}{\sigma^2}e^{-x^2/(2\sigma^2)}, x0x \ge 0. F(x)=1ex2/(2σ2)F(x) = 1 - e^{-x^2/(2\sigma^2)}. Median =σ2ln2= \sigma\sqrt{2\ln 2}, E(X)=σπ/2E(X) = \sigma\sqrt{\pi/2}. Modell for bølgehøyder/vindstyrke.
Weibull
f(x)=αβ(x/β)α1e(x/β)αf(x) = \tfrac{\alpha}{\beta}(x/\beta)^{\alpha-1}e^{-(x/\beta)^\alpha}, x0x \ge 0. F(x)=1e(x/β)αF(x) = 1 - e^{-(x/\beta)^\alpha}. Median =β(ln2)1/α= \beta(\ln 2)^{1/\alpha}. Ved α=1\alpha = 1: eksponensial. Levetidsmodell med aldring.
Weibull → eksponensial

For Weibull er (X/β)α(X/\beta)^\alpha eksponensial med forventning 1 (rate 1). Potenstransformasjonen fjerner formparameteren — grepet i kap. 4.3.

Lognormal
lnXN(μ,σ2)\ln X \sim N(\mu, \sigma^2), så F(x)=Φ ⁣(lnxμσ)F(x) = \Phi\!\left(\tfrac{\ln x - \mu}{\sigma}\right), x>0x > 0. Median =eμ= e^\mu. Modell for inntekt/forbruk med multiplikativ variasjon.
Potenstetthet
f(x)=(k+1)xkf(x) = (k+1)x^k[0,1][0,1] gir F(x)=xk+1F(x) = x^{k+1}, median =(12)1/(k+1)= (\tfrac12)^{1/(k+1)}, E(X)=k+1k+2E(X) = \tfrac{k+1}{k+2}. Enkleste sjanger C — ren normering og integrasjon.
Konvergensbetingelser

Tunge haler har ikke alltid endelige momenter: Pareto har endelig E(X)E(X) bare for θ>1\theta > 1, endelig Var\text{Var} for θ>2\theta > 2. Sjekk at integralet konvergerer før du oppgir et tall.

Ukjent fordelingsnavn — beredskap

Et uvant navn skal ikke skremme. Oppgaven gir tetthet/CDF/MGF. Kjør oppskriften (normer → FFPP → median → momenter). Navnet er bare en etikett; teknikkene er de vanlige.

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 Universitetet i Oslo. Dette er ikke offisielt studiemateriell. Les mer.