时间序列分析上机操作题
8章 时间序列分析练习题参考答案

35.采用几何平均法计算平均发展速度的理由是( )。
A.各年环比发展速度之和等于总速度B.各年环比发展速度之积等于总速度
C.各年环比增长速度之积等于总速度D.各年环比增长速度之和等于总速度
B
36.计算平均发展速度应用几何平均法的目的在于考察( )。
A.最初时期发展水平B.全期发展水平
C.最末时期发展水平D.期中发展水平
B
3.发展速度属于( )
A比例相对数B比较相对数C动态相对数D强度相对数
C
4.计算发展速度的分母是( )
A报告期水平B基期水平C实际水平D计划水平
B
5.某车间月初工人人数资料如下:
月份
1
2
3
4
5
6
7
月初人数(人)
280
284
280
300
302
304
320
则该车间上半年的平均人数约为( )
A 296人B 292人C 295人D 300人
A.0.33倍B.0.50倍C.0.75倍D.2倍
D
16.已知一个数列的环比增长速度分别为3%、5%、8%,则该数列的定基增长速度为()
A.3%×5%×8%B.103%×105%×108%
C.(3%×5%×8%)+1D.(103%×105%×108%)-1
D
17.企业生产的某种产品2002年比2001年增长了8%,2003年比2001年增长了12%,则2003年比2002年增长了()。
时间第一年第二年第三年第四年第五年销售额万元10001100130013501400根据上述资料计算的下列数据正确的有a第二年的环比增长速度二定基增长速度10b第三年的累计增长量二逐期增长量200万元c第四年的定基发展速度为135d第五年增长1绝对值为14万元e第五年增长1绝对值为135万元ace7
8章时间序列分析练习题参考答案

8章时间序列分析练习题参考答案第⼋章时间数列分析⼀、单项选择题1.时间序列与变量数列( )A 都是根据时间顺序排列的B 都是根据变量值⼤⼩排列的C 前者是根据时间顺序排列的,后者是根据变量值⼤⼩排列的D 前者是根据变量值⼤⼩排列的,后者是根据时间顺序排列的 C2.时间序列中,数值⼤⼩与时间长短有直接关系的是( )A 平均数时间序列B 时期序列C 时点序列D 相对数时间序列 B3.发展速度属于( )A ⽐例相对数B ⽐较相对数C 动态相对数D 强度相对数 C4.计算发展速度的分母是( )A 报告期⽔平B 基期⽔平C 实际⽔平D 计划⽔平 B5.某车间⽉初⼯⼈⼈数资料如下:则该车间上半年的平均⼈数约为( )A 296⼈B 292⼈C 295 ⼈D 300⼈ C6.某地区某年9⽉末的⼈⼝数为150万⼈,10⽉末的⼈⼝数为150.2万⼈,该地区10⽉的⼈⼝平均数为( )A 150万⼈B 150.2万⼈C 150.1万⼈D ⽆法确定 C7.由⼀个9项的时间序列可以计算的环⽐发展速度( ) A 有8个 B 有9个 C 有10个 D 有7个 A8.采⽤⼏何平均法计算平均发展速度的依据是( )A 各年环⽐发展速度之积等于总速度B 各年环⽐发展速度之和等于总速度C 各年环⽐增长速度之积等于总速度D 各年环⽐增长速度之和等于总速度 A9.某企业的科技投⼊,2010年⽐2005年增长了58.6%,则该企业2006—2010年间科技投⼊的平均发展速度为( ) A5%6.58 B 5%6.158 C6%6.58 D 6%6.158B10.根据牧区每个⽉初的牲畜存栏数计算全牧区半年的牲畜平均存栏数,采⽤的公式是( ) A 简单平均法 B ⼏何平均法 C 加权序时平均法 D ⾸末折半法 D11.在测定长期趋势的⽅法中,可以形成数学模型的是( )A 时距扩⼤法B 移动平均法C 最⼩平⽅法D 季节指数法12.动态数列中,每个指标数值相加有意义的是()。
时间序列Stata操作题

