第5章-通用標(biāo)準(zhǔn)過程_第1頁
第5章-通用標(biāo)準(zhǔn)過程_第2頁
第5章-通用標(biāo)準(zhǔn)過程_第3頁
第5章-通用標(biāo)準(zhǔn)過程_第4頁
第5章-通用標(biāo)準(zhǔn)過程_第5頁
已閱讀5頁,還剩69頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

第5章通用標(biāo)準(zhǔn)過程5.1

繪制平面圖過程GCHART(CHART)

5.2計(jì)算相關(guān)系數(shù)過程CORR5.3SAS與dBASE間數(shù)據(jù)轉(zhuǎn)換過程DBF5.4簡(jiǎn)單統(tǒng)計(jì)分析過程MEANS5.5二維散點(diǎn)圖繪制過程GPLOT5.6輸出數(shù)據(jù)集過程PRINT5.7

排序過程SORT5.8STANDARD過程格式5.9

數(shù)據(jù)集轉(zhuǎn)置過程TRANSPOSE5.10兩組數(shù)據(jù)均值比較過程TTEST5.11繪制三維圖形過程G3D第5章通用標(biāo)準(zhǔn)過程第5章通用標(biāo)準(zhǔn)過程5.1繪制平面圖過程GCHART(CHART)

GCHART可以繪制條形圖、立體直方圖、餅圖和星圖。用這些圖可以直觀表示兩個(gè)或多個(gè)變量之間的關(guān)系。例5.1有一份單位職工情況數(shù)據(jù),包括姓名(name)、性別(sex)、年齡(age)、學(xué)歷(educa)(分為初(l)、中(m)和高(h)等教育三種情況)、月工資(salary)等。用圖形方式給出職工性別比例、學(xué)歷情況、男女工資比例、學(xué)歷與工資關(guān)系、學(xué)歷與性別等關(guān)系。建立數(shù)據(jù)集程序如下:DATAeg51;INPUTname$sex$agesalaryeduca$;LABELname=“姓名”sex=“性別”age=“年齡”salary="工資"educa="受教育情況";CARDS;李斯男20200初王五女25260初趙柳女28350中史奇男27350高朱巴男30290中劉久男35400中康實(shí)女32410高申山男31410高;第5章通用標(biāo)準(zhǔn)過程頻數(shù)條形圖橫向圖和縱向圖橫向條形圖HBAR例5.1.1頻數(shù)縱向條形圖VBAR例5.1.2a

例5.1.2b百分比縱向條形圖例5.1.3累積頻數(shù)縱向條形圖例5.1.4累積百分比縱向條形圖例5.1.5分組合計(jì)縱向條形圖例5.1.6分組均值縱向條形圖例5.1.7再分組縱向條形圖例5.1.8分組統(tǒng)計(jì)縱向條形圖例5.1.9

立體圖BLOCK例5.1.10

餅圖PIE例5.1.11

