试验设计与统计分析SAS实践教程课件：SAS统计绘图

SAS统计绘图4.1用SAS绘制散点图4.2用SAS绘制折线图4.3用SAS绘制盒须图4.4用SAS绘制曲线图4.5用SAS绘制饼图4.6用SAS绘制柱形图4.7用SAS绘制直方图4.8用SAS绘制误差图4.9用SAS绘制三维网格图4.10用SAS绘制等值线图上机报告

上机目的

掌握SAS绘统计图的编程方法，熟悉程序中的过程、过程选项、语句、语句选项等编程要素和编程格式。学会用统计图展现试验数据、拟合曲线、函数图形、频数和累积频数分布、百分率和累积百分率分布、响应误差、响应面和等值线。

上机内容①采用gplot过程编写SAS程序，绘散点图、折线图和曲线图；②采用gchart过程编写SAS程序，绘饼图、柱形图、直方图和误差图；③采用g3d过程编写SAS程序，绘三维网格图；④采用gcontour过程编写SAS程序，绘等值线图。

4.1用SAS绘制散点图

在执行SAS绘图程序之前，首先需创建或选定用于绘图的SAS数据表。

【例4-1】在某地连续观测13个年份(Order)的冬季积雪时间(X1)、化雪日期(X2)和二化螟第一代成虫发生量(Y)，获得的数据样本如表4-1所示。试创建它的SAS数据表sasuser.chengchong。

(1)用表4-1数据创建Excel数据表chengchong.xls。

(2)将Excel数据表chengchong.xls导入SAS，创建SAS数据表sasuser.chengchong，以备后面的SAS绘图程序调用。4.1.1单变量散点图

单变量散点图，一般按时间顺序或指定顺序展现变量观测值的动态变化。

【例4-2】试编写SAS程序调用数据表sasuser.chengchong(表4-1)绘散点图，展现二化螟第一代成虫发生量(Y)的年份(Order)动态。

(1)采用gplot过程编写绘制散点图的SAS程序如下：

goptionsreset=allftext=swisshtext=1.65; /*设置文本的字体和尺寸*/

symbolV=starH=1.75CV=B;/*设置数据点标记的类型、尺寸和颜色*/

axis1label=(f='宋体''年份编号');/*设置横轴标签的字体和内容*/

axis2label=(A=90f='宋体''成虫发生量'); /*设置纵轴标签的角度、字体和内容*/

procgplotdata=sasuser.chengchong;/*指定gplot过程调用的数据表*/

plotY*Order/noframehaxis=axis1vaxis=axis2;

/*选项去图框、引用纵横轴的设置*/

run;quit;(2)程序的输出结果如图4-1所示。图4-1二化螟成虫发生量的年份动态4.1.2两变量散点图

两变量散点图，一般用于展现一个变量与另一个变量的相关关系。

【例4-3】试编写SAS程序调用数据表sasuser.chengchong绘制散点图，展现二化螟第一代成虫发生量(Y)与积雪时间(X1)的关系。(1)采用gplot过程编写绘制散点图的SAS程序如下：

goptionsreset=allftext=swisshtext=1.65;

symbolV=starH=1.75CV=B;

axis1label=(f='宋体''积雪时间(周)');

axis2label=(A=90f='宋体''成虫发生量');

procgplotdata=sasuser.chengchong;

plotY*X1/noframehaxis=axis1vaxis=axis2;

run;quit;(2)程序的输出结果如图4-2所示。图4-2二化螟成虫发生量与积雪时间的关系4.1.3多变量散点图

多变量散点图，可展现多个变量的时间动态，或多个变量与同一变量的相关关系。

【例4-4】试编写SAS程序调用数据表sasuser.chengchong绘制散点图，展现二化螟第一代成虫发生量(Y)、积雪时间(X1)、化雪日期(X2)的年份动态。(1)采用gplot过程编写绘制散点图的SAS程序如下：

goptionsreset=allftext=swisshtext=1.55;

symbol1V=starH=1.75CV=black;

symbol2V=squareH=1.75CV=B;

symbol3V=hashH=1.75CV=R;

axis1label=(f='宋体''年份编号');

