2
ときの尤度の計算
私はこの事後分布を計算しようとしています： （θ | − ）= ∏んi = 1py私私（1 − p私）1 − y私Σすべてθ 、p私| θΠんi = 1py私私（1 − p私）1 − y私（θ|−）=Π私=1んp私y私（1−p私）1−y私Σすべてθ、p私|θΠ私=1んp私y私（1−p私）1−y私 (\theta|-)=\frac{\prod_{i=1}^{n}p_i^{y_i}(1-p_i)^{1-y_i}}{\sum_{\text{all}\,\theta,p_i|\theta}\prod_{i=1}^{n}p_i^{y_i}(1-p_i)^{1-y_i}} 問題は、ベルヌーイ（p私、y私）ベルヌーイ（p私、y私）\text{Bernoulli}(p_i,y_i)確率の束の積である分子が小さすぎることです。（私のんんnは大きく、約1500です）。 したがって、すべての事後値はすべてθθ\theta0と計算されます（私はRで計算を行っています）。 明確にするために、各y私y私y_iは独自のp私p私p_i、これらのはn yのn要素のp私p私p_iベクトルを作成します。各θには、p iの独自のn要素ベクトルがあります。んんnんんn yyyθθ\thetaんんnp私p私p_iます。 編集：再現例の追加（分子用） p <- sample(seq(0,1,by=0.01), 1500, replace=T) y <- sample(c(0,1), 1500, replace=T) dbern(y, p) # 1500-element vector, each element is < 1 prod(dbern(y, p)) …