刚刚开始学习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)):(30*i),(1+30*(j-1)):(30*j)]
Ajjn<-solve(I-A[(1+30*(j-1)):(30*j),(1+30*(j-1)):(30*j)])
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
运行之后有结果,同时也报错了,不知道啥原因,
请问哪位前辈知道哪里出错了嘛?非常感谢