axis2label=(A=90f='宋体''*成虫发生量□积雪时间

＃化雪日期');

procgplotdata=sasuser.chengchong;

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

run;quit; /*选项overlay使多个散点图叠加*/(2)程序的输出结果如图4-3所示。图4-3成虫发生量、积雪时间、化雪日期的年份动态

4.2用SAS绘制折线图

在执行SAS绘图程序之前，首先需创建或选定用于绘图的SAS数据表。

【例4-5】抽样观测12个高营养玉米杂交种(variety)的产量(X1，t/ha)、穗长(X2，cm)、百粒重(X3，g)、蛋白质含量(X4，%)和百克蛋白赖氨酸含量(X5，%)，获得的数据样本如表4-2所示。试创建它的SAS数据表sasuser.yumi。

(1)用表4-2所示数据创建Excel数据表yumi.xls。

(2)将Excel数据表yumi.xls导入SAS，创建SAS数据表sasuser.yumi，以备后面的SAS绘图程序调用。4.2.1单变量折线图

单变量折线图，可用于单个连续变量观测数据的时序或指定顺序的展现。

【例4-6】试编写SAS程序调用数据表sasuser.yumi绘制折线图，展现12个玉米杂交种(variety)的产量(X1，t/ha)动态。

(1)采用gplot过程编写绘制折线图的SAS程序如下：goptionsreset=allftext=swisshtext=1.95;

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

axis1label=(f='宋体''玉米杂交种编号');

axis2label=(A=90f='宋体''产量(t/ha)');

procsortdata=sasuser.yumi;

byVariety;/*按横坐标变量排序*/

run;

procgplotdata=sasuser.yumi;

plotX1*Variety/noframehaxis=axis1vaxis=axis2;

run;quit;(2)程序的输出结果如图4-4所示。图4-4玉米杂交种的产量动态4.2.2变量相关折线图

变量相关折线图，通常用于展现两连续变量观测数据的相关关系。

【例4-7】试编写SAS程序调用数据表sasuser.yumi绘制折线图，展现12个玉米杂交种上百克蛋白赖氨酸含量(X5)与蛋白质含量(X4)的关系。

(1)采用gplot过程编写绘制折线图的SAS程序如下：

goptionsreset=allftext=swisshtext=1.95;

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

axis1label=(f='宋体''蛋白质含量(%)')order=7.5to10by0.5;

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

procsortdata=sasuser.yumi;

byX4;/*按横坐标变量排序*/

procgplotdata=sasuser.yumi;

plotX5*X4/noframehaxis=axis1vaxis=axis2;

run;quit;(2)程序的输出结果如图4-5所示。图4-5玉米杂交种上百克蛋白赖氨酸含量与蛋白质含量的关系4.2.3多变量折线图

多变量折线图，通常用于展现多个连续变量按时间顺序或指定顺序的动态变化，或多个连续变量与同一连续变量的关系。

【例4-8】试编写SAS程序调用数据表sasuser.yumi绘制折线图，展现玉米杂交种4个测定指标X2、X3、X4和X5的品种动态和相互关系。(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='宋体''玉米杂交种编号');

axis2label=(A=90f='宋体''指标测定值');

legend1value=(f='宋体''产量''穗长''百粒重''蛋白质''百克蛋白赖氨酸');

procsortdata=sasuser.yumi;

byVariety;

procgplotdata=sasuser.yumi;

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

haxis=axis1vaxis=axis2legend=legend1;

run;quit;

(2)程序的输出结果如图4-6所示。图4-6玉米杂交种4个测定指标的品种动态

4.3用SAS绘制盒须图

盒须图(boxandwhiskerplots)以图形方式展现变量重复观测数据的最小值、最大值、中值(50%分位点)和指定百分位点(如5%分位点和95%分位点)，如图4-7所示。图4-7盒须图的框线结构及意义

【例4-9】为考察保鲜处置配合储藏天数(Treat)对梨枣生化性状的影响，分别测定了过氧化物酶活性(POD)、过氧化氢酶活性(CAT)和VC含量(VC)的数据样本，如表4-3所示。试绘制盒须图展现三测定指标的百分位点分布。(1)创建如表4-3所示的单向分组SAS数据表sasuser.lizao。

(2)采用gplot过程编写绘制盒须图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.05;

symbolI=BOXT05BWIDTH=15CI=B;

axis1label=(f='宋体''处理号')offset=(25,15);

axis2label=(A=90f='宋体''指标测定值');

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项等号右端由BOX、T和05三个字符串拼连构成，BOX指定绘盒须图，T指定须线延伸到最大值和最小值线，05指定最小分位点为5%，最大分位点为95%，CI=B指定盒须图线为蓝色。语句axis1的offset=(25,15)项指定第一个盒须图框边界距左轴端25个单位，最后一个盒须图框边界距右轴端15个单位。

(4)程序的输出结果如图4-8、图4-9和图4-10所示。图4-8梨枣储藏保鲜试验POD观测的百分位点分布图4-9梨枣储藏保鲜试验CAT观测的百分位点分布图4-10梨枣储藏保鲜试验VC观测的百分位点分布

4.4用SAS绘制曲线图

曲线图以图形方式展现连续变量的回归函数或已知函数的图形，回归曲线图常附加试验数据点，以观察数据的拟合效果。

4.4.1拟合曲线图

【例4-10】为考察冬枣的多酚氧化酶活性(PPO)和过氧化物酶活性(POD)在某处置下与储藏天数(Days)的相关关系，测定了12种储藏天数(储藏天数0作对照)的数据样本，如表4-4所示。试绘图展现试验点与多项式拟合曲线。

(1)创建如表4-4所示的单向分组SAS数据表sasuser.dongzao。

(2)变量PPO(Y1)和POD(Y2)分别对变量Days(X)做多项式回归，获得如下回归方程(详细过程略)：

(3)采用gplot过程编写绘制试验点和拟合曲线的SAS程序如下：

goptionsreset=allftext=swisshtext=2.05;

symbol01V=hashH=2.25CV=B;

symbol02I=splineCI=R;

axis1label=(f='宋体''储藏天数(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对plot语句中的第1个绘图项(例如PPO*Days)起作用。语句symbol02对plot语句的第2个绘图项(例如E_PPO*Days)起作用。

(5)程序的输出结果如图4-11和图4-12所示。图4-11冬枣PPO活性与储藏天数的试验点和拟合曲线图4-12梨枣POD活性与储藏天数的试验点和拟合曲线4.4.2函数的图形

函数图形通常用于展现两连续变量的函数关系。

【例4-11】已知两变量的关系由下面的函数决定：试编写SAS程序绘出该函数的图形。

(1)采用gplot过程编写绘制函数图形的SAS程序如下：

goptionsreset=allftext=swisshtext=2.05;

symbolI=splineCI=B;

axis1label=(f='宋体''自变量X');

axis2label=(A=90f='宋体''响应Y');

datahanshu;

dox=-30to30by0.1;

y=sin(x)/x;

output;

end;

run;

procgplotdata=hanshu;

ploty*x/noframehaxis=axis1vaxis=axis2;

run;quit;(2)程序的输出结果如图4-13所示。图4-13函数

的图形

4.5用SAS绘制饼图

饼图顾名思义形似“圆饼”，一般用于展示离散或连续变量的百分率分布。

【例4-12】各国(Country)的榛子产量占全球总产量的份额(Output，%)如表4-5所示。试用饼图展现全球榛子产量的百分率分布。

(1)创建绘图用的SAS数据表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)程序的输出结果如图4-14所示。图4-14榛子产量在全球的分布(a)二维饼图(b)三维饼图

4.6用SAS绘制柱形图

柱形图通常用于展现离散变量的频数分布，相邻柱之间有间隙，柱宽无意义。

【例4-13】抽样观测每头蒜的蒜瓣数(garlic)，其频数(frequency)样本如表4-6所示。试创建频数样本的SAS数据表sasuser.suanban1。

(1)使用表4-6的数据创建Excel数据表suanban1.xls。

(2)将Excel数据表suanban1.xls导入SAS，并创建SAS数据表sasuser.suanban1，以备后面的SAS绘图程序调用。

【例4-14】抽样观测每头蒜的蒜瓣数(garlic)，其观测值样本如表4-7所示(表中未列出全部数据，可由表4-6推算)。试创建该样本的SAS数据表sasuser.suanban2。

(1)使用表4-7的数据创建Excel数据表suanban2.xls。

(2)将Excel数据表suanban2.xls导入SAS，创建SAS数据表sasuser.suanban2，以备后面的SAS绘图程序调用。4.6.1离散变量的频数分布图

离散变量的频数分布一般采用柱形图展现。

【例4-15】试编写SAS程序调用频数数据表sasuser.suanban1绘柱形图，展现蒜瓣数观测的频数分布。

(1)采用gchart过程编写绘制柱形图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋体''蒜瓣数');

axis2label=(A=90f='宋体''频数');

procgchartdata=sasuser.suanban1;

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

run;quit;(2)程序的输出结果如图4-15所示。图4-15蒜瓣数观测的频数分布

【例4-16】试编写SAS程序调用观测值数据表sasuser.suanban2绘柱形图，展现蒜瓣数观测的频数分布。

(1)采用gchart过程编写绘制柱形图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f=‘宋体’‘蒜瓣数’);

