試驗(yàn)設(shè)計(jì)與統(tǒng)計(jì)分析SAS實(shí)踐教程課件：SAS統(tǒng)計(jì)繪圖

SAS統(tǒng)計(jì)繪圖4.1用SAS繪制散點(diǎn)圖4.2用SAS繪制折線圖4.3用SAS繪制盒須圖4.4用SAS繪制曲線圖4.5用SAS繪制餅圖4.6用SAS繪制柱形圖4.7用SAS繪制直方圖4.8用SAS繪制誤差圖4.9用SAS繪制三維網(wǎng)格圖4.10用SAS繪制等值線圖上機(jī)報(bào)告

上機(jī)目的

掌握SAS繪統(tǒng)計(jì)圖的編程方法，熟悉程序中的過程、過程選項(xiàng)、語句、語句選項(xiàng)等編程要素和編程格式。學(xué)會(huì)用統(tǒng)計(jì)圖展現(xiàn)試驗(yàn)數(shù)據(jù)、擬合曲線、函數(shù)圖形、頻數(shù)和累積頻數(shù)分布、百分率和累積百分率分布、響應(yīng)誤差、響應(yīng)面和等值線。

上機(jī)內(nèi)容①采用gplot過程編寫SAS程序，繪散點(diǎn)圖、折線圖和曲線圖；②采用gchart過程編寫SAS程序，繪餅圖、柱形圖、直方圖和誤差圖；③采用g3d過程編寫SAS程序，繪三維網(wǎng)格圖；④采用gcontour過程編寫SAS程序，繪等值線圖。

4.1用SAS繪制散點(diǎn)圖

在執(zhí)行SAS繪圖程序之前，首先需創(chuàng)建或選定用于繪圖的SAS數(shù)據(jù)表。

【例4-1】在某地連續(xù)觀測(cè)13個(gè)年份(Order)的冬季積雪時(shí)間(X1)、化雪日期(X2)和二化螟第一代成蟲發(fā)生量(Y)，獲得的數(shù)據(jù)樣本如表4-1所示。試創(chuàng)建它的SAS數(shù)據(jù)表sasuser.chengchong。

(1)用表4-1數(shù)據(jù)創(chuàng)建Excel數(shù)據(jù)表chengchong.xls。

(2)將Excel數(shù)據(jù)表chengchong.xls導(dǎo)入SAS，創(chuàng)建SAS數(shù)據(jù)表sasuser.chengchong，以備后面的SAS繪圖程序調(diào)用。4.1.1單變量散點(diǎn)圖

單變量散點(diǎn)圖，一般按時(shí)間順序或指定順序展現(xiàn)變量觀測(cè)值的動(dòng)態(tài)變化。

【例4-2】試編寫SAS程序調(diào)用數(shù)據(jù)表sasuser.chengchong(表4-1)繪散點(diǎn)圖，展現(xiàn)二化螟第一代成蟲發(fā)生量(Y)的年份(Order)動(dòng)態(tài)。

(1)采用gplot過程編寫繪制散點(diǎn)圖的SAS程序如下：

goptionsreset=allftext=swisshtext=1.65; /*設(shè)置文本的字體和尺寸*/

symbolV=starH=1.75CV=B;/*設(shè)置數(shù)據(jù)點(diǎn)標(biāo)記的類型、尺寸和顏色*/

axis1label=(f='宋體''年份編號(hào)');/*設(shè)置橫軸標(biāo)簽的字體和內(nèi)容*/

axis2label=(A=90f='宋體''成蟲發(fā)生量'); /*設(shè)置縱軸標(biāo)簽的角度、字體和內(nèi)容*/

procgplotdata=sasuser.chengchong;/*指定gplot過程調(diào)用的數(shù)據(jù)表*/

plotY*Order/noframehaxis=axis1vaxis=axis2;

/*選項(xiàng)去圖框、引用縱橫軸的設(shè)置*/

run;quit;(2)程序的輸出結(jié)果如圖4-1所示。圖4-1二化螟成蟲發(fā)生量的年份動(dòng)態(tài)4.1.2兩變量散點(diǎn)圖

兩變量散點(diǎn)圖，一般用于展現(xiàn)一個(gè)變量與另一個(gè)變量的相關(guān)關(guān)系。

【例4-3】試編寫SAS程序調(diào)用數(shù)據(jù)表sasuser.chengchong繪制散點(diǎn)圖，展現(xiàn)二化螟第一代成蟲發(fā)生量(Y)與積雪時(shí)間(X1)的關(guān)系。(1)采用gplot過程編寫繪制散點(diǎn)圖的SAS程序如下：

goptionsreset=allftext=swisshtext=1.65;