星圖STAR例5.1.12圖形類型語句:GCHART過程的一般形式PROCGCHART選項(xiàng)1; 圖形類型語句/選項(xiàng)2;[其它語句;]RUN;(1)選項(xiàng)1DATA=數(shù)據(jù)集名給出畫圖的數(shù)據(jù)集。LPI=aa是常數(shù)。該值給出餅圖與星圖的“圓度”比例,適當(dāng)調(diào)節(jié)此值,可以使屏幕顯示和打印的圖更“圓”一些。a的默認(rèn)值為6,通常a值為(每英寸行數(shù)/每英寸列數(shù))*10。例如,打印機(jī)每英寸為8行12列時(shí),LPI=6.6667。(2)圖形類型語句:VBAR縱向條形圖HBAR橫向條形圖BLOCK立體圖PIE餅圖STAR星圖(3)其他語句: BY變量表:給出分組變量,按該變量取值分層繪制,要求數(shù)據(jù)集已按該變量排序 PATTERNn:定義圖樣 AXISn:定義坐標(biāo)軸,控制坐標(biāo)軸的形狀和顏色 LEGEND:定義圖例(4)選項(xiàng)2顏色:RED(紅),BLACK(黑),WHITE(白),BLUE(藍(lán))GREEN(綠),ORANGE(橙),CYAN(青),GRAY(灰),PINK(粉紅),YELLOW(黃),BROWN(褐)方法:ARROW(左側(cè)或者右側(cè)),INSIDE(內(nèi)側(cè)),NONE(無),OUTSIDE(外側(cè))。VBAR和HBAR的選項(xiàng):ASCENDING以變量的值升序做圖AUTOREF給橫軸畫格線DESCENDING以變量的值降序做圖CAXIS=顏色給出坐標(biāo)軸顏色CFRAME=顏色給出條形圖的框架顏色CTEXT=顏色給出圖中文字的顏色COUTLINE=顏色給出條形圖輪廓線顏色DISCRETE數(shù)值變量取離散值FRAME給條形圖畫框架FREQ=變量給出頻數(shù)變量G100每一組內(nèi)的百分比為100%GAXIS=AXISn分組軸的說明GROUP=變量給定分組變量GSPACE=n給出分組圖圖間距LEGEND=LEGENDn圖例說明MAXIS=AXISn中點(diǎn)軸說明MINOR=n橫軸的最小刻度MIDPOINTS=值表給出中點(diǎn)值及范圍MISSING允許指定的中值變量有缺值NOAXISNOAXES不畫坐標(biāo)軸NOSYMBOLNOLEGEND不給出子組符號(hào)和圖例說明NOZEROS不畫零頻數(shù)中點(diǎn)PATTERNID=SUBGROUP圖樣隨SUBGROUP變化PATTERNID=GROUP圖樣隨GROUP變化PATTERNID=MIDPOINT圖樣隨MIDPOINT變化例5.1.13PATTERNID=BY圖樣隨BY變量變化AXIS=值給出軸的標(biāo)尺號(hào)AXIS=AXISn用第n個(gè)AXIS語句表示軸REF=值在軸上垂線SPACE=n給定條形圖間的距離SUBGROUP=變量以變量值繪制條形圖SUBVAR=變量以變量的和或者均值繪制條形圖TYPE=FREQ圖示頻數(shù)TYPE=CFREQ圖示累積頻數(shù)TYPE=PCT圖示百分比TYPE=CPCT圖示累積百分比TYPE=SUM圖示和TYPE=MEAN圖示平均數(shù)WIDTH=n條形圖寬LEVELS=n給定變量繪制VBAR條數(shù)一般情況下,HBAR條形圖右側(cè)顯示統(tǒng)計(jì)量,用選項(xiàng)NOSTAT則不顯示統(tǒng)計(jì)量。BLOCK的選項(xiàng):BLOCK的選項(xiàng)大部分與VBAR相同,只是描述上以立體圖形替代條形圖即可。BLOCKMAX=n立體圖的最大高度。COUTLINE=顏色立體圖輪廓線顏色。PIE的選項(xiàng):PIE的有些選項(xiàng)與VBAR相近,只是在敘說上將條形圖換成餅圖即可。而以下的選項(xiàng)是PIE有專門定義的。ACROSS=n給出橫向餅圖的個(gè)數(shù)。ANGLE=n第一個(gè)餅圖的起始角度。CFILL=顏色餅圖的顏色。COUTLINE=顏色餅圖輪廓線顏色。DOWN=n給出縱向餅圖的個(gè)數(shù)EXPLODE=數(shù)值表給出餅圖分塊值FILL=X網(wǎng)格線繪圖例5.1.14INVISIBLE=值給出不畫圖的數(shù)MATCHCOLOR文字和餅圖的顏色一致NOGROUPHEADING不輸出子組標(biāo)題NOHEADING不輸出標(biāo)題OTHER=n子組餅圖的大小PERCENT=方法百分比的標(biāo)記方式SLICE=方法中點(diǎn)值的標(biāo)記方式VALUE=方法值的標(biāo)記方式5.2計(jì)算相關(guān)系數(shù)過程CORR第5章通用標(biāo)準(zhǔn)過程CORR過程計(jì)算變量間的相關(guān)系數(shù),包括Pearson積矩相關(guān)系數(shù)和加權(quán)相關(guān)系數(shù),還能產(chǎn)生三個(gè)非參數(shù)的關(guān)聯(lián)測(cè)度(Spearman的秩序相關(guān),Kendall的tau-b,和Hoeffding的相關(guān)性度量D)。該過程也可計(jì)算偏相關(guān)(Pearson的偏相關(guān),Spearman的偏秩相關(guān)和Kendall的偏tau-b)和Cronbach系數(shù)及一些單變量的描述性統(tǒng)計(jì)量。PROCCORR選項(xiàng);[語句;]RUN;(1)選項(xiàng):數(shù)據(jù)集選項(xiàng)

DATA=數(shù)據(jù)集名指定要處理的數(shù)據(jù)集,缺省時(shí)為當(dāng)前數(shù)據(jù)集。

OUTH=數(shù)據(jù)集名產(chǎn)生包含Hoeffding統(tǒng)計(jì)量的數(shù)據(jù)集

OUTP=數(shù)據(jù)集名產(chǎn)生包含Pearson相關(guān)的數(shù)據(jù)集。

OUTK=數(shù)據(jù)集名產(chǎn)生包含Kendall相關(guān)的數(shù)據(jù)集。

OUTS=數(shù)據(jù)集名產(chǎn)生包含Spearman相關(guān)的數(shù)據(jù)集。

相關(guān)類型選項(xiàng)

HOEFFDING計(jì)算并輸出Hoeffding相關(guān)系數(shù)

KENDALL計(jì)算并輸出Kendall相關(guān)系數(shù)

PEARSON計(jì)算并輸出Pearson積矩相關(guān)系數(shù)

SPEARMAN計(jì)算并輸出Spearman相關(guān)系數(shù)默認(rèn)是PEARSON其它選項(xiàng)

VARDEF=參數(shù)參數(shù)是計(jì)算方差和協(xié)方差時(shí)的分母

N觀測(cè)個(gè)數(shù)

DF誤差自由度N-1WGT權(quán)重和

WDF權(quán)重和減1ALPHA計(jì)算Cronbach系數(shù)

COV當(dāng)使用PARTIAL語句時(shí),輸出偏協(xié)方差陣

CSSCP輸出修正平方和及叉積

NOCORR不輸出Pearson相關(guān)

SSCP輸出平方和及叉積

第5章通用標(biāo)準(zhǔn)過程(2)語句:

VAR變量表;(指定列)列出數(shù)據(jù)集中所有參與分析計(jì)算的數(shù)值變量。若無此語句,則數(shù)據(jù)集中的數(shù)值變量都參與運(yùn)算。WITH變量表;(指定行)分析兩組變量間相關(guān)性時(shí),用VAR語句和WITH語句分別列出兩組變量。PARTIAL變量;給出要求計(jì)算Pearson偏相關(guān)或Spearman偏等級(jí)次序相關(guān)或Kendall偏TAU-B的變量。WEIGHT變量;計(jì)算加權(quán)積矩相關(guān)系數(shù)時(shí)給出數(shù)據(jù)集中權(quán)重變量名。FREQ變量;給出數(shù)據(jù)集中代表觀測(cè)頻數(shù)的變量名。當(dāng)規(guī)定FREQ語句時(shí),若某觀測(cè)中FREQ變量的值為n,則輸入數(shù)據(jù)集中該觀測(cè)代表n個(gè)觀測(cè)。

第5章通用標(biāo)準(zhǔn)過程第5章通用標(biāo)準(zhǔn)過程x1x2x3…x10

x1c11c12c13…c110

x2……..

x3

…x10……………...c1010PROCCORR;RUN;x1x2x3…x5

x6c11c12c13…c110

x7……..

x8

x10……………...c1010…PROCCORR;VARX1-X5;WITHX6-X10;RUN;例5.2.1數(shù)據(jù)與說明見例10.1的程序,輸出見表5.12與表5.13。TITLE'法國(guó)12種家庭7種食品年均消費(fèi)的相關(guān)分析';DATAxf;INPUTclass$x1-x7;LABELx1='面包'x2='蔬菜'x3='水果'x4='肉'x5='家禽'x6='牛奶'x7='酒';CARDS;MA23324283541437526247427EM22935593881527567239258CA23727675621948927235433MA34065633411507544324407EM33866083961501558319363CA343884368923451148243341MA45346603671620638414407EM44606994841856762400416CA438578962123661149304282MA56557764231848759495486EM55849955482056893518319CA5515109788726301167561284;PROCCORRPEARSONSPEARMANKENDALLHOEFFDING;PROCCORR;VARx1x4x5x6;WITHx2x3x7;RUN;5.3SAS與dBASE間數(shù)據(jù)轉(zhuǎn)換過程DBF第5章通用標(biāo)準(zhǔn)過程PROCDBF選項(xiàng);RUN;選項(xiàng):DB2=文件邏輯名。文件邏輯名可以是當(dāng)前目錄下dBASEⅡ的以DBF為擴(kuò)展名庫(kù)的文件名(不包括擴(kuò)展名DBF);或者是用FILENAME語句定義的庫(kù)文件的文件邏輯名。本選項(xiàng)指定被轉(zhuǎn)換成SAS數(shù)據(jù)集的dBASE庫(kù)文件。DB3=文件邏輯名。這里的文件邏輯名是指dBASEⅢ庫(kù)文件的文件名或用FILENAME語句定義了的庫(kù)文件的文件邏輯名。其它要求同選項(xiàng)DB2=文件邏輯名。OUT=數(shù)據(jù)集名。該數(shù)據(jù)集名是由dBASE庫(kù)文件轉(zhuǎn)換后的SAS數(shù)據(jù)集名。DATA=數(shù)據(jù)集名。該數(shù)據(jù)集名是將要被轉(zhuǎn)換為dBASE庫(kù)文件的數(shù)據(jù)集名。當(dāng)無此選項(xiàng)時(shí),當(dāng)前數(shù)據(jù)集為被轉(zhuǎn)換數(shù)據(jù)集。說明:(1)dBASE庫(kù)文件轉(zhuǎn)換成SAS數(shù)據(jù)集時(shí),字符型值超過200個(gè)字符時(shí),在SAS數(shù)據(jù)集中只截取200個(gè)字符;邏輯型字段的值轉(zhuǎn)換為字符型值;日期型轉(zhuǎn)換為SAS的日期型;備注型被忽略。(2)SAS轉(zhuǎn)換成dBASE庫(kù)文件時(shí),數(shù)值變量轉(zhuǎn)換為dBASE的長(zhǎng)度為16的數(shù)值型數(shù)據(jù)。第5章通用標(biāo)準(zhǔn)過程 下列程序?qū)⒃搸?kù)文件轉(zhuǎn)換為SAS數(shù)據(jù)集WORK.AA。 FILENAMEbb"c:\myfile\eg5-3.DBF"; PROCDBFOUT=aaDB3=bb; RUN; PROCPRINT; RUN; 同樣,下列程序可以將SAS數(shù)據(jù)集又轉(zhuǎn)換成dBASE的庫(kù)文件: FILENAMEbb"c:\myfile\eg5-3.DBF"; PROCDBFDATA=work.aaDB3=bb; RUN; PROCPRINT; RUN;例5.3.1設(shè)c:\myfile\eg5-3.DBF,其字段名和數(shù)據(jù)如下: NAMESEXAGESALARY 張三m30400.00 李四f28380.00 王五m40500.005.4簡(jiǎn)單統(tǒng)計(jì)分析過程MEANSPROCMEANS選項(xiàng);[語句;]RUN;(1)選項(xiàng):DATA=數(shù)據(jù)集名指定被分析的數(shù)據(jù)集,缺項(xiàng)時(shí)為當(dāng)前數(shù)據(jù)集。MAXDEC=nn為輸出結(jié)果要保留的小數(shù)位,缺省時(shí)n=2。FW=nn為輸出統(tǒng)計(jì)量時(shí)的字段寬度,缺省時(shí)n=12。VARDEF=參數(shù)參數(shù)指出計(jì)算方差時(shí)的分母。DF自由度N-1,默認(rèn)值N觀測(cè)數(shù)WGT權(quán)重之和WDF權(quán)重之和減1第5章通用標(biāo)準(zhǔn)過程MEAN平均值STD標(biāo)準(zhǔn)差MIN最小值MAX最大值RANGE極差SUM總和VAR方差USS未修正的平方和CSS修正的平方和CV變異系數(shù)STDERR標(biāo)準(zhǔn)誤TT值,總體均值為零PRTT值的最大絕對(duì)值的概率SUMWGT加權(quán)變量值之和。(2)語句VAR變量表;

