data <- data.frame(
AA_1 = sample(0:10, 10, replace = TRUE),
AA_2 = sample(0:10, 10, replace = TRUE),
BB_1 = sample(0:10, 10, replace = TRUE),
BB_2 = sample(0:10, 10, replace = TRUE)
)
geneToChoose <- unique(unlist(data))
geneToChoose <- geneToChoose[geneToChoose != 0]
print(geneToChoose)
resultSet <- numeric(0)
while(length(geneToChoose)!=0){
index <- 0
maxNumber <- 0
for (i in 1:nrow(data)) {
maxNumberTemp <- length(intersect(geneToChoose, unique(unlist(data[i,]))))
if(maxNumberTemp > maxNumber){
maxNumber <- maxNumberTemp
index <- i
}
}
geneToChoose <- setdiff(geneToChoose, unique(unlist(data[index,])))
resultSet <- c(resultSet, index)
}
print(resultSet)