axis2label=(A=90f=‘宋体’‘频数’);

procgchartdata=sasuser.suanban2;

vbargarlic/midpoints=8to15maxis=axis1raxis=axis2

width=10noframe;

run;quit;

(2)程序的输出结果如图4-15所示。与【例4-15】相同。4.6.2离散变量的累积频数分布图

【例4-17】试编写调用频数数据表sasuser.suanban1的SAS程序绘柱形图，展现蒜瓣数观测的累积频数分布。

(1)采用gchart过程编写绘制柱形图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋体''蒜瓣数');

axis2label=(A=90f='宋体''累积频数');

procgchartdata=sasuser.suanban1;

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

run;quit;(2)程序的输出结果如图4-16所示。图4-16蒜瓣数观测的累积频数分布

【例4-18】试编写调用观测值数据表sasuser.suanban2的SAS程序绘柱形图，展现蒜瓣数观测的累积频数分布。

(1)采用gchart过程编写绘制柱形图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f=‘宋体’‘蒜瓣数’);

axis2label=(A=90f=‘宋体’‘累积频数’);

procgchartdata=sasuser.suanban2;

vbargarlic/midpoints=8to15type=cfreqmaxis=axis1

raxis=axis2width=10noframe;

run;quit;

(2)程序的输出结果如图4-16所示，与例4-17相同。4.6.3离散变量的百分率分布图

