for(i in 1:(length(h$counts)-1)){
if(h$counts[i]!=0 || h$counts[i+1]!=0){
counts2[j]<-h$counts[i]
- breaks2[j+1]<-h$breaks[i+1];
+ breaks2[j+1]<-h$breaks[i+1]
j<-j+1
}
}
-
+
counts2[j]<-h$counts[length(h$counts)]
breaks2[j+1]<-h$breaks[length(h$breaks)]
h<-merge_empty_bins(h)
}
}
-
+
cat(i, file = outputfile, sep = "\t", append = TRUE)
cat("\t", file = outputfile, append = TRUE)
cat(sum(h$counts), file =outputfile, sep = "\t", append = TRUE)