某學院SAS基礎(chǔ)與金融計算

SAS基礎(chǔ)與金融計算游家興廈門大學經(jīng)濟學院計統(tǒng)系1本節(jié)講授內(nèi)容SAS圖表過程Procchart過程Procgchart過程2統(tǒng)計圖對統(tǒng)計資料圖形化處理就形成統(tǒng)計圖，運用統(tǒng)計圖形能直觀地描述統(tǒng)計分布的特征和變量之間的關(guān)系。統(tǒng)計圖有多種表現(xiàn)形式：1、直方圖2、餅圖3、折線圖4、曲線圖3回顧Procunivariatedata=studentplot;varheight;Run;獲得三個圖:1、莖葉圖2、箱線圖3、正態(tài)概率圖4統(tǒng)計分組的相關(guān)概念根據(jù)統(tǒng)計研究的目的和客觀現(xiàn)象的內(nèi)在特點，按某個標志或幾個標志把被研究的總體劃分為若干個不同性質(zhì)的組，稱為統(tǒng)計分組。分組的性質(zhì)來看：兼有分和合雙重含義。5統(tǒng)計分組的步驟1.確定數(shù)據(jù)的全距（R）：R＝xmax-xmin（即最大變量值與最小變量值之差）2.確定組限、組距和組數(shù)（1）組限分為上限和下限，上限是各組的最大變量值，下限是各組的最小變量值。（2）各組上限與下限之差稱為該組的組距。（3）組數(shù)是指某個變量數(shù)列劃分為多少個組。說明：組數(shù)的多少與組距的大小成反比關(guān)系，組距越大，則可分的組就越少。6分組，應(yīng)注意如下兩個問題：首先，分組時各組數(shù)量界限的確定必須能反映事物質(zhì)的差別——體現(xiàn)量變到質(zhì)變的轉(zhuǎn)化。其次，應(yīng)根據(jù)被研究的現(xiàn)象總體的數(shù)量特征，采用適當?shù)姆纸M形式，確定相宜的組距、組限——主要解決分多少組的問題。7組數(shù)的確定：與組距大小成反比以下是美國學者斯特杰斯經(jīng)驗公式：

n=1+3.322lgN

（n:組數(shù)；N：總體單位數(shù)）

d=R/n=（Xmax-Xmin）/n（n為組數(shù)，d為組距，R為全距）8根據(jù)上述公式，可以得出如下的組數(shù)參考標準：分組組數(shù)參考表N15-2425-4445-8990-179180-359n56789上述數(shù)據(jù)僅供參考，不能生搬硬套，實際分組時要依所研究資料的性質(zhì)而定。9SAS制圖過程SAS統(tǒng)計制圖過程均可以實現(xiàn)對樣本分布特征的圖形表示，一般情況下可以使用的有chart過程、plot過程、gchart過程和gplot過程。前兩個和后兩個只有一個字母“g”（代表graph）的差別，其實它們之間（只差一個字母g的過程之間）的統(tǒng)計描述功能是相同的，區(qū)別僅在于繪制出的圖形的復(fù)雜和美觀程度。10此外，還有一些較為復(fù)雜的制圖命令，如運用G3D過程繪制曲面圖，運用GCONTOUR過程繪制曲面對應(yīng)的等高線圖。例如，我們想繪制一個二維正態(tài)分布曲面的圖形，假設(shè)(X，Y)服從聯(lián)合正態(tài)分布，其均值都是0，方差分別為1和a，相關(guān)系數(shù)為r。11chart（或plot）與gchart（或gplot）chart過程和plot過程繪制的圖形類似于我們用文本字符堆積起來的圖形，只能概括地反映出資料分布的大體形狀，實際上這兩個過程繪制的圖形并不能稱之為圖形，因為它沒有涉及一般意義上圖形的任何一種元素（如顏色、分辨率等）。而gchart過程和gplot過程給出的是真正意義上的圖形，可以用很多的語句和選項來控制圖形的各方面的性質(zhì)和特征。

12Chart（或gchart）與plot（或gplot）