symbolV=starH=1.75CV=B;

axis1label=(f='宋體''積雪時(shí)間(周)');

axis2label=(A=90f='宋體''成蟲發(fā)生量');

procgplotdata=sasuser.chengchong;

plotY*X1/noframehaxis=axis1vaxis=axis2;

run;quit;(2)程序的輸出結(jié)果如圖4-2所示。圖4-2二化螟成蟲發(fā)生量與積雪時(shí)間的關(guān)系4.1.3多變量散點(diǎn)圖

多變量散點(diǎn)圖，可展現(xiàn)多個(gè)變量的時(shí)間動(dòng)態(tài)，或多個(gè)變量與同一變量的相關(guān)關(guān)系。

【例4-4】試編寫SAS程序調(diào)用數(shù)據(jù)表sasuser.chengchong繪制散點(diǎn)圖，展現(xiàn)二化螟第一代成蟲發(fā)生量(Y)、積雪時(shí)間(X1)、化雪日期(X2)的年份動(dòng)態(tài)。(1)采用gplot過程編寫繪制散點(diǎn)圖的SAS程序如下：

goptionsreset=allftext=swisshtext=1.55;

symbol1V=starH=1.75CV=black;

symbol2V=squareH=1.75CV=B;

symbol3V=hashH=1.75CV=R;

axis1label=(f='宋體''年份編號(hào)');

axis2label=(A=90f='宋體''*成蟲發(fā)生量□積雪時(shí)間

?；┤掌?);

procgplotdata=sasuser.chengchong;

plotY*OrderX1*OrderX2*Order /noframeoverlayhaxis=axis1vaxis=axis2;

run;quit; /*選項(xiàng)overlay使多個(gè)散點(diǎn)圖疊加*/(2)程序的輸出結(jié)果如圖4-3所示。圖4-3成蟲發(fā)生量、積雪時(shí)間、化雪日期的年份動(dòng)態(tài)

4.2用SAS繪制折線圖

在執(zhí)行SAS繪圖程序之前，首先需創(chuàng)建或選定用于繪圖的SAS數(shù)據(jù)表。

【例4-5】抽樣觀測(cè)12個(gè)高營養(yǎng)玉米雜交種(variety)的產(chǎn)量(X1，t/ha)、穗長(X2，cm)、百粒重(X3，g)、蛋白質(zhì)含量(X4，%)和百克蛋白賴氨酸含量(X5，%)，獲得的數(shù)據(jù)樣本如表4-2所示。試創(chuàng)建它的SAS數(shù)據(jù)表sasuser.yumi。

(1)用表4-2所示數(shù)據(jù)創(chuàng)建Excel數(shù)據(jù)表yumi.xls。

(2)將Excel數(shù)據(jù)表yumi.xls導(dǎo)入SAS，創(chuàng)建SAS數(shù)據(jù)表sasuser.yumi，以備后面的SAS繪圖程序調(diào)用。4.2.1單變量折線圖

單變量折線圖，可用于單個(gè)連續(xù)變量觀測(cè)數(shù)據(jù)的時(shí)序或指定順序的展現(xiàn)。

【例4-6】試編寫SAS程序調(diào)用數(shù)據(jù)表sasuser.yumi繪制折線圖，展現(xiàn)12個(gè)玉米雜交種(variety)的產(chǎn)量(X1，t/ha)動(dòng)態(tài)。

(1)采用gplot過程編寫繪制折線圖的SAS程序如下：goptionsreset=allftext=swisshtext=1.95;

symbolV=hashH=2.25CV=BLCI=Bi=join;

axis1label=(f='宋體''玉米雜交種編號(hào)');

axis2label=(A=90f='宋體''產(chǎn)量(t/ha)');

procsortdata=sasuser.yumi;

byVariety;/*按橫坐標(biāo)變量排序*/

run;

procgplotdata=sasuser.yumi;

plotX1*Variety/noframehaxis=axis1vaxis=axis2;

run;quit;(2)程序的輸出結(jié)果如圖4-4所示。圖4-4玉米雜交種的產(chǎn)量動(dòng)態(tài)4.2.2變量相關(guān)折線圖

變量相關(guān)折線圖，通常用于展現(xiàn)兩連續(xù)變量觀測(cè)數(shù)據(jù)的相關(guān)關(guān)系。

【例4-7】試編寫SAS程序調(diào)用數(shù)據(jù)表sasuser.yumi繪制折線圖，展現(xiàn)12個(gè)玉米雜交種上百克蛋白賴氨酸含量(X5)與蛋白質(zhì)含量(X4)的關(guān)系。

(1)采用gplot過程編寫繪制折線圖的SAS程序如下：

