SAS基本圖形輸出

SAS/GRAPH模塊的特點：（1）是SAS系統(tǒng)的一部分；（2）提供了許多產生圖形的過程；（3）支持眾多的圖形設備；（4）提供對文本的廣泛控制，包括標題和腳注；（5）提供改善圖形的工具；（6）提供存儲和再現圖形的機制；（7）能夠重疊圖形；（8）能夠在每頁產生多個圖形。SAS/GRAPH模塊可以處理的圖形有：（1）散點圖及折線圖（plots）（2）圖表（charts）（3）地圖（maps）（4）三維圖形（three-dimensionalgraphics）（5）文本縮放及其他（textslides）§1SAS程序的過程步SAS的過程步可看作是已經編寫好的子程序，用戶可以在需要的時候調用它們。過程步的一般形式：PROC過程名<DATA=數據集名><選項>;

過程的專用語句描述;<VAR變量序列;><WHERE條件表達式;><BY變量序列;>

RUN;

說明：過程步都是以PROC語句開始，后面緊跟過程名，用以區(qū)分不同的過程步。尖括號中的內容是可選項，省略時對應的系統(tǒng)默認值如下：處理最近建立的數據集；處理所有變量（對于一個計算過程步處理所有數值變量）；處理數據集中的所有觀測?！?制作散點圖和連線圖散點圖——以數據集中某兩個變量為縱坐標變量和橫坐標變量，每條觀測對應于圖中的一個點，數據集中的多條觀測就在平面上構成一幅散點圖。連線圖——是將散點圖中的各點之間以一定的方式用直線或曲線相連所形成的圖形。散點圖和連線圖的作用：1）揭示一個變量隨另一個變量的變化2）顯示變量間的關系3）對不同的數據進行比較4）顯示數據值的分析實現：調用過程步GPLOT或PLOT。前者是在GRAPH窗口中以精確方式繪圖，后者是在OUTPUT窗口中以非精確方式繪圖。1．圖形選項設置一般形式：GOPTIONS圖形選項;功能：根據需要來進行圖形選項的設置，設置有關圖形輸出的各項參數。

說明：此操作一般在繪圖之前進行。

此語句屬于環(huán)境設置語句?；謴拖到y(tǒng)的缺省設置圖形選項：RESET=ALL|GOPTIONS|GLOBAL|……

說明：GOPTIONS：將所有全局語句（如TITLE等）的設置恢復為缺省狀態(tài)。

GOPTIONS：將所有關于圖形的設置恢復為缺省狀態(tài)。

ALL：同時實現上述兩種功能。

2．基本散點圖一般形式：

PROCGPLOT<DATA=數據集名>;PLOT縱坐標變量*橫坐標變量;RUN;

說明：

1）DATA=數據集名：此選項指定數據集的名稱，要求用二級文件名表示；

2）PLOT語句中指定縱坐標變量和橫坐標變量，要求必須是數值型變量；

3）每條觀測對應圖上的一個點。

例1：假設數據集dst.airqual中保存的是在一段時間內美國47個州的逐月的空氣質量監(jiān)測數據。其中包含下列變量：變量名類型

說明STATE字符型州名，用2個字母的縮寫表示REGION數值型州所在地區(qū)，用1~7的數字表示MONTH數值型月份，用1~12的數字表示AVE_TSP數值型懸浮顆粒物含量的平均值AVE_LEAD數值型鉛含量的平均值取紐約州的數據制作平均懸浮顆粒物含量隨月份變化的散點圖：libnamedst

'd:\dst';goptions

reset=all;proc

gplot

data=dst.airqual;

plotave_tsp*month;

wherestate='NY';run;注意：繪制高分辨率的統(tǒng)計圖非常耗費系統(tǒng)資源，因此當我們提交作圖語句時，SAS系統(tǒng)只是把它們送入程序緩沖區(qū)。只有當用戶打開GRAPH窗口時，程序才真正執(zhí)行，繪出統(tǒng)計圖。然后此過程步處于駐留狀態(tài)，系統(tǒng)可以繼續(xù)執(zhí)行所提交的語句，從而實現了用戶和系統(tǒng)之間的交互操作。結束駐留狀態(tài)的方法有：（1）提交另一個程序步；（2）提交QUIT語句。

