Stata: 如何分組畫regression line

在作regression的時候，分組比較是非常常見的，有時候為了看出兩組或多組之間的差異，最簡單的方式就是用圖表呈現，如下圖所示。要怎麼畫這個圖表呢？Stata有很多方式可以作，這裡提供兩種。
1.png
方法一：
regress inc educ male, beta
predict incfnoi if male==0
predict incmnoi if male==1

twoway (connected incmnoi educ if male==1, lcolor(black) ///
lpattern(dot) msymbol(diamond) msize(large)) ///
(connected incfno educ if male ==0, lcolor(black) ///
lpattern(solid) msymbol(circle) msize(large)), ///
ytitle(Income in thousands) xtitle(Education) ///
legend(order(1 "Men" 2 "Women")) scheme(s2manual)
上面這段程式是從A gentle introduction to Stata 抄來的，這本書的相關介紹可見這篇文章：Stata書籍介紹。出來的圖就向上面的圖一樣。

方法二：

net describe postgr3, from(http://www.ats.ucla.edu/stat/stata/ado/analysis)
net install postgr3.pkg
net describe xi3, from(http://www.ats.ucla.edu/stat/stata/ado/analysis)
net install xi3.pkg
然後作下面的
xi3: regress inc educ male, beta
postgr3 educ, by(male) table
簡簡單單兩行指令就出來了，比上面簡單不少吧？咦？發現有什麼不對嗎？喔，下面的標題好醜，我不想要yhat_, male==0 出現，這對讀者一點意義有沒有。沒關係，點下圖紅色框線處，再作點修改即可。
2.png
先點你要修改的地方，然後到上面箭頭所指的地方將text修改成你想要的即可。
3.png
除了圖以外，剛剛的指令還產生table，可以確定一下有沒有問題。：
4.png
這方法是參考ucla所作的說明：http://www.ats.ucla.edu/stat/stata/ado/analysis/postgr3.htm
最後一種方法是用predxcat 和 predxcon (從findit上面可以找到)。不過我個人覺得這種還滿繁複的，我就沒特別實驗也不再這裡贅述。有興趣嘗試的人就自己試試看了！
至於如何畫interaction，就留待下次有機會再說。

標籤：