R语言绘制结果
数据链接:
1. 数据读取和数据整合:
library(tidyverse)
library(readxl)
data <- read_xlsx("data_1.xlsx",col_names = TRUE,range = "B3:Z1151")
data_extracted<-extract(data,`Taxonomic classification`, into=c('Domain', 'Phylum', 'Genus'),'k__([A-z]*);p__(.*);c__.*g__(.*);', convert=TRUE)
data_extracted[which(data_extracted$Phylum == "",arr.ind = TRUE),]$Phylum="Unknown"
data_extracted[which(data_extracted$Genus == "",arr.ind = TRUE),]$Genus="Unknown"
data_sum<-data_extracted %>% gather(wipe,abundant,-Domain,-Phylum,-Genus)
data的数据结构:
>str(data) ##该命令可用于查看数据结构
Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 1148 obs. of 25 variables:
$ Taxonomic classificat