inst/extdata/Rscript/rscript_chromosome_job.R

In MPIIComputationalEpigenetics/MAGAR: MAGAR: R-package to compute methylation Quantitative Trait Loci (methQTL) from DNA methylation and genotyping data

suppressPackageStartupMessages(library(argparse))
suppressPackageStartupMessages(library(MAGAR))

ap <- ArgumentParser()
ap$add_argument("-m","--methQTL",action="store",help="methQTL object to be used")
ap$add_argument("-j","--json",action="store",help="Configuration JSON file")
ap$add_argument("-c","--chr",action="store",help="chromosome to be used")
ap$add_argument("-u","--covariates",action="store",default=NULL,help="covariates to be included")
ap$add_argument("-p","--p.val",action="store",default=1e-5,help="p-value cutoff")
ap$add_argument("-o","--output",action="store",help="Output directory")
ap$add_argument("-n","--ncores",action="store",default=1,help="Number of cores to be used")
ap$add_argument("-d","--hdf5dir",action="store",help="The HDF5 dump directory")
ap$add_argument("-f","--ffdir",action="store",help="The ff dump directory")
cmd.args <- ap$parse_args()

logger.start(paste("Running on:",Sys.info()["nodename"]))

logger.start("Configuring job")
qtlJSON2options(cmd.args$json)
logger.completed()

logger.start("Loading methQTL object")
meth.qtl <- loadMethQTLInput(cmd.args$methQTL)
logger.completed()

dump.file <- file.path(cmd.args$hdf5dir,cmd.args$chr)
dir.create(dump.file)
setHDF5DumpDir(dump.file)
options(fftempdir=cmd.args$ffdir)

if(!is.null(cmd.args$covariates)){
  logger.start("Reading covariates")
  covs <- unlist(readLines(cmd.args$covariates))
  logger.completed()
}else{
  covs <- NULL
}

p.val <- as.numeric(cmd.args$p.val)
ncores <- as.numeric(cmd.args$ncores)

methQTL.res <- doMethQTLChromosome(meth.qtl,cmd.args$chr,sel.covariates = covs,p.val.cutoff = p.val, out.dir=cmd.args$output,
                                     ncores=ncores)

logger.start("Saving results")
path.save <- file.path(cmd.args$output,paste0("MethQTLResult_",cmd.args$chr))
saveMethQTLResult(methQTL.res,path.save)
logger.completed()