如图所示,if代码一直在出问题,Error in if (K < 2) { : argument is of length zero。
代码应该是在这里出问题的
pdshare.computation <- function(x, nlag) {
cointest <- ca.jo(x, K = nlag, type = "eigen", ecdet = "const",
spec = "transitory")
k <- cointest@lag
vecm <- cajorls(cointest)
varm <- vec2var(cointest)
vma <- Psi(varm)
beta.ort <- as.vector(c(-cointest@V[2,1], cointest@V[1,1]))
alpha.ort <- as.vector(c(-cointest@W[2,1], cointest@W[1,1]))
aa <- bb <- cc <- dd <- 0
for (i in 1:(k-1)) {
aa <- aa + vecm$rlm$coefficients[2*i,1]
bb <- bb + vecm$rlm$coefficients[2*i+1,1]
cc <- cc + vecm$rlm$coefficients[2*i,2]
dd <- dd + vecm$rlm$coefficients[2*i+1,2]
}
但具体是什么问题,一直没弄明白。多谢!


雷达卡





京公网安备 11010802022788号







