用R包 “vars”做B型SVAR模型(A默认为Ik)时,可以出结果;
做AB型SVAR模型时,将A设为Ik,遇到了Error提示:Error in `[<-`(`*tmp*`, pos, i, value = 1) : subscript out of bounds
请问一下,默认设置A为Ik的B型SVAR和手动设置的AB型SVAR有不同的地方吗?
类似的问题在网络上有一个资料:http://afni.nimh.nih.gov/afni/community/board/read.php?1,142479,142507#msg-142507
回答者认为是识别矩阵提供的约束不足,但是,我的模型中提供的约束是足够的——双变量模型K=2,A model的约束有1个,AB model的约束有5个。
具体描述:
出问题的具体的代码如下:
- amat <- matrix(c(1,0,0,1),nrow = 2, ncol = 2)
- bmat <- matrix(c(NA,0,NA,NA),nrow = 2, ncol = 2
- ## 得到VAR模型
- p1ct <- VAR(temp, p, type = "none")
- ## 构建SVAR模型
- sp1ct <- SVAR(p1ct, Amat = amat, Bmat = bmat, max.iter = 1000)
请问一下大家有没有遇到类似的问题,能够解答一下我的困惑,谢谢!!!