goptionsreset=allftext=swisshtext=1.95;

symbolV=hashH=2.25CV=BLCI=Bi=join;

axis1label=(f='宋體''蛋白質(zhì)含量(%)')order=7.5to10by0.5;

axis2label=(A=90f='宋體''百克蛋白賴氨酸含量(%)');

procsortdata=sasuser.yumi;

byX4;/*按橫坐標(biāo)變量排序*/

procgplotdata=sasuser.yumi;

plotX5*X4/noframehaxis=axis1vaxis=axis2;

run;quit;(2)程序的輸出結(jié)果如圖4-5所示。圖4-5玉米雜交種上百克蛋白賴氨酸含量與蛋白質(zhì)含量的關(guān)系4.2.3多變量折線圖

多變量折線圖，通常用于展現(xiàn)多個(gè)連續(xù)變量按時(shí)間順序或指定順序的動(dòng)態(tài)變化，或多個(gè)連續(xù)變量與同一連續(xù)變量的關(guān)系。

【例4-8】試編寫SAS程序調(diào)用數(shù)據(jù)表sasuser.yumi繪制折線圖，展現(xiàn)玉米雜交種4個(gè)測(cè)定指標(biāo)X2、X3、X4和X5的品種動(dòng)態(tài)和相互關(guān)系。(1)采用gplot過程編寫繪制折線圖的SAS程序如下：

goptionsreset=allftext=swisshtext=1.95;

symbol1V=hashH=2.25CV=BLCI=Bi=join;

symbol2V=squareH=2.25CV=BLCI=Bi=join;

symbol3V=starH=2.25CV=BLCI=Bi=join;

symbol4V=circleH=2.25CV=BLCI=Bi=join;

axis1label=(f='宋體''玉米雜交種編號(hào)');

axis2label=(A=90f='宋體''指標(biāo)測(cè)定值');

legend1value=(f='宋體''產(chǎn)量''穗長''百粒重''蛋白質(zhì)''百克蛋白賴氨酸');

procsortdata=sasuser.yumi;

byVariety;

procgplotdata=sasuser.yumi;

plotX1*VarietyX2*VarietyX3*VarietyX4*VarietyX5*Variety/noframeoverlay

haxis=axis1vaxis=axis2legend=legend1;

run;quit;

(2)程序的輸出結(jié)果如圖4-6所示。圖4-6玉米雜交種4個(gè)測(cè)定指標(biāo)的品種動(dòng)態(tài)

4.3用SAS繪制盒須圖

盒須圖(boxandwhiskerplots)以圖形方式展現(xiàn)變量重復(fù)觀測(cè)數(shù)據(jù)的最小值、最大值、中值(50%分位點(diǎn))和指定百分位點(diǎn)(如5%分位點(diǎn)和95%分位點(diǎn))，如圖4-7所示。圖4-7盒須圖的框線結(jié)構(gòu)及意義

【例4-9】為考察保鮮處置配合儲(chǔ)藏天數(shù)(Treat)對(duì)梨棗生化性狀的影響，分別測(cè)定了過氧化物酶活性(POD)、過氧化氫酶活性(CAT)和VC含量(VC)的數(shù)據(jù)樣本，如表4-3所示。試?yán)L制盒須圖展現(xiàn)三測(cè)定指標(biāo)的百分位點(diǎn)分布。(1)創(chuàng)建如表4-3所示的單向分組SAS數(shù)據(jù)表sasuser.lizao。

(2)采用gplot過程編寫繪制盒須圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.05;

symbolI=BOXT05BWIDTH=15CI=B;

axis1label=(f='宋體''處理號(hào)')offset=(25,15);

axis2label=(A=90f='宋體''指標(biāo)測(cè)定值');

procsortdata=sasuser.lizao;

byTreats;

procgplotdata=sasuser.lizao;

plotPOD*Treats/noframeoverlayhaxis=axis1vaxis=axis2;

plotCAT*Treats/noframeoverlayhaxis=axis1vaxis=axis2;

plotVC*Treats/noframeoverlayhaxis=axis1vaxis=axis2;

run;quit;

(3)語句symbol的I=BOXT05項(xiàng)等號(hào)右端由BOX、T和05三個(gè)字符串拼連構(gòu)成，BOX指定繪盒須圖，T指定須線延伸到最大值和最小值線，05指定最小分位點(diǎn)為5%，最大分位點(diǎn)為95%，CI=B指定盒須圖線為藍(lán)色。語句axis1的offset=(25,15)項(xiàng)指定第一個(gè)盒須圖框邊界距左軸端25個(gè)單位，最后一個(gè)盒須圖框邊界距右軸端15個(gè)單位。

