實(shí)驗(yàn)三SAS描述統(tǒng)計(jì)分析

試驗(yàn)三SAS描述統(tǒng)計(jì)分析質(zhì)的初步結(jié)論。據(jù)時(shí)可以給人留下深刻的印象。統(tǒng)計(jì)圖形的種類(lèi)很多，利用SAS可以便利的繪制常用的統(tǒng)計(jì)圖形。試驗(yàn)?zāi)康陌盐帐褂肧AS對(duì)數(shù)據(jù)作描述性統(tǒng)計(jì)分析的方法。把握SAS/GRAPH能，能用SAS的統(tǒng)計(jì)圖形對(duì)數(shù)據(jù)進(jìn)展描述性統(tǒng)計(jì)分析。試驗(yàn)內(nèi)容一、用INSIGHT計(jì)算統(tǒng)計(jì)量、繪制統(tǒng)計(jì)圖形二、用“分析家”計(jì)算統(tǒng)計(jì)量、繪制統(tǒng)計(jì)圖形三、編程實(shí)現(xiàn)描述性統(tǒng)計(jì)〔MEANS、UNIVARIATE、FREQ過(guò)程、編程繪制統(tǒng)計(jì)圖〔GPLOTGCHART〕試驗(yàn)指導(dǎo)一、用INSIGHT計(jì)算統(tǒng)計(jì)量【試驗(yàn)【試驗(yàn)3-1】按性別分別計(jì)算SASHELP.CLASS中身高的均值、標(biāo)準(zhǔn)差、中位數(shù)和其它四分位數(shù)，簡(jiǎn)潔分析學(xué)生身高的狀況。在INSIGHTSolutio〔解決方案Analysi〔分析InteractiveDataAnalysis〔交互式數(shù)據(jù)分析SAS/INSIGHTOpen”對(duì)話(huà)框，在對(duì)話(huà)框中選擇數(shù)據(jù)集：SASHELP.CLASS，單擊“Open〔翻開(kāi)INSIGHT3-1左所示。用DistributionAnalyz〔分析Distribution(Y〔分布Distribution對(duì)話(huà)框。在數(shù)據(jù)集CLASS的變量列表中，選擇height，然后單擊“Y”按鈕，height被選為分析變量，選擇sex，然后單擊“Group”按鈕，sex3-1右所示。 圖3-1 在INSIGHT中計(jì)算統(tǒng)計(jì)量Outpu〔輸出3-。表3-1 描述性統(tǒng)計(jì)選項(xiàng)選項(xiàng)MomentsQuintilesBasicConfidenceIntervalsTestsforLocation

說(shuō)明矩統(tǒng)計(jì)量分位數(shù)根本置信區(qū)間位置檢驗(yàn)

選項(xiàng)FrequencyCountsTestsforNormality

說(shuō)明頻數(shù)統(tǒng)計(jì)尺度的穩(wěn)健估量正態(tài)性檢驗(yàn)BoxPlot/MosaicPloChar3-2所示。圖3-2 描述性統(tǒng)計(jì)選項(xiàng)兩次單擊“OKheight按sex分組的各種矩統(tǒng)計(jì)量〔Moments〕和分位數(shù)Quantile，如圖3-3所示。結(jié)果分析91060.5889低于男63.91，女生的標(biāo)準(zhǔn)差5.01834.9379，說(shuō)明女生的身高較為分散一些。圖3-3分組顯示的各種矩統(tǒng)計(jì)量和分位數(shù)62.5、56.5和64.3，說(shuō)明身高大于或小于62.556.5，而另有四分之一的女生身超群過(guò)了64.3；男生的中位數(shù)、四分之一分位數(shù)和四分之三分位數(shù)分別為64.15、59.0和67.0，說(shuō)明身高大于或小于64.15的男生各占一半，有四分之一的男生身高不超過(guò)，而另有四分之一的男生身超群過(guò)了67.0。二、用“分析家”計(jì)算統(tǒng)計(jì)量【試驗(yàn)二、用“分析家”計(jì)算統(tǒng)計(jì)量【試驗(yàn)3-2500名學(xué)生的“概率統(tǒng)計(jì)學(xué)”課程的考試成績(jī)中，隨機(jī)地603-2所示〔數(shù)據(jù)集為：MYLIB.SY3-2，變量名為F1〕：表3-2 “概率統(tǒng)計(jì)學(xué)”考試成績(jī)6375839145819330728082838176678472588364936375997476959183618285834488726694687888719485827910090838884487280858087766296求出上述數(shù)據(jù)的平均值、方差、極差、中位數(shù)、四分之一分位數(shù)及四分之三分位數(shù)，并分析成績(jī)分布狀況。在“分析家”中翻開(kāi)MYLIB.SY3-2數(shù)據(jù)集Solution〔解決方案Analysi〔分析Analys〔分析家翻開(kāi)“分析家”窗口。Fil〔文件AsSASNameOpen…〔按SAS名稱(chēng)翻開(kāi)的對(duì)話(huà)框中，選擇MYLIB.SY3-2。用SummaryStatisticsStatistic〔統(tǒng)計(jì)Descriptiv〔描述性統(tǒng)計(jì)SummaryStatistics〔匯總統(tǒng)計(jì)量SummaryStatisticsF1Analysis”3-4左所示。圖3-4 選定分析變量與統(tǒng)計(jì)量單擊“Statistics”按鈕，翻開(kāi)“SummaryStatistics：Statistics”對(duì)話(huà)框。對(duì)話(huà)框中列3-3。3-3可計(jì)算統(tǒng)計(jì)量Mean均值Numberofobservations觀(guān)測(cè)數(shù)Standarddeviation標(biāo)準(zhǔn)差Numberofmissingvalues缺失數(shù)Standarderror標(biāo)準(zhǔn)誤Skewness偏度Variance方差Kurtosis峰度MinimumMaximumRange極差Student”sttCoefficientofvariationt檢驗(yàn)值p值變異系數(shù)Sum和Correctedsumofsquares校正平方和Median中位數(shù)Uncorrectedsumofsquares未校平方和選擇輸出統(tǒng)計(jì)量平均值、方差、極差，如圖3-4右所示。兩次單擊“OK”按鈕后，得3-5所示。圖3-5 描述性統(tǒng)計(jì)量用Distributions菜單項(xiàng)計(jì)算分位點(diǎn)在“分析家”窗口選擇主菜單“Statistic〔統(tǒng)計(jì)”→“Descriptiv〔描述性統(tǒng)計(jì)”→按鈕，選定分析變量。單擊“OK”按鈕，即可得到關(guān)于變量F13-6所示。圖3-6 矩統(tǒng)計(jì)量和根本統(tǒng)計(jì)測(cè)度另外，還有位置檢驗(yàn)、分位數(shù)以及極值觀(guān)測(cè)〔5個(gè)觀(guān)測(cè)，如圖3-7所示。圖3-7 位置檢驗(yàn)、分位數(shù)以及極值觀(guān)測(cè)結(jié)果分析學(xué)生成績(jī)的均值、方差和極差分別為77.72、194.07和70，方差和極差較大，說(shuō)明成績(jī)的分布很分散。中位數(shù)、四分之一分位數(shù)及四分之三分位數(shù)分別為81、72和86，說(shuō)明成績(jī)大于或小于81的學(xué)生各占一半，有四分之一的學(xué)生成績(jī)不超過(guò)72；而有四分之一的學(xué)生成86，從均值和中位數(shù)可以看出有一大半學(xué)生的成績(jī)超過(guò)了平均數(shù)。三、編程實(shí)現(xiàn)描述性統(tǒng)計(jì)三、編程實(shí)現(xiàn)描述性統(tǒng)計(jì)【試驗(yàn)3-3】編程實(shí)現(xiàn)描述性統(tǒng)計(jì)。1.1.用FREQ編制SASUSER.CLASS中age(年齡)的頻數(shù)表，代碼如下：procfreqdata=SASUSER.CLASS;tablesage;RUN; 圖3-8 頻數(shù)表2.用NEANS2.用NEANS按性別分別計(jì)算SASUSER.CLASS中體重的均值、標(biāo)準(zhǔn)差、極差、中位數(shù)、四分之一分位數(shù)和四分之三分位數(shù)，代碼如下：proc means data=SASUSER.CLASS mean std range median q1 var weight;class run;3-9所示。圖3-9 使用class語(yǔ)句分類(lèi)計(jì)算均值、標(biāo)準(zhǔn)差、極差、中位數(shù)和分位數(shù)或者proc sort data=SASUSER.CLASS;bysex;run;proc means data=SASUSER.CLASS meanstdrangemedianq1q3;var bysex;run;3-10所示。說(shuō)明:1)通常使用BY語(yǔ)句時(shí)要求數(shù)據(jù)集須按BY變量排序，使用CLASS語(yǔ)句無(wú)此要求。2)meansBYBY變量的每個(gè)值分別供給一個(gè)表，使用CLASS語(yǔ)句則將全部結(jié)果排列在一個(gè)表之中。3.用UNIVARIATE圖3-10 3.用UNIVARIATE計(jì)算SASUSER.CLASSage(年齡)的眾數(shù)、2.5%97.5%分位數(shù)。proc univariate data=SASUSER.CLASS;var age;output out=pct pctlpre=p pctlpts=2.5 run;proc print run;3-11所示，age(年齡)12，2.5%11，97.5%16。 圖3-11 眾數(shù)和分位數(shù)說(shuō)明:UNIVARIATE過(guò)程中的output語(yǔ)句表示將UNIVARIATE產(chǎn)生的統(tǒng)計(jì)量輸出到“out=”后定義的數(shù)據(jù)集中，本例數(shù)據(jù)集名稱(chēng)為pct，output語(yǔ)句中的選項(xiàng)pctlpre用于指定數(shù)pctlpts2.5%97.5%分位數(shù)。四、用INSIGHT四、用INSIGHT繪制統(tǒng)計(jì)圖形3-4】爭(zhēng)論Sasuser.gpagpa分?jǐn)?shù)的分布：畫(huà)出gpa的直方圖、盒形圖；通過(guò)直方圖、盒形圖簡(jiǎn)述不同性別的gpa數(shù)據(jù)的特點(diǎn)。INSIGHT中翻開(kāi)數(shù)據(jù)集Sasuser.gpa。1.1.繪制直方圖選擇菜單“Analyze〔分析”→“Chart(Y)〔條形圖Histogram/BarChart(Y)”對(duì)話(huà)框。gpagpa變量，然后單擊“Y”按鈕，gpa3-12所示。單擊“OK”按鈕，即可得到直方圖，如圖3-13左所示。單擊圖形框左下角處的三角形，在翻開(kāi)的菜單中alue〔值，即可在圖中顯示頻數(shù)值，如圖3-13右所示。

