在作regression的時候,分組比較是非常常見的,有時候為了看出兩組或多組之間的差異,最簡單的方式就是用圖表呈現,如下圖所示。要怎麼畫這個圖表呢?Stata有很多方式可以作,這裡提供兩種。
方法一:
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(<a href="http://www.ats.ucla.edu/stat/stata/ado/analysis">http://www.ats.ucla.edu/stat/stata/ado/analysis</a>)
net install postgr3.pkg
net describe xi3, from(<a href="http://www.ats.ucla.edu/stat/stata/ado/analysis">http://www.ats.ucla.edu/stat/stata/ado/analysis</a>)
net install xi3.pkg
然後作下面的
xi3: regress inc educ male, beta
postgr3 educ, by(male) table
簡簡單單兩行指令就出來了,比上面簡單不少吧?咦?發現有什麼不對嗎?喔,下面的標題好醜,我不想要yhat_, male==0 出現,這對讀者一點意義有沒有。沒關係,點下圖紅色框線處,再作點修改即可。
先點你要修改的地方,然後到上面箭頭所指的地方將text修改成你想要的即可。
除了圖以外,剛剛的指令還產生table,可以確定一下有沒有問題。:
這方法是參考ucla所作的說明:http://www.ats.ucla.edu/stat/stata/ado/analysis/postgr3.htm
最後一種方法是用predxcat 和 predxcon (從findit上面可以找到)。不過我個人覺得這種還滿繁複的,我就沒特別實驗也不再這裡贅述。有興趣嘗試的人就自己試試看了!
至於如何畫interaction,就留待下次有機會再說。