(4)程序的輸出結(jié)果如圖4-8、圖4-9和圖4-10所示。圖4-8梨棗儲(chǔ)藏保鮮試驗(yàn)POD觀測(cè)的百分位點(diǎn)分布圖4-9梨棗儲(chǔ)藏保鮮試驗(yàn)CAT觀測(cè)的百分位點(diǎn)分布圖4-10梨棗儲(chǔ)藏保鮮試驗(yàn)VC觀測(cè)的百分位點(diǎn)分布

4.4用SAS繪制曲線圖

曲線圖以圖形方式展現(xiàn)連續(xù)變量的回歸函數(shù)或已知函數(shù)的圖形，回歸曲線圖常附加試驗(yàn)數(shù)據(jù)點(diǎn)，以觀察數(shù)據(jù)的擬合效果。

4.4.1擬合曲線圖

【例4-10】為考察冬棗的多酚氧化酶活性(PPO)和過氧化物酶活性(POD)在某處置下與儲(chǔ)藏天數(shù)(Days)的相關(guān)關(guān)系，測(cè)定了12種儲(chǔ)藏天數(shù)(儲(chǔ)藏天數(shù)0作對(duì)照)的數(shù)據(jù)樣本，如表4-4所示。試?yán)L圖展現(xiàn)試驗(yàn)點(diǎn)與多項(xiàng)式擬合曲線。

(1)創(chuàng)建如表4-4所示的單向分組SAS數(shù)據(jù)表sasuser.dongzao。

(2)變量PPO(Y1)和POD(Y2)分別對(duì)變量Days(X)做多項(xiàng)式回歸，獲得如下回歸方程(詳細(xì)過程略)：

(3)采用gplot過程編寫繪制試驗(yàn)點(diǎn)和擬合曲線的SAS程序如下：

goptionsreset=allftext=swisshtext=2.05;

symbol01V=hashH=2.25CV=B;

symbol02I=splineCI=R;

axis1label=(f='宋體''儲(chǔ)藏天數(shù)(d)')offset=(5,5);

axis2label=(A=90f='宋體''PPO值');

axis3label=(A=90f='宋體''POD值');

dataaa;

dodays=0to88by1;

E_PPO=0.886552+0.003550*days*days-.00007425*days*days*days+0.00000056*days*days*days*days;E_POD=0.0663663+0.0307833*days-.0007755*days*days+0.00000504*days*days*days;

output;

end;

run;

databb;

setsasuser.dongzaoaa;

run;

procsortdata=bb;

byDays;

procgplotdata=bb;

plotPPO*DaysE_PPO*Days/noframeoverlayhaxis=axis1vaxis=axis2;

plotPOD*DaysE_POD*Days/noframeoverlayhaxis=axis1vaxis=axis3;

run;quit;

(4)語句symbol01對(duì)plot語句中的第1個(gè)繪圖項(xiàng)(例如PPO*Days)起作用。語句symbol02對(duì)plot語句的第2個(gè)繪圖項(xiàng)(例如E_PPO*Days)起作用。

(5)程序的輸出結(jié)果如圖4-11和圖4-12所示。圖4-11冬棗PPO活性與儲(chǔ)藏天數(shù)的試驗(yàn)點(diǎn)和擬合曲線圖4-12梨棗POD活性與儲(chǔ)藏天數(shù)的試驗(yàn)點(diǎn)和擬合曲線4.4.2函數(shù)的圖形

函數(shù)圖形通常用于展現(xiàn)兩連續(xù)變量的函數(shù)關(guān)系。

【例4-11】已知兩變量的關(guān)系由下面的函數(shù)決定：試編寫SAS程序繪出該函數(shù)的圖形。

(1)采用gplot過程編寫繪制函數(shù)圖形的SAS程序如下：

goptionsreset=allftext=swisshtext=2.05;

symbolI=splineCI=B;

axis1label=(f='宋體''自變量X');

axis2label=(A=90f='宋體''響應(yīng)Y');

datahanshu;

dox=-30to30by0.1;

y=sin(x)/x;

output;

end;

run;

procgplotdata=hanshu;

ploty*x/noframehaxis=axis1vaxis=axis2;

run;quit;(2)程序的輸出結(jié)果如圖4-13所示。圖4-13函數(shù)

的圖形

4.5用SAS繪制餅圖

餅圖顧名思義形似“圓餅”，一般用于展示離散或連續(xù)變量的百分率分布。

【例4-12】各國(Country)的榛子產(chǎn)量占全球總產(chǎn)量的份額(Output，%)如表4-5所示。試用餅圖展現(xiàn)全球榛子產(chǎn)量的百分率分布。

