■どう進めようか悩んでいた作業を、bysortコマンドとSystem Variablesの_nを使って、階差をとればよいだけだと気がついた。分かってしまえば何ということはないのだが。


■ダミー変数Xを、基準カテゴリーを指定して入れたいとき(例えば3番目)、主効果だけの場合はib3.Xとすればよい。しかし、変数Yとの交互作用項を入れたいときの場合は、ib3.X*ib2.Yというように入れられないのがやや不便だと思った。X#Yで入れればよいわけだが、読み取るのに微妙に時間がかかる。


■来週くらいからやろうと思っている作業。

global xvars "x1 x2 x3 x4 x5"

logit yvar1 $xvars
predict pi, p
gen ipw = 1
replace ipw =1/pi if yvar1 == 1
replace ipw=1/(1-pi) if yvar1 == 0
sum ipw

regress yvar2 yvar1 [pweight=ipw]