Stanで体重の推移をみつめてみた(状態空間モデル)

この記事は Stan Advent Calendar 2016 のエントリーです。 Stanいいですね。いいですね,Stan. 私はもっぱらRから触るrstanユーザーなんですが,いろいろなところから触れる間口の広さがい … 続きを読む

不倫報道の間隔を推定する

2016年は不倫報道が多いですね!全くどうでもいいんですが,これほど多いとデータとしてみなし,なんかの分析をしたくなる。 そこでデータをつくって分析してみます。 [crayon-59ca6689b865699777606 … 続きを読む

カテゴリー: R, stan

Hijiyama.R #5でプレゼンしてきました

こんかいのHijiyama.RはRをつかって文書作成,R markdown祭りという裏テーマがありまして,それに乗っかって話をしてきました。 内容は,心理統計の課題レポートを課すとき,乱数でデータ作るんですけど,思い通り … 続きを読む

習作;ラベルスイッチング問題をなんとかする関数

前回の続き。Stanコードに対数尤度を出す部分を追加。。。 (どう見てもカッコ悪い。スマートに書き直したいものだ。ていうか,そもそもこれで合ってるのかな?) で,鎖が複数になるとラベルスイッチング問題が生じるんだけど,こ … 続きを読む

カテゴリー: R, stan

習作;混合正規分布モデル(Gaussian Mixture Model) by rstan

Stanのマニュアルに載っている簡単なモデルなんだけど,最後の「各要素の所属確率」を算出するのになぜか手間取った。softmax関数を使うときの型の問題でした。ちぇ。 Stanコードはこの通り。generated qua … 続きを読む

カテゴリー: R, stan

塾合宿のおもひで

お盆休みを満喫しております。   お休みに入る前,8/8-8/11の二泊三日で魁!ベイズ塾夏合宿2016に参加してきました。今年から新塾生も加入し,総勢19名というソコソコの集団サイズ。残念ながら,最近の若手は … 続きを読む

カテゴリー: R, stan

毒くらわば皿まで

rstan 2.10が来てますね。 2.9とちがって,式の代入が< -から=にすることが推奨されているようで,そうなっちゃうと今までのコードも全部書き換えないといけない。 とはいえこれから先はこうなるってんなら,早 … 続きを読む

読書会を始めます

学内で,読書会を始めます。 これをみんなで読んでいきます。 Q&Aを読み終わって,自分も何か始めないではいられなかった。半ば強引に,ゼミ生を誘い込みました。 私の考えでは,ベイズ統計学になると,もっと情報仮説が増 … 続きを読む

ベイズ塾例会にて発表してきました

ベイズ塾のミーティングがありました。 一人一つ,担当の「分布」を決めて発表するというお話。私はコーシー分布を担当。たまたま先日読んだ本にコーシー分布のことが書いてあったので,話が膨らんでよかったわ。 資料を公開しておきま … 続きを読む