Variance Function Regressionをmlexpで解く

Western and Bruceが示しているコードではなく，尤度関数を最大化する方法をとる．これはstataのmlexpを用いた方法である．このブログを参考にする．ただし，ローカルマクロに関してはwindowsで利用可能なように書き換えた． sysuse auto, clear local x =…

2016-05-15

Stataでラグを取る――Lを使う，[_n-1]は使わない

Stataのページにも[_n-1]を使えばよいと書いてあるが，パネルの場合，[_n-1]を使うと大変なことになるので注意． use http://www.stata-press.com/data/r14/nlswork . list idcode year ttl_exp in 1/12, clean idcode year ttl_exp 1. 1 70 1.083333 2. 1 7…

2016-03-16

日本社会心理学会　春の方法論セミナー　統計的因果推論への招待　メモ

大塚淳 Causation, Prediction, and Search 林岳彦『交尾行動の新しい理解』周辺尤度，リスクー社会心理学，記述・断りの追求・予想・介入の違い（背反ではない），内的・外的妥当性星野崇宏コンジョイント分析部分効用に分解因果効果それ単独の効果…

2016-03-07

太郎丸本 3章のデータをSASで

/*データ入力*/ data datach3; infile cards; input score; cards; 85 84 71 48 16 64 94 74 69 92 ; proc print;run;proc means data=datach3; run; Data datach3; set datach3; zscore=(score-69.7)/23.4143830; run; /*手で標準化*/ Data datach3;set dat…