3．散點圖和連線圖的修飾

點的符號、大小、顏色；點間連線、連線方式、線型、寬度；坐標軸：包括刻度、標簽、顏色等；圖例的配置與修飾。

3.1SYMBOL語句

一般形式：

SYMBOL<n>選項;

功能：用于控制圖中點和連線的顯示方式。說明：n是SYMBOL語句序號，取值范圍：1~99。系統(tǒng)缺省值為1。取消SYMBOL語句的設置

不加選項的SYMBOL語句SYMBOL<n>;

通過恢復系統(tǒng)缺省設置的圖形選項語句取消前面對所有SYMBOL語句的設置。

GOPTIONSRESET=SYMBOL;

性質

全局性：可以出現在過程步內或過程步之外，一旦定義則一直保持到重新定義或退出SAS系統(tǒng)；可加性：一個選項的定義不影響其他選項。

1）設置點

VALUE|V=符號：設置圖中點的表示符號

符號名稱符號表示PLUS+（缺省方式）XXSTAR*SQUARE□DIAMOND

TRIANGLE△HASH＃PAW::POINT.DOT●CIRCLE○NONE無符號COLOR|C=顏色：設置圖中點和連線的顏色。CV=顏色：設置圖中點的顏色。BLACK黑色（缺省方式）RED紅色GREEN綠色BLUE藍色YELLOW黃色CYAN青色MAGENTA洋紅PINK粉紅ORANGE橙色BROWN棕色GREY灰色HEIGHT|H=n<單位>：設置圖中點的大小。CELL單元，SAS默認單位CM厘米PCT屏幕畫面大小的百分數

PT象素IN英寸

POINTLABEL：在圖中每個點的旁邊標注縱坐標變量的值。例2：對于例1我們作如下修飾：goptions

reset=all;symbol

v=diamondcv=redh=1.5

pointlabel;proc

gplot

data=dst.airqual;

plotave_tsp*month;

wherestate='NY';run;2）設置連線

選項：INTEPOL|I=連線方式：設置相鄰點之間的連線方式。名稱連線方式NONE不連線（缺省方式）JOIN直線SPLINE光滑的樣條曲線NEEDLE通過點的并且垂直于橫坐標軸的直線SMnn繪制不經過散點的樣條曲線。nn曲線光滑性與擬合度的折衷，取值：00~99。RL<CLInn>線性回歸直線，CLInn：另外繪制預測值的置信區(qū)間曲線，nn：置信區(qū)間大小。RQ二次曲線RC三次曲線CI=顏色：設置連線的顏色。

LINE|L=n：設置連線的線型。n為線型號碼，取值范圍：0~46。其中常用的有：0（空）、1（實線，缺省值）、2（由點組成的虛線）……。WIDTH|W=n：設置連線的寬度。n為寬度值。系統(tǒng)缺省值為1。例3：對例1的修飾改為如下形式：goptions

reset=all;symbol

v=diamondh=2

cv=redi=joinci=bluel=2

w=2

pointlabel;proc

gplot

data=dst.airqual;

plotave_tsp*month;

wherestate='NY';run;

說明：連線是對數據（子）集中的每一對表示相鄰觀測的點之間連線的，所以在繪圖前要首先處理好觀測在數據集中的順序。最好是首先將數據集按橫坐標變量進行排序。例4：對sashelp.class作圖，顯示變量身高（Height）和體重（Weight）之間的關系。goptions

reset=all;symbol

v=squarei=join;proc

gplot

data=sashelp.class;

plotweight*height;run;對數據集sashelp.class按變量height排序后描散點圖。goptions

reset=all;symbol

v=squarei=join;Procsortdata=sashelp.class;Byheight;proc

gplot

data=sashelp.class;

plotweight*height;run;3.2設置標題和腳注

設置標題的一般格式：

