# 清除工作空间
rm(list = ls())# 设置工作目录
setwd("C:\\Users\\Administrator\\Desktop\\machine learning\\Heatmap\\CAZy")# 读取数据
otu <- read.table("TPM.CAZy.raw.txt", sep = "\t", header = TRUE, row.names = 1)
lda <- read.table("1.txt", sep = "\t")# 假设 lda 文件中只有一列，提取这一列作为匹配对象
lda_items <- lda[[1]] # 提取第一列# 根据 lda 的内容匹配 otu 的行名，保留匹配到的行
matched_rows <- otu[rownames(otu) %in% lda_items, ]# 将匹配结果输出到一个新的 txt 文件
write.table(matched_rows, file = "matched_otu.txt", sep = "\t", row.names = TRUE, col.names = NA, quote = FALSE)cat("匹配完成，结果已保存为 matched_otu.txt 文件。\n")
