Bagaimana plot pdf ini

Oct 08 2020

Membiarkan $x\sim N(\mu,\sigma^2)$ dan $y=\max\{e^{x-1}-\kappa,0\}$. Bagaimana cara mendapatkan dan memplot pdf dari y? Secara khusus, akan ada massa probabilitas pada$y=0$, yang merupakan titik yang tidak dapat dibedakan. Saya mencoba

pdfY[y_] := \[Piecewise] {
   {1/(y + κ) f[Log[y + κ] + 1], y > 0},
   {N[Integrate[f[x], {x, -10, Log[κ] + 1}]], y = 0},
   {0, y < 0}
  }

tapi tidak bekerja dengan baik.

Jawaban

wuyudi Oct 08 2020 at 10:23
TransformedDistribution[Max[E^(x-1)-κ,0],x\[Distributed] NormalDistribution[μ,σ]]//PDF[#,y]&
1/2 DiracDelta[y+κ] (1+Erf[(1-μ+Log[y+κ])/(Sqrt[2] σ)]) (-1+UnitStep[-y]) (UnitStep[-y] (-1+UnitStep[-κ])-UnitStep[-κ])-DiracDelta[y] (1+Erf[(1-μ+Log[κ])/(Sqrt[2] σ)]) (-1+UnitStep[-κ])-(E^(-((1-μ+Log[y+κ])^2/(2 σ^2))) (-1+UnitStep[-y]) (1+UnitStep[-y] (-1+UnitStep[-y-κ]) (-1+UnitStep[-κ])-UnitStep[-y-κ] UnitStep[-κ]))/(Sqrt[2 π] (y+κ) σ)

lalu Anda bisa memplotnya.