TITLE<n><選項>‘字符串’<選項>‘字符串’……;設置腳注的一般格式：

FOOTNOTE<n><選項>‘字符串’<選項>‘字符串’……;

功能：設置在圖上顯示的標題和腳注。

選項：

FONT|F=字體：設置字符串的字體。COLOR|C

=

顏色：設置字符串的顏色。HEIGHT|H

=

n<單位>：設置字符的高度。說明：

1）TITLE和FOOTNOTE語句都是全局語句，性質與SYMBOL語句相同。

2）可選項<n>標明了TITLE或FOOTNOTE語句的序號，n的取值范圍1~99，缺省值為1。

3）希望在標題或腳注位置顯示的信息在TITLE或FOOTNOTE語句中用字符串的形式給出。在字符串前面的選項是對字符串顯示性質的控制，如果不給出，則系統(tǒng)以缺省方式顯示字符串。4）SAS系統(tǒng)提供的2類字體：硬字庫：由操作系統(tǒng)提供的字庫；軟字庫：由SAS系統(tǒng)提供的字庫。注意：軟字庫的字體名直接輸入即可，而硬字庫的字體名兩端需加引號。

獲得想用的字體名：1）在命令框中輸入FONTLIST并提交；2）在屏幕上出現的SelectFont（選擇字體）對話框中選擇；

軟字庫的所有字體點擊可獲得硬字庫中的字體名稱例5：在例3中加入title語句：goptions

reset=all;symbol

v=diamondh=2

cv=redi=joinci=bluew=2

pointlabel;titlef='隸書'c=green'紐約州每月的懸浮顆粒物平均值';title2h=1.2

'2002';proc

gplot

data=dst.airqual;

plotave_tsp*month;

wherestate='NY';run;

3.3設置坐標軸

AXIS語句的一般形式：

AXIS<n>選項;

功能：設置坐標軸的顯示形式，包括：坐標軸在畫面中的位置和長度，變量值在軸上的出現次序、尺度、刻度的標注，相關文字的字體、顏色、布局。

說明：可選項<n>標明了AXIS語句的序號，n的取值范圍1~99，缺省值為1。AXIS語句也是全局語句。取消設置：

要取消第n號坐標軸的設置

AXIS<n>;要取消所有坐標軸的設置

GOPTIONSRESET=AXIS;

選項：COLOR|C=顏色：設置坐標軸及相關文字的顏色。LABEL=(字符屬性‘字符串’)|NONE：設置軸的標簽及其屬性。說明：軸的標簽就是軸的說明信息，可在括號中的字符串給出。系統(tǒng)缺省設置是顯示對應變量的標簽或名字。VALUE=(字符屬性‘字符串’)|NONE：設置軸的主刻度的描述字符及其屬性。

ORDER=(值列表)：設置主刻度位置依次出現的數值。MINNOR|N=n|NONE：設置在主刻度之間加入的次刻度數n。標簽和主刻度的字符屬性：C=顏色：設置字符顏色；F=

字體：設置字符字體；H=n<單位>：設置字符高度；R=

角度：設置單個字符的旋轉角度，單位是度，逆時針方向為正。A=角度：設置整個字符串的旋轉角度。J=LEFT|CENTER|RIGHT：設置字符串的對齊方式。例6：設置如下坐標軸：axis1

label=(f=complexc=blueh=3pct)c=magentaw=3

minor=none;axis2

label=(a=-90

r=90

f=complexc=blueh=3pct'PartPerMillion')c=magentaw=3;3.4PLOT語句中的選項

PLOT語句的一般形式：

PLOT縱坐標變量*橫坐標變量

/選項;

選項：FRAME|NOFRAME：設置圖中是否顯示邊框。AUTOHREF|AUTOVREF：自動在圖中添加經過主刻度的水平/垂直參考線。NOAXIS：取消坐標軸以及與坐標軸相關的圖形元素。CAXIS=

顏色：設置坐標軸的顏色。CTEXT=

顏色：設置坐標軸旁字符的顏色。HAXIS|VAXIS=AXIS<n>|值列舉：設置圖中的橫坐標軸或縱坐標軸。

