File: Figure_for_sequence_iteration.R

package info (click to toggle)
pyfastx 2.3.0-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 1,336 kB
  • sloc: ansic: 4,826; python: 1,816; sh: 505; perl: 66; makefile: 31
file content (48 lines) | stat: -rw-r--r-- 1,276 bytes parent folder | download | duplicates (3)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
library(ggplot2)
library(patchwork)

df1 <- read.table("matrix_fasta_sequence_iteration.tsv", header = TRUE, sep = "\t")

#df2 <- df1[which(df1$tool != 'pyfasta'),]
#df3 <- df1[which(df1$tool == 'pyfasta'),]

p1 <-
  ggplot(df1, aes(
    x = reorder(genome, size),
    time,
    colour = tool,
    group = tool
  )) +
  geom_point(aes(size = memory)) +
  geom_line() +
  ylim(0,500) +
  facet_grid(rows = vars(tool)) +
  theme_bw() +
  theme(axis.text.x = element_text(angle = 30, hjust = 1),
        axis.title.x = element_blank()) +
  ylab("Elapsed time (s)") +
  guides(colour = guide_legend(title = "Tools"),
         size = guide_legend(title = "Memory (GB)"))

df2 <- read.table("matrix_fastq_sequence_iteration.tsv", header = TRUE, sep = "\t")
p2 <-
  ggplot(df2, aes(
    x = reorder(file, fsize),
    time,
    colour = tool,
    group = tool
  )) +
  geom_point(aes(size = memory)) +
  geom_line() +
  ylim(0,2500) +
  facet_grid(rows = vars(tool)) +
  theme_bw() +
  theme(axis.text.x = element_text(angle = 30, hjust = 1),
        axis.title.x = element_blank()) +
  ylab("Elapsed time (s)") +
  guides(colour = guide_legend(title = "Tools"),
         size = guide_legend(title = "Memory (GB)"))

pp <- p1 | p2
pp <- pp + plot_annotation(tag_levels = 'A')
show(pp)