bug fix: GRCh38 in annotateRanges

gagneurlab · Aug 14, 2020 · eb45f63 · eb45f63
1 parent 5db7597
commit eb45f63
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/R/annotationOfRanges.R b/R/annotationOfRanges.R
@@ -55,6 +55,11 @@ annotateRanges <- function(fds, feature="hgnc_symbol", featureName=feature,
     # check input
     stopifnot(is(fds, "FraserDataSet"))
     if(length(fds) == 0) return(fds)
+
+    # useEnsembl only understands GRCh=37 or GRCh=NULL (uses 38 then)
+    if(GRCh == 38){
+        GRCh <- NULL
+    }
 
     if(is.null(ensembl)){
         tryCatch({