fix verbose multiclass, add header to location data

guigolab · Nov 17, 2021 · 9525e63 · 9525e63
1 parent 0ac689c
commit 9525e63
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 1 deletion.
diff --git a/R/sqtl.seeker.R b/R/sqtl.seeker.R
@@ -131,7 +131,7 @@ sqtl.seeker <- function(tre.df, genotype.f, gene.loc, covariates = NULL,
                 covariates <- covariates[com.samples, , drop = FALSE]
                 multiclass <- apply(covariates, 2, function(x){length(table(x)) > 1})
                 covariates <- covariates[, multiclass, drop = FALSE]
-                if (verbose){
+                if (verbose & any(!multiclass)){
                   message("\t", "Covariates removed due to only one value: ", 
                           paste(names(multiclass)[!multiclass], collapse = ", "))
                 }

diff --git a/data/genes.bed b/data/genes.bed
@@ -1,3 +1,4 @@
+chr	start	end	geneId
 19	38826415	38861589	ENSG00000099338.14
 19	52901102	52921657	ENSG00000167555.8
 19	52956829	53016982	ENSG00000258405.3