《应用时间序列分析(第四版)》王燕编著中国人民大学出版社第四章习题71974年1月至1994年12月,某地胡椒价格数据如下:(21行*12列)1102 1151 1093 1118 1168 1118 1085 1135 1138 1135 1235 1301 1283 1250 1210 1135 1085 1060 1102 1151 1127 1226 1217 1215 1250 1210 1268 1402 1486 1534 1567 1585 1717 2002 2086 2059 1250 1210 1268 1402 1486 1534 1567 1585 1717 2002 2086 2059 2425 2326 2176 2121 2000 2000 1850 1640 1700 1925 1850 1830 1850 1790 1700 1700 1750 1775 1925 2000 1975 1940 1889 1881 2000 2024 1900 1750 1649 1601 1625 1609 1649 1640 1640 1620 1590 1526 1451 1424 1424 1329 1199 1179 1285 1349 1265 1299 1373 1440 1451 1376 1325 1261 1199 1219 1250 1274 1365 1424 1420 1385 1321 1235 1215 1310 1319 1319 1279 1481 1956 2165 2125 2087 1895 1840 1874 1863 1836 1894 2105 2159 2131 2029 2270 2411 2652 3294 3360 3686 3593 3482 3615 3963 4328 4309 4336 4382 4326 4009 4000 4070 4200 4278 4435 4772 4812 4908 4857 4865 4711 4640 4877 4902 4884 4833 4903 4963 4804 4679 4810 4571 4250 3850 3775 3357 2946 2342 1994 2420 2464 2763 2993 3108 2729 2525 2457 2136 2272 2175 2100 2068 1955 1950 1969 2025 1726 1579 1768 1766 1621 1692 1634 1750 1620 1515 1508 1525 1502 1374 1212 1198 1107 1052 1069 1050 1098 1150 1126 1200 1193 1058 1043 1026 980 976 1000 1210 1264 1150 1117 1188 1100 1040 1028 1113 1154 1350 1722 1616 1525 1403 1497 1522 1550 1575 1538 1650 1800 1933 2219 2606 2563 2433 1检验序列的平稳性(Stata 语句). drop B-T . generate n=_n . rename A price. tsset ntime variable:n, 1 to 252 delta: 1 unit . tsline price=>p r i c e{price}的时序图由时序图观测得price 变化落差很大,该序列不平稳...。
时间序列分析上机操作题

20.1971年9月—1993年6月澳大利亚季度常住人口变动(单位:千人)情况如下表。
问题:(1)判断该序列的平稳性与纯随机性。
(2)选择适当模型拟合该序列的发展。
(3)绘制该序列拟合及未来5年预测序列图。
针对问题一:将以下程序输入SAS编辑窗口,然后运行后可得图1.data example3_1;inputx@@;time=_n_;cards;63.26ﻩ7.95ﻩ5.8 49.5ﻩ50.255.4ﻩ49.9 45.348.1 61.755.2ﻩ 53.149.5ﻩ59.9ﻩ30.4ﻩ30.6ﻩ33.8 42.135.8ﻩ28.4ﻩ44.1ﻩ32.9ﻩ45.5 36.639.5 49.8 48.8 29 37.33ﻩ4.2 47.637ﻩ.339ﻩ.2 47.6 43.9ﻩ4951.2ﻩ60.8 67ﻩ 48.9 65.4ﻩ65.467.6 62.555.1 49.6ﻩ57.3 47.345.544.5ﻩ 4847.9 49.1 48.859.451ﻩ.651.4 60.9 60.9 55.8 58.662.1ﻩ64ﻩ 60.3 64.6ﻩ7179.459.983.4 75.4 80.255ﻩ.9 58.5 65.269.55ﻩ9.1 21.5 62.5 170 ﻩ-47.462.2ﻩ 60ﻩ33.135ﻩ.343.4ﻩ42.758ﻩ.434ﻩ.4;procgplotdata=example3_1;plotx*time=1;symbol1c=red I=join v=star;run;图1该序列的时序图由图1可读出:除图中170和-47.4这两个异常数据外,该时序图显示澳大利亚季度常住人口变动一般在在60附近随机波动,没有明显的趋势或周期,基本可视为平稳序列。
再接着输入以下程序运行后可输出五方面的信息。
具体见表1-表5.procarima data=example3_1;identifyVar=x nlag=8;run;表1 分析变量的描述性统计从表1可读出分析变量的名称、该序列的均值;标准差及观察值的个数(样本容量)。
时间序列分析上机操作题

