SVAR模型制作过程

1、设置月度数据 MONTHLY>start date： 2008M01>end date 2018M08一，数据的季节调整(利用 x-12进行季节性调整)由于在建模时所选取的是宏观经济的月度数据，而月度数据容易受到季节因素的影响，从而掩盖经济运行的客观规律，因此我们采用 Census X13 (功能时最强大的) 调整方法对各个变量数据进行季节性调整。分别记做CPI' FOOD、HOUSE、M2、VMI'。时间序列按照时间次序排列的随机变量序列，任何时间序列经过合理的函数变换后都可以被认为由几个部分叠加而成。三个部分：趋势部分(T),季节部分(S)和随机噪声部分(I)。

2、常见的时间序列都是等间隔排列的。时间序列调整各部分构成的基本模型Xt= T+ Tt+ It对任何时刻有，E (It) =0, Var (It) =c 2加法模型 Xt= T *Tt* It对任何时刻有，E (It) =1, Var (It)=异 加法模型(1) 判定一个数据序列究竟适合乘法模型还是加法模型，可考察其趋势变化持性及季节变化的波动幅度。(2) 所谓季节调整就是按照上述两种模型将经济时间序列进行分解，去掉季节项的序列成为调过序列。对于时间序列而言是否存在整体趋势？如果是，趋势是显示持续存在还是显示将随时间而消逝？对于时间序列而言是否显示季节性变化？如果是，那么这种季节的波动是随时间而

3、加剧还是持续稳定存在？I这阿候脚珮可以看到旳同序列團了 1TSHS h H9_ rrta-yN A3-kuYg Ttf UK d-HS 卡 r uos mi Ls-y8 4 Ms -MJT* 31it I ISZS-我们看到：JW列3示整嗥上严超却，曲用列11摄讨司向增加。二卄趋势似乎碍样熱即为註性理势，上序沁 有一个明异托淳特彳即年度高点.在+二月*季节娈化显林1升序列而噌民的時輒耒明岸乘:云李节樟型而 不定加尿季节樽军。对于时间序列的分解模型主要有加法模型和乘法模型。加法模型适用于T、S C相互独立的情形。乘法模型适用于 T、S C相关的情形。由于时间序列分解的四大要素一般都存在相 互影响

4、，因此大多数的经济数据都采用乘法模型进行季节性分解。第一步：双击进行季节性调整的变量组 CP, proc >Seasonal Adjustment>x-12殛肚栈型加法粒Mtasntl Filltr 恤诃口 12血udJ : Li £ "iii ti vFlit *r刃敷加法模世呆巽于节洞整序列Y-SA怎线千节同于¥心卜區毎鉛斡稲坏対了 ¥.TCI轴pemEil Sridi i# SfttBue u-ue 1 |y7i nal *E-«45Qnilly djusttd wer 厂 fi j；*i£Q£iJil f

