Load libraries
library(tidyverse)
Import data
df<-read_delim(file="filtering.csv", delim=",")
##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## Sample = col_character(),
## Reads = col_double(),
## NoHost = col_double(),
## Filtered = col_double()
## )
head(df)
## # A tibble: 6 x 4
## Sample Reads NoHost Filtered
## <chr> <dbl> <dbl> <dbl>
## 1 Sample13 1000000 989840 938220
## 2 Sample22 1000000 758380 678168
## 3 Sample25 1000000 972564 913842
## 4 Sample30 1000000 944690 879807
## 5 Sample31 1000000 983343 901878
## 6 Sample3 1000000 989295 915797
Plot reads
df_gg1<-df %>%
gather(group,reads, -Sample)
head(df_gg1)
## # A tibble: 6 x 3
## Sample group reads
## <chr> <chr> <dbl>
## 1 Sample13 Reads 1000000
## 2 Sample22 Reads 1000000
## 3 Sample25 Reads 1000000
## 4 Sample30 Reads 1000000
## 5 Sample31 Reads 1000000
## 6 Sample3 Reads 1000000
ggplot(df_gg1, aes(x=Sample, y=reads)) +
geom_point(aes(color=group, shape=group), size=2, alpha=0.8) +
theme(axis.text.x = element_text(angle=45, hjust=1))