问题遇到的现象和发生背景
刚刚开始学习R 不久,在对大矩阵A分块运行的时候遇到这么一个问题,代码运行老是报错,不知道错误出在哪里,以下是程序:
问题相关代码,请勿粘贴截图
i<-1 #第一行处理
m<-matrix(0,nrow=30,ncol=30)
I<-diag(30)
for (j in 1:31) {
if (j==i) Int1_<-matrix(0,nrow=30,ncol=1)
else {
Aij<-A[(1+30*(i-1)):(30i),(1+30(j-1)):(30j)]
Ajjn<-solve(I-A[(1+30(j-1)):(30j),(1+30(j-1)):(30j)])
Int1_<-Aij%%Ajjn%%Y_pro[(1+30(j-1)):(30*j)]
}
Int1_<-as.matrix(Int1_)
m[,j]<-Int1_
}
运行结果及报错内容
运行后报错:
Error in [<-
(*tmp*
, , j, value = Int1_) : subscript out of bounds
运行之后有结果,同时也报错了,不知道啥原因,
请问大家知道哪里出错了嘛?非常感谢