習作;混合正規分布モデル(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とちがって,式の代入が< -から=にすることが推奨されているようで,そうなっちゃうと今までのコードも全部書き換えないといけない。 とはいえこれから先はこうなるってんなら,早 … “毒くらわば皿まで”の続きを読む

Lesson10.Bayesian Modeling 2

#2016.01.07現在,Rpubsが不調なため(?),こちらにアップします。 #2016.01.22 Rpubsが復旧したので,そちらにリンクを貼るようにします。なお,前に載せてあった記事にはstanコードにバグがあ … “Lesson10.Bayesian Modeling 2″の続きを読む

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

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

2群の平均の差について,あるいはGitHubデビューの件

基礎からのベイズ統計学: ハミルトニアンモンテカルロ法による実践的入門に打ち震えたという話から,学部生の統計の授業でもベイズを紹介したくなった。 ということで,とりあえずt検定をベイズ流にはこうやりますよ,というコードを … “2群の平均の差について,あるいはGitHubデビューの件”の続きを読む

MCMCでモコモコしちゃう

MCMCって大量にサンプルを発生させる技ですが,あー計算しているなー,発生させているなー,というのが目に見えてわかるようになったら楽しいと思うんです。イメージもしやすいだろうし。 ということで,Rでアニメーション画像を作 … “MCMCでモコモコしちゃう”の続きを読む

ベイズ塾で発表してきました

ベイズ塾,MCMCワークショップ第二弾で,RStanの導入について話してきました。 まあ午後の岡田先生のお話がメインなので,その前座です。 そのなかで,世界一簡単なrstanコードについて紹介したのですが,岡田先生の方か … “ベイズ塾で発表してきました”の続きを読む