plot_regions.R 355 Bytes
Newer Older
Beatrice Tan's avatar
Beatrice Tan committed
1
2
3
4
5
6
7
8
9
10
11
12
13
#!/usr/bin/env Rscript
args = commandArgs(trailingOnly=TRUE)

source("https://bioconductor.org/biocLite.R")
biocLite("karyoploteR")

bed = read.table(args[1], sep="\t", header=TRUE)
#Read output plot files?
#chrom as "chr1"
kp <- plotKaryotype(chromosomes=c(chrom))
#kpDataBackground(kp)
kpAxis(kp)
kpRect(kp, chr=chrom, x0=start, x1=end, y0=0.2, y1=0.4)