我喜欢提供更多信息。我从 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
}
然后我加载gtools
和heirpart
.
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:找不到对象'输出'