-2

我喜欢提供更多信息。我从 4 个基因开始,并使用以下脚本定义了它们的开启和关闭条件:

Proli <- function (DISC1, GSK3B, DIXDC1, CTNNB1){     
  inputs <- permutations(2,4,v=c(0,1),repeats.allowed=TRUE);
  if (DISC1 == 1 && GSK3B == 0 && DIXDC1 == 1 && CTNNB1 == 1){
    Proliferation <- "TRUE"  
  }
  else
  {
    Proliferation <- "FALSE"
  }
  Proliferation
}

然后我加载gtoolsheirpart.

inputs <- permutations(2,26,v=c(0,1),repeats.allowed=TRUE)
len <- length[inputs]
for(i in 1:len) 
{
  if(inputs[i,1] == 1 && inputs[i,2] == 0 && inputs[i,3] == 1 && inputs[i,4] == 1){
    #Constructing  a Truth Table
    output[i] <- 1
  }
  else{
    output[i] <- 0
  }
}

现在我得到错误输出[i] <- 0:找不到对象'输出'

4

1 回答 1

-1

#只需添加到您的脚本之上

len <- lenght(object on with you do your loop)

或者

len <- nrow(object on with you do your loop)



for(i in 1:len) {
        if(inputs[i,1] == 1 && inputs[i,2] == 0 && inputs[i,3] == 1 && inputs[i,4] == 1){
            #Constructing  a Truth Table
            output[i] <- 1
        }
        else{
            output[i] <- 0
        }
    }
于 2016-03-01T08:16:54.300 回答