(1)創(chuàng)建繪圖用的SAS數(shù)據(jù)表sasuser.zhenzi，如表4-5所示。

(2)采用gchart過程編寫繪制餅圖的SAS程序如下：

goptionsreset=allftext='宋體'htext=2.25;

procgchartdata=sasuser.zhenzi;

piecountry/noheadingtype=pctfreq=output;/*繪二維餅圖*/

pie3dcountry/noheadingtype=pctfreq=output;/*繪三維餅圖*/

run;quit;

(3)程序的輸出結(jié)果如圖4-14所示。圖4-14榛子產(chǎn)量在全球的分布(a)二維餅圖(b)三維餅圖

4.6用SAS繪制柱形圖

柱形圖通常用于展現(xiàn)離散變量的頻數(shù)分布，相鄰柱之間有間隙，柱寬無意義。

【例4-13】抽樣觀測(cè)每頭蒜的蒜瓣數(shù)(garlic)，其頻數(shù)(frequency)樣本如表4-6所示。試創(chuàng)建頻數(shù)樣本的SAS數(shù)據(jù)表sasuser.suanban1。

(1)使用表4-6的數(shù)據(jù)創(chuàng)建Excel數(shù)據(jù)表suanban1.xls。

(2)將Excel數(shù)據(jù)表suanban1.xls導(dǎo)入SAS，并創(chuàng)建SAS數(shù)據(jù)表sasuser.suanban1，以備后面的SAS繪圖程序調(diào)用。

【例4-14】抽樣觀測(cè)每頭蒜的蒜瓣數(shù)(garlic)，其觀測(cè)值樣本如表4-7所示(表中未列出全部數(shù)據(jù)，可由表4-6推算)。試創(chuàng)建該樣本的SAS數(shù)據(jù)表sasuser.suanban2。

(1)使用表4-7的數(shù)據(jù)創(chuàng)建Excel數(shù)據(jù)表suanban2.xls。

(2)將Excel數(shù)據(jù)表suanban2.xls導(dǎo)入SAS，創(chuàng)建SAS數(shù)據(jù)表sasuser.suanban2，以備后面的SAS繪圖程序調(diào)用。4.6.1離散變量的頻數(shù)分布圖

離散變量的頻數(shù)分布一般采用柱形圖展現(xiàn)。

【例4-15】試編寫SAS程序調(diào)用頻數(shù)數(shù)據(jù)表sasuser.suanban1繪柱形圖，展現(xiàn)蒜瓣數(shù)觀測(cè)的頻數(shù)分布。

(1)采用gchart過程編寫繪制柱形圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋體''蒜瓣數(shù)');

axis2label=(A=90f='宋體''頻數(shù)');

procgchartdata=sasuser.suanban1;

vbargarlic/freq=frequencymaxis=axis1raxis=axis2width=10noframe;

run;quit;(2)程序的輸出結(jié)果如圖4-15所示。圖4-15蒜瓣數(shù)觀測(cè)的頻數(shù)分布

【例4-16】試編寫SAS程序調(diào)用觀測(cè)值數(shù)據(jù)表sasuser.suanban2繪柱形圖，展現(xiàn)蒜瓣數(shù)觀測(cè)的頻數(shù)分布。

(1)采用gchart過程編寫繪制柱形圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f=‘宋體’‘蒜瓣數(shù)’);

axis2label=(A=90f=‘宋體’‘頻數(shù)’);

procgchartdata=sasuser.suanban2;

vbargarlic/midpoints=8to15maxis=axis1raxis=axis2

width=10noframe;

run;quit;

(2)程序的輸出結(jié)果如圖4-15所示。與【例4-15】相同。4.6.2離散變量的累積頻數(shù)分布圖

【例4-17】試編寫調(diào)用頻數(shù)數(shù)據(jù)表sasuser.suanban1的SAS程序繪柱形圖，展現(xiàn)蒜瓣數(shù)觀測(cè)的累積頻數(shù)分布。

(1)采用gchart過程編寫繪制柱形圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋體''蒜瓣數(shù)');

axis2label=(A=90f='宋體''累積頻數(shù)');

procgchartdata=sasuser.suanban1;

vbargarlic/freq=frequencytype=cfreqmaxis=axis1raxis=axis2width=10noframe;

run;quit;(2)程序的輸出結(jié)果如圖4-16所示。圖4-16蒜瓣數(shù)觀測(cè)的累積頻數(shù)分布

【例4-18】試編寫調(diào)用觀測(cè)值數(shù)據(jù)表sasuser.suanban2的SAS程序繪柱形圖，展現(xiàn)蒜瓣數(shù)觀測(cè)的累積頻數(shù)分布。