圖3-12 話(huà)框假設(shè)在“Histogram/BarChart(Y)”對(duì)話(huà)框中選定gpa變量作為分析變量，然后選擇sexGrou〔分組O”后可將gpa分?jǐn)?shù)按男女學(xué)生分別畫(huà)出3-14所示。 圖3-13 在INSIGHT中繪制直方圖2.繪制盒形圖圖3-14 不同性別gpa2.繪制盒形圖Analyz〔分析BoxPlot/MosaicPlot(Y〔馬賽克圖“BoxPlot/MosaicPlot(Y)”對(duì)話(huà)框。在數(shù)據(jù)集gpa的變量列表中，選擇gpa變量，然后單擊“Y”按鈕，gpa3-15左所示。Outpu〔輸出YAxisVertica〔垂直的Y軸”如圖3-15右所示。單擊“OK”按鈕返回對(duì)話(huà)框。 圖3-15 “BoxPlot/MosaicPlot(Y)”對(duì)話(huà)框單擊“OK”按鈕，即可得到盒形圖，如圖3-16左所示。單擊圖形框左下角處的三角形，在翻開(kāi)的菜單中選擇“Mean〔均值3-16右所示。圖3-16 在INSIGHT中繪制盒形圖按(1)中方法翻開(kāi)“BoxPlot/MosaicPlot(Y)”對(duì)話(huà)框，選擇gpa變量，單擊“Y”按鈕，再選擇sex變量，然后單擊“X”按鈕，如圖3-17左所示，可以將男女生gpa分?jǐn)?shù)的兩個(gè)盒形圖畫(huà)在同一個(gè)坐標(biāo)系下，如圖3-17右所示，便于查看男女生gpa分?jǐn)?shù)的分布狀況。圖3-17 畫(huà)在一個(gè)坐標(biāo)系下單擊圖3-17MeanValue在盒形圖上標(biāo)出男女生gpa分?jǐn)?shù)的均值、中位數(shù)、四分之一分位數(shù)和四分之三分位數(shù)等。圖3-18 顯示均值、中位數(shù)、四分之一分位數(shù)和四分之三分位數(shù)等3.結(jié)果分析3-133-16可以看到gpa分?jǐn)?shù)在均值兩邊的分布是不對(duì)稱(chēng)gpa分?jǐn)?shù)超過(guò)了平均值。這局部高分學(xué)生的分?jǐn)?shù)比較集中，由于盒形圖中均值線(xiàn)右邊局部較短，而其他低分局部的學(xué)生成績(jī)較為分散，gpagpa分?jǐn)?shù)是極端值，它不具有代表性。gpa分?jǐn)?shù)直方圖〔3-14〕和盒形圖〔3-17右〕可以看到，男女生gpa分?jǐn)?shù)的分布狀況類(lèi)似，都略微左偏。兩條中位數(shù)線(xiàn)位置沒(méi)有明顯區(qū)分，但女生中位數(shù)線(xiàn)左方局部更長(zhǎng)一些，說(shuō)明女生的gpa分?jǐn)?shù)低分局部更分散一些。五、用“分析家”繪制統(tǒng)計(jì)圖形五、用“分析家”繪制統(tǒng)計(jì)圖形3-5】某高校對(duì)學(xué)生逃課的狀況進(jìn)展抽樣調(diào)查的結(jié)果如表3-4所示：3-4學(xué)生逃課狀況表年級(jí)性別逃課的狀況年級(jí)性別逃課的狀況2023級(jí)男從不逃課2023級(jí)男從不逃課2023級(jí)男間或逃課2023級(jí)男從不逃課2023級(jí)男間或逃課2023級(jí)男從不逃課2023級(jí)男間或逃課2023級(jí)男從不逃課2023級(jí)男間或逃課2023級(jí)女從不逃課2023級(jí)男間或逃課2023級(jí)女從不逃課2023級(jí)男間或逃課2023級(jí)女從不逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)男間或逃課2023級(jí)女間或逃課2023級(jí)男間或逃課2023級(jí)女間或逃課2023級(jí)男間或逃課2023級(jí)女間或逃課2023級(jí)男間或逃課2023級(jí)女間或逃課2023級(jí)男間或逃課2023級(jí)女間或逃課2023級(jí)男間或逃課2023級(jí)女間或逃課2023級(jí)男間或逃課2023級(jí)女間或逃課2023級(jí)男間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)女間或逃課2023級(jí)男常常逃課2023級(jí)男常常逃課2023級(jí)女常常逃課2023級(jí)男常常逃課2023級(jí)女常常逃課2023級(jí)男常常逃課2023級(jí)女常常逃課2023級(jí)男常常逃課2023級(jí)女常常逃課2023級(jí)男常常逃課2023級(jí)女常常逃課2023級(jí)女常常逃課2023級(jí)女常常逃課2023級(jí)女常常逃課2023級(jí)女常常逃課2023級(jí)女常常逃課利用統(tǒng)計(jì)圖形從整體上對(duì)逃課狀況進(jìn)展分析。利用統(tǒng)計(jì)圖形對(duì)兩個(gè)年級(jí)逃課狀況進(jìn)展比較分析。1.數(shù)據(jù)預(yù)處理1.數(shù)據(jù)預(yù)處理由于“分析家”對(duì)漢字處理效果不是很好，用變量grade、sexcase分別代表年級(jí)、2023202320232023”分別NSA首先將表3-4處理成Excel〔sy4_2.xl3-19所示。sy4_2.xls導(dǎo)入為SAS數(shù)據(jù)集Mylib.sy4_2?；蛟凇胺治黾摇敝兄苯臃_(kāi)sy4_2.xls后存為數(shù)據(jù)集Mylib.sy4_2。以下操作總假定已經(jīng)在“分析家”中翻開(kāi)數(shù)據(jù)集2.利用統(tǒng)計(jì)圖形從整體上對(duì)逃課狀況進(jìn)展分析Mylib.sy4_2。2.利用統(tǒng)計(jì)圖形從整體上對(duì)逃課狀況進(jìn)展分析繪制三維餅圖

