時(shí)間序列分析之試驗(yàn)一2021/6/281SAS簡(jiǎn)介及數(shù)據(jù)集的建立一、SAS簡(jiǎn)介2021/6/282SAS(Statistical

Analysis

System)是一個(gè)管理數(shù)據(jù)、分析數(shù)據(jù)和打印各種報(bào)告的大型組合統(tǒng)計(jì)軟件系統(tǒng)。最早由美國(guó)北卡羅納州州立大學(xué)的兩名教授研發(fā)。1976年創(chuàng)建了SAS研究所，正式推出SAS軟件，期初限于統(tǒng)計(jì)分析，現(xiàn)在是用于決策支持的大型集成信息系統(tǒng)，但統(tǒng)計(jì)分析功能仍是其重要組成部分和核心功能。1、SAS版本2021/6/283目前最新版本SAS

9.3，要近20G。一般個(gè)人用版本SAS

8.1(2)(300M多)（低版本對(duì)中文支持不好），由數(shù)十個(gè)模塊構(gòu)成，功能包括數(shù)據(jù)訪問、數(shù)據(jù)儲(chǔ)存及管理、應(yīng)用開發(fā)、圖形處理、數(shù)據(jù)分析、報(bào)告編制、運(yùn)籌學(xué)方法、計(jì)量經(jīng)濟(jì)學(xué)與預(yù)測(cè)等。2、SAS的主要模塊2021/6/284SAS系統(tǒng)是多模塊的集成系統(tǒng)，主要有：Base

SAS是SAS系統(tǒng)的核心-負(fù)責(zé)數(shù)據(jù)管理，交互應(yīng)

用環(huán)境管理，用戶語言處理及其它軟件產(chǎn)品的調(diào)用。SAS/STAT—統(tǒng)計(jì)分析方法軟件包，是國(guó)際統(tǒng)計(jì)分析領(lǐng)域的標(biāo)準(zhǔn)軟件包。SAS/GRAPH—統(tǒng)計(jì)圖形軟件包。（一般統(tǒng)計(jì)學(xué)處理用上面三個(gè)模塊）其它還有很多專用模塊，如：SAS/OR--運(yùn)籌學(xué)方法軟件包。SAS/QC--全面質(zhì)量管理方法軟件包?！?、SAS/ETS介紹2021/6/285SAS/ETS

(Econometric

&

Time

Series)是專門進(jìn)行計(jì)量經(jīng)濟(jì)與時(shí)間序列分析的軟件。具有全球一流的數(shù)據(jù)倉(cāng)庫功能，因此在進(jìn)行

海量數(shù)據(jù)的時(shí)間序列分析時(shí)具有優(yōu)勢(shì)。時(shí)間序列分析還可采用Eviews

(EconometricsViews)，稱為計(jì)量經(jīng)濟(jì)學(xué)軟件包。4、SAS的運(yùn)用2021/6/286SAS的運(yùn)用以程序來完成，其程序的特點(diǎn)：1、程序由語句組成，語句以分號(hào)(;)結(jié)束；2、程序由“兩步”實(shí)現(xiàn)，即數(shù)據(jù)步（以data開始）和過程步（以proc開始），提交時(shí)程序都由run語句結(jié)束。3、數(shù)據(jù)步的功能是創(chuàng)建數(shù)據(jù)集，常用語句：DATA語句，INPUT語句，CARDS語句，賦值語句，SET

語句。4、過程步的主要功能是對(duì)數(shù)據(jù)集作分析、報(bào)表，語句非常多。如print，var，where，title，sum，means等。;SAS操作界面

SAS的啟動(dòng)與退出；

菜單欄

工具欄

窗口：程序編輯窗口(ProgramEditor)運(yùn)行記錄窗口(Log);

結(jié)果輸出窗口(Output)瀏覽器窗查看、管理文件日志窗口F6顯示sas程序運(yùn)行進(jìn)程及信息輸出窗口F7瀏覽sas程序運(yùn)行結(jié)果編輯窗口F5Sas程序編輯結(jié)果窗口2021/6/287查看、管理SAS運(yùn)行結(jié)果5、SAS的界面二、SAS數(shù)據(jù)集的創(chuàng)建2021/6/2881、SAS數(shù)據(jù)集類型SAS數(shù)據(jù)文件（SASdatafiles）--包括描述部分及數(shù)據(jù)部分；SAS數(shù)據(jù)視窗（SASdataviews）--只有描述部分，無數(shù)據(jù)部分（含映射關(guān)系），能使SAS過程訪問。二、SAS數(shù)據(jù)集的創(chuàng)建2021/6/2892、SAS數(shù)據(jù)集命名三部分：庫標(biāo)志.數(shù)據(jù)集名.集類型庫標(biāo)記─SAS數(shù)據(jù)庫的邏輯名字;數(shù)據(jù)集名字─SAS數(shù)據(jù)集的名字;成員類型─即擴(kuò)展名，指DATA或view。使用時(shí)一般不必給出,默認(rèn)DATA。如

