跳到主要内容

電気通信大学 情報理工学研究科 情報・ネットワーク工学専攻 2025年8月実施 選択問題 確率統計

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

XBin(m1,q)X\sim\operatorname{Bin}(m_1,q)YBin(m2,q)Y\sim\operatorname{Bin}(m_2,q) を独立とする。XX の積率母関数、期待値・分散、qq の最尤推定量を求めよ。さらに QBeta(α,β)Q\sim\operatorname{Beta}(\alpha,\beta) の下で、X=xX=x の観測後分布と YY の予測分布、および S=X+YS=X+Y の条件付き分布を求めよ。

题目描述

给定参数相同的两个独立二项分布,求矩生成函数、期望与方差、极大似然估计;再以 Beta 分布为先验,求后验参数、后验预测分布及和的条件分布。

Kai

(1)

MX(t)=x=0m1(m1x)(qet)x(1q)m1x=(qet+1q)m1.\begin{aligned} M_X(t) &=\sum_{x=0}^{m_1}\binom{m_1}{x}(qe^t)^x(1-q)^{m_1-x}\\ &=\boxed{(qe^t+1-q)^{m_1}}. \end{aligned}

(2)

MX(0)M_X'(0)MX(0)M_X''(0) より、

E[X]=m1q,V[X]=m1q(1q).\boxed{E[X]=m_1q},\qquad \boxed{V[X]=m_1q(1-q)}.

(3)

L(q)=(m1x)qx(1q)m1xL(q)=\binom{m_1}{x}q^x(1-q)^{m_1-x}

とすると、

ddqlogL(q)=xqm1x1q=0.\frac{d}{dq}\log L(q)=\frac{x}{q}-\frac{m_1-x}{1-q}=0.

したがって、

q^=xm1.\boxed{\widehat q=\frac{x}{m_1}}.

(4)

Bayes の定理より、

fQX=x(q)qx+α1(1q)m1x+β1.f_{Q\mid X=x}(q)\propto q^{x+\alpha-1}(1-q)^{m_1-x+\beta-1}.

よって、

α=x+α,β=m1x+β.\boxed{\alpha'=x+\alpha},\qquad \boxed{\beta'=m_1-x+\beta}.

(5)

P(Y=yX=x)=01(m2y)qy(1q)m2yqα1(1q)β1B(α,β)dq=(m2y)B(y+x+α, m2y+m1x+β)B(x+α, m1x+β).\begin{aligned} P(Y=y\mid X=x) &=\int_0^1\binom{m_2}{y}q^y(1-q)^{m_2-y} \frac{q^{\alpha'-1}(1-q)^{\beta'-1}}{B(\alpha',\beta')}\,dq\\ &=\boxed{ \binom{m_2}{y} \frac{B(y+x+\alpha,\ m_2-y+m_1-x+\beta)} {B(x+\alpha,\ m_1-x+\beta)}}. \end{aligned}

(6)

Vandermonde の恒等式を用いると、

P(S=sQ=q)=k(m1k)(m2sk)qs(1q)m1+m2s=(m1+m2s)qs(1q)m1+m2s.\begin{aligned} P(S=s\mid Q=q) &=\sum_k\binom{m_1}{k}\binom{m_2}{s-k} q^s(1-q)^{m_1+m_2-s}\\ &=\binom{m_1+m_2}{s}q^s(1-q)^{m_1+m_2-s}. \end{aligned}

したがって、

SQ=qBin(m1+m2,q).\boxed{S\mid Q=q\sim\operatorname{Bin}(m_1+m_2,q)}.