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

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

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

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

カテゴリー: R, stan

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

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

キングオブコント2016の評価

キングオブコント2016が終わりまして,案の定「あいつの評価がおかしい」とかいう話が出てきたわけですが,実際はどうだったかというのをデータから見てみようと。 スコアはネットにも転がっていたし,少ないので手打ちのデータにし … 続きを読む

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

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

Kandai.Rで発表してきました

6/18,第一回Kandai.Rで発表してきました。母校で行われる記念すべき第一回目に,初心者講習を引き受けさせていただきました。 Rとベイズは私の最後の恋です。Rを広めるためであれば,初心者講習をずっとやり続けることも … 続きを読む