WWW.SAMPLE1兩類：永久數(shù)據(jù)集和臨時(shí)數(shù)據(jù)集。一個(gè)SAS數(shù)據(jù)集是臨時(shí)或永久，取決于該數(shù)據(jù)集所附屬的SAS數(shù)據(jù)庫是臨時(shí)的或永久的。如果數(shù)據(jù)集名字只有第二部分或庫標(biāo)記為WORK時(shí)，SAS數(shù)據(jù)集是臨時(shí)的；（work庫可不寫出）而用Libname語句定義庫名，并用這個(gè)庫名作為SAS數(shù)據(jù)集名的第一部分時(shí)，SAS數(shù)據(jù)集是永久的。SAS自動(dòng)指定兩個(gè)庫標(biāo)記：SASUSER（永久）和WORK（臨時(shí)）。2021/6/28103、SAS的存儲(chǔ)方式DATA

example1；//標(biāo)志數(shù)據(jù)步開始；并定義或打開數(shù)據(jù)集名INPUT

變量名表；//為數(shù)據(jù)定義變量名；確定變量格式和讀入方式；讀入指定數(shù)據(jù)列。三種方式：（1）自由格式input name

$

age

height@@;注：$--標(biāo)示name為字符型變量；@@連續(xù)讀入記錄，否則只讀一次。2021/6/28114、數(shù)據(jù)步（data）常用語句（2）列方式Cards;//讀入數(shù)據(jù)列表王軍男40.5孫雪女18.9//每行一條記錄，一組變量值；若用

@@可以一行放多條記錄//前面數(shù)據(jù)最好與CARDS對(duì)齊。；//以；結(jié)束，一般要另起一行；input name

$1-20

sex

$22-23

age

25-26;注：變量排列位置固定。（3）格式化方式input name

$10.

sex

$2.

age

5.1;注：變量寬度位數(shù)格式固定。格式最后要有.2021/6/2812賦值語句;2021/6/2813//利用現(xiàn)有變量產(chǎn)生新變量如y=AGE+10;

X=SQRT(AGE);SET

語句//調(diào)用數(shù)據(jù)集或縱向連接數(shù)據(jù)集如

DATA

CCC;SETAAA;//打開數(shù)據(jù)集CCC,調(diào)入數(shù)據(jù)集AAA的數(shù)據(jù)如

DATA

CCC;SET

AAA+BBB;INFILE

語句//打開數(shù)據(jù)集CCC,調(diào)入數(shù)據(jù)集AAA+BBB的數(shù)據(jù)//讀入外部文件數(shù)據(jù)如

data

abc；infile

“d:\練習(xí)\abc.TXT”;proc

print//查看數(shù)據(jù)集語句如proc

print

data=example1；【實(shí)例1】創(chuàng)建臨時(shí)數(shù)據(jù)集2021/6/2814data

example1_1;//輸入格式//輸出格式input

time

monyy7.

price;format

time

monYY5.

;cards;Jan2005

101Feb2005

82Mar2005

66Apr2005

35May2005

31Jun2005

7;proc

print

data=example1_1;Run;創(chuàng)建永久數(shù)據(jù)集sasuser永久數(shù)據(jù)庫（系統(tǒng)默認(rèn)）data

sasuser.example1;用libname建立自己的永久數(shù)據(jù)庫

libname

文件夾名‘盤符:\上層文件夾\文件夾名’;//建庫，庫名即文件夾名，文件夾要先建好。data

文件夾名.example1;//將數(shù)據(jù)集存入庫中。2021/6/2815三、時(shí)間序列數(shù)據(jù)集的處理2021/6/28161、間隔函數(shù)—自動(dòng)產(chǎn)生等時(shí)間間隔的時(shí)間數(shù)據(jù)time=intnx('month','01jan2005'd,_n_-1);第一參數(shù)：等時(shí)間間隔，常用的：month,day,week,quarter,year第二參數(shù)：起始參照時(shí)間，如'01jan2005‘d第三參數(shù)：_n_k,用于調(diào)整開始觀測(cè)指針。k為整數(shù)，k取正值，指針由參照時(shí)間向未來(不包括參照時(shí)間)撥k期，k取負(fù)值，指針由參照時(shí)間向過去(包括參照時(shí)間)撥k期。(注意：在DATA步內(nèi)用有效)data

example1_2;input

price

;time=intnx('month','01jan2005'd,_n_-1);format

time

monyy.;cards;3.413.453.423.533.45;proc

print

data=

example1_2

;run;

結(jié)果見P132021/6/2817【實(shí)例2】間隔函數(shù)應(yīng)用2、序列變換2021/6/2818對(duì)觀察值序列{xt}，需分析其某個(gè)函數(shù)變換值序列時(shí)。如要分析對(duì)數(shù)序列{lnxt}，在建立數(shù)據(jù)集時(shí)，可以通過簡(jiǎn)單的賦值命令實(shí)現(xiàn)：【實(shí)例3】data

example1_3;input

price

;logprice=log(price);time=intnx('month','01jan2005'd,_n_-1);format

time

monyy.;cards;3.413.453.423.533.45;proc

print

data=example1_3;

run;

結(jié)果見P143、子集2021/6/2819只需分析時(shí)間序列的部分序列值時(shí)，可在DATA步中建立一個(gè)子集【實(shí)例4】data

example1_4;set

example1_3;

\\數(shù)據(jù)集1_4從1_3中產(chǎn)生(子集)keep

time

logprice;\\只保留兩變量where

time>='01mar2005'd;proc

print

data=

example1_4;run;4、缺失值插值2021/6/2820觀察值序列有缺失值時(shí)，可用EXPAND過程，用插值方法補(bǔ)全?！緦?shí)例5】dataexample1_5;Input

price

;time=intnx('month','01jan2005'd,_n_-1);format

time

date.;ca