R代码 | 检查H的列中的基因名是否也存在于另一个数据框（C_unique）的行名中

最新推荐文章于 2024-10-07 00:05:15 发布

2301_77397699

最新推荐文章于 2024-10-07 00:05:15 发布

阅读量106

点赞数 1

文章标签： r语言

本文链接：https://blog.csdn.net/2301_77397699/article/details/140525032

版权

# 假设C_unique是一个数据框，其行名包含基因名  
row_names_C_unique <- rownames(C_unique)  

# 初始化一个空向量来存储不在C_unique行名中的基因名  
genes_not_in_C_unique <- vector("character")  

# 遍历H的每一列  
for (col in seq_along(H)) {  
  # 获取当前列的非NA基因名  
  current_col_genes <- H[[col]][!is.na(H[[col]])]  
  
  # 检查这些基因名是否不在C_unique的行名中  
  not_in_C_unique <- current_col_genes[!current_col_genes %in% row_names_C_unique]  
  
  # 如果找到这样的基因名，将它们添加到结果向量中  
  if (length(not_in_C_unique) > 0) {  
    genes_not_in_C_unique <- c(genes_not_in_C_unique, not_in_C_unique)  
  }  
}  

# 去除结果中的重复项（如果需要）  
genes_not_in_C_unique <- unique(genes_not_in_C_unique)  

# 输出结果  
print(genes_not_in_C_unique)