列入變量表中的數(shù)據(jù)集變量將被MEANS過程進(jìn)行分析。CLASS變量;

將數(shù)據(jù)集中的全體觀測(cè)以CLASS語句中的變量之值進(jìn)行分組分析處理。FREQ變量;

如果數(shù)據(jù)集中的某一變量的值是各觀測(cè)的頻數(shù),那么要用該語句指出該變量。OUTPUTOUT=數(shù)據(jù)集名;

將MEANS過程的結(jié)果輸出給定數(shù)據(jù)集名的數(shù)據(jù)集中去。

patient數(shù)據(jù)集包含4個(gè)變量(性別sex,年齡age,身高h(yuǎn)eight,體重weight),共237個(gè)觀測(cè)。不加參數(shù),默認(rèn)輸出每個(gè)數(shù)值變量非缺失值的頻數(shù)、均值、標(biāo)準(zhǔn)差、最大值和最小值。

procmeansdata=patient;run;使用class參數(shù),對(duì)變量進(jìn)行分組統(tǒng)計(jì),同時(shí)可以指定需要匯總統(tǒng)計(jì)的值。procmeansdata=patientnrangesumvarmeannmissmedian;classsex;run;使用output命令,指定輸出結(jié)果至數(shù)據(jù)集。procmeansdata=patient;classsex;varheight;outputout=htstatn=sum=mean=nmiss=median=/autoname;run;procprintdata=htstat;run;autoname是指定sas自動(dòng)命名匯總變量名稱。_TYPE_=0:全體。_TYPE_=1:按sex分類。若class后面有兩個(gè)變量,如sexage,則_TYPE_會(huì)出現(xiàn)4個(gè)值。0-整體,1-sex,2-age,3-sex*age第5章通用標(biāo)準(zhǔn)過程利用MEANS過程進(jìn)行成對(duì)數(shù)據(jù)比較 數(shù)據(jù)源于南京農(nóng)業(yè)大學(xué)主編《田間試驗(yàn)和統(tǒng)計(jì)方法》第77頁表4.4,程序?yàn)椋?TITLE“成對(duì)數(shù)據(jù)比較”; DATAa; INPUTgroup$x1x2; x3=x1-x2; LABELgroup=“參加者”x1=“用藥前”x2=“用藥后”x3=“x1-x2”; CARDS; 11025 21312 3814 4315 5512 62027 7618 ; 輸出結(jié)果:

表5.16成對(duì)數(shù)椐比較MEANS過程分析變量:x3=x1-x2均值標(biāo)準(zhǔn)誤差t值Pr>|t|----------------------------------------------------8.28571431.9965957-4.150.0060---------------------------------------------------PROCMEANSMEANSTDERRTPRT;VARx3;RUN;5.5二維散點(diǎn)圖繪制過程GPLOT(PLOT)

第5章通用標(biāo)準(zhǔn)過程 在一個(gè)數(shù)據(jù)集中,指定二變量作為縱橫軸,繪制高質(zhì)量的平面散點(diǎn)圖。例如在主成分分析后,用第一主成分和第二主成分繪制散點(diǎn)圖等。例5.5.1運(yùn)行下列程序:OPTIONSPS=25;DATAa;INPUTno$xyz;CARDS;a-4160b-397c-2412d-1115e0016f1115g2412h397i4160;PROCGPLOT;PLOTy*x;RUN;圖5.12PLOTy*x圖5.13PLOTx*y將程序中的PLOTy*x;換成PLOTx*y;將產(chǎn)生圖5.13的圖形。圖5.14PLOTy*x=“#”將PLOTy*x;換為PLOTy*x=“#”;后的圖形為圖5.14。將PLOTy*x;換為PLOTy*x=“#”z*x=“a”/OVERLAY;后的圖形為圖5.15。圖5.15PLOTy*x=“#”z*x=“a”/OVERLAYPROCGPLOT選項(xiàng)1;PLOT圖形說明1圖形說明2.../選項(xiàng)2;[其他語句;]RUN;OPTIONSPS=mLS=n;改變頁面的大小,m和n是正整數(shù),m為指定的每頁行數(shù),n為每行的字符數(shù)(64<n<256,15<n<32767),在SAS輸出中,默認(rèn)一頁為21行78列。

