LISの可処分所得データからジニ係数を計算する必要があります。LISトレーニングドキュメントによると、これを行うためのStataコードは次のとおりです。
di "** INCOME DISTRIBUTION II – Exercise 13 **"
program define bottop
qui sum ey [w=hweight*d4]
replace ey = .01*r(mean) if ey<.01*r(mean)
qui sum dpi [w=hweight*d4], de
replace ey = (10*r(p50)/(d4^.5)) if dpi>10*r(p50)
end
foreach file in $us00h $fi00h {
display "`file'"
use hweight d4 dpi if (!mi(dpi) & !(dpi==0)) using "`file'", clear
gen ey=dpi/(d4^0.5)
bottop
ineqdeco ey [w=hweight*d4]
}
トレーニングドキュメントからこのコードをコピーして貼り付けただけです。スニペット
qui sum ey [w=hweight*d4]
replace ey=0.01*r(mean) if ey<0.01*r(mean)
と
qui sum dpi [w=hweight*d4], de
replace ey=(10*r(p50)/(d4^0.5)) if dpi>10*r(p50)
それぞれボトムコーディングとトップコーディングです。
このコードを実行しようとしたところ、変数hweightが見つかりませんでした。誰かがLISでhweightの新しい名前が何であるか知っていますか?または、他の方法でこの行き詰まりを克服する方法を誰かが提案できますか?
私はstataに精通していますが、このコードの洗練度は私の知識を超えています。
とても有難い。