任务目标:批量绘制每个RNA文库reads比对情况的饼图;
任务流程: 数据预处理 和 图样式处理 + 循环出图
library(RColorBrewer);library(ggforce);set.seed(123);;library(ggplot2);library(dplyr);library(tidyverse)
###数据集概况
加载进来的的数据集是按行记录了每个文库的reads比对信息,其中比对类别存在列向量,绘图注意数据格式转换,绘制一个文库的饼图需要提取数据集的一行来进行处理。
###数据处理脚本
- 第一步,绘图数据格式转换
data.set <- read.table("*.tsv") #加载数据集
df <- data.set[1,] #提取数据集的一行进行脚本测试
df %>% data.frame() %>% t() %>% data.frame() %>% tibble::rownames_to_column(var = "group") %>% filter_all( any_vars(grepl("Reads", .)) ) %>%
dplyr::rename(labels_1 = 2) %>% filter(!grepl("Mapped.to.Genome", group)) %>%
mutate(perc = as.numeric(sub("%", "", labels_1))/100) %>%
mutate(labels = scales::percent(perc),labels = paste( gsub("\\."," ",x = group) ,labels,sep = "-")) -> t
> t
group labels_1 perc