Mercurial > repos > dktanwar > test_csaw_1

changeset 3:ce3ad612a104 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author dktanwar
date Mon, 18 Dec 2017 11:20:06 -0500
parents ed929b85b051
children 609b5a2ae0e8
files csaw.R
diffstat 1 files changed, 4 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- a/csaw.R	Mon Dec 18 10:58:39 2017 -0500
+++ b/csaw.R	Mon Dec 18 11:20:06 2017 -0500
@@ -15,6 +15,7 @@
 library("stringr")
 library("data.table")
 library("getopt")
+library("Rsamtools")
 
 
 options(stringAsfactors = FALSE, useFancyQuotes = FALSE)
@@ -29,6 +30,9 @@
 # Separate multiple input files into a list of individual files
 files <- unlist(strsplit(inputs, ','))
 
+# Index bamfiles
+indexBam(files = files)
+
 # Create windows and count reads in them ----
 Sys.time()
 windows <- windowCounts(files, spacing=150, width=200, bin=F)