时间序列分析上机操作题20.1971年9月—1993年6月澳大利亚季度常住人口变动(单位:千人)情况如下表。
问题:(1)判断该序列的平稳性与纯随机性。
(2)选择适当模型拟合该序列的发展。
(3)绘制该序列拟合及未来5年预测序列图。
针对问题一:将以下程序输入SAS编辑窗口,然后运行后可得图1.data example3_1;input x@@;time=_n_;cards;63.2 67.9 55.8 49.5 50.2 55.4 49.9 45.3 48.1 61.7 55.2 53.1 49.5 59.9 30.6 30.4 33.8 42.1 35.8 28.4 32.9 44.1 45.5 36.6 39.5 49.8 48.8 29 37.3 34.2 47.6 37.3 39.2 47.6 43.9 49 51.2 60.8 67 48.9 65.4 65.4 67.6 62.5 55.1 49.6 57.3 47.3 45.5 44.5 48 47.9 49.1 48.8 59.4 51.6 51.4 60.9 60.9 55.8 58.6 62.1 64 60.3 64.6 71 79.4 59.9 83.4 75.4 80.2 55.9 58.5 65.2 69.5 59.1 21.5 62.5 170 -47.4 62.2 60 33.1 35.3 43.4 42.7 58.4 34.4;proc gplot data=example3_1;plot x*time=1;symbol1c=red I=join v=star;run;图1 该序列的时序图由图1可读出:除图中170和-47.4这两个异常数据外,该时序图显示澳大利亚季度常住人口变动一般在在60附近随机波动,没有明显的趋势或周期,基本可视为平稳序列。
再接着输入以下程序运行后可输出五方面的信息。
具体见表1-表5.proc arima data= example3_1;identify Var=x nlag=8;run;表1 分析变量的描述性统计从表1可读出分析变量的名称、该序列的均值;标准差及观察值的个数(样本容量)。
时间序列上机操作

时间序列上机操作第⼀讲:基本的eviews操作(包括试验⼀和实验⼆)实验⼀ EVIEWS中时间序列相关函数操作【实验⽬的】熟悉Eviews的操作:菜单⽅式,命令⽅式;练习并掌握与时间序列分析相关的函数操作。
【实验内容】⼀、EViews软件的常⽤菜单⽅式和命令⽅式;⼆、各种常⽤差分函数表达式;三、时间序列的⾃相关和偏⾃相关图与函数;【实验步骤】⼀、EViews软件的常⽤菜单⽅式和命令⽅式;㈠创建⼯作⽂件⒈菜单⽅式启动EViews软件之后,进⼊EViews主窗⼝在主菜单上依次点击File/New/Workfile,即选择新建对象的类型为⼯作⽂件,将弹出⼀个对话框,由⽤户选择数据的时间频率(frequency)、起始期和终⽌期。
选择时间频率为Annual(年度),再分别点击起始期栏(Start date)和终⽌期栏(End date),输⼊相应的⽇期,然后点击OK按钮,将在EViews软件的主显⽰窗⼝显⽰相应的⼯作⽂件窗⼝。
⼯作⽂件窗⼝是EViews的⼦窗⼝,⼯作⽂件⼀开始其中就包含了两个对象,⼀个是系数向量C(保存估计系数⽤),另⼀个是残差序列RESID(实际值与拟合值之差)。
⒉命令⽅式在EViews软件的命令窗⼝中直接键⼊CREATE命令,也可以建⽴⼯作⽂件。
命令格式为:CREATE 时间频率类型起始期终⽌期则菜单⽅式过程可写为:CREATE A 1985 1998㈡输⼊Y、X的数据⒈DATA命令⽅式在EViews软件的命令窗⼝键⼊DATA命令,命令格式为:DATA <序列名1> <序列名2>…<序列名n>本例中可在命令窗⼝键⼊如下命令:DATA Y X⒉⿏标图形界⾯⽅式在EViews软件主窗⼝或⼯作⽂件窗⼝点击Objects/New Object,对象类型选择Series,并给定序列名,⼀次只能创建⼀个新序列。
再从⼯作⽂件⽬录中选取并双击所创建的新序列就可以展⽰该对象,选择Edit+/-,进⼊编辑状态,输⼊数据。
第五章时间序列分析习题