【例4-19】试编写调用频数数据表sasuser.suanban1的SAS程序绘柱形图，展现蒜瓣数观测的百分率分布。

(1)采用gchart过程编写绘制柱形图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f=‘宋体’‘蒜瓣数’);

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

procgchartdata=sasuser.suanban1;

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

run;quit;

(2)程序的输出结果如图4-17所示。图4-17蒜瓣数观测的百分率分布

【例4-20】试编写调用观测值数据表sasuser.suanban2的SAS程序绘制柱形图，展现蒜瓣数观测的百分率分布。

(1)采用gchart过程编写绘制柱形图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋体''蒜瓣数');

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

procgchartdata=sasuser.suanban2;

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

run;quit;

(2)程序的输出结果如图4-17所示。与例4-19相同。4.6.4离散变量的累积百分率分布图

【例4-21】试编写调用频数数据表sasuser.suanban1的SAS程序绘制柱形图，展现蒜瓣数观测的累积百分率分布。

(1)采用gchart过程编写绘制柱形图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋体''蒜瓣数');

axis2label=(A=90f='宋体''累积百分率');

procgchartdata=sasuser.suanban1;

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

run;quit;

(2)程序的输出结果如图4-18所示。图4-18蒜瓣数观测的累积百分率分布

【例4-22】试编写调用观测值数据表sasuser.suanban2的SAS程序绘柱形图，展现蒜瓣数观测的累积百分率分布。

(1)采用gchart过程编写绘制柱形图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.25;

patternV=E;

axis1label=(f='宋体''蒜瓣数');

axis2label=(A=90f='宋体''累积百分率');

procgchartdata=sasuser.suanban2;

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

run;quit;

(2)程序的输出结果如图4-18所示，与例4-21相同。

4.7用SAS绘制直方图

直方图通常用于展示连续变量的频数和累积频数分布、百分率和累积百分率分布。

【例4-23】在来自某果园的10箱串枝红杏中各抽测10个果重(weight)数据，共100个数据，如表4-8所示。试用这些观测创建SAS数据表sasuser.czh01。

(1)将表4-8中的100个数据以表4-9所示的格式输入Excel，输入时不用考虑表4-8数据的排列和顺序(一种记录而已)，输入完毕存盘为Excel数据表czh01.xls。

(2)将数据表czh01.xls导入SAS，创建SAS数据表sasuser.czh01，以备后面程序调用。

【例4-24】利用freq过程或手工对SAS数据表sasuser.czh01做分组频数统计，结果如表4-10所示。试将表4-10创建为SAS数据表sasuser.czh02，包括组下限(group_low)、组上限(group_up)、组中值(median)和组频数(frequency)。

SAS解决问题的过程如下：