5、BCtvrf JSFJ 厂 Fin匚 4r±-ela ( T匚f" Ii n&Lr e gul ax cocnpone-fit. I. 防u理终不規则闵了 ¥皿IT" t<Htbi-+ 厂Kali d-ay/ tTadinj dla- £« 1 ors掘日空射日混今冈丁臨區 第二步：用Eviews软件进行季节调整的操作步骤：1,准备一个用于调整的时间序列（GDP）（注意：序列需同口径（当月或当季）、不变价、足够长）2，在Eviews中建立工作文件，导入序列数据3，序列图形分析（1）观察序列中的是否有季节性（2）是否有离群

6、值或问题值（3）序列的趋势变动（是加法还是乘法模型）（加法模型主要适用于呈线性增长的数据序列，或者是围绕某一个中指波动的数据序列，如pmi数据序列）（乘法模型主要适用于呈指数级数增长的序列，如GDP、工业增加值，投资数据的名义值、实际值及物价的指数序列等。）（对数加法模型主要适用于同比增速呈线性增长的数据序列，如GDP、工业、投资及cpi的同比增速数据；伪加法模型则主要是对某些非负时间序列 进行季节调整，他们具有这样的性质：在每一年中的相同月份出现接 近与0的正值，在这些月份含有接近于0的季节因子，受这些小因子的影响，季节调整结果将出现偏差。在一年的特定时期，农产品产量 就是这样的数据序列）C

7、pi, vmi为对数加法模型，（4）必要时还要分析谱图和自相关、偏相关图4，季节调整参数设定（1）季节调整选择项（模型分解方法、季节虑子、调整后的序列变量名）a. 勾选 x11 method 中的 multiplicative , seasonal filter 中的 auto x12 defaultb. Component series to save 选择 final seasonal factor （_SF）Trend Filter 选择 Auto （X12 fefat）（2）ARIMA模型参数（序列是否需要做转换、ARIMA说明）（主要是做预测用）（3）交易节假日设定（西方模式，不适合

8、中国模式）（4）离群值设定（5）模型诊断（选上）5，执行季节调整6，查看季节调整后的结果7分析季节调整的结果诊断报告主要查看M1-M11、以及Q统计量有没有通过检验如果诊断报告不好，返回第 4步8，导出数据，在EXCEL中计算环比增长率在建立SVAR模型时，需要考虑变量序列的平稳性，这就要求在建模前需要对变量进行 平稳性检验，如果变量序列是平稳的，那么可以直接进行SVAR模型的构建，但是如果变量为非平稳序列那么需要对变量序列进行平稳性处理，常用的方法是做差分和取对数，如若变量序列满足同阶单整，则可以进行协整检验，如若各个变量序列满足协整检验，具有长期的均衡关系，则可以建立SVAR模型。PROC

9、>Seas onal Adjustme nt>Ce nsus X12Sensonal adjustment（季节调整选择设定），ARMIA Option , Trading Day/Holiday （交易日、 节假日设定），Outliers （离群值设定），diagnostics （诊断）。做的比较粗糙一点：（1）打开变量列，proc>x-13>method>x-11>additive（加法）（2）Output>seas on ally adjusted一，对各变量序列的平稳性检验（ADF检验）原因:模型要求所需的变量数据为平稳序列。（1）单位根检验单

10、位根检验是检验数据的平稳性，或是说单整阶数。引用高人的回答：滞后阶数的问题。最佳滞后阶数主要根据AIC SC准则判定，当你选择好检验方式，确定好常数项、趋势项选择后，在lagged differences栏里可以从0开始尝试，最大可以尝试 到7。你一个个打开去观察，看哪个滞后阶数使得结论最下方一栏中的AIC和SC值最小，那么该滞后阶数则为最佳滞后阶数。单位根是否应该包括常数项和趋势项可以通过观察序列图确定，通过 Quick-graph-line操作观察你的数据，若数据随时间变化有明显的上升或 下降趋势，则有趋势项，若围绕0值上下波动，则没有趋势项；其二，关于是否包括常数项有两种观点，一种是其截

11、距为非零值，则取常数项，另一种是序列均值不为零则取常数项。使得t大于1%, 5%，10%条件小的值步骤：第一：利用图形确定常数项和趋势项Quick>series statistic>u nit root test其中：检验对象Level （水平序列），1st differenee （阶差分序列），2st differenee （二阶差分序列） 检验附加项Intercept （常数项，漂移项），trend and intercept （趋势项和漂移项），none （无附加项） Lag length （之后长度） lagged differencesAutomatic selectio

12、n （系统自动选择之后长度）AIC SIC 等。User specified （用户自己选择）第二，确定滞后项方法一是在User specified （用户自选模式）中选择从0开始慢慢增加，看下面的AIC与sic的大小，最后 AIC与sic最小时，就是滞后项数。方法二是在Automatic selection中选择AIC模式，可以把最大滞后项数选大一点（7或者以上），软件会自动选择 AIC最小时的项，即为滞后项。D（x（-1）为滞后1项。（3）Joha nsen检验（视单整情况而定）Johansen检验的关键是有同阶单整可以进行协整检验。非同阶单整可不需要进行Joha nsen 检验。协整检验

13、是两个或多个变量之间具有长期的稳定关系。但变量协整的必要条件是他们之间时同阶单整，也就是说在进行协整检验之前进行单位根检验。根据SIMS（ 1990）的研究结果，只有在变量序列之间存在长期的均衡关系即协整关 系时，VAR模型才能避免出现错误识别，才能通过最小二乘法得到一致估计。（4）建立VAR模型（不断重复直至模型通过三项检验：稳定性，滞后阶数正确，外 生变量与生变量明晰）第一步估计 var 模型，Objects>New object/Var选择 VAR type 为：unrestrictedEn doge nous Variables :生变量（d（vmi_d11）差分）（有生变量为

14、1,有外生变量为 0） Exogenous Variables:外生变量估计系数的标准差（圆括号中）及 t-统计量（方括号中）d（cpi_d11） d（food_d11） house_d11 d（m1_d11） d（vmi_d11）不断改变 En doge nous Variables 中（1,?）? =1,2,3比较结果最下面的 AIC与SC DE值越小越好，最后确定VAR模型的滞后阶数。（注意：1，其实在初始设置VAR模型的时候可以任意设置为（1,?）（后面检验的时候才会确认？的滞后阶数是什么）。（1,1）自己2，默认为全体变量为生变量（后面检验的时候可以确定哪些是外生变 量）。第二步检验

15、所估计的VAR模型（三个检验）1 , VAR的滞后阶数检验在 VAR工作表中 VIEWS>lag structure>lag length criteria （填写最大阶数） 软件将会用“ *”给出某个AIC或者SC准则的最小值。（滞后阶数越小 越好）。oar LagOrjorSfliedion CriteriaEnd&gmous wiawes o(cpldi 1)d(fcx3Djdi 1)housejm 1 EMgenous v,Hnabl#9, CDale-10flQ? 16 Time- 03-20Sample:jSkKJl 20ldkl06Unclkpded 口 bs

16、E制讣仃百 115LagLaglLRFP'ESCHQ0-1109536初歸血la.nno16H4B4Z19 77912Il ;i; -I;.'157 20CM3& 6790117 76071A 0*62r-0K750959,9036631170CS17J62&2216-91069T18 14780书竝.酬対5 E.90695_2fl.270£71E EI6SDMOUA48.mia'2B.462M* <Uy嗣曲Q，1S.«0S621844220B16 9»8«7227030 74374tT?777720C162

17、®18 B1211-9Cf.M5315 2112739 B3 &4 217S15?2143S091'3U7 BE'5.937421.D61B743 M91117.30164221Q51519 BECBlB-&&E.910T2? blEMO/皓9勢1U _孤 422 ajMiJ 吳：12* indcs lag crdtrbrflhii ailetioRLRtiaO moddfiG d LR 忧出 sta1i£Uc 抽曲也门温 E 补 lEvalFPE Final predction vroi科U AjfikeinTixmsiicwi c

18、merwiSC Schh3.!z|irrnrrn3ticin cnlenonHflr Hannai-Cmnn inhrra1icn mlenon2,的稳定性检验(AR根小于1 ,在单位圆才能满足脉冲分析及方差分解所 需条件)。VIEWS>lag structure>AR ROOTS TABLE/ GRAPH 3,Gra nger 检验VIEWS>lag structure>Pairwise Gran ger Causality Testsffrn | rrvt WW¥U| I FTII111 FMirrW 1 口 曜 I | 匸共111 rrfi b >

19、WI4 J J pm¥|LUIHTiv'AR Granger c：3U"irt)fl3fd( I：他證 rwity 钟 aid resii Dsrts WDSIrtB Ti- ?ri?Same 0 20OSW1 1SMO0In dudsd j L 皐代L1 歐 126Dependent'/a riahlF DiCP? C-1)iduaeiaChiMPtodDlFODO_D11)3 W737fl1OOT3SHOUSEJ011f JUO307312 02Z-D(M1 JJ1D:辺110 0014D5UI_D1I10 0UT111C 82：3All3S327604

20、QCIIKHD.<FIj 'miFv .I'Jtil. DiFOQDDVFEidudedChi-sqdfProkD(CH Dt10 033QM1HnueE_DH101745 IT JOE90974E81OOCM4doui-DHjC QM4汰1;驱/All117.214 B640001BD«pMd«fitvariAbl# HUUBL_J1*ErdijiieaChi sqPTDbDCCFI cm2 31910 12W3，建立的简约式 VAR(?) ?为滞后阶数的模型输出样式VIEW>REPRESPENTATION(5) 在构建成功VAR模型后，为了验

21、证扰动项之间是否存在同期相关关系，可用残差的同期相关矩阵来描述，可以利用这个模型进行预测即下一步的分析，为了 验证扰动项之间是否存在同期相关关系，可以利用残差的同期相关矩阵来描述。在构建的 VAR窗口中：VIEWS>Residuals>correlation matrixAR1 OOQOODO3W?3O-0 39869®DL3C(M1_-D. 355301 000 QOQDLQ9(G0FJ0.3696980151449l.OOCWOQPfinl IJanie Freeze E&tirrite SEdti impuite Res rd &RMHiiM Cul

22、l hilh >«i I4dtl i>RR DLOglLP" DLOgDPHI埠JW詁JH!切Ff4D4«ii4u atie-LriN4vflCt> Of"，D 00 a 187 (DUU24t| ix Mnn IIf.l：f: pm 3i|triBriidi|i1电*1 Ti-tpurn斜PJCFIHe-HU JljJL4il N| i j11 wiglOOJtltID 001WI20Q521-11 JIM3 014)331和曲刃|C 1«4眄L Ffe'-rw WiiiweiEhpr«4ai|si n-u

23、s0 I43OII)|T他叭I i I -IAfcncfilt： SVAR5a(l：SJJ_3_5叵| Van vafjoiPioc ObjiKt从表中可以看判实肛利睾弘实龄Ml ffldln（1i方程和实陆切户的伽如方程的残琴 功壬间存在的冋期相关耶数比越高，进一加衰明实靂利事、实际贷币供箱量M）和实际的进一步表明可以利用同期的影响来构建SVAR模型。（5）在已构建的VAR模型上构建SVAR模型第一步：实施约束识别条件为k（k-1） 12个，识别约束条件可以是短期约束条件，也可以长期约束条 件。短期约束意味着脉冲响应函数随着时间的变化将会消失，（对Do进行影响）而长期约束意味着对响应变量未来

24、的值有一个长期的影响。（更像是累计影响如?=o?不能同时施加长期与短期约束。短期约束是基于 A-B型SVAR莫型（Aet=But），长期约束基于脉冲响应的累积响应函数。（1 ）短期约束可识别条件：AB型SVAR模型至少需要2k2-k（k+1）/2个约束可识别条件一般假设结构新息ut有单位方差，因此通常对矩阵B的约束为对角阵（约束个数为k2-k）或者单位矩阵（约束个数为K2）,以致获得冲击的标准偏差A矩阵主对角元素一般设为 1 （约束个数为k）在矩阵B为单位阵情况下，对A矩阵的约束相当于对Co矩阵施加约束，即对变量间同期相关关系的约束，如有三个生变量税收（1）,政府支出（2）,产即矩阵C0中C2

25、3=0,0,比如第2个生变量出（3）,根据经济理论当期产出不会影响当期政府支出, 在约束时当B为单位阵时，直接写成 a23=约束矩阵中未知元素定义为NA（2 ）长期约束建立包括长期响应矩阵屮模块，约束处填写对第1个结构冲击的长期影响为0,则长期响应矩阵模块中第2行第1列约束为0,其他类同，无约束的填写 NA施加在当期的约束就是短期约束，（3）为了简便起见应按如下进行SVAR短期约束条件的设立1，AB型SVAR模型至少需要 2k2-k（ k+1） 12个约束（AB型的特点 是，可以明确建立系统各个生变量的当期结构关系，并且可以 直观地分析标准正交随机扰动项对系统产生冲击后的影响情况， 即et就是

26、所谓的“标准正交随机扰动项”，因为它的组成元素之 间互相正交（即互相独立），并且其方差-协方差为单位阵）2，若约束矩阵B为单位阵，此时约束个数为K2个3，若约束矩阵A为主对角元素为1，约束个数为K4，再加上经济原理上，使得在矩阵A中至少增加2k2-k（ k+1）/2 -（ k2+k）个 0 约束5，构造的约束按照 012=0或者C21=0来进行。1，EVIEWS在计算过程中同时限定了矩阵A必须为单位矩阵，对于n变量的SVAR这实质上又给出了n2个限制条件。所以，当在EVIEWS中设立长期约束条件时，实际上对矩阵C的约束条件，只要有 2n2-n（n+1）/2 -n2=n2-n （n+1）/2-n

27、2个就满足了 SVAR模型的可识别条件。2，长期约束，实质上就是要限定短期条件下的矩阵A和B与长期条件下的矩阵C之间的关系。3，在长期约束过对矩阵 C中的元素加以限制，然后利用这些限制条件以及C与矩阵A,B的关系估计出矩阵 A与B的系数。因此，在给定一个限制约 束条件的矩阵C后，EViews部算法会给出相应的 SVAR模型A、B矩阵中 的系数，而无论如何限制矩阵C，EViews给出的结果中矩阵 A总是单位矩阵。所以，当在 EViews中设立长期约束条件时，实际上对矩阵 C的约 束条件，只要有【2n2-n （n+1） /2】-n2=n 2-n （n+1）/2个就满足了 SVAR 模型的可识别条件

28、。4，在进行短期约束时，可以根据经济理论或者Gran ger来得到约束条件。5，构架的短期约束以格兰杰因果检验为主6，最好不要把B矩阵设为单位阵7 长期约束条件多用于货币政策的分析，情况较为复杂，在其他领域应用 较少。8，Svar模型根据其建模特点，主要分为3种类型：K型，c型和AB型，其中AB其中型是最通常的类型，而K型、C型都可视为 AB型的特殊形式。 如果模型中的 A矩阵为单位矩阵，则AB模型就转化为 C模型;如果AB模 型中的B矩阵为单位矩阵，则模型就转化为 K模型。9，在EVIEWS中利用约束条件生成矩阵A与BObjects/new object/matrix -vector-coe

29、f 填写矩阵名称 A在命令面板中可以如下输入：A（1,1）=1,A（1,2）=NA等，表示在 A矩阵中第1行第2列中设置为1，和在第1行 第2列或者矩阵命令 matrix （ 5,5） xdata创建一个5*5的矩阵xdata。第二步：估计SVAR从 VAR对象窗口 的菜单中选择procs>estimate structural factorizationSvar OPTIONS的对话框中，击中 MATRIX按钮和short -run pattern ，并在相应 的编辑框中填入模板矩阵的名字。e仁-c(1)*e2 -c(2) *e3 -c(3)*e4-c (4) *e5+1*u1 e2=

30、-c(5)*e1 -c(6) *e4 -c(7)*e5+1*u2 e3=-c(8)*e1 -c(9)*e4 -c(10)*e5+1*u3 e4=-c(11)*e1 -c(12)*e3+1*u4 e5=-c(13)*e1 -c(14)*e2 -c(15)*e3+1*u5如上设置的短期约束条件（多增加两个好一点）d(cpi_d11) d(food_d11) house_d11 d(m1_d11) d(vmi_d11)结果如图所示约束条件如下:里面不包含 M1不影响当期食物价格,c24M1不影响当期房屋价格 C3441,4 "nahodtMiiMjni mljiii 附轲 vtHiHEQQHnHA1QICK»D i uimo i04«M0HAi： HXKKM'HJUQCUC H44M4第三步：分析A,脉冲响应函数分析d(cpi_d11) d(food_d11) house_d11 d(m1_d11) d(vmi_d11) 主界面 view>impulse response 。在impulses产生冲击的变量处填写在response处填写观测其脉冲响应的变量图片输岀，点击输岀结果