(Stata)在Marginal effect plot裡面放置 histogram plot的code

在Marginal effect plot裡面放置 histogram plot的code

1) loop

2) regression

3) margins and marginsplot

local mili militaryexpGDP arm_total 

foreach i of local mili{

    xtreg  `i' l.`i' c.v2x_polyarchy##c.lagGDPgrowth lngdppc wardummy lngdppc,fe

margins, at(v2x_polyarchy = (0.009 (0.06) 0.9) lagGDPgrowth=(-6.142) ) at( v2x_polyarchy = (0.009 (0.06) 0.9) lagGDPgrowth=(6.142) ) force

marginsplot,  xlabel(0.009 (0.06) 0.9) addplot(histogram v2x_polyarchy, freq yaxis(2) yscale(alt) fcolor(%15) lc(black%50)) name(gr143`i') title(Predictive margins for `i' with 95% Cls)

}

*freq: Stata的histogram plot預設值是density，指的是probability的density distribution。但是這個數值有時候不太好理解，所以建議更換成freq，樣本的次數(count)。

*  yaxis(2): 表示你希望有兩個y軸；一個是marginal plot的predicted value，另一個是histogram的density

* fcolor(%15): 由於這張圖的重點還是在marginal effect，我們只是添加上histogram plot來讓大家瞭解說我這個資料分佈的狀況。所以用fcolor()這個command來調降顏色。%15是把histogram 的顏色變得比較淺，這樣才不會阻礙大家辨別marginal effect。%100顏色最深、透明度為0，%0顏色最淺，近乎看不到、透明度為100%。如果你把透明度調到0%，記得加上線條的顏色，用lc(black%20)。

*讓stata output看起來像是用R做的: plotplain, ppt講解版