第五章时间序列分析习题第五章时间序列分析习题⼀、填空题1.时间序列有两个组成要素:⼀是,⼆是。
2.在⼀个时间序列中,最早出现的数值称为,最晚出现的数值称为。
3.时间序列可以分为时间序列、时间序列和时间序列三种。
其中是最基本的序列。
4.绝对数时间序列可以分为和两种,其中,序列中不同时间的数值相加有实际意义的是序列,不同时间的数值相加没有实际意义的是序列。
5.已知某油⽥1995年的原油总产量为200万吨,2000年的原油总产量是459万吨,则“九五”计划期间该油⽥原油总产量年平均增长速度的算式为。
6.发展速度由于采⽤的基期不同,分为和两种,它们之间的关系可以表达为。
7.设i=1,2,3,…,n,a i为第i个时期经济⽔平,则a i/a0是发展速度,a i/a i-1是发展速度。
8.计算平均发展速度的常⽤⽅法有⽅程式法和.9.某产品产量1995年⽐1990年增长了105%,2000年⽐1990年增长了306.8%,则该产品2000年⽐1995增长速度的算式是。
10.如果移动时间长度适当,采⽤移动平均法能有效地消除循环变动和。
11.时间序列的波动可分解为长期趋势变动、、循环变动和不规则变动。
12.⽤最⼩⼆乘法测定长期趋势,采⽤的标准⽅程组是。
⼆、单项选择题1.时间序列与变量数列( )A都是根据时间顺序排列的B都是根据变量值⼤⼩排列的C前者是根据时间顺序排列的,后者是根据变量值⼤⼩排列的D前者是根据变量值⼤⼩排列的,后者是根据时间顺序排列的2.时间序列中,数值⼤⼩与时间长短有直接关系的是( )A平均数时间序列B时期序列C时点序列D相对数时间序列3.发展速度属于( )A⽐例相对数B⽐较相对数C动态相对数D强度相对数4.计算发展速度的分母是( )A报告期⽔平B基期⽔平C实际⽔平D计划⽔平5.某车间⽉初⼯⼈⼈数资料如下:A 296⼈B 292⼈C 295 ⼈D 300⼈6.某地区某年9⽉末的⼈⼝数为150万⼈,10⽉末的⼈⼝数为150.2万⼈,该地区10⽉的⼈⼝平均数为( )A150万⼈ B150.2万⼈ C150.1万⼈ D ⽆法确定 7.由⼀个9项的时间序列可以计算的环⽐发展速度( ) A 有8个 B 有9个 C 有10个 D 有7个 8.采⽤⼏何平均法计算平均发展速度的依据是( )A 各年环⽐发展速度之积等于总速度B 各年环⽐发展速度之和等于总速度C 各年环⽐增长速度之积等于总速度D 各年环⽐增长速度之和等于总速度9.某企业的科技投,3,2000年⽐1995年增长了58.6%,则该企业1996—2000年间科技投⼊的平均发展速度为( ) A5%6.58 B5%6.158 C6%6.58 D6%6.15810.根据牧区每个⽉初的牲畜存栏数计算全牧区半年的牲畜平均存栏数,采⽤的公式是( )A 简单平均法B ⼏何平均法C 加权序时平均法D ⾸末折半法 11.在测定长期趋势的⽅法中,可以形成数学模型的是( )A 时距扩⼤法B 移动平均法C 最⼩平⽅法D 季节指数法三、多项选择题1.对于时间序列,下列说法正确的有( )A 序列是按数值⼤⼩顺序排列的B 序列是按时间顺序排列的C 序列中的数值都有可加性D 序列是进⾏动态分析的基础E 编制时应注意数值间的可⽐性 2.时点序列的特点有( )A 数值⼤⼩与间隔长短有关B 数值⼤⼩与间隔长短⽆关C 数值相加有实际意义D 数值相加没有实际意义E 数值是连续登记得到的3.下列说法正确的有( )A 平均增长速度⼤于平均发展速度B 平均增长速度⼩于平均发展速度C 平均增长速度=平均发展速度-1D 平均发展速度=平均增长速度-1E 平均发展速度×平均增长速度=14.下列计算增长速度的公式正确的有( )A 增长速度=%100?基期⽔平增长量 B 增长速度=%100?报告期⽔平增长量C 增长速度= 发展速度—100%D 增长速度=%100?-基期⽔平基期⽔平报告期⽔平E 增长速度= %100?基期⽔平报告期⽔平5.采⽤⼏何平均法计算平均发展速度的公式有( ) A 1 231201-?=n n a a a a a a a a nx B 0a a nx n =C 1a a nx n = D R n x = E nx x ∑=A 第⼆年的环⽐增长速度⼆定基增长速度=10%B 第三年的累计增长量⼆逐期增长量=200万元C 第四年的定基发展速度为135%D 第五年增长1%绝对值为14万元E 第五年增长1%绝对值为13.5万元7.下列关系正确的有( )A 环⽐发展速度的连乘积等于相应的定基发展速度B 定基发展速度的连乘积等于相应的环⽐发展速度C 环⽐增长速度的连乘积等于相应的定基增长速度D 环⽐发展速度的连乘积等于相应的定基增长速度E 平均增长速度=平均发展速度-1 8.测定长期趋势的⽅法主要有( )A 时距扩⼤法B ⽅程法C 最⼩平⽅法D 移动平均法E ⼏何平均法 9.关于季节变动的测定,下列说法正确的是( ) A ⽬的在于掌握事物变动的季节周期性 B 常⽤的⽅法是按⽉(季)平均法C 需要计算季节⽐率D 按⽉计算的季节⽐率之和应等于400%E 季节⽐率越⼤,说明事物的变动越处于淡季 10.时间序列的可⽐性原则主要指( )A 时间长度要⼀致B 经济内容要⼀致C 计算⽅法要⼀致D 总体范围要⼀致E 计算价格和单位要⼀致四、判断题1.时间序列中的发展⽔平都是统计绝对数。
统计学习题答案 第9章 时间序列分析