(1)采用gchart過程編寫繪制柱形圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f=‘宋體’‘蒜瓣數(shù)’);

axis2label=(A=90f=‘宋體’‘累積頻數(shù)’);

procgchartdata=sasuser.suanban2;

vbargarlic/midpoints=8to15type=cfreqmaxis=axis1

raxis=axis2width=10noframe;

run;quit;

(2)程序的輸出結(jié)果如圖4-16所示，與例4-17相同。4.6.3離散變量的百分率分布圖

【例4-19】試編寫調(diào)用頻數(shù)數(shù)據(jù)表sasuser.suanban1的SAS程序繪柱形圖，展現(xiàn)蒜瓣數(shù)觀測(cè)的百分率分布。

(1)采用gchart過程編寫繪制柱形圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f=‘宋體’‘蒜瓣數(shù)’);

axis2label=(A=90f=‘宋體’‘百分率’);

procgchartdata=sasuser.suanban1;

vbargarlic/freq=frequencytype=pctmaxis=axis1raxis=axis2width=10noframe;

run;quit;

(2)程序的輸出結(jié)果如圖4-17所示。圖4-17蒜瓣數(shù)觀測(cè)的百分率分布

【例4-20】試編寫調(diào)用觀測(cè)值數(shù)據(jù)表sasuser.suanban2的SAS程序繪制柱形圖，展現(xiàn)蒜瓣數(shù)觀測(cè)的百分率分布。

(1)采用gchart過程編寫繪制柱形圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋體''蒜瓣數(shù)');

axis2label=(A=90f='宋體''百分率');

procgchartdata=sasuser.suanban2;

vbargarlic/midpoints=8to15type=pctmaxis=axis1raxis=axis2width=10noframe;

run;quit;

(2)程序的輸出結(jié)果如圖4-17所示。與例4-19相同。4.6.4離散變量的累積百分率分布圖

【例4-21】試編寫調(diào)用頻數(shù)數(shù)據(jù)表sasuser.suanban1的SAS程序繪制柱形圖，展現(xiàn)蒜瓣數(shù)觀測(cè)的累積百分率分布。

(1)采用gchart過程編寫繪制柱形圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋體''蒜瓣數(shù)');

axis2label=(A=90f='宋體''累積百分率');

procgchartdata=sasuser.suanban1;

vbargarlic/freq=frequencytype=cpctmaxis=axis1raxis=axis2width=10noframe;

run;quit;

(2)程序的輸出結(jié)果如圖4-18所示。圖4-18蒜瓣數(shù)觀測(cè)的累積百分率分布

【例4-22】試編寫調(diào)用觀測(cè)值數(shù)據(jù)表sasuser.suanban2的SAS程序繪柱形圖，展現(xiàn)蒜瓣數(shù)觀測(cè)的累積百分率分布。

(1)采用gchart過程編寫繪制柱形圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋體''蒜瓣數(shù)');

axis2label=(A=90f='宋體''累積百分率');

procgchartdata=sasuser.suanban2;

vbargarlic/midpoints=8to15type=cpctmaxis=axis1raxis=axis2width=10noframe;

run;quit;

(2)程序的輸出結(jié)果如圖4-18所示，與例4-21相同。

4.7用SAS繪制直方圖

直方圖通常用于展示連續(xù)變量的頻數(shù)和累積頻數(shù)分布、百分率和累積百分率分布。

【例4-23】在來自某果園的10箱串枝紅杏中各抽測(cè)10個(gè)果重(weight)數(shù)據(jù)，共100個(gè)數(shù)據(jù)，如表4-8所示。試用這些觀測(cè)創(chuàng)建SAS數(shù)據(jù)表sasuser.czh01。

(1)將表4-8中的100個(gè)數(shù)據(jù)以表4-9所示的格式輸入Excel，輸入時(shí)不用考慮表4-8數(shù)據(jù)的排列和順序(一種記錄而已)，輸入完畢存盤為Excel數(shù)據(jù)表czh01.xls。

(2)將數(shù)據(jù)表czh01.xls導(dǎo)入SAS，創(chuàng)建SAS數(shù)據(jù)表sasuser.czh01，以備后面程序調(diào)用。

【例4-24】利用freq過程或手工對(duì)SAS數(shù)據(jù)表sasuser.czh01做分組頻數(shù)統(tǒng)計(jì)，結(jié)果如表4-10所示。試將表4-10創(chuàng)建為SAS數(shù)據(jù)表sasuser.czh02，包括組下限(group_low)、組上限(group_up)、組中值(median)和組頻數(shù)(frequency)。

SAS解決問題的過程如下：