chart和gchart過程可以繪制出的圖形主要有條形圖（包括橫條和豎條）、圓圖、環(huán)形圖和星形圖等；plot和gplot過程通常用一個記錄中的兩個變量值表示點的坐標來繪制圖形，如散點圖和線圖等。13gchart過程的一般格式

Procgchart選項列表；圖形關(guān)鍵詞變量名稱/選項列表；

Run；此過程格式簡單，復(fù)雜的地方在于圖形關(guān)鍵字（每個圖形關(guān)鍵字對應(yīng)一種圖形類型）所引導(dǎo)的語句，它是控制圖形類型及圖形要素的地方，涉及到眾多的關(guān)鍵字和選項。14gchart過程可以使用的圖形關(guān)鍵字及其所繪制的圖形類型見下表：

15圖形關(guān)鍵字后的變量名用以指定進行圖形描述時的分組變量，可以是數(shù)值型的（此時以各組的組中值為分組的標志），也可以是字符型的。選項列表中比較重要的有：（1）type=統(tǒng)計量關(guān)鍵字，表示圖形對變量（sumvar所指定的變量）的哪一種統(tǒng)計量進行描述，比如頻數(shù)（freq）、均數(shù)（mean）、總計（sum）、頻數(shù)百分比（percentage）等；16（2）subgroup=變量名（分組變量），指定要進行分組（各組內(nèi)再分組）的變量；（3）sumvar=變量名（數(shù)值變量），指定要進行統(tǒng)計計算的變量，也就是“type=統(tǒng)計量關(guān)鍵字”選項中統(tǒng)計量的計算所依據(jù)的變量。17下表為長沙市某大學160名正常成年女子2006年體檢資料中的血清甘油三酯的測量結(jié)果（mmol/L）試以直方圖來表示上述資料的分布特征。

18例1：單變量數(shù)據(jù)第一步：建立數(shù)據(jù)集；datatest;inputx@@;cards;0.910.88 1.41 0.96 1.48 1.46 0.91 1.101.261.69 1.14 1.24 0.98 0.68 0.83 1.77……………1.491.02 1.17 0.99 0.61 1.33 0.79 0.951.051.65 1.20 1.24 0.84 1.40 0.72 1.06;Run;雙尾符，它表示為了執(zhí)行下一個input語句，保持指針在這個輸入的記錄行上。比較有加雙尾符和沒加雙尾符的區(qū)別。19第二步三，進行三初步的三統(tǒng)計分三析，確三定分組三的個數(shù)三。Pro三cu三niv三ari三ate三da三ta=三tes三t;var三x;Run三;20所有160個觀測三值中，三最大的三為1.7三7，最小三的為0.5三1，如果三以0.1為組距三，將所三有觀測三值歸入13個組。如果根三據(jù)前面三的斯特三杰斯經(jīng)三驗公式三，則應(yīng)三將所有三觀測值三歸入8組。The三ME三ANS三Pr三oce三dur三eAna三lys三is三Var三iab三le三:xN三Me三an三S三td三Dev三Min三imu三m三Ma三xim三um160三1三.13三12三0.三294三5三0三.51三00三1.7三70021第三步三，畫圖三。pro三cg三cha三rt三dat三a=t三est三;vba三rx三/typ三e=fr三eqlev三els=13ins三ide=fr三eqspa三ce=0wid三th=6;run三;Vba三r語句表三示我們?nèi)龅娜秦Q立三的條形三圖，其三后跟的三是作圖三所要依三據(jù)的分三類變量三，這里三我們以三變量x來做直三方圖。三斜杠后三的內(nèi)容三是對所三做條形三圖進行三控制的三選項，三本段程三序所涉三及的各三選項的三含義和三功能見三下表。22各選項三的含義三和功能23Pro三cg三cha三rt三dat三a=t三est三;vba三rx三/typ三e=fr三eqlev三els=13sum三var=xins三ide=fr三eqspa三ce=0wid三th=6;Run三;24Pro三cg三cha三rt三dat三a=t三est三;vba三rx三/typ三e=fr三eqlev三els三=8sum三var=xins三ide=fr三eqspa三ce=0wid三th=6;Run三;25Pro三cg三cha三rt三dat三a=t三est三;vba三rx三/typ三e=fr三eqlev三els=13sum三var=xins三ide三=me三anspa三ce=0wid三th=6;Run三;26Pro三cg三cha三rt三dat三a=t三est三;vba三rx三/typ三e=fr三eqlev三els=13sum三var=xins三ide=me三anspa三ce=三5wid三th=6;Run三;27Pro三cg三cha三rt三dat三a=t三est三;vba三rx三/typ三e=fr三eqlev三els=13sum三var=xins三ide=me三anspa三ce=0wid三th=三9;Run三;28比較Pro三cg三cha三rt與Pro三cc三har三t;比較vba三r和hba三r;Vba三r用來列三出用戶三想要畫三的垂直三條形圖三的那些三變量，三每個圖三表變量三對應(yīng)一三張垂直三條形圖三；Hba三r對列出三的每個三圖變量三生成水三平的條三形圖。比較直三方圖與三餅圖。29Pro三cg三cha三rt三dat三a=t三est三;hba三rx/typ三e=fr三eqlev三els=13三;Run三;30Pro三cg三cha三rt三dat三a=t三est三;piex/typ三e=fr三eqlev三els=13三;Run三;31Pro三cg三cha三rt三dat三a=t三est三;piex/typ三e=me三ansum三var=xlev三els=13三;Run三;32Pro三cg三cha三rt三dat三a=t三est三;piex/typ三e=me三ansum三var=xlev三els=5;Run三;33例2：

