Browse code

handled case when GmapSnp obj is created with path to VCF, but no genome specified

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

Cory Barr authored on 21/08/2012 23:49:18
Showing 1 changed files

... ...
@@ -47,6 +47,10 @@ setReplaceMethod("snps", c("GmapSnpDirectory", "VCF"),
47 47
 setReplaceMethod("snps", c("GmapSnpDirectory", "character"),
48 48
                  function(x, name, genome, which, ..., value)
49 49
                  {
50
+                   if (missing(genome)) {
51
+                     stop("Please supply the \"genome\" argument")
52
+                   }
53
+
50 54
                    param <- ScanVcfParam(fixed = "ALT", info = NA, geno = NA)
51 55
                    if (!missing(which)) # FIXME: waiting for vcfWhich<-
52 56
                      param@which <- as(which, "RangesList")