習作;混合正規分布モデル(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を広めるためであれば,初心者講習をずっとやり続けることも … 続きを読む

RstudioとパッケージとGithubと

大事な午前中がまるまる潰れてしまった。 Rstudioでパッケージを作るのが楽になったのはわかる。RstudioでGithubのバージョンコントロールができることも知ってる。ちょっとできる。Githubで文章を加筆修正し … 続きを読む

因子軸の回転の可視化

「因子軸を回転って結局どういうイメージですか」 という質問をされたのですね。 探索的因子分析の後は,解釈をしやすくするために軸を回転しますよ,と説明しているのだけど,そのイメージがわかりにくいということなので,Rでプロッ … 続きを読む

初めてのスクレイピング

スクレイピングとぞいふもの へー,webからデータを取ってくることをスクレイピングっていうんだ,というレベルから始めたのですがね。 毎年,データ演習の授業では色々なデータを使うのだけど,分布が偏ってるのもあって,適度に関 … 続きを読む

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

先週末に行われたキングオブコント2015。リアルタイムで観て,昨日また録画したものを見直して。 いやー,なかなかいい戦いをしていますね。 ところで今回は評価方法が変わって,5人の審査員が100点ずつ持ち寄って500点満点 … 続きを読む