圖3-19 數(shù)據(jù)表sy4_2.xlsGraph〔圖形PieChar〔餅圖PieChar變量列表中的case，單擊“Chart3-D3-20左。 圖3-20 總體狀況三維餅圖再單擊“Options”按鈕，進(jìn)入“Options”窗口。單擊“SliceValues”選項(xiàng)卡，選中“Percent”單項(xiàng)選擇按鈕，如圖3-20右所示，兩次單擊“OK”按鈕，得到標(biāo)有百分比的立體三維餅圖，如圖3-21所示。繪制三維條形圖Graph〔圖形BarChart〔條形圖Vercita〔垂直BarChart”對(duì)話(huà)框，選擇變量列表中的case，單擊“Chart3-D選3-22上所示。

圖3-21 總體狀況三維餅圖單擊“Options”按鈕，進(jìn)入“Options”窗口，單擊“Statistics”選項(xiàng)卡，選中“DisplayStatistics”復(fù)選框，如圖3-22下所示。兩次單擊“OK”按鈕，得到標(biāo)有頻數(shù)的立體三維條3-23所示。結(jié)果分析

圖3-22 設(shè)置三維條形圖3-213-238%26%，間或逃課的占66，這說(shuō)明絕大多數(shù)學(xué)生逃過(guò)課，眾數(shù)為“間或逃課況，有一小局部常常逃課，這是值得我們留意的。3.利用統(tǒng)計(jì)圖形對(duì)兩個(gè)年級(jí)逃課狀況進(jìn)展比較分析圖3.利用統(tǒng)計(jì)圖形對(duì)兩個(gè)年級(jí)逃課狀況進(jìn)展比較分析分別繪制兩個(gè)年級(jí)的二維餅圖Graph〔圖形PieChar〔餅圖PieCharcase，單擊“Chart”按鈕，選定分析變量，選中2-D選項(xiàng)，再單擊“Options”按鈕，翻開(kāi)“Options”對(duì)話(huà)框。Percent”單項(xiàng)選擇OKVariables”按GRADE選為分組3-24所示。

