我想根据下面给出的代码填充矩阵
num_feats = 4726
num_obsvs = 500
mod_size = 35
num_causl = 35
W<-matrix(0,num_feats,num_feats)
for(i in 2:(as.integer(num_feats/mod_size)-1))
w[(i-1)*mod_size:(i)*mod_size,(i-1)*mod_size:(i)*mod_size]<-matrix(1,ncol= mod_size,nrow = mod_size)
if (i!=as.integer(num_feats/mod_size-1))
{
w[(i+1)*mod_size-1,(i+1)*mod_size]<-1
w[(i+1)*mod_size,(i+1)*mod_size-1]<-1
}
但是它给了我这个错误:
w [(i-1)* mod_size:(i)* mod_size,(i-1)* mod_size:(i)*中的错误 :要替换的项目数不是替换的倍数 长度
我尝试调整大小,但现在开始!我该如何解决?