R语言中用做脉冲响应时遇到Error in solve.default(B) : 解决个案

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

正在用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包作者的方法算出来了。大家可以参考一下。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Default Solve fault Error 脉冲响应 Error