第四章 其他圖形繪制

第四章其他圖形繪制4.1繪制曲線標(biāo)繪圖和連線標(biāo)繪圖4.1.1繪制曲線標(biāo)繪圖曲線標(biāo)繪圖的命令語(yǔ)法如下：[twoway]linevarlist[if][in][,options]其中，twoway可以省略，varlist包括了參加繪圖的變量名稱，選項(xiàng)options包括了connect_options、axis_choice_options和各種其他的twoway_options，例如控制標(biāo)題、圖例、坐標(biāo)軸、增加線的選項(xiàng)，所有這些選項(xiàng)都在3.1和3.2節(jié)中詳細(xì)介紹了?！纠?-1】使用系統(tǒng)自帶的數(shù)據(jù)集auto.dta中的變量mpg和weight，通過(guò)散點(diǎn)圖和曲線圖的復(fù)合來(lái)繪制mpg對(duì)weight一元線性回歸的95%的預(yù)測(cè)置信區(qū)間?！纠?-2】使用系統(tǒng)數(shù)據(jù)集uslifeexp.dta，該數(shù)據(jù)集已經(jīng)多次使用，這里不再重復(fù)介紹。要求讀者使用這個(gè)數(shù)據(jù)集完成如下要求的繪圖：繪制美國(guó)全體人口預(yù)期壽命隨時(shí)間變化的曲線標(biāo)繪圖。要求在同一圖形中繪制如下四條曲線：美國(guó)男性預(yù)期壽命隨時(shí)間變化的曲線標(biāo)繪圖美國(guó)女性預(yù)期壽命隨時(shí)間變化的曲線標(biāo)繪圖美國(guó)男女預(yù)期壽命差異隨時(shí)間變化的曲線標(biāo)繪圖美國(guó)男女預(yù)期壽命差異對(duì)時(shí)間的一次擬合線4.1.2繪制連線標(biāo)繪圖

在曲線標(biāo)繪圖中，數(shù)據(jù)點(diǎn)是看不見(jiàn)的，只能看到連線，但是下面命令的繪制結(jié)果，圖中的數(shù)據(jù)點(diǎn)是由散點(diǎn)圖加號(hào)加以標(biāo)記的，這就是連線標(biāo)繪圖，其命令如下：twowayconnectedvarlist[if][in][weight][,scatter_options]這個(gè)命令的選項(xiàng)與散點(diǎn)圖繪制命令scatter完全一致，這里不再贅述。需要說(shuō)明的是連線標(biāo)繪圖與散點(diǎn)圖本質(zhì)上是一樣的，下面兩條命令的結(jié)果是一直的：.twowayconnectedyvarxvar.scatteryvarxvar,connect(l)同樣的，下面兩條命令的結(jié)果是也是一致的：.scatteryvarxvar.twowayconnectedyvarxvar,connect(none)這說(shuō)明連線標(biāo)繪圖完全可以和默認(rèn)情況下的散點(diǎn)圖相互轉(zhuǎn)換。4.2繪制擬合圖形

4.2.1繪制一次擬合圖形一次擬合圖形的繪制實(shí)際上分兩步，首先Stata使用yvar為因變量，xvar為自變量進(jìn)行一元線性回歸，然后得到y(tǒng)var的擬合值比如說(shuō)是hat，然后使用hat對(duì)xvar做曲線標(biāo)繪圖，同時(shí)復(fù)合原始數(shù)據(jù)的散點(diǎn)圖。一次擬合的命令如下：twowaylfityvarxvar[if][in][weight][,options]其中，yvar是作為y軸變量的變量名，xvar是作為x軸變量的變量名【例4-3】使用數(shù)據(jù)集auto.dta中的mpg和weight兩個(gè)變量，繪制如下要求的圖形；（1）繪制復(fù)合圖形，包括mpg和weight的散點(diǎn)圖以及mpg和weight的一次擬合圖。（2）根據(jù)分類變量foreign繪制復(fù)合圖形，包括mpg和weight的散點(diǎn)圖以及mpg和weight的一次擬合圖。4.2.2繪制二次擬合圖形二次擬合圖形的繪制實(shí)際上分兩步，首先Stata使用yvar為因變量，xvar和xvar的平方為自變量進(jìn)行二元線性回歸，然后得到y(tǒng)var的擬合值比如說(shuō)是hat，然后使用hat對(duì)xvar做曲線標(biāo)繪圖，同時(shí)復(fù)合原始數(shù)據(jù)的散點(diǎn)圖。二次擬合的命令如下：

