gff2svg/gffdiff.R at master · OpenScienceRender/gff2svg · GitHub

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
#! /usr/bin/Rscript

args <- commandArgs()
url <- args[-(1:match("--args", args))]

## url <- "https://raw.githubusercontent.com/OpenScienceRender/gff2svg/master/data/my.gffdiff"

stopifnot(require("rtracklayer"))
stopifnot(require("ggbio"))

download.file(url, destfile = "gffdiff.tmp", method = "wget")

gffls0 <- readLines("gffdiff.tmp")
gffls <- paste0(dirname(url), "/", gffls0)

sapply(gffls, function(f)
       download.file(f, dest = basename(f), method = "wget"))

gffl <- lapply(gffls0, function(f) {
    gff <- import(f)
    mcols(gff)$file <- f
    gff
})

gff <- do.call(c, gffl)

diff <- compare(gffl[[1]], gffl[[2]]) == 0
different <- rep(diff, 2)

svg("output.svg")
autoplot(gff, aes(fill = different), facets = file ~ .) +
    theme(legend.position="none")
dev.off()

unlink(c("gffdiff.tmp", gffls0))