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

最新推荐文章于 2024-03-19 20:03:54 发布

管理学phd

最新推荐文章于 2024-03-19 20:03:54 发布

阅读量3.3k

点赞数 14

分类专栏：应用时间序列分析文章标签：数据分析

本文链接：https://blog.csdn.net/m0_72318954/article/details/130851151

版权

应用时间序列分析专栏收录该内容

4 篇文章 6 订阅

订阅专栏

一、1-9题

二、10-15题

三、16题

3.1 题目

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）时序图

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

（2）单位根检验

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

四、17题

4.1 题目

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;