stanとRでMCMCの結果をいじりたい。同時確立を計算してみる。

stanでMCMCやって、2つ以上の確立をそれぞれ別々に計算していたとする。その時、それらが同時に成立する確立を計算する（連言命題の確立を計算）方法をメモ。
要するにそれぞれの確立の積の平均を計算してやる。下記のように…（なお、自分用メモなので表現は不適切かも）

mean(fit_data$p1*fit_data$p2*fit_data$p3)

ここでfit_data$p1とかp2とかp3とかはMCMCで計算した結果から抽出したパラメータ。