跳到主要内容

早稲田大学 創造理工学研究科 経営システム工学専攻 2018年7月実施 情報数理応用 問題6

Author

祭音Myyura

Description

統計学、情報理論、パターン認識、機械学習に関する次の問いに答えよ。

  1. P(Xi=1θ)=θP(X_i=1\mid\theta)=\theta のベルヌーイ情報源に一様事前分布を仮定する。長さ nn の列で1が yy 回観測されたとき、θ\theta の事後密度、平均、最頻値を求めよ。
  2. 同じベルヌーイモデルで 2yn22\leq y\leq n-2 とする。尤度関数を示し、概形と最尤推定量を求めよ。
  3. kk-means 法の概要、長所、短所を説明せよ。
  4. 多次元データを分類するナイーブベイズ法の概要、長所、短所を説明せよ。
  5. 自然言語テキストに適用される形態素解析を説明せよ。
  6. 記号 a1,,aMa_1,\ldots,a_M が確率 p1,,pMp_1,\ldots,p_M で独立同分布に生成される情報源 XX のエントロピーを示し、その意味を複数の観点から説明せよ。

Kai

[小問 1]

一様事前分布は Beta(1,1)\operatorname{Beta}(1,1) である。尤度は θy(1θ)ny\theta^y(1-\theta)^{n-y} に比例するので、共役性より

θxnBeta(y+1,ny+1).\theta\mid x^n\sim\operatorname{Beta}(y+1,n-y+1).

したがって事後密度は

p(θxn)=Γ(n+2)Γ(y+1)Γ(ny+1)θy(1θ)ny,0θ1.\boxed{ p(\theta\mid x^n) =\frac{\Gamma(n+2)}{\Gamma(y+1)\Gamma(n-y+1)} \theta^y(1-\theta)^{n-y},\quad0\leq\theta\leq1 }.

事後平均と、0<y<n0<y<n の場合の最頻値は

E[θxn]=y+1n+2,θmode=yn.\boxed{E[\theta\mid x^n]=\frac{y+1}{n+2}},\qquad \boxed{\theta_{\mathrm{mode}}=\frac yn}.

y=0y=0 なら最頻値は0、y=ny=n なら1である。

[小問 2]

観測された具体的な列に対する尤度は

L(θ)=θy(1θ)ny,qquad0θ1.\boxed{L(\theta)=\theta^y(1-\theta)^{n-y}},qquad0\leq\theta\leq1.

2yn22\leq y\leq n-2 なので両端で0となり、区間内部の1点まで増加してから減少する。対数尤度を微分すると

ddθlogL(θ)=yθny1θ.\frac{d}{d\theta}\log L(\theta) =\frac y\theta-\frac{n-y}{1-\theta}.

これを0とおけば

θ^ML=yn.\boxed{\hat\theta_{\mathrm{ML}}=\frac yn}.

[小問 3]

kk-means 法は、各標本を最も近いクラスタ中心へ割り当てる操作と、各クラスタ中心を所属標本の平均へ更新する操作を、割当てが変化しなくなるまで交互に行う。クラスタ内平方和を局所的に最小化する方法である。

長所は実装が簡単で計算量が比較的小さく、大規模な数値データにも適用しやすいことである。短所は kk を事前に指定する必要があり、初期値によって局所解が変わること、外れ値・尺度に敏感で、球状かつ同程度の大きさでないクラスタを捉えにくいことである。

[小問 4]

特徴ベクトル x=(x1,,xd)\boldsymbol x=(x_1,\ldots,x_d) とクラス CkC_k に対し、条件付き独立

P(xCk)=j=1dP(xjCk)P(\boldsymbol x\mid C_k)=\prod_{j=1}^dP(x_j\mid C_k)

を仮定し、

C^=argmaxkP(Ck)jP(xjCk)\hat C=\arg\max_kP(C_k)\prod_jP(x_j\mid C_k)

で分類する。少量データでも学習と予測が高速で、高次元・疎なテキスト特徴に強い。一方、特徴間の依存を無視するため仮定が大きく破れると精度や確率校正が悪化し、未観測カテゴリには平滑化が必要となる。

[小問 5]

形態素解析は、空白で単語境界が明示されない日本語などの文を、意味を担う最小単位である形態素へ分割し、各形態素に原形、品詞、活用形、読みなどを付与する処理である。辞書と接続コストを用いる動的計画法や統計モデルにより、複数の分割・品詞候補から文全体として最も妥当な系列を選ぶ。検索、構文解析、機械翻訳などの前処理となる。

[小問 6]

対数の底を2とすれば、エントロピーは

H(X)=i=1Mpilog2pi\boxed{H(X)=-\sum_{i=1}^Mp_i\log_2p_i}

で、単位は bit である。

  • 観測前の不確実性、または観測によって平均的に得られる情報量を表す。
  • 一意復号可能な2進符号の平均符号長の理論的下限を与える。
  • 分布が一様なとき最大 log2M\log_2M、1記号に集中すると0となるため、記号の多様性・予測困難性を表す。
  • 独立な記号列ではエントロピーが加法的となり、典型列の個数がおよそ 2nH(X)2^{nH(X)} になるため、圧縮可能率を決める。