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