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 列。