圖3-24 Variable對(duì)話(huà)框O3-25所示。圖3-25 兩個(gè)年級(jí)二維餅圖分別繪制兩個(gè)年級(jí)的二維條形圖Graph〔圖形BarChar〔條形圖Vercita〔垂直“VercitalBarChart”對(duì)話(huà)框，選擇變量列表中的case，單擊“Chart”按鈕，選定分析變量，2-Dgrade變量，單擊“GroupBy3-26左所示。單擊“Options”按鈕，進(jìn)入“Options”窗口，單擊“Statistics”選項(xiàng)卡，選中“DisplayStatisticsOK”按鈕，得到標(biāo)有頻數(shù)的兩個(gè)年級(jí)的二維條形圖，如圖3-26右所示。 結(jié)果分析

圖3-26 兩個(gè)年級(jí)三維條形圖3-253-26右看出兩個(gè)年級(jí)學(xué)生逃課的狀況的差異，主要在兩個(gè)極端，即“從不逃課”和“常常逃課”上。2023202320%，而“從不逃課”12%20232023級(jí)嚴(yán)峻得多。4.利用統(tǒng)計(jì)圖形對(duì)男女生逃課狀況進(jìn)展比較分析兩個(gè)年級(jí)的逃課狀況中，眾數(shù)均為“間或逃課2023級(jí)的異眾比例為38，2023級(jí)的異眾比例為30%，這說(shuō)明“間或逃課”并不能很好地代表2023級(jí)逃課的狀況，相比之20234.利用統(tǒng)計(jì)圖形對(duì)男女生逃課狀況進(jìn)展比較分析分別繪制男女生的餅圖Graph〔圖形PieChar〔餅圖PieChar”對(duì)話(huà)框，選擇變量列表中的case，單擊“Chart”按鈕，選定分析變量，選中3-D選項(xiàng)，再單擊“Options”按鈕，進(jìn)入“Options”窗口，單擊“SliceValues”選項(xiàng)卡，選中“Percent”單項(xiàng)選擇按鈕，單擊“OK”返回，再單擊“Variable”選項(xiàng)Sex選為分組變量，3-27所示。兩次單擊“OK”分別得到標(biāo)有百分3-28所示。

