metabolite <- read.csv("D:/isoquinoline/benzylisoquinoline.csv")
rna_seq_tpm_normal <- read.csv(paste0("/xuyu_data/TPM.csv"),stringsAsFactors = F,header = T)
colnames(rna_seq_tpm_normal) <- c("gene","Capsule","Petal","Stem","Fine Root","Leaf","Tap Root")
rna_seq_tpm_normal <- rna_seq_tpm_normal[rna_seq_tpm_normal[,1] %in% metabolite[,1],]
row.names(rna_seq_tpm_normal) <- rna_seq_tpm_normal[,1]
rna_seq_tpm_normal <- rna_seq_tpm_normal[,-1]
library(pheatmap)
pheatmap(as.matrix(rna_seq_tpm_normal),scale = "row",cluster_row = FALSE,cluster_cols = TRUE)
pheatmap(log2(as.matrix(rna_seq_tpm_normal)+1),scale = "row",cluster_row = TRUE,cluster_cols = TRUE)