例6：對于例1中的數據集，繪制如下連線圖：goptions

reset=all;symbol1

v=triangleh=1.5

i=joinc=greenw=2;titlec=blue'NewYorkSuspendedParticleAverage';title2c=blueh=1.2

'2002';axis1

label=(f=complexc=blueh=3pct)c=magentawidth=3

minor=none;axis2

label=(a=-90

r=90

f=complexc=blueh=3pct'PartsPerMillion')c=magentawidth=3;proc

gplot

data=dst.airqual;

plotave_tsp*month/haxis=axis1vaxis=axis2noframe;

wherestate='NY';run;4．多條連線的圖形繪制

4.1多幅圖形的繪制

在一個過程步中繪制多幅圖形，可以采用以下幾種方法：（1）使用BY語句對數據集中的觀測分組；（2）在PLOT語句中給出多個縱橫軸變量的組合；；（3）在過程步中使用多個PLOT語句。

說明：在PROCGPLOT過程中給出可選項UNIFORM可使不同圖中同一變量的坐標軸尺度相同。

例7：分別繪制紐約州每月平均懸浮顆粒物和含鉛量的連線圖。goptions

reset=all;proc

gplot

data=dst.airqual;

wherestate='NY';

plotave_tsp*monthave_lead*month/caxis=bluectext=blue;

symbol1

v=trianglei=join;

titlef=swissic=red'NewYorkSuspendedParticleAverage';run;

4.2將多條連線圖繪制在同一幅畫面上

1）在PLOT語句中給出選項OVERLAY；

2）在同一個過程步中使用語句PLOT、PLOT2繪制多條連線；

3）按第三個變量的值對數據集進行分組作圖。1）選項OVERLAY例8：將紐約州每月平均懸浮顆粒物和含鉛量的連線繪制在同一幅圖中。goptions

reset=all;proc

gplot

data=dst.airqual;

wherestate='NY';

plotave_tsp*monthave_lead*month/overlaylegendhminor=0

vaxis=axis1;

symbol1

v=diamondi=joinc=bluel=1

w=3;

symbol2

v=stari=joinc=redl=2

w=3;

axis1

label=(a=-90

r=90

f=swiss

'PARTSPERMILLION');

titlef=swissic=magentah=6pct'DelawareMonthlyPollutantAverages';run;說明：1）在同一幅畫面中坐標軸尺度相同，因此無需UNIFORM選項；

2）為了區(qū)分圖中的不同連線，每條連線的顯示形式應該不同。我們可以通過SYMBOL語句設置連線的參數。系統(tǒng)默認每條SYMBOL語句設置的連線依次使用一次。

3）PLOT語句中的選項LEGEND的功能是在圖中加入圖例以說明每條連線的含義。

2）PLOT2語句一般形式：

PLOT2縱坐標變量*橫坐標變量

……/<選項>;功能：為縱坐標變量在圖的右側另外設置一個縱坐標軸，它與用PLOT語句在圖的左側設置的縱坐標軸可以用不同的AXIS語句設置。例9：修改例8，為變量ave_lead在圖的右側建立一個縱坐標軸。goptions

reset=all;proc

gplot

data=dst.airqual;

wherestate='NY';

plotave_tsp*month/overlaylegendhminor=0;

plot2ave_lead*month/overlaylegend;

symbol1

v=diamondi=joinc=bluel=1

w=3

h=1.5;

symbol2

v=stari=joinc=redl=2

w=3

h=1.5;

titlef=zapfc=blueh=6pct'DelawareMonthlyPollutantAverages';run;3）按照第三個變量分類作圖

一般形式：PLOT縱坐標變量

*橫坐標變量

=第三變量

</選項>;

功能：繪圖時將數據集中的觀測按第三個變量的值分類（第三變量值相等的觀測分為一類），然后分別繪制不同的連線并置于同一幅圖中。

例10：將紐約、特拉威和新澤西三個州的逐月平均懸浮顆粒物的連線圖畫在同一幅畫面中。goptions

reset=all;proc

gplot

