正在用R写毕业论文,参考了清华大学出版社,许启发等人编著的《R软件及其在金融定量分析中的应用》,过程中遇到代码出现Error in solve.default(B) : Lapack routine dgesv: system is exactly singular: U[1,1] = 0
的问题,看了很多国外论坛都没有看懂,最后看的是vars包的作者原版论文,他里面在做脉冲响应时和许启发那本书在设置LR和SR的方法上有所区别。
许启发版:
LR<-matrix(NA,nrow=4,ncol=4)
SR<-matrix(NA,nrow=4,ncol=4)
SR[1,2:4]<-c(0,0,0)
SR[2,1]<-0
SR[2,3:4]<-c(0,0)
SR[3,1:2]<-c(0,0)
SR[3,4]<-0
SR[4,1:3]<-c(0,0,0)
用这样设置的后果就是出现
Error in solve.default(B) :
Lapack routine dgesv: system is exactly singular: U[1,1] = 0
这好像是说矩阵正奇异。。。什么鬼不懂。而vars包的作者用的是另外一种设置矩阵的方法
SR[4,2]<-0
LR[1,2:4]<-0
LR[2:4,4]<-0
虽然我感觉许启发版本比较有道理,但怎么运行都做不出来,最后用的vars包作者的方法算出来了。大家可以参考一下。