sas學(xué)習(xí)系列17.繪圖ⅲ-散點(diǎn)圖,時(shí)間序列圖,擬合曲線圖

17.繪圖出一散點(diǎn)圖，時(shí)間序列圖，擬合回歸圖三、散點(diǎn)圖以兩個(gè)定量變量的數(shù)值在二維坐標(biāo)系中描點(diǎn)，生成的圖形就是散點(diǎn)圖。散點(diǎn)圖能直觀地顯示出因變量隨自變量變化而變化的大致趨勢，據(jù)此可以選擇合適的函數(shù)對數(shù)據(jù)點(diǎn)進(jìn)行擬合?；菊Z法：PROCSGPLOT;SCATTER洎變量Y=因變量/options;可選項(xiàng)：DATALABEL=variable-name指定該變量的值作為各個(gè)數(shù)據(jù)點(diǎn)的標(biāo)簽，若不指定變量名，則使用Y變量的值；GROUP=variablename指定分組變量，按其值對數(shù)據(jù)分組分別處理；NOMISSINGGROUP指定不包含分組變量缺省值的觀測；TRANSPARENCY-^指定透明度。例1例1關(guān)于鳥類的數(shù)據(jù)(),變量包括名稱、種類(S-鳴禽，R-猛禽)、身長(cm)、翼展(cm),注意每行有多組數(shù)據(jù):讀入數(shù)據(jù)，繪制散點(diǎn)圖，要求按類型分組。代碼：datawings;infile'c:\MyRawData\';inputName$12.Type$LengthWingspan@@;run;*PlotWingspanbyLength;procformat;value$birdtype'S'='Songbirds''R'='Raptors';run;procsgplotdata=wings;scatterX=WingspanY=Length/ GROUPType;formatType$birdtype.;title'ComparisonofWingspanvs.Length' ;run;運(yùn)行結(jié)果：四、時(shí)間序列圖時(shí)間序列圖，是以時(shí)間為X軸,變量為Y軸的一種圖(數(shù)據(jù)點(diǎn)之間依次用線段連接)，觀察變量是否隨時(shí)間變化而呈現(xiàn)某種趨勢?；菊Z法：PROCSGPLOT;SERIESX=^間變量Y=因變量/options;注：X變量必須事先排好序?？蛇x項(xiàng)：CURVELABEL='textstring'給曲線添加標(biāo)簽，若不指定標(biāo)簽文本，則使用 Y變量；DATALABEL=variable-name指定該變量的值作為各個(gè)數(shù)據(jù)點(diǎn)的標(biāo)簽，若不指定變量名，則使用Y變量的值；MARKERS給每個(gè)數(shù)據(jù)點(diǎn)加上標(biāo)簽；GROUP=variablename指定分組變量，按其值對數(shù)據(jù)分組分別處理；NOMISSINGGROUP指定不包含分組變量缺省值的觀測；TRANSPARENCY—^指定透明度。例2某作家記錄了某天的用電量數(shù)據(jù)( )，變量包括時(shí)間(24個(gè)時(shí)刻)、用電量(千瓦時(shí))：讀取數(shù)據(jù)，繪制用電量隨時(shí)間變化的曲線圖，畫一條用電量=的基準(zhǔn)線。代碼：dataelectricity;infile'c:\MyRawData\';inputTimekWh@@;run;*Plottemperaturesbytime;procsgplotdata=electricity;seriesX=TimeY=kWh/MARKERSREFLINE/LABEL=('kh')TRANSPARENCY=;xaxisvalues=(0TO24BY1);yaxislabel='Electricity(kWh)';title'HourlyUseofElectricity';run;運(yùn)行結(jié)果：例32013-2014年股票代碼為sh600000的股票數(shù)據(jù)如下：（部分）（）讀入數(shù)據(jù)，日期只選擇2014年，只保留變量date,open（開盤價(jià)），high（最高價(jià)），low（最低價(jià)），close（收盤價(jià)）。按日期從小到大排序，并繪制該股票2014年的開盤價(jià)和收盤價(jià)的變化趨勢圖。代碼：procimportdatafile""='C:\MyRawData\'""DBMSCSVOUT=sh600000REPLACE;getnames=YES;datarow=2;run;datash600000_2014;setsh600000;keepdateopenhighlowclose;where'1Jan2014'D<=date<='31Dec2014'D;procsortdata=sh600000_2014;bydate;run;procprintdata=sh600000_2014;title 'sh600000股票2014年后數(shù)據(jù)’；run;procsgplotdata=sh600000_2014;series X=date Y= open/ name='kaipan'legendlabel='開盤價(jià)’；series X=date Y= close/ name='shoupan'legendlabel ='收盤價(jià);;title '股票（sh600000）2014年開盤價(jià)、收盤價(jià)的變化趨勢’;keylegend'kaipan''shoupan'/LOCATION=INSIDEPOSITION=TOP;xaxislabel='日期'；yaxislabel='股票價(jià)格（元）';run;運(yùn)行結(jié)果：（部分）程序說明：csv文件是逗號分隔的文本文檔，不是Excel表格；"/name="語句是給曲線起名，后面“keylegend”語句使用它。五、擬合（回歸）曲線圖散點(diǎn)圖只能大致顯示兩個(gè)變量之間的關(guān)系，進(jìn)一步還可以繪制一條擬合曲線（應(yīng)用最小二乘法原理，讓誤差的平方和最?。?。這種擬合曲線大概只能為后續(xù)真正做數(shù)據(jù)擬合提供參考信息， 要擬合數(shù)據(jù)是用REGGLMNLIN過程步。用SGPLOTt程步可以繪制擬合曲線圖?；菊Z法：PROCSGPLOTdata嗷據(jù)集；擬合方式名X=自變量Y=因變量/options;可選擬合方式有3種：REG-一線性或曲線回歸，；LOES——局部加權(quán)回歸（loess）;PBSPLINE——懲罰B樣條(B樣條擬合，加入懲罰函數(shù)避免邊界附近“過擬合”)注：RE用以帶可選參數(shù) ”/degree=n"，指定以n次多項(xiàng)式擬合。例如，regx=heighty=weight/degree=3;可選項(xiàng)：ALPHA=p設(shè)定置信限：［0,1］，默認(rèn)為(95%!信)⑵CLI給個(gè)體預(yù)測值加上預(yù)測限(只限RE函PBSPLINE;CLM給均值預(yù)測值加上置信限；CURVELABEL='textstring'為擬合曲線添加標(biāo)注(默認(rèn)是用Y變量)；GROUP=variablename指定分組變量，按其值對數(shù)據(jù)分組分別處理；NOLEGCLINOLEGCLMNOLEGFIT不顯示CLI帶、CLM^＞擬合曲線的圖例說明；⑺NOMARKERS不對數(shù)據(jù)點(diǎn)做標(biāo)記；(8)CLMTRANSPARENCY=n指定置信限的透明度。）,變量包括年例4奧運(yùn)會(huì)男子1500米跑的數(shù)據(jù)（）,變量包括年份、時(shí)間（秒），注意每行有多組數(shù)據(jù):讀入數(shù)據(jù)，用REGf口LOESS以合方式繪制時(shí)間值的擬合曲線。代碼:dataOlympic1500;infile'C:\MyRawData\';inputYearMen@@;run;procsgplotdata=Olympic1500;LOESSX=YearY=Men/NOMA