第3章SAS過(guò)程初步本章介紹SAS過(guò)程步的一般用法和常用語(yǔ)句的含義，以及幾個(gè)簡(jiǎn)單報(bào)表、分析、繪圖過(guò)程§3.1SAS過(guò)程用法及常用語(yǔ)句SAS過(guò)程用法 過(guò)程名DATA=輸入數(shù)據(jù)集選項(xiàng);過(guò)程語(yǔ)句/ 過(guò)程語(yǔ)句 選項(xiàng)PROC語(yǔ)句的選項(xiàng)是可選的，用來(lái)規(guī)定過(guò)程運(yùn)行的一些設(shè)置，如果有多個(gè)選項(xiàng)用空格RUNRUN語(yǔ)句而在下一個(gè)過(guò)程步或數(shù)據(jù)步的開(kāi)始處結(jié)束，另外PROC語(yǔ)句之后、結(jié)束之前可以有若干SAS過(guò)程步有些是對(duì)數(shù)據(jù)集作某種變換（SORT過(guò)程對(duì)數(shù)據(jù)集排序OUTPUT窗口（高精度SAS過(guò)程步常用語(yǔ)句SAS過(guò)程步中常見(jiàn)的語(yǔ)句，更具體的用法可以在以后實(shí)際用到一、VAR語(yǔ)句VAR語(yǔ)句在很多過(guò)程中用來(lái)指定分析變量。在VAR后面給出變量 變量名1變量名 ?變量名變量名列表可以使用省略的形式，如X1-X3，MATH--等。VAR用法例如varmathMODEL語(yǔ)句 ＝/選項(xiàng)modelmath=三、BYCLASSprocprintdata=c9501; TheSAS 98 45當(dāng)于因變量。而在另一些過(guò)程（如MEANS）中，CLASS語(yǔ)句作用與BY語(yǔ)句類似，可以指定四、OUTPUTOUTPUT語(yǔ)句指定輸出結(jié)果存放的數(shù)據(jù)集。不同過(guò)程中把輸出結(jié)果存入數(shù)據(jù)集的方法各有不同，OUTPUT語(yǔ)句是用得最多的一種，其一般格式為：OUTPUTOUT＝輸出數(shù)據(jù)集名關(guān)鍵字＝＝OUT＝給出了要生成的結(jié)果數(shù)據(jù)集的名字，用“關(guān)鍵字＝變量名”的方式指定了輸procmeansvaroutputout=resultn=nmean=meanmathvar=varmath;procprintdata=result;五、FREQ語(yǔ)句和WEIGHT完全相同的重復(fù)觀測(cè)。FREQ變量只取整數(shù)值。如freqWEIGHT語(yǔ)句指定一個(gè)權(quán)重變量，在某些允許的過(guò)程中代表權(quán)重，其值與觀測(cè)對(duì)六、ID語(yǔ)句但是，如果數(shù)據(jù)集中有一個(gè)變量可以用來(lái)區(qū)分觀測(cè)（如人名、省市名，就可以用ID語(yǔ)句id指定用變量NAME的值來(lái)標(biāo)識(shí)觀測(cè)。七、WHEREwheremath>=60and八、LABELFORMAT 語(yǔ)句為變量指定一個(gè)，很多過(guò)程可以使用這樣的。其格式 變量名＝' 變量名＝' procprintdata=sasuser.c9501label;idname;varmathlabelname=''math='數(shù)學(xué)成績(jī)'='語(yǔ)文成績(jī)';procprintdata=sasuser.c9501;formatmath5.15.1;LABEL語(yǔ)句規(guī)定變量的，用LENGTH語(yǔ)句規(guī)定變量的存貯長(zhǎng)度，用ATTRIB語(yǔ)句同§3.2PRINT過(guò)程制作列表報(bào)告基本用法proc了列出一個(gè)指定的數(shù)據(jù)集，在PROC語(yǔ)句中使用DATA＝選項(xiàng)指定要列表的輸入數(shù)據(jù)集procprintdata=sasuser.gpa;在過(guò)程內(nèi)使用VAR語(yǔ)句可以指定要列出的變量并指定順序procprintdata=c9501;varname;TheSAS 3男109男5110女可以在PROCPRINT語(yǔ)句中加入NOOBS選項(xiàng)，如：procprintdata=c9501noobs;使用WHERE可以指定只列出滿足條件的觀測(cè)。比如，procprintwherenamein('',''); MATH AVG1男4男使用中文列標(biāo)題PRINT語(yǔ)句中加LABEL選項(xiàng)，例如：procprintdata=c9501noobslabel;varname mathavg;labelname='' =''math='數(shù)學(xué)'='語(yǔ)文'avg='平均分';TheSAS 男女男男女3.2.3標(biāo)題及全程語(yǔ)句我們從上面的輸出結(jié)果看到，在每頁(yè)輸出結(jié)果上面有一行標(biāo)題，內(nèi)容為“TheSAS '標(biāo)題內(nèi)容title'951則輸出結(jié)果的標(biāo)題為“951班成績(jī)表”。要注意的是，TITLE語(yǔ)句是一個(gè)所謂的“全程”用另一個(gè)同樣的全程語(yǔ)句來(lái)修改它。比如，我們?cè)谏厦嬗肨ITLE語(yǔ)句指定了一個(gè)標(biāo)題，那么，這個(gè)標(biāo)題的作用將持續(xù)下去，雖然下一個(gè)過(guò)程沒(méi)有用TITLE語(yǔ)句指定標(biāo)題它也會(huì)出現(xiàn)在下一個(gè)過(guò)程的輸出中，例如在上面用TITLEC9501數(shù)據(jù)集的列表輸出加了標(biāo)題procmeansdata=sasuser.gpa;用全程語(yǔ)句FOOTNOTE可以為輸出加腳注，如：footnote'第三章例新的腳注，或用空FOOTNOTE語(yǔ)句取消腳注為止。optionsnopagenodatelinesize=78NOPAGE表示輸出不顯示頁(yè)號(hào)（PAGE則規(guī)定顯示頁(yè)號(hào)）NODATE表示不在每，LINESIZE=78個(gè)字符，PAGESIZE=60規(guī)定輸出每頁(yè)為60行，不足時(shí)用空行補(bǔ)齊。用BY語(yǔ)句分組處理使用帶有BY語(yǔ)句的過(guò)程之前一般用SORT過(guò)程對(duì)數(shù)據(jù)集按照分類變量排序。例如：procsortdata=c9501; procprintdata=c9501; AVG 92 98 AVG 89 80 計(jì)算總計(jì)和小計(jì)datainputname$amount;2015王思明2050;中指定要求和的變量名AMOUNT：procprintdata=bkmoneynoobs;sumamount; 20

