Bayes estimator

deejayosamu·2025년 8월 3일

통계 기본 개념

목록 보기
15/20

분포의 모수가 상수가 아니라 모수도 어떤 분포를 따른다고 봄.
=> f(x,θ)=f(xθ)h(θ)f(x,\theta) = f(x|\theta)h(\theta)

  • Joint pdf of (X,θ)(\underline{X},\theta)

    g(x1,...,xn,θ)=if(xiθ)h(θ)=L(θ;x)h(θ)g(x_1,...,x_n,\theta) = \prod_i f(x_i|\theta)h(\theta)=L(\theta;\underline{x})h(\theta)

  • Marginal pdf of X\underline{X}

    gX(x)=g(x,θ)dθg_X(\underline{x}) = \int g(\underline{x},\theta) d\theta

  • Posterior distriburion of θ\theta given the sample

    k(θx)=g(x,θ)gX(x)k(\theta|\underline{x}) = \frac{g(\underline{x},\theta)}{g_X(\underline{x})}

ex1) XiθiidPoisson(θ),θGamma(α,β) (α,β:known)X_i|\theta \overset{iid}{\sim} Poisson(\theta), \theta \sim Gamma(\alpha,\beta) \space (\alpha,\beta:known)
bayes-ex1

Bayesian point estimation

  • δ(x)\delta(\underline{x}): decision function
  • l(θ,δ(x))l(\theta,\delta(\underline{x})): loss function
    ex) (θδ(x))2,θδ(x),(θ/δ(x)1)2(\theta - \delta(\underline{x}))^2,|\theta - \delta(\underline{x})|,(\theta/\delta(\underline{x}) - 1)^2
  • Bayes estimator minimizes the posterior expected loss
    E[l(θ,δ(x))X=x]=l(θ,δ(x))k(θx)dθE[l(\theta,\delta(\underline{x}))|\underline{X}=\underline{x}]=\int l(\theta,\delta(\underline{x}))k(\theta|\underline{x}) d\theta

ex1) XiθiidPoisson(θ),θGamma(α,β) (α,β:known)X_i|\theta \overset{iid}{\sim} Poisson(\theta), \theta \sim Gamma(\alpha,\beta) \space (\alpha,\beta:known)
θxGamma(ixi+α,βnβ+1)When l(θ,δ(x))=(θδ(x))2, θ^=E(θx)θ^B=E(θx)=β(ixi+α)nβ+1=ixi+αn+1β=nn+1βx+1/βn+1βαβ\theta| \underline{x} \sim Gamma(\sum_i x_i + \alpha, \frac{\beta}{n \beta + 1}) \\ When \space l(\theta, \delta(\underline{x})) = (\theta - \delta(\underline{x}))^2, \space \hat{\theta}=E(\theta | \underline{x}) \\ \hat{\theta}^B = E(\theta | \underline{x}) = \frac{\beta(\sum_i x_i + \alpha)}{n \beta + 1} = \frac{\sum_i x_i + \alpha}{n + \frac{1}{\beta}} = \frac{n}{n + \frac{1}{\beta}} \overline{x} + \frac{1/ \beta}{n + \frac{1}{\beta}} \alpha \beta : weighted sum of MLE & mean of prior distritbution

ex2) Xiθiidb(1,θ),θBeta(α,β)(α,β:known)X_i|\theta \overset{iid}{\sim} b(1,\theta),\theta \sim Beta(\alpha,\beta)(\alpha,\beta:known)
bayes-ex2

✔︎ conjugate prior
: If the posterior pdfs is in the same family of distributions as the prior, the prior is called conjugate prior

ex)
① When XPoisson(θ)X \sim Poisson(\theta), Gamma is conjugate prior
② When Xb(n,θ)X \sim b(n,\theta), Beta is conjugate prior
③ When XN(μ,1)X \sim N(\mu,1), Normal is conjugate prior

0개의 댓글