Agent-based modelのチュートリアル

MesaというPythonライブラリの中に所得分布の生成モデル（Dragulescu & Yakovenko 2000）を題材にしたチュートリアルが含まれている．なかなかよくできているのだが，学部生向けとしてはちょっとコードが複雑すぎるので授業用資料として簡略化してみた*1． …

2017-08-17

ggplot2による事後分布の視覚化

絶対忘れるので自分用にメモしておく．目的 Stanで推定した結果を分布のまま直接比較したい． ggplot2用のデータフレームの作り方流れはこんな感じ． Stanで推定した結果をfitに記録（CSVで記録すると結果的に遅かった．sample_file からのサンプル列の取…

2017-06-28

Stanによる推定例：ベルヌーイ分布のパラメータ

連勝記録から勝率を推定するある棋士が29連勝した，というデータから，勝利確率の事後分布をMCMCで計算してみる．統計モデルの仮定は次の通りである．毎回の勝ち負けは独立にベルヌーイ分布に従っているベルヌーイ分布のパラメータpの事前分布は範囲[0,1]…

2017-06-20

StanからWAICを計算するためのメモ

StanからWAICを計算するためのメモ豊田本（2017）を参考に，Stanの結果からWAICを計算する方法のメモ．テキスト198ページ記載のコードは略記のため，そのまま書き写して実行してもStanが動きません*1．公開されたコードを参考に確率分布のパラメータを適…

2017-06-13

TeXstudioショートカットキーのカスタマイズ方法

オリジナルのショートカットキーを設定する方法 TeXstudio便利ですね．私は長い間秀丸+マクロを使っていましたが，TeX Live導入後は TeXwokrsとTeXstudioを併用しています．さて，TeXstudioでオリジナルのショートカットキーを作成したいと思うことが時々…

2017-03-06

tikz-decoration

TikZのdecorationがうまくいかないので検索してみると，次のような解決策がweb上にあった．マニュアルに記載の \draw[decorate,decoration={coil,segment length=4pt}] (0,1) -- (3,1); \draw[decorate,decoration={coil,aspect=0}] (0,.5) -- (3,.5); を実…