如何在R中实现stata里的varlmar，varwle以及vargranger检验？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

R的vars包里没找到与stata相对应的检验，查看stata的说明文件后发现varwle和vargranger都是用Wald test。尝试用restrict拟合restricted model，然后调用lmtest包里的waldtest，得出的Chisq都比stata的结果低10%左右，p-value也不大相同，求助原因。
另外stata的varlmar采用的是Johansen的LM test，R的vars包里的serial.test没有包含此检验，请问如何实现？
万分感谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：Granger Grange Anger range Stata 如何

帮顶

这问题有意思

需要请联系QQ：1298997509

下面是协整的Johansen的LM test程序，看看没有么帮助：

johansen.test<-
function(x, L = 2)
{
        x <- ts(x)
        n <- nrow(x)
        p <- ncol(x)
        Ly <- lag(x[, 1], -1)
        D <- diff(x[, 1])
        for(i in 1:p) {
                if(i > 1) {
                        D <- ts.intersect(D, diff(x[, i]))
                        Ly <- ts.intersect(Ly, lag(x[, i], -1))
                }
                if(L > 0)
                        for(j in 1:L)
                                D <- ts.intersect(D, lag(diff(x[, i]),  - j))
        }
        iys <- 1 + (L + 1) * (0:(p - 1))
        Y <- D[, iys]
        X <- D[,  - iys]
        Ly <- ts.intersect(Ly, D)[, 1:p]
        ZD <- lm(Y ~ X)$resid
        ZL <- lm(Ly ~ X)$resid
        df <- nrow(X) - ncol(X) - 1
        S00 <- crossprod(ZD)/df
        S11 <- crossprod(ZL)/df
        S01 <- crossprod(ZD, ZL)/df
        M <- solve(S11) %*% t(S01) %*% solve(S00) %*% S01
        eigen(M)$values   #返回特征值
}

关于该检验的部分理论见http://en.wikipedia.org/wiki/Johansen_test

grangertest: grangertest is simply a Wald test comparing the unrestricted model—in which y is explained by the lags (up to order order) of y and x—and the restricted model—in which y is only explained by the lags of y.
causality:Two causality tests are implemented. The first is a F-type Granger-causality test and the second is a Wald-type test that is characterized by testing for nonzero correlation between the error processes of the cause and effect variables.

感谢分享。。。

hubifeng？ 发表于 2014-9-2 23:12
grangertest: grangertest is simply a Wald test comparing the unrestricted model—in which y is expla ...

是说wald type test是检验相关性是否不为零，是这个意思吗？
在eviews里怎么实现呢？需要先做var吗?
请指教，多谢啦~