Browse code

Check that genome argument is provided in GsnapParam. No way to get a safe default for this one.

git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/gmapR@71082 bc3139a8-67e5-0310-9ffc-ced21a209358

Michael Lawrence authored on 09/11/2012 12:54:49
Showing 1 changed files

... ...
@@ -154,6 +154,8 @@ GsnapParam <- function(genome, unique_only = FALSE,
154 154
                        split_output = !unique_only,
155 155
                        novelsplicing = FALSE, splicing = NULL, 
156 156
                        nthreads = 1L, part = NULL, batch = "2", ...) {
157
+  if (missing(genome))
158
+    stop("The 'genome' must be specified (should be coercible to GmapGenome)")
157 159
   args <- formals(sys.function())
158 160
   params <- mget(names(args), environment())
159 161
   params$unique_only <- NULL