时间序列第五章作业.docVIP

  1. 1、本文档共13页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
时间序列第五章作业

第五章SAS作业 问题:2203 2360 2254 2165 2024 2078 2214 2292 2207 2119 2119 2137 2132 1955 1785 1747 1818 1909 1958 1892 1919 1853 1868 1991 2111 2119 1991 1859 1856 1924 1892 1916 1968 1928 1898 1850 1841 1824 1823 1843 1880 1968 2029 1996 1933 1805 1713 1726 1752 1795 1717 1648 1512 1338 1383 1344 1384 1484 1597 1686 1707 1640 1611 1632 1775 1850 1809 1653 1648 1665 1627 1791 1、选择恰当模型,拟合该序列的发展; 2、利用拟合模型预测1938-1945年英国绵羊的数量; 3、按照书本相应例题的格式完成问题,并附上SAS程序。 解: (1) 时序图显示,序列具有长期趋势,对序列进行1阶差分▽Xt=Xt-Xt-1,观察差分后序列▽Xt的时序图。 时序图显示长期趋势信息基本被差分运算提取充分,考察差分后序列的自相关图和偏自相关图。 自相关图显示延迟3阶后自相关系数基本在2倍标准差范围内,因此认为该序列为平稳序列。自相关图表现出拖尾现象,偏自相关图表现出3阶结尾现象,且自相关图中2阶自相关系数在2倍标准差范围内,所以考虑构造疏系数模型AR(1,3)。 残差自相关检验结果显示延迟6期后P值都大于0.05,因此认为残差为白噪声序列,即拟合模型显著有效。 参数估计结果显示两参数P值都小于0.05,都显著有效。则拟合的AR(1,3)模型为 ▽Xt=0.32196▽Xt-1 – 0.37616▽Xt-3 + εt (2) 利用拟合模型对1938-1945年英国绵羊的数量进行预测结果如上图所示, 预测图为 (3)SAS程序为 data a; input x@@; dif1=dif(x); t=1867+_n_-1; format time year4.; cards; 2203 2360 2254 2165 2024 2078 2214 2292 2207 2119 2119 2137 2132 1955 1785 1747 1818 1909 1958 1892 1919 1853 1868 1991 2111 2119 1991 1859 1856 1924 1892 1916 1968 1928 1898 1850 1841 1824 1823 1843 1880 1968 2029 1996 1933 1805 1713 1726 1752 1795 1717 1648 1512 1338 1383 1344 1384 1484 1597 1686 1707 1640 1611 1632 1775 1850 1809 1653 1648 1665 1627 1791 ; run; proc gplot data=a; plot x*t dif1*t; symbol c=black i=join v=dot; proc arima; identify var=x(1) ; estimate p=(1 3) noint; forecast lead=7 id=t out=out; proc gplot data=out; plot x*t=1 forecast*t=2 l95*t=3 u95*t=3/overlay; symbol1 c=black i=none v=star; symbol2 c=red

您可能关注的文档

文档评论(0)

zhanghc + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档