(1)對(duì)表4-9的數(shù)據(jù)進(jìn)行頻數(shù)統(tǒng)計(jì)：確定分組數(shù)9、組距3.5、第1組下限43.75，然后確定所有組限，統(tǒng)計(jì)各個(gè)組包含的觀測(cè)個(gè)數(shù)(組頻數(shù))并計(jì)算組中值，結(jié)果如表4-10所示。

(2)將表4-10數(shù)據(jù)按本身格式創(chuàng)建為Excel數(shù)據(jù)表czh02.xls。

(3)將數(shù)據(jù)表czh02.xls導(dǎo)入SAS，創(chuàng)建SAS數(shù)據(jù)表sasuser.czh02，以備后面的SAS繪圖程序調(diào)用。4.7.1連續(xù)變量的頻數(shù)分布圖

【例4-25】試編寫調(diào)用觀測(cè)數(shù)據(jù)表sasuser.czh01的SAS程序繪制直方圖，以展現(xiàn)果重觀測(cè)的頻數(shù)分布。

SAS解決問題的過程如下：

(1)采用gchart過程編寫繪制直方圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋體''串枝紅杏果重(g)');

axis2label=(A=90f='宋體''頻數(shù)');

procgchartdata=sasuser.czh01;

vbarweight/width=10space=0axis=0to35by5midpoints=45.5to73.5by3.5

maxis=axis1raxis=axis2noframe;

run;quit;(2)程序的輸出結(jié)果如圖4-19所示。圖4-19串枝紅杏果重觀測(cè)的頻數(shù)分布

【例4-26】試編寫調(diào)用頻數(shù)統(tǒng)計(jì)數(shù)據(jù)表sasuser.czh02的SAS程序繪制直方圖，以展現(xiàn)果重觀測(cè)的頻數(shù)分布。

(1)采用gchart過程編寫繪制直方圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋體''串枝紅杏果重(g)');

axis2label=(A=90f='宋體''頻數(shù)');

procgchartdata=sasuser.czh02;

vbarmedian/width=10space=0freq=frequencyaxis=0to35by5

midpoints=45.5to73.5by3.5maxis=axis1raxis=axis2noframe;

run;quit;

(2)程序的輸出結(jié)果如圖4-19所示，與例4-25相同。4.7.2連續(xù)變量的累積頻數(shù)分布圖

【例4-27】試編寫調(diào)用觀測(cè)值數(shù)據(jù)表sasuser.czh01的SAS程序繪制直方圖，以展現(xiàn)果重觀測(cè)的累積頻數(shù)分布。

(1)采用gchart過程編寫繪制直方圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋體''串枝紅杏果重(g)');

axis2label=(A=90f='宋體''累積頻數(shù)');

procgchartdata=sasuser.czh01;

vbarweight/width=10space=0type=cfreqmidpoints=45.5to73.5by3.5

maxis=axis1raxis=axis2noframe;

run;quit;(2)程序的輸出結(jié)果如圖4-20所示。圖4-20串枝紅杏果重觀測(cè)的累積頻數(shù)分布

【例4-28】試編寫調(diào)用頻數(shù)統(tǒng)計(jì)數(shù)據(jù)表sasuser.czh02的SAS程序繪制直方圖，以展現(xiàn)果重觀測(cè)的累積頻數(shù)分布。

(1)采用gchart過程編寫繪制直方圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋體''串枝紅杏果重(g)');

axis2label=(A=90f='宋體''累積頻數(shù)');

procgchartdata=sasuser.czh02;

vbarmedian/width=10space=0type=cfreqfreq=frequency

midpoints=45.5to73.5by3.5maxis=axis1raxis=axis2noframe;

run;quit;

(2)程序的輸出結(jié)果如圖4-20所示，與例4-27相同。4.7.3連續(xù)變量的百分率分布圖

【例4-29】試編寫調(diào)用觀測(cè)值數(shù)據(jù)表sasuser.czh01的SAS程序繪制直方圖，以展現(xiàn)果重觀測(cè)的百分率分布。

(1)采用gchart過程編寫繪制直方圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋體''串枝紅杏果重(g)');

axis2label=(A=90f='宋體''百分率');

procgchartdata=sasuser.czh01;

vbarweight/width=10space=0type=pctaxis=0to35by5

midpoints=45.5to73.5by3.5maxis=axis1raxis=axis2noframe;

run;quit;(2)程序的輸出結(jié)果如圖4-21所示。圖4-21串枝紅杏果重觀測(cè)的百分率分布

【例4-30】試編寫調(diào)用頻數(shù)統(tǒng)計(jì)數(shù)據(jù)表sasuser.czh02的SAS程序繪制直方圖，以展現(xiàn)果重觀測(cè)的百分率分布。