data=dst.airqual;

wherestatein('NY'

'DE'

'NJ');

plotave_tsp*month=state/hminor=0;

symbol1

v=diamondi=joinc=bluel=1

w=2

h=1.5;

symbol2

v=stari=joinc=redl=2

w=2

h=1.5;

symbol3

v=squarei=joinc=greenl=3

w=2

h=1.5;

titlef=zapfc=blueh=5pct'NewYork,Delaware,NewJerseyTSPAverage';run;說明：（1）在按第三個變量分類作圖時，系統(tǒng)缺省的提供圖例。（2）由于新澤西州后幾個月的數據沒有提供，所以在圖中只顯示了前幾個月的連線。

4.3LEGEND語句

主要用于對圖例進行設置的全局語句。一般形式：

LEGENDn

選項;

功能：設置圖例。說明：n為LEGEND語句的序號，取值范圍1~99，缺省值為1。

恢復所有LEGEND語句的缺省設置：

GOPTIONSRESET=LEGEND;

在圖中顯示我們已設置的圖例：

PLOT縱坐標變量*橫坐標變量/LEGEND=LEGENDn;選項：

ACROSS=n：將圖例元素排為包含n個元素的行；

DOWN=n：將圖例元素排為包含n個元素的列；POSITION=(bottom|middle|topleft|center|rightoutside|inside)：設置圖例在畫面上的位置，畫橫線的是缺省值；

FRAME：設置圖例區(qū)加邊框；CFRAME=顏色：設置圖例區(qū)邊框的顏色；CSHADOW=顏色：設置圖例區(qū)邊框陰影的顏色；LABEL=(字符屬性‘字符串’)：設置圖例標簽；VALUE=(字符屬性1‘字符串1’……)：設置圖例中的每條連線的說明標簽。例11：修改例10中的圖例。goptions

reset=all;proc

gplot

data=dst.airqual;

wherestatein('NY'

'DE'

'NJ');

plotave_tsp*month=state/hminor=0

legend=legend1;

symbol1

v=diamondi=joinc=bluel=1

w=2

h=1.5;

symbol2

v=stari=joinc=redl=2

w=2

h=1.5;

symbol3

v=squarei=joinc=greenl=3

w=2

h=1.5;

legend1

down=3

position=(toprightinside)frame

cshadow=blacklabel=(f=duplexh=1.5)value=(f=duplex);

titlef=zapfc=blueh=5pct'NewYork,Delaware,NewJerseyTSPAverage';run;§2制作柱狀圖和餅圖SAS系統(tǒng)提供GCHART過程繪制二維或三維的柱狀圖或餅圖，用于表現某種分類下某個變量的統(tǒng)計量的匯總信息。

1．GCHART過程步

1.1一般形式PROCGCHARTDATA=數據集名;

圖形名

作圖變量

</選項>;

RUN;1.2

定義要素

1）定義圖表的形式；

2）定義一個作圖變量（ChratVariable）用于決定圖表的自變量，也稱分類變量；

3）定義一個分析變量用于計算各種統(tǒng)計量。1.3

圖表的形式

圖形形式圖形名自變量的值對應的圖形元素分析變量的值水平柱狀圖HBAR一個柱柱長三維水平柱狀圖HBAR3D一個柱柱長垂直柱狀圖VBAR一個柱柱高三維垂直柱狀圖VBAR3D一個柱柱高三維平鋪柱狀圖BLOCK一個柱柱高餅圖PIE餅的一角扇形的角度三維餅圖PIE3D餅的一角扇形的角度星形圖STAR一條射線射線長度2.基本匯總圖

系統(tǒng)根據作圖（分類）變量的值將觀測分類，再用柱狀圖或餅圖展示匯總結果。說明：1）作圖變量：字符型：變量的一個值是一類；數值型：系統(tǒng)認為值是連續(xù)的，會自動將值的變化范圍分為若干等長區(qū)間。并以區(qū)間的中點作標簽。2）不指明分析變量，系統(tǒng)默認的統(tǒng)計量是頻數（Freq）；若指明分析變量，系統(tǒng)默認的統(tǒng)計量是總和（Sum）。3）可添加TITLE、FOOTNOTE等語句修飾圖。例：假設數據集dst.airqual中保存的是在一段時間內美國47個州的逐月的空氣質量監(jiān)測數據。其中包含下列變量：變量名類型