(1)選項(xiàng)1 DATA=數(shù)據(jù)集名指定繪圖數(shù)據(jù)集(2)圖形說明縱軸變量*橫軸變量[=變量或“單字符”]例如:(1)y*x/*y為縱軸,x為橫軸,圖形用默認(rèn)字符“+"形成。*/(2)x*y/*x為縱軸,y為橫軸,圖形用默認(rèn)字符“+"形成。*/(3)y*x=“#”/*y為縱軸,x為橫軸。圖形符號(hào)是“”*/(4)y*x=“#”z*x=“a”/*同一圖內(nèi)畫兩個(gè)散點(diǎn)圖,第一個(gè)圖用字符“”形成,第二個(gè)圖用字母“+”形成*/第5章通用標(biāo)準(zhǔn)過程

(3)選項(xiàng)2

VAXIS=值表給出縱軸的刻度標(biāo)記信息

HAXIS=值表給出橫軸的刻度標(biāo)記信息

值表可以是:用空格分隔的多個(gè)升序排列的數(shù)值。首刻度TO末刻度BY間隔。例如:PLOTy*x/VAXIS=5TO50BY5OVERLAY

將PLOT后的多個(gè)圖形說明給定的散點(diǎn)圖繪制于一個(gè)坐標(biāo)系下。

CHREF=顏色給出水平參考線的顏色CVREF=顏色給出垂直參考線的顏色HREF=值給出水平參考線的位置LHREF=值給出線型,“值”的范圍是1至46VREF=值給出垂直參考線的位置LVREF=值給出線型,“值”的范圍是1至46CAXI=顏色給出坐標(biāo)軸的顏色CFRAME=顏色給出所圍區(qū)域的顏色CTEXT=顏色給出文字的顏色FRAME畫框架VMINOR=n給出主標(biāo)記之間的次標(biāo)記的數(shù)VZERO垂直軸上的第一個(gè)標(biāo)記是0VREVERSE垂直軸上的次序顛倒(4)其他語句BY變量表;系統(tǒng)依“變量表”中的變量取值進(jìn)行分組顯示圖形。SYMBOLn選項(xiàng)3;給出繪圖的符號(hào)、顏色及各點(diǎn)之間的連接方式,n為1到255的正整數(shù),給定SYMBOL語句編號(hào)。

“選項(xiàng)3”取值:

C=顏色給出繪圖符號(hào)和連線顏色CV=顏色給出繪圖符號(hào)顏色CI=顏色給出連線顏色W=數(shù)值給出連線的寬度L=n給出連線的線型,1≤n≤46V=字符任一鍵盤字符H=數(shù)值給出字符的高F=字型給出字符的字型I=NONE點(diǎn)間不連線I=JOIN點(diǎn)間連直線I=SPLINE點(diǎn)間連曲線I=RXXXXXXX曲線擬合繪圖的選項(xiàng),第一個(gè)X取值為L(zhǎng)、Q和C,分別表示線性回歸擬合、二次回歸擬合和三次回歸擬合;第二個(gè)X的取值為零或非零分別表示有無截距;第三、第四和第五個(gè)X的取值為每個(gè)預(yù)測(cè)值的置信限CLI,或者為平均預(yù)測(cè)值的置信限;第六和第七個(gè)X的取值可以是90、95和99,表示不同置信水平的值。

例5.5.2以下程序段繪制曲線圖(圖5.16)PROCGPLOT;PLOTy*x;SYMBOLCI=BLUEI=SPLINE;RUN;圖5.16PLOT繪制曲線例5.5.3繪制Sinx在0到2π內(nèi)的圖形 OPTIONSPS=30; DATAa; DOx=0TO8*ATAN(1)BY0.05; y=SIN(x); OUTPUT; END; PROCGPLOT; PLOTy*x=“*”; RUN;輸出如圖5.17圖5.17y=SIN(x)例5.5.4goptionsreset=globalgunit=pctbordercback=whitecolors=(blackbluered)ftitle=swissbftext=swisshtitle=6htext=3;datacitytemp;inputmonthfarencity$;cards;140.5Raleigh112.2Minn152.1Phoenix...moredatalines...1241.2Raleigh1218.6Minn1252.5Phoenix;title1'AverageMonthlyTemperature';footnote1j=l'Source:1984AmericanExpress';footnote2j=l'AppointmentBook'j=r'GR21N08(a)';symbol1interpol=joinvalue=dotheight=3;/*將三個(gè)城市十二個(gè)月的氣溫圖顯示在同一幅圖中*/procgplotdata=citytemp;plotfaren*month=city/hminor=0;run;/*將三個(gè)城市十二個(gè)月的氣溫圖顯示在三幅圖中*/procgplotdata=citytemp;plotfaren*month;bycity;run;

plotfaren*month=cityPROCPLOT選項(xiàng)1;PLOT圖形說明1圖形說明2.../選項(xiàng)2;RUN;(1)選項(xiàng)1 DATA=數(shù)據(jù)集名指定繪圖數(shù)據(jù)集 UNIFORM分組畫圖時(shí),各組都用同一刻度VPERCENT=數(shù)值用給出的數(shù)值控制圖形在頁面上縱向的分布

