The One with ...

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

2016-01-01から1年間の記事一覧

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値…