R循环错误：下标超出范围

Question

我试图做一个循环来探索模型组合。 因此，我正在执行以下操作：

# k=ncol(tmp1)-1
k = 5  ## example k without knowing "tmp1"

#number of models
M=2^k
#a matrix to store indicators/switch for predictors. Need two loops to do this (for loops)
incl_pred<-matrix(0,M,4)
tmpm=M
i=1
for(i in 1:k){
  tmpm=tmpm/(2)
  incl_pred[1:tmpm,i]=1
}
for(i in 2:k){
  tmpm=M
  for(j in i:k){
    incl_pred[((tmpm/2)+1):tmpm,j]=incl_pred[1:(tmpm/2),j]
    tmpm=tmpm/2
  }
}

但这给了我两个错误：

[<- （ *tmp* ，1：tmpm，i，value = 1）中的错误：下标超出范围

和

incl_pred [1：（tmpm / 2），j]中的错误：下标超出范围。

你能帮我解决吗？

谢谢！

Answer 1

您应该设置：

incl_pred<-matrix(0,M,k)

因为你的内循环：

  for(j in i:k){
     incl_pred[((tmpm/2)+1):tmpm,j]=incl_pred[1:(tmpm/2),j]

正在访问incl_pred k列。