圖3-27 Variable對(duì)話(huà)框 分別繪制男女生的三維條形圖

圖3-28 男女生三維餅圖Graph〔圖形BarChar〔條形圖Vercita〔垂直“VercitalBarChart”對(duì)話(huà)框，選擇變量列表中的case，單擊“Chart”按鈕，選定分析變量，再選擇sex變量，單擊“GroupBy”按鈕，選定分組變量。3-D選項(xiàng)，再單擊“Options”按鈕，進(jìn)入“Options”窗口，單擊“Statistics”選項(xiàng)卡，選中“DisplayStatistics”按鈕，兩次單擊“OK”按鈕，得到標(biāo)有頻數(shù)的男女生的兩個(gè)3-29所示。結(jié)果分析3-283-29。男女生均以“間或58%和74%，42%，女生的異眾比例為26%，這說(shuō)明用“間或逃課”代表男生逃課的明白，這是由于男生走“兩個(gè)極端”的不少。她們逃課的狀況。圖3-29 男女生三維條形圖六、編程繪制統(tǒng)計(jì)圖六、編程繪制統(tǒng)計(jì)圖使用GCHART3-6】使用GCHART3-5中有關(guān)圖形?！矆D形略〕goptionsftext=“宋體“htitle=1;title”總體狀況三維餅圖”;proc gchart data=Mylib.sy4_2;pie3d case/type=PCT;run;title”總體狀況三維條形圖”;procgchartdataMylib.sy4_2;vbar3d case/outside=FREQ;run;title”兩個(gè)年級(jí)二維餅圖”;procgchartdata=Mylib.sy4_2;pie case/group=grade type=PCT;run;title”兩個(gè)年級(jí)二維條形圖”;procgchartdataMylib.sy4_2;vbar case/group=grade outside=FREQ;run;title”男女生三維餅圖”;procgchartdata=Mylib.sy4_2;pie3d case/group=sex type=PCT;run;title”男女生三維條形圖”;procgchartdata=Mylib.sy4_2;vbar3d case/group=sex outside=FREQ;run;title””;使用GPLOT3-7GPLOT過(guò)程繪制SASUSER.AIRcodatetime的散點(diǎn)圖以及winddatetime的連線(xiàn)圖。goptionsftext=“宋體“htitle=1;proc gplot data=sasuser.air;plotco*datetime;title”codatetime的散點(diǎn)圖”;symbolv=doti=nonecv=orangeproc gplot data=sasuser.air;plotwind*datetime;title”winddatetime的連線(xiàn)圖”;symbolv=stari=joincv=redc