twowayqfityvarxvar[if][in][weight][,options]其中，yvar是作為y軸變量的變量名，xvar是作為x軸變量的變量名?！纠?-4】使用數(shù)據(jù)集auto.dta中的mpg和weight兩個(gè)變量，繪制如下要求的圖形；（1）繪制復(fù)合圖形，包括mpg和weight的散點(diǎn)圖以及mpg和weight的二次擬合圖。（2）根據(jù)分類變量foreign繪制復(fù)合圖形，包括mpg和weight的散點(diǎn)圖以及mpg和weight的二次擬合圖。4.2.3繪制lowess擬合圖形graphtwowaylowess具有簡(jiǎn)潔的優(yōu)點(diǎn)，并遵循graphtwoway這一族命令的習(xí)慣語(yǔ)法。lowess命令的語(yǔ)法如下：twowaylowessyvarxvar[if][in][,options]【例4-5】使用數(shù)據(jù)集auto.dta中的mpg和weight兩個(gè)變量，繪制如下要求的圖形；（1）繪制復(fù)合圖形，包括mpg和weight的散點(diǎn)圖、mpg和weight的一次擬合圖以及l(fā)owess修勻擬合圖。（2）根據(jù)分類變量foreign繪制復(fù)合圖形，包括mpg和weight的散點(diǎn)圖、mpg和weight的一次擬合圖以及l(fā)owess修勻擬合圖。4.3繪制條形圖雖然條形圖的圖形相對(duì)比較容易閱讀，但是它所包含的信息缺失非常龐雜的。它可以顯示總多的描述性統(tǒng)計(jì)量，同時(shí)，同一個(gè)條形圖中可以顯示多個(gè)變量的統(tǒng)計(jì)量。條形圖顯示的描述性統(tǒng)計(jì)量見(jiàn)表4-4。表4-4條形圖顯示的描述性統(tǒng)計(jì)量統(tǒng)計(jì)量Stata寫(xiě)法均值mean中位數(shù)median和sum計(jì)數(shù)count標(biāo)準(zhǔn)差sd最大值max最小值min各種百分位數(shù)p1……p99Stata可以繪制的條形圖種類繁多，既顯示了條形圖的廣泛用途，也說(shuō)明了Stata制圖功能的強(qiáng)大。Stata繪制條形圖的命令是：graphbaryvars[if][in][weight][,options]graphhbaryvars[if][in][weight][,options]前者用于繪制縱向的條形圖，在縱向條形圖中，y軸是數(shù)值型的，x軸是分類變量。。后者用于繪制橫向條形圖。表4-5條形圖選項(xiàng)分類選項(xiàng)（options）含義group_options用于設(shè)定分類的組別yvar_options用于設(shè)定繪制條形圖的變量lookofbar_options設(shè)定條形圖的外觀legending_options設(shè)定條形圖的圖例axis_options設(shè)定y軸的刻度標(biāo)識(shí)title_and_other_options設(shè)置添加標(biāo)題等【例4-6】數(shù)據(jù)集citytemp.dta是Stata系統(tǒng)自帶的數(shù)據(jù)集，該數(shù)據(jù)集共有6個(gè)變量和956個(gè)觀測(cè)值，其中主要變量如表4-13所示，部分?jǐn)?shù)據(jù)羅列在了表4-14中：表4-13數(shù)據(jù)集citytemp.dta變量情況變量名存儲(chǔ)類型顯示類型變量含義divisionint%8.0g調(diào)查區(qū)域劃分regionint%8.0g地理區(qū)域heatddint%8.0g最熱天數(shù)coolddint%8.0g最冷的天數(shù)tempjanfloat%9.0g一月份氣溫tempjulyfloat%9.0g七月份氣溫表4-14數(shù)據(jù)集citytemp.dta（部分?jǐn)?shù)據(jù)）regioncountrypopgrowthlexpgnppcsafewaterN.A.Canada1.2791917099N.A.Cuba0.77693N.A.DominicanRepublic2.171177071N.A.ElSalvador1.569185055N.A.Guatemala2.664164067N.A.Haiti25441028N.A.Honduras36974065N.A.Jamaica175174070N.A.Mexico1.972384083N.A.Nicaragua2.868189681N.A.Panama1.974299084N.A.PuertoRico17697N.A.TrinidadandTobago