R Advent Calendar 2013の11日担当記事です。 日付を一日勘違いしてしまい一日遅れになってしまいました。すみません。 テーマは、スプライン曲線で自由に設定した確率分布を使ってベイズ推定をやってみよう、というものです。ガウス分布もベータ分布も出てきません。え、共益分布?なにそれ、尤度関数にマッチする事前分布を設定すればパラメータの更新だけで事後分布出てきますよだって?ノンノン、本記事はそんな小難しいこと知らなくてもベイズ推定のキモは試せるぜ、という話です。 さて能書きはこれくらいにして本題。 R言語でスプライン曲線 飛び飛びの点を用意してsmooth.splineというの関数に与えると、スプライン補間した関数(に相当するオブジェクト)を返してくれます。それに任意のxの点を与えると対応するyの点を返してくれます。 例えば xy<- rbind( c(0,10), c(20,
![スプライン曲線で簡単お手軽ベイズ推定 - アドファイブ日記(ミラー版)](https://cdn-ak-scissors.b.st-hatena.com/image/square/24a473742f7135d7bb860fdf7dd76b5b629fa746/height=288;version=1;width=512/http%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fi%2Fisobe1978%2F20131213%2F20131213004626.png)