案の定，2週間ほど前にやった計算をすっかり忘れて，同じコードを書くところであった． ブログにメモを残しておいたおかげで，思い出した． というわけで gmcから直接対数正規分布のパラメータを指定して，乱数を発生 そのデータからパラメータをMCMC推定す…

Refreshに関するメモ Mathematica内蔵関数Refreshの引数である，UpdateIntervalは更新を制御できたりできなかったりして，いまいち仕組みが分からない．そこでいろいろと実験してみた． Manipulate[ Refresh[Random[], UpdateInterval -> If[go, 0.105, Infi…

Dynamic Models of Segregation (Schelling 1971)をリアルタイムに表現するコードを更新したのでそのメモを残しておきます． 今回はNetLogo風に，エージェントの幸福度をリアルタイムに計算してプロットする部分を追加してみました． Mathematica ノートブッ…

ノートPC SSD換装メモ PCがモッサリしてきたけど，買い換えると再設定が面倒．そんなときOSごとクローンしてSSD換装するわけですが， 数年に1回しかやらないので，手順をすぐに忘れてしまう．そこで，次に作業するとき参照するために記録を残しておきます． …

Kruschke氏のblogでは,対数正規分布のデータを生成するのに，一度正規乱数作ってから基準化して指数化するというやや煩雑な手順を用いている． doingbayesiandataanalysis.blogspot.jp 最初にコードを見たとき，BUGSのクセとあいまって，何をやっているのか…

校正をしてて気づいた点をメモしておきます． 図は完全にreproducibleなコードを残しておくべきである． 例えば，Mathematicaで独自関数をplotする場合， パラメータの数値なども図と一緒に出力するようなコードを書いておく． Module[{plot1, point = 18, t…

前回の続き． そもそも分布を導出する理論モデルのパラメータを階層ベイズMCMCで推定できるのか，という確認をやってみる． まずデータとして，理論モデルからのサンプルをつくる． 今回はJAGSをRから呼び出すので, コードは全部Rで済むように書き換える． d…

ずっと気になっている（しかしマイナな）問題が解決しないので，そのメモを記録しておきます． 所得分布の生成モデルは，投資行動の反復から対数正規分布を導出する数理モデルである．仮定は次のとおりである（Hamada 2003; 2004; 2016）. 人々が最初に初期…