2. .libPaths(c("/home/.opt/sysoft/R-4.0.2/lib64/R/library", "/home/shwzhao/R/x86_64-pc-linux-gnu-library/4.0"))
    4. library(tidyverse)
    5. library(pheatmap)
    3. TPM <- read_tsv("Quercus_acutissima.pep.tpm2.SY.tsv")
    4. ERFfam <- read_tsv("ERF.2.tsv", col_names = c("Geneids", "subFam"))
    6. TPM2 <- TPM %>%
    7.   pivot_longer(-Geneids) %>%
    8.   mutate(name = str_replace(name, "-.*", "")) %>%
    9.   group_by(Geneids, name) %>%
    10.   mutate(value = sum(value)/3) %>%
    11.   distinct() %>%
    12.   pivot_wider(names_from = "name", values_from = "value")
    14. ERFfamTPM <- ERFfam %>%
    15.   left_join(TPM2, by = "Geneids") %>%
    16.   mutate(subFam = factor(subFam, levels = c("I", "II", "III", "IV", "V", "VI", "VII", "VIII", "IX", "X"))) %>%
    17.   arrange(subFam) %>%
    18.   select(-subFam) %>%
    19.   column_to_rownames("Geneids") %>%
    20.   as.matrix()
    22. annotation_row <- ERFfam %>%
    23.   column_to_rownames("Geneids")
    25. pdf("test.2.pdf", width = 8, height = 15)
    26. pheatmap(ERFfamTPM,
    27.          scale = "row",
    28.          cluster_rows = F,
    29.          cluster_cols = F,
    30.          # display_numbers = TRUE,
    31.          # number_format = "%1.0f",
    32.          # gaps_row = c(15, 30, 45),
    33.          # gaps_col = c(3, 6, 9, 12, 15),
    34.          annotation_row = annotation_row,
    35.          # annotation_col = annotation_col
    36. )
    37. dev.off()