分三組數(shù)據(jù)下表是三財政系三問卷設(shè)三計資料三匯總34第一步三，建立三數(shù)據(jù)集三；Dat三as三cor三e;inp三ut三sub三jec三t$三se三x$三nu三m1三num三2;Car三ds;………；Run三;Dat三as三cor三e;set三sc三ore三;lab三elsub三jec三t='專業(yè)'s三ex=三'性別'n三um1三='人數(shù)'n三um2三='題量';Run三;35Gro三up語句Pro三cgch三art;vba三rsub三jec三t/gro三up=se三xsum三var=nu三m2typ三e=mea三n;run;Num三2表示題三量。上述語三句意思三是：先三按性別三分兩組三，接著三在每組三中按專三業(yè)進一三步分三三組，最三后再將三每組設(shè)三計的題三量的平三均數(shù)列三示出來三。36Pro三cgch三art;vba三rsub三jec三t/gro三up=se三xtyp三e=fre三q;run;Num三2表示題三量。上述語三句意思三是：先三按性別三分組，三接著在三每組中三按專業(yè)三進一步三分組，三最后再三將每組三包含的三單位數(shù)三列示出三來。37Pro三cgch三art;vba三rsex三/gro三up=su三bje三ctsum三var=nu三m2typ三e=mea三n;run;Num三2表示題三量。上述語三句意思三是：先三按專業(yè)三分三組三，接著三在每組三中按性三別分別三分兩組三，最后三再將每三組設(shè)計三的題量三的平均三數(shù)列示三出來。38By語句Pro三cgch三art;vba三rsub三jec三t/sum三var=nu三m2typ三e=mea三n;bysex三;run;上述語三句將得三到由by變量定三義的幾三個觀測三組分別三畫出的三圖表。39Sub三gro三up語句Pro三cgch三art;vba三rsex三/sub三gro三up=su三bje三cttyp三e=fr三eq三;run;上述語三句意思三是：先三按性別三分兩組三，接著三在每組三中按專三業(yè)進一三步分三三組，最三后再將三每組包三含的單三位數(shù)列三示出來三。40Pro三cgch三art;vba三rsub三jec三t/sub三gro三up=se三xtyp三e=fr三eq三;run;上述語三句意思三是：先三按專業(yè)三分三組三，接著三在每組三中按專三業(yè)進一三步分兩三組，最三后再將三每組包三含的單三位數(shù)列三示出來三。41Blo三ck語句該語名三列出用三戶想要三畫的塊三形圖的三變量。Pro三cgch三art;blo三cksex三/gro三up=su三bje三ctsum三