【应用时间序列分析】第三章课后题

目录

一、1-9题

二、10-15题

三、16题

3.1 题目

3.2 代码

3.3结果分析

(1)时序图

(2)单位根检验

四、17题

4.1 题目

4.2 代码

4.3 结果分析

(1)时序图

(2)单位根检验


一、1-9题

【应用时间序列分析】第三章课后题_第1张图片

 二、10-15题

【应用时间序列分析】第三章课后题_第2张图片

三、16题

3.1 题目

【应用时间序列分析】第三章课后题_第3张图片

【应用时间序列分析】第三章课后题_第4张图片 

3.2 代码

data Mysas.data2;                                                                                                                                                                                                                                               
input year death;                                                                                                                                                                                                                                              
cards;                                                                                                                                                                                                                                                          
1915      0.5215052                                                                                                                                                                                                                                             
1916      0.4248284                                                                                                                                                                                                                                             
1917      0.4250311                                                                                                                                                                                                                                             
1918      0.4771938                                                                                                                                                                                                                                                     
(数据太多,只放前四行)                                                                                                                                                                                                                                                 
;        
run;                                                                                                                                                                                                                                                       
proc gplot data=Mysas.data2;                                                                                                                                                                                                                                    
plot death*year;                                                                                                                                                                                                                                                  
symbol c=black v=star i=join;                                                                                                                                                                                                                                   
run;                                                                                                                                                                                                                                                                   
proc arima data=Mysas.data2;                                                                                                                                                                                                                                    
identify var=death stationarity=(adf);                                                                                                                                                                                                                          
run;
proc arima data=Mysas.data2;                                                                                                                                                                                                                                    
identify var=death nlag=24 outcov=out1;                                                                                                                                                                                                                         
run;                                                                                                                                                                                                                                                            
symbol c=red i=needle v=none;                                                                                                                                                                                                                                   
proc gplot data=out1;                                                                                                                                                                                                                                           
plot corr*lag partcorr * lag;                                                                                                                                                                                                                                   
run;

3.3结果分析

(1)时序图

【应用时间序列分析】第三章课后题_第5张图片

        由时序图可以看出,序列在0.5上下波动,无周期性和趋势性,直观上看是平稳序列。

(2)单位根检验

【应用时间序列分析】第三章课后题_第6张图片

【应用时间序列分析】第三章课后题_第7张图片

 

        由单位根检验P值知,序列非平稳,这从自相关系数未立即衰减到零也可以得到验证。

四、17题

4.1 题目

【应用时间序列分析】第三章课后题_第8张图片

 4.2 代码

data Mysas.data2;                                                                                                                                                                                                                                               
input year shuiwei;                                                                                                                                                                                                                                              
cards;                                                                                                                                                                                                                                                         
1860 83.3
1861 83.5
1862 83.2
1863 82.6                                                                                                                                                                                                                                               
(数据太多,只放前四行)                                                                                                                                                                                                                                                 
;        
run;                                                                                                                                                                                                                                                       
proc gplot data=Mysas.data2;                                                                                                                                                                                                                                    
plot shuiwei*year;                                                                                                                                                                                                                                                  
symbol c=black v=star i=join;                                                                                                                                                                                                                                   
run;                                                                                                                                                                                                                                                            
proc arima data=Mysas.data2;                                                                                                                                                                                                                                    
identify var=shuiwei stationarity=(adf);                                                                                                                                                                                                                          
run;
proc arima data=Mysas.data2;                                                                                                                                                                                                                                    
identify var=shuiwei nlag=24 outcov=out1;                                                                                                                                                                                                                         
run;                                                                                                                                                                                                                                                            
symbol c=red i=needle v=none;                                                                                                                                                                                                                                   
proc gplot data=out1;                                                                                                                                                                                                                                           
plot corr*lag partcorr * lag;                                                                                                                                                                                                                                   
run;

4.3 结果分析

(1)时序图

【应用时间序列分析】第三章课后题_第9张图片

         由时序图可以看出,序列无明显周期性和趋势性,因而直观上认为该序列平稳。

(2)单位根检验

【应用时间序列分析】第三章课后题_第10张图片

【应用时间序列分析】第三章课后题_第11张图片 

        由单位根检验结果知,序列非平稳,这从自相关系数短时间内未衰减到零也可以得到验证。 


个人见解,还行各位读者批评指正!

你可能感兴趣的:(应用时间序列分析,数据分析)