必要な人がいると思ったので、MCMCによる、平均値の分布範囲と、重回帰の切片と係数の推定範囲の計算のスクリプトを作っておきました。平均値の方はネットで拾ったものをパクったものです。重回帰はS先生に教えてもらいました。

両方とも動きます。解説はとりあえず意味が分かって使えるということを目指したものです。細かい解説は将来加えます。