例如:50是輸出兩個(gè)大小相同的圖形;502525是輸出三個(gè)圖形,一個(gè)占二分之一頁面,另兩個(gè)各占四分之一;330是輸出一個(gè)圖形,占三分之一頁面。HPERCENT=數(shù)值用給出的數(shù)值控制圖形在頁面上橫向的分布NOMISS含有缺值的觀測(cè)不參與計(jì)算(2)圖形說明縱軸變量*橫軸變量[=變量或“單字符”]例如:(1)y*x/*y為縱軸,x為橫軸,圖形用默認(rèn)字符“A"形成。*/(2)x*y/*x為縱軸,y為橫軸,圖形用默認(rèn)字符“A"形成。*/(3)y*x=“#”/*y為縱軸,x為橫軸。圖形符號(hào)是“#”*/(4)y*x=“#”z*x=“a”/*同一圖內(nèi)畫兩個(gè)散點(diǎn)圖,第一個(gè)圖用字符“#”形成,第二個(gè)圖用字母“a”形成*/(3)選項(xiàng)2OVERLAY

將PLOT后的多個(gè)圖形說明給定的散點(diǎn)圖繪制于一個(gè)坐標(biāo)系下。CONTOUR=n

其中n為1到10的正整數(shù)。選取本選項(xiàng),可以繪制二元函數(shù)的平面等高值圖,n稱為灰度級(jí),也叫分辨率。n值越大,圖的層次越多。(4)其他語句BY變量表;系統(tǒng)依“變量表”中的變量取值進(jìn)行分組顯示圖形。

例5.5.5 OPTIONSPS=60; DATAa; DOx=0TO500BY5; DOy=0TO500BY10; z=50+x/105*x**2/10000+y/105*y**2/10000+4*x*y/1000;OUTPUT; END; END; PROCPLOT; PLOTy*x=z/CONTOUR=10; RUN;圖5.18等高值圖ProcGcontour把三維的關(guān)系表現(xiàn)在二維的坐標(biāo)軸上,線或區(qū)域代表Z對(duì)應(yīng)的在平面坐標(biāo)軸上(x,y)的值。(x,y)代表湖的面積,而Z代表湖的深度,不同的深度值用唯一的一種顏色的線條表示。配合symbol語句可以更多樣化的顯示圖形。等高線繪制過程GCONTOUR下列程序的數(shù)據(jù)集與例5.5.5相同(圖5.18),繪圖過程用GCONTOUR,輸出如圖5.19。DATAa; DOx=0TO500BY5; DOy=0TO500BY10; z=50+x/10-5*x**2/10000+y/10-5*y**2/10000+4*x*y/10000; OUTPUT; END; END;PROCGCONTOUR;PLOTy*x=z;RUN;等高線繪制過程GCONTOUR將上例的繪圖過程步用以下程序段替換,輸出結(jié)果如圖5.20。PROCGCONTOUR;PLOTy*x=z/PATTERNCOUTLINE=GRAYCTEXT=GREEN;RUN;圖5.19GCONTOUR繪制等高線圖圖5.20使用PATTERN的GCONTOUR繪圖5.6輸出數(shù)據(jù)集過程PRINT

將一個(gè)已經(jīng)產(chǎn)生了的數(shù)據(jù)集輸出到OUTPUT窗口,可以輸出全部變量或部分變量的全部觀測(cè)值,也可以輸出數(shù)值變量的總和。PROCPRINT選項(xiàng);[語句;]RUN;(1)選項(xiàng)DATA=數(shù)據(jù)集名指定輸出數(shù)據(jù)集名。若無此選項(xiàng),則輸出當(dāng)前數(shù)據(jù)集。LABEL

在輸出的第一行中,用變量的標(biāo)號(hào)代替變量名。DOUBLE隔行輸出。NOOBS不輸出觀測(cè)的序號(hào)。N在輸出的最后一行給出觀測(cè)的個(gè)數(shù)。(2)語句

VAR變量表;只輸出變量表中給定的變量

SUM變量表;輸出變量表中所列變量的和第5章通用標(biāo)準(zhǔn)過程例5.6.1用例5.1的數(shù)據(jù)步程序,加上如下過程步:

PROCPRINTNLABELNOOBS;SUMagesalary;RUN;則輸出為:姓名性別年齡工資受教育情況李斯男20200初史奇男20350高王老五女25260初趙柳女28350中朱巴男30290中申山男31410高康實(shí)女32410高劉久男35400中========2212670N=85.7排序過程SORTPROCSORT選項(xiàng);BY選項(xiàng)1變量1...選項(xiàng)n變量i;RUN;(1)選項(xiàng)DATA=數(shù)據(jù)集名指定被排序的數(shù)據(jù)集。OUT=數(shù)據(jù)集名給出排序后的新的數(shù)據(jù)集名,若無此項(xiàng),則排序后的數(shù)據(jù)集用原數(shù)據(jù)集名。(2)BY語句;用來指定排序變量,選項(xiàng)n是DESCENDING,當(dāng)選此項(xiàng)時(shí),其后的排序變量將以從大到小的順序進(jìn)行比較排序,若無此選項(xiàng),則以升序的方式排序。第5章通用標(biāo)準(zhǔn)過程

例5.7.1將例5.1中的數(shù)據(jù)步,填上如下過程步:

