The One with ...

思いついたことや作業メモなどを公開しています

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

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

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

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

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

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

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); を実…

JAGSでWAICを計算するためのメモ

試しに計算してみたのですが,まだ細部を理解できていないので,検証用にメモを残しておきます. 豊田本2017(『実践ベイズモデリング』)にStanからWAICを計算する手順が解説されていたので,それを参考にして JAGSから作ったサンプルでWAICを計算してみま…

TikZ

TikZで簡単かつ綺麗な図がかけるので,これがあれば研究の分析・報告・文章化は全てTeX+beamerだけで事足りるでしょう. 使わないとすぐに忘れそうなので,ポイントだけメモしておきます. 簡易マニュアルでも推奨しているように,最初にグリッドを書いてお…

Texworks ショートカットキーのカスタマイズ

Texworksのショートカットキーの設定を編集すると,任意のTeXコマンド入力を登録できます. 例えばTeXソース編集中に 図 + tab と入力するだけで, \begin{figure}[h] \begin{center} \includegraphics[scale=0.5]{ .pdf} \caption{•} \label{american} \end…

Income Segregation

www.jstor.org この論文を読んでいて,Income Segregation と Race Segregation を組み合わせて表現するコードを思いついたので,メモを残しておきます. 以前から,シェリングのsegregation modelにエージェントの所得情報を別のレイヤーで表示させる,とい…

研究会用メモ

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

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

Dynamic Models of Segregation

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

SSD換装メモ

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

JAGSによる対数正規分布のパラメータ推定

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

再現可能な図

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

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

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

順番を決めるコード

新学期がはじまるので,演習報告担当者の順番決め用のあみだくじコード(Mathematica)を再掲しておきます. (* 動的表現 サイズ固定*) Manipulate[ Dynamic[Refresh[RandomSample[Range[n]], UpdateInterval -> x]], {{x, Infinity, "シャッフル"}, {0 -> "…

探索的分析の自動化実験

内容 データフレームを突っ込むと,データに含まれる全変数の中から1つを順番に応答変数として選び,残りの変数に回するような分析を自動的にやってくれるコードを書いてみた. 「本格的な分析の前に,変数間の関係を探索的・網羅的に知りたい」場合に使うこ…

mathjax実験

はてぶろがマークダウンに対応しているのを知ったとき,へえそう,くらいにしか思ってなかった. が... いやいや,とっくにMathjax使えるようになっとるやん. とか とか. まあ,Rmarkdownやそのものより使いづらいけど,少なくとも見た目はRpubsより綺麗や…

p値の分布を計算する

人工的に作成したデータを分析すると,統計の理解が捗ります.「p値ちっちゃい,やったー」みたいな意味のない分析するより,こっちのほうが統計初学者には勉強になるのではないか,と思います. 次のような例を考えてみましょう. 仮説検定を繰り返してp値…