(1)对表4-9的数据进行频数统计：确定分组数9、组距3.5、第1组下限43.75，然后确定所有组限，统计各个组包含的观测个数(组频数)并计算组中值，结果如表4-10所示。

(2)将表4-10数据按本身格式创建为Excel数据表czh02.xls。

(3)将数据表czh02.xls导入SAS，创建SAS数据表sasuser.czh02，以备后面的SAS绘图程序调用。4.7.1连续变量的频数分布图

【例4-25】试编写调用观测数据表sasuser.czh01的SAS程序绘制直方图，以展现果重观测的频数分布。

SAS解决问题的过程如下：

(1)采用gchart过程编写绘制直方图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋体''串枝红杏果重(g)');

axis2label=(A=90f='宋体''频数');

procgchartdata=sasuser.czh01;

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

maxis=axis1raxis=axis2noframe;

run;quit;(2)程序的输出结果如图4-19所示。图4-19串枝红杏果重观测的频数分布

【例4-26】试编写调用频数统计数据表sasuser.czh02的SAS程序绘制直方图，以展现果重观测的频数分布。

(1)采用gchart过程编写绘制直方图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋体''串枝红杏果重(g)');

axis2label=(A=90f='宋体''频数');

procgchartdata=sasuser.czh02;

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

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

run;quit;

(2)程序的输出结果如图4-19所示，与例4-25相同。4.7.2连续变量的累积频数分布图

【例4-27】试编写调用观测值数据表sasuser.czh01的SAS程序绘制直方图，以展现果重观测的累积频数分布。

(1)采用gchart过程编写绘制直方图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋体''串枝红杏果重(g)');

axis2label=(A=90f='宋体''累积频数');

procgchartdata=sasuser.czh01;

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

maxis=axis1raxis=axis2noframe;

run;quit;(2)程序的输出结果如图4-20所示。图4-20串枝红杏果重观测的累积频数分布

【例4-28】试编写调用频数统计数据表sasuser.czh02的SAS程序绘制直方图，以展现果重观测的累积频数分布。

(1)采用gchart过程编写绘制直方图的SAS程序如下：

goptionsreset=allftext=swisshtext=2.0;

patternV=E;

axis1label=(f='宋体''串枝红杏果重(g)');

axis2label=(A=90f='宋体''累积频数');

procgchartdata=sasuser.czh02;

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

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

run;quit;

(2)程序的输出结果如图4-20所示，与例4-27相同。4.7.3连续变量的百分率分布图

【例4-29】试编写调用观测值数据表sasuser.czh01的SAS程序绘制直方图，以展现果重观测的百分率分布。

(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)程序的输出结果如图4-21所示。图4-21串枝红杏果重观测的百分率分布

【例4-30】试编写调用频数统计数据表sasuser.czh02的SAS程序绘制直方图，以展现果重观测的百分率分布。

(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连续变量的累积百分率分布图

【例4-31】试编写调用观测值数据表sasuser.czh01的SAS程序绘制直方图，以展现果重观测的累积百分率分布。

(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)程序的输出结果如图4-22所示。图4-22串枝红杏果重观测的累积百分率分布

【例4-32】试编写调用频数统计数据表sasuser.czh02的SAS程序绘制直方图，以展现果重观测的累积百分率分布。

(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绘制误差图

柱形图或直方图附加误差条称做误差图。误差图通常用于展示变量重复抽样的观测值分布特征，一般以重复抽样观测的均值绘数据点，缺省以均值的95%置信区间绘误差条。

【例4-33】在某榛园观测了4个榛树品种(Treats)的土壤氮含量(N)、土壤磷含量(P)、土壤钾含量(K)和叶片净光合速率(Pn)，如表4-11所示。试绘制误差图展现各个观测指标的品种动态(均值和置信区间)。

(1)使用表4-11的数据创建Excel数据表errorbar.xls。

(2)将Excel数据表errorbar.xls导入SAS，创建SAS数据表sasuser.errorbar，以备后面的SAS程序调用。(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)程序的输出结果如图4-23所示。图4-23榛树土壤养分及光合速率的品种动态(a)氮含量

(b) P含量图4-23榛树土壤养分及光合速率的品种动态(c) K含量

(d)净光合速率

4.9用SAS绘制三维网格图

网格图(网状曲面图)可在三维空间坐标系中展现函数Z = f (X,Y)的图形(连续变量的响应面)，即一个响应变量对两个自变量的回归曲面。

网格图的基本原理是，先