【正文】 中用于畫點(diǎn)的作圖字符。 Plot 縱坐標(biāo)變量 Y*橫坐標(biāo)變量 X?? /選項(xiàng)列表 。 PROC PLOT 過(guò)程是用來(lái)畫易生成的低分 辨 率的圖形，輸出在 OUTPUT 窗口。其中 ， 圖形大小的兩個(gè)選項(xiàng)較為常用： ? VPCT=百分比列表 —— 規(guī)定產(chǎn)生圖形在垂直方向占一頁(yè)的百分比。例如，輸出圖中有一個(gè) Z符號(hào)，表示在這點(diǎn)上共有 26 條觀測(cè)的點(diǎn)或更多點(diǎn)。此例中每個(gè)散點(diǎn)不僅用加號(hào)（ +）表示，而且還標(biāo)出了圖中每一個(gè)散點(diǎn)在數(shù)據(jù)集的同一條觀測(cè)中這個(gè)標(biāo)記變量的值。 Run 。 VPOS 或 HPOS=軸長(zhǎng)度，規(guī)定水平或垂直軸上打印輸出的數(shù)目。程序如下： Libname Study “D:/Sasdata/Mydir”。d by 98 vaxis=100 200 400 800 1600 href=39。Shanghai StockIndex39。從上面的分析我們可以看到，水平和垂直參考線提供了一種看得見的暗示。 另 外 ， 我 們 還 可 以 通 過(guò) 選 擇 菜 單 命 令 Globals/Options/Global options 或Globals/Graph/File/Print…/SAS Graph Driver 來(lái)設(shè)定圖形輸出設(shè)備。最主要的區(qū)別是，在 PROC GPLOT 過(guò)程中，可使用 SYMBOLn 語(yǔ)句和 AXISn 語(yǔ)句。 ? C=定義顏色 —— 例如 ， c=red。 4. 用 AXISn 語(yǔ)句定義坐標(biāo)軸 在 PROC GPLOT 中，可以通過(guò) PLOT 語(yǔ)句的選項(xiàng)、 LABEL 語(yǔ)句和 FORMAT 語(yǔ)句來(lái)修飾坐標(biāo)軸。 Run 。 Symboln 選項(xiàng)列表 。要生成最高 價(jià)－ 最低 價(jià)－ 收盤價(jià)圖，原來(lái) INDEXSH 數(shù)據(jù)集中的數(shù)據(jù)排列格式不符合要求，需要用數(shù)據(jù)步修改 INDEXSH 為 INDEXSH2，將一條觀測(cè)中最高 價(jià) HIGH、上海財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息管理系 IS/SHUFE Page 8 of 28 最低 價(jià) LOW、收盤價(jià) CLOSE 三個(gè)數(shù)據(jù)，轉(zhuǎn)換成三條觀測(cè)中 VALUE 變量的數(shù)據(jù)，且三條觀測(cè)的日期 DATE 不變。 value=close。 symbol2 i=hiloc。d to 39。 Run 。 為了便于理解散布圖在回歸模型中的應(yīng)用，我們先用已知的線性方程來(lái)虛擬構(gòu)造一組樣本數(shù)據(jù)，畫出散布圖，從圖 中各散點(diǎn)形成的趨勢(shì)看，是否符合線性關(guān)系，然后決定用線性回歸模型求回歸方程 。 weight=40+60*ranuni(1) 。 plot height*weight predict*weight / overlay。最后 ， 用 PROC GPLOT 過(guò)程把實(shí)際高度對(duì)重量及預(yù)測(cè)高度對(duì)重量?jī)蓮埳⒉紙D重疊放在一張圖上，通常兩張圖應(yīng)該用不同的作圖字符來(lái)表示散點(diǎn)，實(shí)際高度對(duì)重量用無(wú)任何連線的小圓點(diǎn)表示 ，預(yù)測(cè)高度對(duì)重量用平滑曲線連接散點(diǎn)。 plot height*weight height*weight /vaxis=axis1 haxis=axis2 overlay。) order=(40 to 100 by 10) 。 程序的運(yùn)行結(jié)果顯示在 GRAPH窗口中， 如圖 。用時(shí)間序列的數(shù)據(jù)建立模型時(shí)，可以去掉序列的 確定分量，也可以對(duì)確定部分建立一個(gè)明確的模型。一般地，通過(guò)計(jì)算數(shù)據(jù)的移動(dòng)平均可實(shí)現(xiàn)數(shù)據(jù)的平滑。 Set 。 retain ewma2 。Price39。 Run 。例如 ， 線性規(guī)劃中的產(chǎn)品混合問(wèn)題，經(jīng)濟(jì)學(xué)的一個(gè)核心問(wèn)題是在有限的資源情況下，如何合理的分配資源，使效益最大化。 do x=0 to 9 by 。 plot y*x=z /contour 。 plot y*x=z / nlevels=20 autolabel xticknum=10 yticknum=10 nolegend 。 用等值線和可行域方法研究?jī)蓚€(gè)產(chǎn)品 X 和 Y 的混合問(wèn)題，還可以推廣到求解 n個(gè)產(chǎn)品的混合問(wèn)題。也就是說(shuō)， GCHART 過(guò)程具有統(tǒng)計(jì)和顯示圖表的兩大功能。 Run 。 我們以 VBAR 圖形要求語(yǔ)句為例來(lái)說(shuō)明各個(gè)圖形要求語(yǔ)句的基本使用情況。 2. PROC GCHART 過(guò)程的三要素 由 PROC GCHART 過(guò)程產(chǎn)生圖形要決定以下三個(gè)要素： （ 1） 選擇圖形的表示方法 下面的幾種圖形至少要選一種： ? VBAR 語(yǔ)句 —— 繪制垂直條形圖或垂直直方圖。統(tǒng)計(jì)量的類型有： ? TYPE=FREQ—— 統(tǒng)計(jì)圖形變量的各個(gè)給定值或間隔的頻數(shù)。把每一 個(gè)數(shù)字值作為圖形的一個(gè)分開的條形或線段。 ? SUMVAR=變量 —— 指定用于計(jì)算總和或均數(shù)的變量。 ? 分組條形 圖 —— 是指由每?jī)蓷l或以上條組成一組的條形圖，組與組之間有大間隔，每組內(nèi)條形之間無(wú)間隙或小間隙。設(shè)定從 INDEXSH數(shù)據(jù)集中選擇 1991 年 1 月 1 日到 1999 年 12 月 31 日的觀測(cè)數(shù)據(jù)作為統(tǒng)計(jì)計(jì)算對(duì)象。 set 。 Value m_day 1=39。 5=39。 9=39。 。Increase by an average of percent a month39。 PROC GCHART 過(guò)程確定對(duì)新數(shù)據(jù)集INDEXSH2 繪制垂直條形圖，圖形變量為 MONTH，由于 MONTH 是數(shù)字變量，故要選用DISCRETE 選項(xiàng)，確定它為數(shù)字型離散變量。事實(shí)上，數(shù)據(jù)集 CLASS 的數(shù)據(jù)庫(kù)結(jié)構(gòu)不滿足數(shù)據(jù)庫(kù)設(shè)計(jì)的 1NF（第一范式），也同樣不滿足 PROC GCHART 過(guò)程對(duì)制作此圖的要求。 COURSE=39。 TEST=TEST2 。 Title 39。新數(shù)據(jù)集 GCLASS，用 PROC PRINT 打印過(guò)程輸出，如圖 。 Hbar NAME /Subgroup=COURSE sumvar=TEST sum mean ascending ref=180 270 。 SUM和 MEAN 選項(xiàng)指定對(duì)相同 NAME 值的 TEST 統(tǒng)計(jì)變量求和和求平均值。 Proc gchart data=。 程序如下： Goptions reset=global gunit=pct cback=white border htitle=6 htext=3 ftext=swissb colors=(back) 。 Title 39。 ? EXPLODE=‘ Wangwu’ ‘zhaoliu’ —— 把在列表中列出的 NAME 變量值所對(duì)應(yīng)的扇形分離出去，可以從視覺上起到強(qiáng)調(diào)這一塊扇形的目的。 程序的運(yùn)行結(jié)果顯示在 GRAPH 窗口中，如圖 所示。Jan39。May39。Sep39。 Proc gchart data= 。 。總體上看，以十月份和四 月份為分界線，十一、十二、一、二、三月份的成交量之和遠(yuǎn)遠(yuǎn)小于五、六、七、八、九月份的成交量之和，值得引起注意，究其原因可能與我國(guó)的企業(yè)體制、財(cái)務(wù)制度、銀行借貸、資金結(jié)構(gòu)和資金來(lái)源等規(guī)律性的操作有內(nèi)在關(guān)系，也就是說(shuō)在上述這些因素未根本改變前，這種成交量分布的情況還將繼續(xù)。 程序的運(yùn)行結(jié)果顯示在 GRAPH 窗口中 ，如圖 所示。Shanghai StockIndex39。Nov39。Jul39。Mar39。 程序如下： Goptions reset=global gunit=pct cback=white border htitle=6 htext=3 ftext=swissb colors=(back) 。 ? CFILL=BLACK—— 設(shè)置圖中所有文本為黑色，還可以是其他顏色。 PROC GCHART 過(guò)程的PIE語(yǔ)句，首先要指定圖形變量為 NAME，統(tǒng)計(jì)變量為 TEST，統(tǒng)計(jì)類型為求平均值 MEAN。Wangwu39。 四、 圓餅圖 圓餅圖是以圓的整體面積代表被研究現(xiàn)象的總體，按各構(gòu)成部分占總體比重的大小把原面積分割成若干扇形，用以表示現(xiàn)象的部分對(duì)總體的比例關(guān)系的統(tǒng)計(jì)圖。 三、 塊形圖 塊形圖與分組條形圖和分段條形圖有點(diǎn)類似 ，最主要的區(qū)別僅在于塊形圖以三維立體方式顯示條形。 Run 。 圖 符合分組條形圖的 GCLASS 數(shù)據(jù)集 上海財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息管理系 IS/SHUFE Page 23 of 28 3. 分段條形圖 例如，如果我們想要比較每個(gè)學(xué)生的考試總分和平均分，還要比較構(gòu)成總分的各科成績(jī)的大小，最恰當(dāng)?shù)姆椒ň褪鞘褂梅侄螚l形圖來(lái)表示。 Title 39。 。 OUTPUT 。 程序如下： Goptions reset=global gunit=pct cback=white border htitle=6 htext=3 ftext=swissb colors=(back) 。其中五月份是全年上漲最大的月份，七月份下跌風(fēng)險(xiǎn)較大，八月份反彈機(jī)會(huì)較大。 。 title 39。 11=39。 7=39。 3=39。01jan202039。然后調(diào)用 PROC GCHART 圖表過(guò)程確定 MONTH為圖形變量，周漲跌 CHANGE_為統(tǒng)計(jì)變量，統(tǒng)計(jì)類型為求平均值，這樣將對(duì)所有相同月份的 CHANGE%值求平均，并用條形長(zhǎng)短表示。由于每個(gè)條圖由若干段組成，因此分段條形圖又可稱作分量條形圖。條形圖可分成水平條形圖（又稱帶形圖）和垂直條形圖（又稱柱形圖），分別用 HBAR和 VBAR 語(yǔ)句實(shí)現(xiàn)。 ? SUBGROUP=變量列表 —— 將條形或線段按照指定變量的值分成段。 ? TYPE=SUM—— 統(tǒng)計(jì)圖形變量所有值的總和。 ? STAR 語(yǔ)句 —— 繪制星形圖。垂直條形圖的縱坐標(biāo)軸用來(lái)表示頻數(shù)、累計(jì)頻數(shù)、百分?jǐn)?shù)、累計(jì)百分?jǐn)?shù)、總和或平均值。由于每個(gè)圖形要求語(yǔ)句的選項(xiàng)眾多，將它們的選項(xiàng)分成兩部分：通用選項(xiàng)列表和專用選項(xiàng)列表。 圖 用 GCONTOUR 過(guò)程繪制的利潤(rùn)等值線圖 上海財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息管理系 IS/SHUFE Page 17 of 28 BLOCK 變量列表 /通用選項(xiàng)列表 BLOCK 專用選項(xiàng)列表 。在這里我們將只介紹 GCHART 過(guò)程。 Run 。每一層用一個(gè)不同的作圖字符表示。 end。設(shè) f1(x)=7X+1/X，f2(y)=9Y+1/Y，那么： Z=2+7XX2+9YY2 圖 對(duì)時(shí)間序列的數(shù)據(jù)用移動(dòng)平均估計(jì)趨勢(shì) 上海財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息管理系 IS/SHUFE Page 15 of 28 用上面的這個(gè)函數(shù)生成一個(gè)包含三個(gè)變量 X、 Y、 Z的利潤(rùn)數(shù)據(jù)集 PROFIT，然后用 PLOT語(yǔ)句的 CONTOUR 選項(xiàng)就可以繪制出等值線圖。調(diào)用 PROC GPLOT 過(guò)程繪制了收盤價(jià)和十日簡(jiǎn)單移動(dòng)平均的圖形，其中對(duì)縱軸（ Price）定義上海財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息管理系 IS/SHUFE Page 14 of 28 的 AXIS1 語(yǔ)句中加上了 LOGBASE=10（或 PI或 e）的選項(xiàng)，就表示該軸是以 10 為底的對(duì)數(shù)軸。) 。 plot close*date ma10*date /vaxis=axis1 haxis=axis2 overlay 。 retain ewma1 。然后繪制 原始序列和移動(dòng)平均的圖形。 直接使用一些股票市場(chǎng)的時(shí)間序列數(shù)據(jù)繪制折線圖，全部趨勢(shì)可能難以直觀識(shí)別。廣義上說(shuō)，模型可劃分為回歸模型與時(shí)間序列模型。 程序說(shuō)明：對(duì)樣本數(shù)據(jù)集 不先調(diào)用統(tǒng)計(jì)計(jì)算過(guò)程求回歸預(yù)測(cè)值，而是直接繪制實(shí)際身高對(duì)體重的兩張重疊散點(diǎn)圖。Height39。我們可以確認(rèn)回歸模型的估計(jì)式為： HEIGHT=*WEIGHT+ SAS 系統(tǒng)的許多圖形過(guò)程都帶有統(tǒng)計(jì)計(jì)算，許多統(tǒng)計(jì)分析過(guò)程都帶有圖形輸出。 程序說(shuō)明：數(shù)據(jù)步 DATA 自動(dòng)生成 200 條觀測(cè)， SAS 系統(tǒng)提供了兩個(gè)均勻分 布的隨機(jī)函數(shù) RANUNI（）和 UNIFORM（）產(chǎn)生 0～ 1之間的隨機(jī)數(shù)，用隨機(jī)函數(shù)構(gòu)成的 40+60*ranuni(1)表達(dá)式，將使 WEIGHT 變量的取值范圍為 40～ 100公斤之間的隨機(jī)值。 Proc reg data=。假設(shè)我們需要繪制回歸直線及置信帶等 ， 可以先畫出高度對(duì)體重的散布圖進(jìn)行觀察和分析，然后利用回歸過(guò)程建立高度（因變 量）和體重（自變量）的線性方程式后，進(jìn)一步畫出高度預(yù)測(cè)值對(duì)體重的散布圖 ， 進(jìn)行比較和分析。它是以散點(diǎn)的分布反映變量之間相關(guān)情況的統(tǒng)計(jì)圖。Price39。 axis