岡田(2014)のベイズファクターの計算をStanでやってみた件。 これについては，次のスライドに詳しい。 ただ，途中のStanコードがちょっとわかりにくいのと，計算をstan内部でgenerated quantitie… 続きを読む 習作；ベイズファクター
カテゴリー: stan
Stanで体重の推移をみつめてみた(状態空間モデル)
この記事は Stan Advent Calendar 2016 のエントリーです。 Stanいいですね。いいですね，Stan. 私はもっぱらRから触るrstanユーザーなんですが，いろいろなところから触れる間口の広さがい… 続きを読む Stanで体重の推移をみつめてみた(状態空間モデル)
不倫報道の間隔を推定する
2016年は不倫報道が多いですね！全くどうでもいいんですが，これほど多いとデータとしてみなし，なんかの分析をしたくなる。 そこでデータをつくって分析してみます。 dataset <- as.data.frame( … 続きを読む 不倫報道の間隔を推定する
Hijiyama.R #5でプレゼンしてきました
こんかいのHijiyama.RはRをつかって文書作成，R markdown祭りという裏テーマがありまして，それに乗っかって話をしてきました。 内容は，心理統計の課題レポートを課すとき，乱数でデータ作るんですけど，思い通り… 続きを読む Hijiyama.R #5でプレゼンしてきました
習作；ラベルスイッチング問題をなんとかする関数
前回の続き。Stanコードに対数尤度を出す部分を追加。。。 (どう見てもカッコ悪い。スマートに書き直したいものだ。ていうか，そもそもこれで合ってるのかな？) [code] data{ int<lower =2>… 続きを読む 習作；ラベルスイッチング問題をなんとかする関数
習作；混合正規分布モデル(Gaussian Mixture Model) by rstan
Stanのマニュアルに載っている簡単なモデルなんだけど，最後の「各要素の所属確率」を算出するのになぜか手間取った。softmax関数を使うときの型の問題でした。ちぇ。 Stanコードはこの通り。generated qua… 続きを読む 習作；混合正規分布モデル(Gaussian Mixture Model) by rstan
塾合宿のおもひで
お盆休みを満喫しております。 お休みに入る前，8/8-8/11の二泊三日で魁！ベイズ塾夏合宿2016に参加してきました。今年から新塾生も加入し，総勢19名というソコソコの集団サイズ。残念ながら，最近の若手は… 続きを読む 塾合宿のおもひで
毒くらわば皿まで
rstan 2.10が来てますね。 2.9とちがって，式の代入が< -から=にすることが推奨されているようで，そうなっちゃうと今までのコードも全部書き換えないといけない。 とはいえこれから先はこうなるってんなら，早… 続きを読む 毒くらわば皿まで
読書会を始めます
学内で，読書会を始めます。 [amazonjs asin=”4254122144″ locale=”JP” title=”はじめての 統計データ分析 ―ベイズ的… 続きを読む 読書会を始めます
Lesson10.Bayesian Modeling 2
#2016.01.07現在，Rpubsが不調なため(?)，こちらにアップします。 #2016.01.22 Rpubsが復旧したので，そちらにリンクを貼るようにします。なお，前に載せてあった記事にはstanコードにバグがあ… 続きを読む Lesson10.Bayesian Modeling 2