線形回帰モデル
、 ε 〜N (0 、Iは、σ 2)Y= Xβ+ ϵε 〜N(0 、私σ2)
、 X ∈ R N × Pと β ∈ R PY∈RnX∈Rn×pβ∈Rp
ϵ=Y−XββL2
最小二乗
(x1,y1),...,(xn,yn)xip
βˆLS=argminβ||ϵ||2=argminβ||Y−Xβ||2=argminβ∑i=1n(yi−xiβ)2
最尤
β
L(Y|X,β)=∏i=1nf(yi|xi,β)
f(yi|xi,β)σ2
L(Y|X,β)=∏i=1n12πσ2−−−−√e−(yi−xiβ)22σ2
一般に、可能性を処理する場合、続行する前に数学的にログを取る方が簡単です(積は合計になり、指数はなくなります)。そうしましょう。
logL(Y|X,β)=∑i=1nlog(12πσ2−−−−√)−(yi−xiβ)22σ2
ββ
βˆMLE=argmaxβ∑i=1n−(yi−xiβ)22σ2
β
βˆMLE=argminβ∑i=1n(yi−xiβ)2=βˆLS
これが機能するためには、特定のモデルの仮定(誤差項の正規性、0平均、一定の分散)を作成する必要があったことを思い出してください。これにより、特定の条件下で最小二乗がMLEと同等になります。詳細については、こちらとこちらをご覧ください。
完全を期すために、ソリューションは次のように記述できることに注意してください。
β=(XTX)−1XTy