check file isn't empty

MyersGroup · Mar 10, 2019 · 486d3c1 · 486d3c1
1 parent 3cdd301
commit 486d3c1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/scripts/functions.R b/scripts/functions.R
@@ -85,7 +85,7 @@ compute_likelihood = function(r1, r2, g, alpha1, alpha2, beta, return_slim=TRUE)
 
 get_frag_overlap_counts <- function(posfiles = c(posfileA, posfileB, posfileG), infiles = c(infile1, infile2, infile3), windowfilepath=windowfile, chromosomes=chrs, bedtools=btpath, recalculate_coverage=TRUE){
 
-  if(all(file.exists(infiles)) & !recalculate_coverage){
+  if(all(file.exists(infiles)) & all(file.size(infiles)>1) & !recalculate_coverage){
     print("Skipping recalculation of coverage")
   }else{
     print("Calculating coverage")