(1)采用gchart過程編寫繪制直方圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋體''串枝紅杏果重(g)');

axis2label=(A=90f='宋體''百分率');

procgchartdata=sasuser.czh02;

vbarmedian/width=10space=0type=pctfreq=frequencyaxis=0to35by5midpoints=45.5to73.5by3.5maxis=axis1raxis=axis2noframe;

run;quit;4.7.4連續(xù)變量的累積百分率分布圖

【例4-31】試編寫調(diào)用觀測(cè)值數(shù)據(jù)表sasuser.czh01的SAS程序繪制直方圖，以展現(xiàn)果重觀測(cè)的累積百分率分布。

(1)采用gchart過程編寫繪制直方圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋體''串枝紅杏果重(g)');

axis2label=(A=90f='宋體''累積百分率');

procgchartdata=sasuser.czh01;

vbarweight/width=10space=0type=cpctmidpoints=45.5to73.5by3.5

maxis=axis1raxis=axis2noframe;

run;quit;(2)程序的輸出結(jié)果如圖4-22所示。圖4-22串枝紅杏果重觀測(cè)的累積百分率分布

【例4-32】試編寫調(diào)用頻數(shù)統(tǒng)計(jì)數(shù)據(jù)表sasuser.czh02的SAS程序繪制直方圖，以展現(xiàn)果重觀測(cè)的累積百分率分布。

(1)采用gchart過程編寫繪制直方圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋體''串枝紅杏果重(g)');

axis2label=(A=90f='宋體''累積百分率');

procgchartdata=sasuser.czh02;

vbarmedian/width=10space=0type=cpctfreq=frequency

midpoints=45.5to73.5by3.5maxis=axis1raxis=axis2noframe;

run;quit;4.8用SAS繪制誤差圖

柱形圖或直方圖附加誤差條稱做誤差圖。誤差圖通常用于展示變量重復(fù)抽樣的觀測(cè)值分布特征，一般以重復(fù)抽樣觀測(cè)的均值繪數(shù)據(jù)點(diǎn)，缺省以均值的95%置信區(qū)間繪誤差條。

【例4-33】在某榛園觀測(cè)了4個(gè)榛樹品種(Treats)的土壤氮含量(N)、土壤磷含量(P)、土壤鉀含量(K)和葉片凈光合速率(Pn)，如表4-11所示。試?yán)L制誤差圖展現(xiàn)各個(gè)觀測(cè)指標(biāo)的品種動(dòng)態(tài)(均值和置信區(qū)間)。

(1)使用表4-11的數(shù)據(jù)創(chuàng)建Excel數(shù)據(jù)表errorbar.xls。

(2)將Excel數(shù)據(jù)表errorbar.xls導(dǎo)入SAS，創(chuàng)建SAS數(shù)據(jù)表sasuser.errorbar，以備后面的SAS程序調(diào)用。(3)采用gchart過程編寫繪制誤差圖的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋體''處理');

axis2label=(A=90f='宋體''N含量(%)');

axis3label=(A=90f='宋體''P含量(%)');

axis4label=(A=90f='宋體''K含量(%)');

axis5label=(A=90f='宋體''光合速率Pn');

procgchartdata=sasuser.errorbar;

vbarTreats/sumvar=Ntype=meanwidth=15maxis=axis1raxis=axis2

noframeerrorbar=bothcoutline=black;

vbarTreats/sumvar=Ptype=meanwidth=15maxis=axis1raxis=axis3

noframeerrorbar=bothcoutline=black;

vbarTreats/sumvar=Ktype=meanwidth=15maxis=axis1raxis=axis4

noframeerrorbar=bothcoutline=black;

vbarTreats/sumvar=Pntype=meanwidth=15maxis=axis1raxis=axis5

noframeerrorbar=bothcoutline=black;

run;quit;

(4)程序的輸出結(jié)果如圖4-23所示。圖4-23榛樹土壤養(yǎng)分及光合速率的品種動(dòng)態(tài)(a)氮含量

(b)?P含量圖4-23榛樹土壤養(yǎng)分及光合速率的品種動(dòng)態(tài)(c)?K含量

(d)凈光合速率

4.9用SAS繪制三維網(wǎng)格圖

網(wǎng)格圖(網(wǎng)狀曲面圖)可在三維空間坐標(biāo)系中展現(xiàn)函數(shù)Z?=?f?(X,Y)的圖形(連續(xù)變量的響應(yīng)面)，即一個(gè)響應(yīng)變量對(duì)兩個(gè)自變量的回歸曲面。

網(wǎng)格圖的基本原理是，先