dismiss_macs2_extractor.R

# Copyright (C) 2014-2015  Umar Niazi
#   
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
# 
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
# 
# You should have received a copy of the GNU General Public License
# along with this program.  If not, see <http://www.gnu.org/licenses/>

# File: dismiss_macs2_extractor.R
# Auth: umn@aber.ac.uk, uhkniazi@gmail.com
# DESC: main script called from command line
#       1 - takes excel file name and bam file names are command line arguments
#       arg1 = excel file name
#       arg2 = bam file name (index file for bam should be in same directory)
#       arg3 = p (for paired end data) or s (for single end data)
# Date: 24/10/2014


# get command line arguments
args = commandArgs(trailingOnly = TRUE)
csHelp = 'Run script with: \nRscript dismiss_macs2_extractor.R MACS2_result.xls File.bam [p or s]
the bam index file should be in same directory with same name and .bai extension e.g. File.bam.bai'

# check argument size
if (args[1] == 'h' || args[1] == 'help' || length(args) != 3) stop(paste(csHelp))

# variable setting and some error checking
macs = args[1]
bam = args[2]
paired = NA
if (args[3] == 'p') { paired = TRUE
  } else if (args[3] == 's') paired = FALSE
if (is.na(paired)) stop(paste(csHelp))

# header file with functions
source(file='dismiss_header.R')

oGRdismiss = f_oGRMACS2toGRanges(macs)
oGRdismiss = f_oGRSeparateStrands(oGRdismiss, bam, paired)

# save output to granges object 
name.out = paste('oGRdismiss_', make.names(date()), '.rds', sep='')
save(oGRdismiss, file = name.out)

# save output as gff3 files
f_WriteGff3(oGRdismiss)