使用R语言进行协整关系检验协整检验是为了检验非平稳序列的因果关系,协整检验是解决伪回归为问题的重要方法。首先回归伪回归例子:伪回归Spurious regression伪回归方程的拟合优度、显著性水平等指标都很好,但是其残差序列是一个非平稳序列,拟合一个伪回归:#调用相关R包library(lmtest)library(tseries)#模拟序列set.seed(123456)e1=rnorm(500)e2=rnorm(500)trd=1:500y1=0.8*trd+cumsum(e1)y2=0.6*trd+cumsum(e2)sr.reg=lm(y1~y2)#提取回归残差error=residuals(sr.reg)#作残差散点图plot(error, main="Plot of error")#对残差进行单位根检验adf.test(error)## Dickey-Fuller = -2.548, Lag order = 7, p-value = 0.3463## alternative hypothesis: stationary#伪回归结果,相关参数都显著summary(sr.reg)## Residuals:## Min 1Q Median 3Q Max## -30.654 -11.526 0.359 11.142 31.006## Coefficients:## Estimate Std. Error t value Pr(>|t|)## (Intercept) -29.32697 1.36716 -21.4 <2e-16 ***## y2 1.44079 0.00752 191.6 <2e-16 ***## Residual standard error: 13.7 on 498 degrees of freedom## Multiple R-squared: 0.987, Adjusted R-squared: 0.987## F-statistic: 3.67e+04 on 1 and 498 DF, p-value: <2e-16dwtest(sr.reg)## DW = 0.0172, p-value <2.2e-16恩格尔-格兰杰检验Engle-Granger第一步:建立两变量(y1,y2)的回归方程,第二部:对该回归方程的残差(resid)进行单位根检验其中,原假设两变量不存在协整关系,备择假设是两变量存在协整关系。利用最小二乘法对回归方程进行估计,从回归方程中提取残差进行检验。set.seed(123456)e1=rnorm(100)e2=rnorm(100)y1=cumsum(e1)y2=0.6*y1+e2# (伪)回归模型lr.reg=lm(y2~y1)error=residuals(lr.reg)adf.test(error)## Dickey-Fuller = -3.988, Lag order = 4, p-value = 0.01262## alternative hypothesis: stationaryerror.lagged=error[-c(99,100)]#建立误差修正模型ECM.REGdy1=diff(y1)dy2=diff(y2)diff.dat=data.frame(embed(cbind(dy1, dy2),2))#emed表示嵌入时间序列dy1,dy2到diff.datcolnames(diff.dat)=c("dy1","dy2","dy1.1","dy2.1")ecm.reg=lm(dy2~error.lagged+dy1.1+dy2.1, data=diff.dat)summary(ecm.reg)## Residuals:## Min 1Q Median 3Q Max## -2.959 -0.544 0.137 0.711 2.307## Coefficients:## Estimate Std. Error t value Pr(>|t|)## (Intercept) 0.0034 0.1036 0.03 0.97## error.lagged -0.9688 0.1585 -6.11 2.2e-08 ***## dy1.1 0.8086 0.1120 7.22 1.4e-10 ***## dy2.1 -1.0589 0.1084 -9.77 5.6e-16 ***## Residual standard error: 1.03 on 94 degrees of freedom## Multiple R-squared: 0.546, Adjusted R-squared: 0.532## F-statistic: 37.7 on 3 and 94 DF, p-value: 4.24e-16par(mfrow=c(2,2))plot(ecm.reg)Johansen-Juselius(JJ)协整检验法,该方法是一种用向量自回归(VAR)模型进行检验的方法,适用于对多重一阶单整I(1)序列进行协整检验。JJ检验有两种:特征值轨迹检验和最大特征值检验。我们可以调用urca包中的ca.jo命令完成这两种检验。其语法:ca.jo(x, type = c("eigen", "trace"), ecdet = c("none", "const", "trend"), K = 2,spec=c("longrun", "transitory"), season = NULL, dumvar = NULL)其中:x为矩阵形式数据框;type用来设置检验方法;ecdet用于设置模型形式:none表示不带截距项,const表示带常数截距项,trend表示带趋势项。K表示自回归序列的滞后阶数;spec表示向量误差修正模型反映的序列间的长期或短期关系;season表示季节效应;dumvar表示哑变量设置。set.seed(12345)e1=rnorm(250,0,0.5)e2=rnorm(250,0,0.5)e3=rnorm(250,0,0.5)#模拟没有移动平均的向量自回归序列;u1.ar1=arima.sim(model=list(ar=0.75), innov=e1, n=250)u2.ar1=arima.sim(model=list(ar=0.3), innov=e2, n=250)y3=cumsum(e3)y1=0.8*y3+u1.ar1y2=-0.3*y3+u2.ar1#合并y1,y2,y3构成进行JJ检验的数据库;y.mat=data.frame(y1, y2, y3)#调用urca包中cajo命令对向量自回归序列进行JJ协整检验vecm=ca.jo(y.mat)jo.results=summary(vecm)#cajorls命令可以得到限制协整阶数的向量误差修正模型的最小二乘法回归结果vecm.r2=cajorls(vecm, r=2)vecm.r2## Call:lm(formula = substitute(form1), data = data.mat)## Coefficients:## y1.d y2.d y3.d## ect1 -0.33129 0.06461 0.01268## ect2 0.09447 -0.70938 -0.00916## constant 0.16837 -0.02702 0.02526## y1.dl1-0.22768 0.02701 0.06816## y2.dl1 0.14445 -0.71561 0.04049## y3.dl1 0.12347 -0.29083 -0.07525## $beta## ect1 ect2## y1.l2 1.000e+00 0.0000## y2.l2 -3.402e-18 1.0000## y3.l2 -7.329e-01 0.2952数据框是由多个不同模式的不同变量组成的,而变量分为连续型变量、名义型变量、有序型变量,其中名义型变量和有序型变量在R中被称为因子。数据框的每一列并不都是因子。只有名义型和有序型的列才是因子。你可以使用函数str(数据框)来显示此数据框的信息,会显示哪个变量是因子,哪个是不是。
使用R语言进行协整关系检验
给您推荐相同类型的内容:
"=="在c语言中表示什么意思
1、c语言中=!是两个运算符摆放在一起,“=”是赋值运算符,“!”逻辑非运算符;2、逻辑操作符一般将其操作数视为条件表达式,返回结果为Bool类型;3、首先对其操作数求值,若结果为0,则条件为假(false),否则为真(true);4js基础篇:if else和else if 的使用
1.if和else if的使用条件: if可以单独使用,也可以和else块一起使用。 ②else 必须和if配合使用(或者和if else if 配合使用) ③else if必须和if配合使用,最后可跟接else判断 elspython基本语法
python语法如下:1、Python标识符在Python里,标识符有字母、数字、下划线组成。在Python中,所有标识符可以包括英文、数字以及下划线(_),但不能以数字开头。Python中的标识符是区分大小写的。以下划线开头的标识jquery.superslide怎么做手风琴图
以下是全部代码,请自行修改图片及js路径 SuperSlide - 实用焦点图 *{margin:0padding:0list-style:none}body{ background:#ffffont:normal 12px22px 宋体css怎么做斜线
代码已经测试,兼容主流浏览器:<style>.line{width:100pxheight:25pxborder-bottom:1px solid #333filter: progid:DXImageTransCSS 文字渐变的几种方式
css文件渐变虽然兼容性比较差,但是用在移动端和chrome中还是没有问题的。 实现文件渐变的方法有两种 效果如下 -webkit-background-clip W3C支持的属性说明 但是并没有text 属性,所以这个只能java中的参数
try,catch是一个语句对,catch用来捕获try块中抛出的异常exception,而Exception是Java中所有异常的基类,它实现了Thorwable接口,两者都能够通过throw来抛出。至于e,它代表的是本次捕获到的异常对象js基础篇:if else和else if 的使用
1.if和else if的使用条件: if可以单独使用,也可以和else块一起使用。 ②else 必须和if配合使用(或者和if else if 配合使用) ③else if必须和if配合使用,最后可跟接else判断 els请分析一下asp.netjsppythonphpruby哪个更有前途?
Python最有前途,我很喜欢用。Python并不是一种新兴的语言,1991年就已经出现。你列举的这五种言语的共性是都用于Web开发,Web现在最牛的是Google,Python就是由Google支持的语言,Python的创始人现在就在Gocss如何让div横向排列
可以使用弹性布局,display: flex 就可以让div横排,默认布局轴方向就是横向,justify-content: space-between可以是div平均分布在水平轴上。如果是导航菜单可以使用ul无需列表来制作同时使用CSS中的CSS 文字渐变的几种方式
css文件渐变虽然兼容性比较差,但是用在移动端和chrome中还是没有问题的。 实现文件渐变的方法有两种 效果如下 -webkit-background-clip W3C支持的属性说明 但是并没有text 属性,所以这个只能JAVA中反射是什么
JAVA中反射是动态获取信息以及动态调用对象方法的一种反射机制。Java反射就是在运行状态中,对于任意一个类,都能够知道这个类的所有属性和方法;对于任意一个对象,都能够调用它的任意方法和属性;并且能改变它的属性。而这也是Java被视为动态js中如何动态向枚举中添加元素
1、动态创建元素时,将事件触发方式加入到元素创建中例如:document.getElementById('id').innerHTML='<input type="button"如何用HTML做树状的目录
可以通过jQuery的treeview插件实现,基本实现代码如下:<link rel="stylesheet" href="..jquery.treeview.css" type=&quJS中typeof的用法
js是一门弱语言,它在声明变量时无需确定变量的类型,js在运行时会自动判断。那么如何判断一个变量的类型呢,js提供了typeof运算符,用来检测一个变量的类型。 typeof是一个运算符,有2种使用方式:typego语言能做什么?
很多朋友可能知道Go语言的优势在哪,却不知道Go语言适合用于哪些地方。 1、 Go语言作为服务器编程语言,很适合处理日志、数据打包、虚拟机处理、文件系统、分布式系统、数据库代理等;网络编程方面。Go语言广泛应用于Web应用、API应js中如何动态向枚举中添加元素
1、动态创建元素时,将事件触发方式加入到元素创建中例如:document.getElementById('id').innerHTML='<input type="button"HTML列表有哪些类型呢?
HTML 列表(List)可以将若干条相关的内容整理起来,让内容看起来更加有条理。在列表内您可以放置文本、图像、链接等,也可以在一个列表中定义另一个列表(列表嵌套)。HTML 为我们提供了三种不同形式的列表:有序列表,使用 &l如何在windows下安装GIT
本文在Windows7下测试成功。安装和设置Git下载Git for Windows,采用默认安装,安装完成后就可以在本地使用Git了。但要将内容放到Github上,必须先在Github网站上注册个账户,然后在本机使用Git创建SSH Kehtml语言背景图片拉伸代码是什么?
分类:电脑网络 >>程序设计 >>其他编程语言解析:在引用图片时使用<body bgproperties="fixed" backgrouR语言绘图——数据可视化ggplot2 介绍和主要的参数
R 有几种用于制作图形的系统,但 ggplot2 是最优雅和最通用的系统之一。与大多数其他图形包不同,ggplot2 具有基于图形语法的底层语法,它允许您通过组合独立组件来组合图形。如果想要更加了解ggplot2,请阅读 ggplot2:用HTML做的照片墙怎么用手机打开
输入QQ号码和密码,点击登录登录成功后点击正下方“动态”找到并点击“好友动态”进入好友动态,点击正上方的背景图片进入背景商城,找到并点击“照片墙”进入照片墙,点击右上角的“添加”,你可以选择从手机里或者从空间相册里添加照片Microsoft求助R语言如何运算permutation test
原始数据HAPLOTYPE CASEFREQUENCE CONTROLFREQUENCE 10.412027 0.421346 20.206045 0.233228 30.349419 0.296320 40.032509CSS如何嵌入到HTML中?
CSS通过内联、外联等样式嵌入到HTML中。1、使用样式表的处理指令语句在HTML文档的开头部分写一个关于样式表的指令处理语句复制代码代码如下:<?xml-stylesheet type="textcss&qJAVA能在WINDOWSXP中运行吗?
java可以在几乎现在的所有的操作系统中运行,java程序是运行在java虚拟机之上的,所以不管是什么样的操作系统只要安装了java的虚拟机(jvm)就可以运行java程序。具体可以参考sun的网站:www.sun.com ,也可以在sun怎样用css写出圆形边框
1、首先打开sublime text编辑器,新建一个html文件,里面写入一个p标签:2、然后设置p标签的样式,这里先设置一个边框,然后设置圆角边框,主要使用CSS3属性border-radius属性定义圆角效果。其中的数值为参数lengt小米RUBY电脑如何设置锁屏密码
你好。MIUI系统并没有直接设置开机密码的选项。你可以通过设置屏幕密码和开机锁屏来实现。屏幕密码:系统设置-安全和隐私-屏幕安全保护开机锁屏:系统设置-开发者选项-直接进入系统-关闭。只能去服务器更改了。更改步骤如下:1. 修改 安装目录csgo语音按键是哪个
csgo语音按键是:Csgo游戏中默认的语音按键是【K】按键。长按住【k】的时候你就可以跟队友交流了。设置语音的话我们首先需要打开电脑中的【csgo】。然后在游戏界面中找到【设置】,然后我们点击进去。接着在设置中找到【键盘鼠标】。然后在【有关运用JQ为CSS添加样式
JQ为CSS添加样式,可以使用jQuery CSS 操作 的css() 方法定义和用法css() 方法返回或设置匹配的元素的一个或多个样式属性。设置 CSS 属性设置所有匹配元素的指定 CSS 属性。$(selector).css(nameJS中typeof的用法
js是一门弱语言,它在声明变量时无需确定变量的类型,js在运行时会自动判断。那么如何判断一个变量的类型呢,js提供了typeof运算符,用来检测一个变量的类型。 typeof是一个运算符,有2种使用方式:type