說明STATE字符型州名，用2個字母的縮寫表示REGION數值型州所在地區(qū)，用1~7的數字表示MONTH數值型月份，用1~12的數字表示AVE_TSP數值型懸浮顆粒物含量的平均值AVE_LEAD數值型鉛含量的平均值例1：用垂直柱狀圖顯示不同區(qū)域的監(jiān)測次數goptions

reset=all;proc

gchart

data=dst.airqual;

whereregion=1;

vbarstate;titlec=redh=5pctf=swissl

'RegionalContrctFrequencies';run;例2：用數值變量作分類變量，以垂直柱狀圖顯示紐約州不同含鉛量的監(jiān)測次數goptions

reset=all;proc

gchart

data=dst.airqual;

wherestate='NY';

vbarave_lead;

titlec=redh=5pctf=swissl

'RegionalContrctFrequencies';run;3.與分類有關的選項說明：在CHART語句的選項中設置。LEVELS=n：對數值型變量設置其分類數n；MIDPOINTS=值列：對字符型分類變量依次列出每類的值，對數值型分類變量依次列出每個區(qū)間的中點；格式：值列形式：值1值2……初值to終值by步長DISCRETE：對于每個數值型分類變量取值作為一類，適用于離散型數值變量；ASCENDING|DESCENDING：按統(tǒng)計量值的升序或降序排列各柱或餅的各角。例：在下面的過程步中分別以下列vbar語句代替原有的vbar語句，可以看到不同效果。goptions

reset=all;proc

gchart

data=dst.airqual;vbarave_lead;run;vbarave_lead/levels=2;vbarmonth/discrete;vbarmonth/midpoints=1to12by3;vbarstate/midpoints='NY'

'DE'

'NJ';4.選擇分類變量和統(tǒng)計量說明：在CHART語句的選項中設置。SUMVAR=變量名：設置分析變量；TYPE=FREQ|CFREQ|PERCENT|CPERCENT|MEAN|SUM：設置統(tǒng)計量。FREQ：頻數；CFREQ：累積頻數；PERCENT：百分數；CPERCENT：累積百分數；MEAN

：均值；SUM：總和。例：以水平柱狀圖顯示各地區(qū)懸浮顆粒物含量的平均值的平均值goptions

reset=all;proc

gchart

data=dst.airqual;

hbarregion/sumvar=ave_tsp

type=mean;run;5.顯示統(tǒng)計量說明：在CHART語句的選項中設置。NOSTAT：在圖中不顯示統(tǒng)計量；OUTSIDE=統(tǒng)計量：在柱的外部顯示；INSIDE=統(tǒng)計量：在柱的內部顯示；例：以三維垂直柱狀圖顯示各地區(qū)懸浮顆粒物含量的平均值的平均值，并標注頻數和累計頻數。proc

gchart

data=dst.airqual;

vbar3dregion/sumvar=ave_tsp

type=mean

inside=freqoutside=cfreq;run;6.柱狀圖的修飾6.1在CHART語句的選項中設置。1）坐標軸NOAXIS：不顯示與坐標軸有關的內容；CAXIS=顏色：設置坐標軸的顏色；CTEXT=顏色：設置軸旁文字的顏色；MAXIS=AXISn：設置分類變量的坐標軸；RAXIS=AXISn：設置統(tǒng)計量的坐標軸；MINOR=n：統(tǒng)計量坐標軸的主刻度之間加入的次刻度數；AUTOREF：自動顯示參照線，垂直于統(tǒng)計量坐標軸。2）邊框NOFRAME：不顯示邊框；CFRAME=顏色：設置邊框的顏色。3）柱寬WIDTH=n：設置柱寬（單位：CELL）SPACE=n：設置柱間寬度（單位：CELL）6.2