PROCSORTOUT=eg51new;BYagesalary;RUN;輸出結(jié)果為:姓名性別年齡工資受教育情況李斯男20200初史奇男20350高王老五女25260初趙柳女28350中朱巴男30290中申山男31410高康實(shí)女32410高劉久男35400中5.8數(shù)據(jù)標(biāo)準(zhǔn)化過程STANDARD第5章通用標(biāo)準(zhǔn)過程PROCSTANDARD選項(xiàng);[語句;]RUN;(1)選項(xiàng)

DATA=數(shù)據(jù)集名指定被標(biāo)準(zhǔn)化的數(shù)據(jù)集名

OUT=數(shù)據(jù)集名指定一個(gè)數(shù)據(jù)集名,其內(nèi)包含的數(shù)據(jù)為標(biāo)準(zhǔn)化了的變量值及不要求標(biāo)準(zhǔn)化的變量值。

VARDEF=參數(shù)給定方差計(jì)算時(shí)分母選取的方式,可選參數(shù):DF自由度N-1,默認(rèn)值N觀測(cè)數(shù)WGT權(quán)重之和WDF權(quán)重之和減1 M=mm為指定的標(biāo)準(zhǔn)化均值,通常的標(biāo)準(zhǔn)化中M取0。

S=ss為指定的標(biāo)準(zhǔn)化標(biāo)準(zhǔn)差,通常的標(biāo)準(zhǔn)化中S取1。若無M=和S=兩個(gè)選項(xiàng)時(shí),標(biāo)準(zhǔn)化后的數(shù)據(jù)與原數(shù)據(jù)一樣。

REPLACE有此選項(xiàng)時(shí),缺值數(shù)據(jù)將用均值代替。(2)語句

VAR變量表;出現(xiàn)在變量表中的變量都將被標(biāo)準(zhǔn)化。若無此語句時(shí),數(shù)據(jù)集中的全部數(shù)值變量將被標(biāo)準(zhǔn)化。

FREQ變量;變量是數(shù)據(jù)集中的頻數(shù)變量。

WEIGHT變量;變量是數(shù)據(jù)集中的權(quán)重變量。例5.8.1

將下列數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化并產(chǎn)生新的數(shù)據(jù)集B。XYZ0.1318000.3717000.222500程序?yàn)椋?DATAa; INPUTxyz; CARDS; 0.131800 0.371700 0.222500 ; PROCSTANDARDOUT=bM=0S=1; RUN; PROCPRINTDATA=b; RUN;

輸出結(jié)果如下:NAMEMEANSTDNX0.20.13Y42.64575131113Z2000435.889894353OBSXYZ

1-1.00000-0.37796-0.4588321.000001.13389-0.688253-0.00000-0.755931.147085.9數(shù)據(jù)集轉(zhuǎn)置過程TRANSPOSE第5章通用標(biāo)準(zhǔn)過程PROCTRANSPOSE選項(xiàng);[語句;]RUN;(1)選項(xiàng)

DATA=數(shù)據(jù)集名指定被轉(zhuǎn)置的數(shù)據(jù)集名,無此選項(xiàng)時(shí),當(dāng)前數(shù)據(jù)集被轉(zhuǎn)置。

PREFIX=SAS名數(shù)據(jù)集被轉(zhuǎn)置后,原數(shù)據(jù)集的變量變成新數(shù)據(jù)集的觀測(cè),而原來的觀測(cè)變成新數(shù)據(jù)集的變量,那么,在新數(shù)據(jù)集中的變量名是什么呢?本選項(xiàng)中的“SAS”名就是為新變量名提供一個(gè)變量名的前綴,不同的變量用前綴“SAS名”后接自然數(shù)來命名。OUT=數(shù)據(jù)集名給定轉(zhuǎn)置后的新數(shù)據(jù)集名,若無此項(xiàng),則新數(shù)據(jù)集名為DATAn。NAME=SAS名SAS名將作為轉(zhuǎn)置后的數(shù)據(jù)集的一個(gè)變量名,其值為原數(shù)據(jù)集中的變量名,若無此選項(xiàng)時(shí),原數(shù)據(jù)集的變量名構(gòu)成的新數(shù)據(jù)集中變量名為_NAME_的變量。(2)語句VAR變量表;

變量表中給出原數(shù)據(jù)集中參與轉(zhuǎn)置的變量,若無此語句,則原數(shù)據(jù)集中的全部數(shù)值型變量將被轉(zhuǎn)置。例5.9.1下列程序?qū)崿F(xiàn)了由數(shù)據(jù)行中的數(shù)據(jù)產(chǎn)生的數(shù)據(jù)集的轉(zhuǎn)置,并且輸出原數(shù)據(jù)集和轉(zhuǎn)置后的數(shù)據(jù)集,程序如下:

DATAa;INPUTA$XYZ;CARDS;a123b456c789d011;PROCPRINT;PROCTRANSPOSEOUT=b;RUN;PROCPRINTDATA=b;RUN;

輸出結(jié)果:

原數(shù)據(jù)集OBSAXYZ1a1232b4563c7894d011轉(zhuǎn)置后數(shù)據(jù)集OBS_NAME_COL1COL2COL3COL41X14702Y25813Z3691

proctransposename=ovnameprefix=newout=b;

varAXYZ;