BY語(yǔ)句SUM語(yǔ)句就可以既計(jì)算總和也計(jì)算分組小計(jì)。比如，我們除了要計(jì)算學(xué)生課外書(shū)總支出外還想分男、計(jì)算總支出，可以用下面的程序。注意由于數(shù)據(jù)集BKMONEY中沒(méi)有的信息，我們用了帶MERGE語(yǔ)句的數(shù)據(jù)步來(lái)橫向合并C9501和BKMONEY兩個(gè)數(shù)據(jù)集。procsortdata=c9501;byname;procsortdata=bkmoney;byname;datamergec9501bkmoney;byname;procsortdata=c9501new; procprintdata=c9501new; AVG 92 98 50 AVG 80 §3.3用TABULATE過(guò)程作統(tǒng)計(jì)量表格

PRINT過(guò)程可以制作列表，它列出所有觀測(cè)。當(dāng)觀測(cè)個(gè)數(shù)很多時(shí)，這樣的列表意義不 可以作出很復(fù)雜的表，其一般 頁(yè)維說(shuō)明，行維說(shuō)明，說(shuō)明/選項(xiàng);分別計(jì)算。VAR語(yǔ)句給出分析變量，分析變量必須是數(shù)值型（區(qū)間）變量，要計(jì)算的統(tǒng)計(jì)proctabulatedata=c9501new;class;vartable,amount; | | 用的值區(qū)分，AMOUNT在，它畫在列標(biāo)題中。占全部總和的百分比。例如，我們可以用如下程序求男、的數(shù)學(xué)、語(yǔ)文成績(jī)平均值及proctabulatedata=c9501new;class;varmathtable,(math)*(meanstd); | | | | | | 字。把上面的例子中的TABLE語(yǔ)句改成： all,(math)*(mean | | | | | 量和分析變量以構(gòu)成分類計(jì)算分析變量的統(tǒng)計(jì)量。例如，上面的TABLE語(yǔ)句寫成table all)*math*(mean | KEYLABEL關(guān)鍵字＝'proctabulatedata=c9501new;class;varmathtable( all),(math)*(meanstd);keylabelmean='平均值std=all='總計(jì)';label=''math='數(shù)學(xué)'='語(yǔ)文'; 數(shù) 語(yǔ) | | | | |總 §3.4SORT過(guò)程排序SASBYSORT過(guò)程據(jù)集C9501的觀測(cè)（行）按排序（分類可以用此程序：procsortdata=c9501; procsortdata=c9501; descendingavg;結(jié)果可以用procprint;run;顯示出來(lái)?！?.5MEANS過(guò)程、UNIVARIATE過(guò)程、FREQ過(guò)程計(jì)算均值、標(biāo)準(zhǔn)差等數(shù)字特征，而FREQ過(guò)程對(duì)離散變量計(jì)算取值頻數(shù)分布。procmeansdata=c9501;varmath;NStd55 如果使用UNIVARIATE過(guò)程則可以計(jì)算較多的統(tǒng)計(jì)量，例如我們分析中的變量GPA，可以用如procunivariatedata=sasuser.gpa;vargpa;Univariate CollegeGradePointAverageNSumStd-StdNum^=Num>Sgn100%6675%50%25%0%Q3-它顯示的各統(tǒng)計(jì)量已在1.3.7中作了介紹。其中T:Mean=0是用來(lái)檢驗(yàn)變量的平均值是否為0的我們想了解C9501中的分布情況，可以用：procfreqdata=c9501; Frequency

男33女25百分比。也可以對(duì)區(qū)間變量使用FREQ過(guò)程列出頻數(shù)分布，如tables“Tables|FrequencyTable”得到?！?.6用CORR過(guò)程計(jì)算相關(guān)系數(shù)proccorrdata=sasuser.gpa;varhsmhsshse;Correlation3'VAR'Variables:HSM SimpleStatisticsNStdSimple HighSchoolMathHighSchoolScienceHighSchoolEnglishPearsonCorrelationCoefficients/Prob>|R|underHo:Rho=0/N=HighSchoolMathHighSchoolScienceHighSchoolEnglish系數(shù)矩陣。比如，HSMHSS0.57569，這個(gè)相關(guān)系數(shù)為零的顯著性概相關(guān)系數(shù)，第二個(gè)數(shù)是檢驗(yàn)這兩個(gè)變量之間相關(guān)系數(shù)為0的檢驗(yàn)的顯著性概率值（p值§3.7用SAS/GRAPH繪圖SAS可以把存貯在SAS數(shù)據(jù)集中的數(shù)據(jù)以圖形的方式形象直觀地顯示出來(lái)。在SAS/GRAPH模塊的支持下，SAS可以作散點(diǎn)圖、曲線圖、直方圖、扇面圖、三維曲面圖、散點(diǎn)圖和曲線圖SASM的散點(diǎn)圖，只要用procgplotdata=sasuser.gpa;symboli=nonev=star;plotsatv*satm;結(jié)果顯示了一個(gè)GRAPHICS窗口，繪出了以SATV為縱軸、以SATM為橫軸的散點(diǎn)圖（1GPLOTPLOT語(yǔ)句指定繪圖用的變量。SYMBOL語(yǔ)SYMBOL2，SYMBOL3帶編號(hào)的相當(dāng)于SYMBOL1為了繪制連線，只要在procgplot;

1SATVSATM2symboli=joinv=star;plotco*datetime;符號(hào)可以在SYMBOL語(yǔ)句中用PLOT語(yǔ)句中指定多個(gè)因變量（自變量一般應(yīng)為同一，procgplotdata=sasuser.air;symbol1i=joinv=noneline=1;symbol2i=joinv=noneline=2plotco*datetime=1so2*datetime=2/overlay;SYMBOLSYMBOL語(yǔ)句指定LINE=1，表示SYMBOL語(yǔ)句指定了LINE=2表示線型為虛P(yáng)LOTSYMBOL語(yǔ)句的規(guī)定來(lái)畫，n對(duì)應(yīng)于SYMBOL語(yǔ)句的序號(hào)。SYMBOL語(yǔ)句的I=選項(xiàng)還可以取SPLINE表示在散點(diǎn)間連接連樣條曲線，

395procgplotdata=sasuser.gpa;symboli=rlcli95v=star;plotsatv*satm; |Graphics|GPLOT直方圖和扇形圖GCHART過(guò)程繪制直方圖、扇形例如要繪SASUSER.GPAGPA的分procgchartdata=sasuser.gpa;vbargpa;INSIGHT中繪制的直方圖有所不INSIGHT中橫軸標(biāo)的是區(qū)間的端點(diǎn)值?？梢岳L制分組的直方圖，例如按procgchartdata=sasuser.gpa;vbargpa/group=;結(jié)果見(jiàn)圖5在GCHART中用PIE語(yǔ)句繪制表示頻procgchartdata=sasuser.gpa; 結(jié)果見(jiàn)圖6。如果想顯示百分比值，只要在PIE語(yǔ)句中加入TYPE=PERCENT選項(xiàng)， /type=percent;。繪制三維直方圖。例如，在SASUSER.HOUSES有各種房子的情況，

4GPA56procgchart;blockstyle/group=bedrooms;圖形見(jiàn)圖7線

7我們可以用G3D過(guò)程繪制曲面圖形，用GCONTOUR繪制曲面的等高線圖。例如，我們想繪制一個(gè)二維正態(tài)分布曲面的圖形，假設(shè)(X，Y)服從聯(lián)合正態(tài)分布，11的為f(x,y)

ax2y22r2a(1r2

2a(1r2 dox=-3to3bydoy=-3*a2to3*a2by0.3*a2; (a*x*x+y*y-2*r*a2*x*y));keepxyz;procg3ddata=dnorm2;plotx*y=z;8GCONTOUR過(guò)程可以繪制曲面對(duì)應(yīng)的procgcontourdata=dnorm2;plotx*y=z;見(jiàn)圖9

89圖形的調(diào)整與輸出各繪圖過(guò)程中都指定了豐富的選項(xiàng)來(lái)調(diào)整圖形，讀者可以參考有關(guān)資料或查系統(tǒng)的幫TTLE