第9章 时间序列分析——练习题●1. 某汽车制造厂2003年产量为30万辆。
(1)若规定2004—2006年年递增率不低于6%,其后年递增率不低于5%,2008年该厂汽车产量将达到多少?(2)若规定2013年汽车产量在2003年的基础上翻一番,而2004年的增长速度可望达到7.8%,问以后9年应以怎样的速度增长才能达到预定目标?(3)若规定2013年汽车产量在2003年的基础上翻一番,并要求每年保持7.4%的增长速度,问能提前多少时间达到预定目标?解:设i 年的环比发展水平为x i ,则由已知得:x 2003=30, (1)又知:320042005200620032004200516%x x x x x x ≥+(),2200720082006200715%x x x x ≥+(),求x 2008由上得32200820072008200320032007(16%)(15%)x x x x x x =≥++ 即为3220081.061.0530x ≥,从而2008年该厂汽车产量将达到 得 x 2008≥30× 31.06×21.05= 30×1.3131 = 39.393(万辆) 从而按假定计算,2008年该厂汽车产量将达到39.393万辆以上。
(2)规定201320032x x =,20042003x x =1+7.8%由上得=107.11%==可知,2004年以后9年应以7.11%的速度增长,才能达到2013年汽车产量在2003年的基础上翻一番的目标。
(3)设:按每年7.4%的增长速度n 年可翻一番, 则有 201320031.0742na a == 所以 1.074log 20.30103log 29.70939log1.0740.031004n ====(年)可知,按每年保持7.4%的增长速度,约9.71年汽车产量可达到在2003年基础上翻一番的预定目标。
原规定翻一番的时间从2003年到2013年为10年,故按每年保持7.4%的增长速度,能提前0.29年即3个月另14天达到翻一番的预定目标。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
时间序列分析上机操作题-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN(2)选择适当模型拟合该序列的发展。
(3)绘制该序列拟合及未来5年预测序列图。
针对问题一:将以下程序输入SAS编辑窗口,然后运行后可得图1. data example3_1;input x@@;time=_n_;cards;2949674864 71170 60;proc gplot data=example3_1;plot x*time=1;symbol1c=red I=join v=star;run;图1 该序列的时序图由图1可读出:除图中170和这两个异常数据外,该时序图显示澳大利亚季度常住人口变动一般在在60附近随机波动,没有明显的趋势或周期,基本可视为平稳序列。
再接着输入以下程序运行后可输出五方面的信息。
具体见表1-表5.proc arima data= example3_1;identify Var=x nlag=8;run;表1 分析变量的描述性统计从表1可读出分析变量的名称、该序列的均值;标准差及观察值的个数(样本容量)。
表2 样本自相关图由表2可知:样本自相图延迟3阶之后,自相关系数都落入2倍标准差范围以内,而且自相关系数向零衰减的速度非常快,故可以认为该序列平稳。
表3 样本自相关系数该图从左到右输出的信息分别为:延迟阶数、逆自相关系数值和逆自相关图。
表4 样本偏自相关图该图从左到右输出信息是:延迟阶数、偏自相关系数值和偏自相关图。
表5 纯随机性检验结果由上表可知在延迟阶数为6阶时,LB检验统计量的P值很小,所以可以断定该序列属于非白噪声序列。
针对问题二:将IDENTIFY命令中增加一个可选命令MINIC,运行以下程序可得到表6.表6 IDENTIFY命令输出的最小信息量结果通过上表可知:在自相关延迟阶数小于等于5,移动平均延迟阶数也小于等于5的所有ARMA(p,q)模型中,BIC信息量相对最小的是ARMA(1,3)模型。
进行参数估计,输入以下命令,运行可得到表7—表10estimate p=1q=3;run;表7 ESTIMATE命令输出的位置参数估计结果表8 ESTIMATE命令输出的拟合统计量的值表9 ESTIMATE命令输出的系数相关阵表10 ESTIMATE命令输出的残差自相关检验结果拟合模型的具体形式如表11所示。
表11 ESTIMATE命令输出的拟合模型形式针对问题三:对拟合好的模型进行短期预测。
输入以下命令,运行可得表12和图2.forecast lead=5id=time out=results;run;proc gplot data=results;plot x*time=1 forecast*time=2 l95*time=3 u95*time=3/overlay;symbol1c=black i=none v=star;symbol2c=red i=join v=none;symbol3c=green i=join v=none l=32;run;表12 forecast命令输出的预测结果图2 拟合效果图5416755196563005748258796602666146562828 6465365994672076620765859672956917270499 7253874542763687853480671829928522987177 8921190859924209371794974962599754298705 100072101654103008104357105851107507109300111026 112704114333115823117171118517119850121121122389 123626124761125786126743127627128453129227129988 130756131448132129132802采用SAS软件运行下列程序:data example5_1;input x@@;t=_n_;cards;54167 55196 56300 57482 58796 60266 61465 6282864653 65994 67207 66207 65859 67295 69172 7049972538 74542 76368 78534 80671 82992 85229 8717789211 90859 92420 93717 94974 96259 97542 98705100072 101654 103008 104357 105851 107507 109300 111026112704 114333 115823 117171 118517 119850 121121 122389123626 124761 125786 126743 127627 128453 129227 129988130756 131448 132129 132802;proc gplot;plot x*t=1;symbol1i=join v=none c=blavk;run;图3 该序列的时序图通过时序图可以得知,该序列有明显的线性递增趋势,故用线性回归模型来拟合。
在接着在编辑窗口输入以下命令,运行程序:proc autoreg data=example5_1;model x=t;run;表12 AUTOREG过程输出线性拟合结果通过该表可得知:(1)因变量的名称,本例中因变量为x。
(2)普通最小二乘统计量,误差平方和、均方误差、SBC信息量、回归模型的R^2、DW统计量、误差平方和的自由度、均方根误差、AIC信息量、包括自回归误差过程在内的整体模型R^2。
(3)参数估计量。
该部分从左到右输出的信息分别是:变量名、自由度、估计值、估计值的标准差、t值以及统计量的t值的近似概率P值。
对于进行5期预测,再接着输入以下命令运行:proc forecast data=example5_1 method=stepar trend=2 lead=5out=out outfull outtest=est;id t;var x;proc gplot data=out;plot x*t= _type_ / href=2008;symbol1i=none v=star c=black;symbol2i=join v=none c=red;symbol3i=join v=none c=green l=2;symbol4i=join v=none c=green l=2;run;表13 FORECAST过程OUT命令输出数据集图示该表有四个变量:时间变量,类型变量,预测时期标示变量,序列值变量。
表14 FORECAST过程OUTSET命令输出数据集图示此表可以查看预测过程中相关参数及拟合效果。
这些信息分为三部分:(1)关于序列的基本信息。
序列样本个数、非缺失数据个数、拟合模型自由度、残差标准差。
(2)关玉预测模型的参数估计信息。
线性模型的常数估计值、线性模型的斜率、残差自回归的参数估计值。
(3)拟合优度统计量信息。
图4 FORECAST过程预测效果图589561640656727697640599568577553582 600566653673742716660617583587565598 62861866705770736678639604611594634 658622709722782756702653615621602635 677635736755811798735697661667645688 713667762784837817767722681687660698 717696775796858826783740701706677711 734690785805871845801764725723690734 750707807824886859819783740747711751(2)使用X-11方法,确定该序列的趋势。
针对问题一:运行以下程序可得到该序列的时序图,见图5。
data example4_3;input x@@;time=intnx ('month','01jan1962'd, _n_-1);format time data;cards;589 561 640 656 727 697 640 599 568 577 553 582 600 566 653 673 742 716 660 617 583 587 565 598 628 618 688 705 770 736 678 639 604 611 594 634 658 622 709 722 782 756 702 653 615 621 602 635 677 635 736 755 811 798 735 697 661 667 645 688 713 667 762 784 837 817 767 722 681 687 660 698 717 696 775 796 858 826 783 740 701 706 677 711 734 690 785 805 871 845 801 764 725 723 690 734 750 707 807 824 886 859 819 783 740 747 711 751 ;proc gplot data=example4_3;plot x*time=1;symbol1c=red I=join v=star;run;图5 1962-1970年平均每头奶牛的月度产奶量的时序图通过时序图,我们可以发现1962-1970年平均每头奶牛的月度奶产量随着月度的变动有着非常明显的规律变化,此外该序列有线性递增趋势,故此时序图具有“季节”效应。
针对问题二:采用x-11过程。
在编辑窗口输入以下命令,然后运行后可得到以下几个表和图。
data example4_3;input x@@;t=intnx ('monthly','1jan1962'd, _n_-1);cards;589 561 640 656 727 697 640 599 568 577 553 582 600 566 653 673 742 716 660 617 583 587 565 598 628 618 688 705 770 736 678 639 604 611 594 634 658 622 709 722 782 756 702 653 615 621 602 635 677 635 736 755 811 798 735 697 661 667 645 688 713 667 762 784 837 817 767 722 681 687 660 698 717 696 775 796 858 826 783 740 701 706 677 711 734 690 785 805 871 845 801 764 725 723 690 734 750 707 807 824 886 859 819 783 740 747 711 751 ;proc x11 data=example4_3;monthly date=t;var x;output out=out b1=x d10=season d11=adjusted d12=trend d13=irr;data out;set out;estimate=trend*saeson/100;proc gplot data=out;plot x*t=1 estimate*t=2/overlay;plot adjusted*t=1trend *t=1irr*t=1;symbol1c=black i=join v=star;symbol2c=red i=join v=none w=2l=3;run;消除季节趋势,得到调整后的序列图,见图6。