run;Obsovnamenew1new2new3new41Aabcd2X14703Y25814Z36915.10兩組數(shù)據(jù)均值比較過程TTESTPROCTTEST選項(xiàng)CLASS變量;[其它語句;]RUN;(1)選項(xiàng)DATA=數(shù)據(jù)集名指定要分析的數(shù)據(jù)集名,若無此選項(xiàng),當(dāng)前數(shù)據(jù)集將被處理。(2)CLASS變量;

本語句是必須的,指定區(qū)分兩個(gè)樣本觀察的變量,該變量取值只能是兩個(gè)不同的值,可以是字符型,也可以是數(shù)值型。(3)其它語句VAR語句,本語句是可選的,指明參與分析的變量,若無此語句,數(shù)據(jù)集中全部數(shù)值變量將被分析。輸出結(jié)果說明

(1)被分析變量名(Variable)(2)分類變量及其水平(no)(3)觀察個(gè)數(shù)(不含缺項(xiàng))(N)(4)均值(Mean)(5)標(biāo)準(zhǔn)差(StdDev)(6)標(biāo)準(zhǔn)誤差(StdErr)(7)檢驗(yàn)兩組數(shù)據(jù)均值相等假設(shè)的T近似值(Pooled,tValue)(8)近似自由度(Satterthwaite,DF)(9)大于t絕對(duì)值的雙側(cè)顯著性概率(Satterthwaite,Pr>|t|)(10)檢驗(yàn)兩組數(shù)據(jù)均值相等假設(shè)的T值(Satterthwaite,tValue)(11)自由度(PooledDF)(12)大于t絕對(duì)值的雙側(cè)顯著概率(Pooled,Pr>|t|)(13)F'統(tǒng)計(jì)量(FoldedF,FValue)(14)每組中的自由度(FoldedF,DenDF)(15)大于F'值的概率(FoldedF,Pr>F)

例5.10.1數(shù)據(jù)源于南京農(nóng)業(yè)大學(xué)主編<<田間試驗(yàn)和統(tǒng)計(jì)方法>>一書74頁表4.2。程序如下: DATAa; INPUTno$x@@; CARDS; x1800x2900 x1840x2880 x1870x2890 x1920x2890 x1850x2840 ; PROCTTEST; CLASSno; VARx; RUN;TheTTESTProcedureStatisticsLowerCLUpperCLLowerCLUpperCLVariablenoNMeanMeanMeanStdDevStdDevStdDevStdErrxx1801.45856910.5526.32143.932126.2419.647xx25850.88880909.1214.05123.45267.39110.488xDiff(1-2)-75.36-2427.35723.78535.21467.46122.271T-TestsVariableMethodVariancesDFtValuePr>|t|

xPooledEqual8-1.080.3126xSatterthwaiteUnequal6.11-1.080.3219EqualityofVariancesVariableMethodNumDFDenDFFValuePr>FxFoldedF443.510.25155.11繪制三維圖形過程G3DG3D過程能夠繪制網(wǎng)格三維曲面和散點(diǎn)三維曲面等三維圖形。PROCG3D選項(xiàng)1;PLOTY*X=Z/選項(xiàng)2;SCATTERY*X=Z/選項(xiàng)3;[其它語句;]RUN;在G3D過程中,PLOT語句和SCATTER語句至少出現(xiàn)其中之一。PLOT用于繪制網(wǎng)格三維曲面圖;SCATTER用于繪制散點(diǎn)三維曲面圖。Y*X=Z:繪圖說明。其中Y、X和Z是三個(gè)變量,分別表示三維坐標(biāo)系下的Y軸、X軸和Z軸。(1)選項(xiàng)1與GPLOT過程的“選項(xiàng)1”相同。(2)選項(xiàng)2

ANNODATE=數(shù)據(jù)集名給定一個(gè)注釋數(shù)據(jù)集。CBOTTOM=顏色給出曲面底面顏色。CTOP=顏色給出曲面上面顏色。ROTATE=角度值列表給出一個(gè)或多個(gè)繞Z軸旋轉(zhuǎn)的角度值。SIDE給曲面畫側(cè)斷面。TILT=角度值列表給出一個(gè)或多個(gè)相對(duì)Y軸的傾斜角度值。XYTYPE=n當(dāng)n=1時(shí),畫圖的線與X軸平行;當(dāng)n=2時(shí),畫圖的線與Y軸平行;當(dāng)n=3時(shí),用網(wǎng)格畫圖且網(wǎng)格線分別與X軸和Y軸平行。CAXIS=顏色給出軸的顏色。CTEXT=顏色給出文字的顏色。GRID在軸的刻度上畫參考線。NOAXIS不畫軸及刻度等。NOLABEL不畫刻度等。XTICKNUM=nn≥2,默認(rèn)值是4。給出X軸上刻度數(shù)。YTICKNUM=nn≥2,默認(rèn)值是4。給出Y軸上刻度數(shù)。ZTICKNUM=nn≥2,默認(rèn)值是4。給出Z軸上刻度數(shù)。ZMAX=數(shù)值給定變量Z的最大值。ZMIN=數(shù)值給定變量Z的最小值。(3)選項(xiàng)3對(duì)于選項(xiàng):ANNOTATE、ROTATE、TILT、CAXIS、CTEXT、GRID、NOAXIS、NOLABEL、XTICKNUM、YTICKNUM、ZMAX、ZMIN、ZTICKNUM的意

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論