《stata操作介绍之时间序列-》由会员分享，可在线阅读，更多相关《stata操作介绍之时间序列-（38页珍藏版）》请在金锄头文库上搜索。
1、stata操作介绍之时间序列 平滑分析自相关分析ARIMA模型 平滑分析 大多数时间序列数据都具有上下起伏的波动 对于时间序列数据的识别十分困难 平滑分析可以把数据分为两个部分 一部分逐渐发生变化 便于分析处理 另一部分则含有突变的成份 时间平滑 用某时刻及其前后若干时刻的值进行加权平均 所得值作为该时刻的替代值以滤去小扰动的方法 平滑分析 该部分以1983年1月到7月Milford城镇的自来水消费量为例 文件 ch51 dta 导入数据 use C Users Administrator Desktop 时间序列 数据 ch51 dta clear 平滑分析 描述性统计 describe 平
2、滑分析 生成日期变量 一 gendate mdy month day year listin1 5 平滑分析 设置时间 二 tssetdate format d listin1 5 平滑分析 趋势图 graphtwolinewaterdate ylabel 300 100 900 平滑分析 生成移动平均值 1 genwater3 water n 1 water n water n 1 3 平滑分析 生成移动平均值 2 tssmoothmawater5 water window 212 注 tssmooth 表示移动平均值平滑 加权或不加权 window 212 表示使用该值的前两个值 该值与该
3、值的后两个值进行平均计算 平滑分析 趋势图 graphtwolinewater5date clwidth thick linewaterdate clwidth thin clpattern solid 平滑分析 波动幅度 gench water water5 listin1 5 平滑分析 波动幅度 graphtwolinechdate 平滑分析 滞后变量 导入数据 use C Users Administrator Desktop 时间序列 数据 ch52 dta clear描述性统计 describe时间设置 tssetyear yearly 平滑分析 滞后变量 时间平滑 tssmooth
4、mafylln fylltemp window 212 平滑分析 滞后变量 趋势图 graphtwospikefylltempyear base 1 67 yline 1 67 linefyllnyear clpattern solid 平滑分析 滞后变量 生成n阶滞后变量的两种方法 genwNAO n wNAO n 1 genwNAO n Ln wNAO注 第二种方法中的 Ln 表示Lag n 平滑分析 滞后变量 生成一阶滞后变量 genwNAO 1 wNAO n 1 生成二阶滞后变量 genwNAO 2 L2 wNAO 平滑分析 滞后变量 平滑分析 滞后变量 三种滞后回归方法 regfyl
5、ltempwNAOwNAO 1wNAO 2ifyear 1970 year 1970 year 1970 year 1990 平滑分析 滞后变量 回归结果 相关分析 自相关系数是对变量自身与其滞后变量之间相关关系的估计 偏相关系数是在消除其他变量影响的条件下 所计算的某两变量之间的相关系数 交叉相关是分析两个时间序列之间的关系 注 本部分继续使用ch52 dta数据 自相关分析 自相关表 corrgramfylltemp lags 9 说明 大多数P值都小于0 05 故认为fylltemp具有显著的自相关性 相关关系或偏相关关系越强 相应的线条越长 自相关分析 自相关图 acfylltemp
6、lags 9 注 阴影部分是95 的置信区间 自相关分析 偏相关图 pacfylltemp lags 9 yline 0 ciopts bstyle outline 注 ciopts bstyle outline 表示将偏相关图中阴影部分改为矩阵区域 自相关分析 交叉相关图 xcorrwNAOfylltempifyear 1970 year 1990 lags 7 xlabel 9 1 9 grid 说明 由图可知 lag为0时 交叉相关性最强 线条最长 且为负 自相关分析 交叉相关表 xcorrwNAOfylltempifyear 1970 year 1990 lags 7 table AR
7、IMA模型 时间序列中的自相关集成移动平均模型 autoregressiveintegratedmovingaverage简称ARIMA 是指将非平稳时间序列转化为平稳时间序列 然后将因变量仅对它的滞后值以及随机误差项的现值和滞后值进行回归所建立的模型 ARIMA模型根据原序列是否平稳以及回归中所含部分的不同 包括移动平均过程 MA 自回归过程 AR 自回归移动平均过程 ARMA 以及ARIMA过程 ARIMA模型 ARIMA模型操作步骤 1 对变量进行检验 若变量具有稳定性 进行第二步 否则 不能使用ARIMA模型 2 做出变量的自相关图 3 根据变量的自相关图 选择合适的模型 4 根据选定
8、的模型进行分析 并检验系数是否显著 若有的系数不显著 所选择的模型可能存在问题 若所有系数都显著 进行第五步 5 检验残差是否具有自相关性 若残差具有自相关性 则所选择的模型存在问题 若残差不具有自相关性 则所选择的模型是合适的 ARIMA模型 三种单位根检验方法 pperron检验 pperronfylltemp lag 3 说明 P值为0 0003 小于0 05 故拒绝原假设 说明该变量不满足稳定性检验 ARIMA模型 dfuller检验 dfullerfylltemp lag 3 说明 P值为0 089 大于0 05 故不能拒绝原假设 说明该变量满足稳定性检验 ARIMA模型 dfgls
9、检验 dfglsfylltemp maxlag 3 notrend说明 P值为0 0 小于0 05 故拒绝原假设 说明该变量不满足稳定性检验 虽然pperron检验和dfgls检验拒绝了变量fylltemp具有稳定性的假设 但是dfuller检验不能拒绝原假设 还是可以认为该变量具有稳定性 ARIMA模型 自相关表 corrgramfylltemp lags 4 说明 该变量的自相关关系随着滞后期的增加而减少 偏自相关关系在一期自后滞后消失 故适合模型AR 1 来分析该变量 ARIMA模型 AR 1 arimafylltemp arima 1 0 0 nolog说明 由上图可知 常数项与一期滞后变量系数都是统计显著的 卡方检验也显著 故上图可得到模型为 在得出结论之前 还需对残差进行检验 检验残差是否存在自相关性 在进行检验之前需要生成残差 ARIMA模型 生成残差 predictfylres res corrgramfylres lags 9 说明 由图可知 统计量对应的P值为0 6574 不能拒绝原假设 即认为残差不存在自相关 因此 认为变量fylltemp使用AR 1 进行分析是合适的
