Mercurial > repos > davidvanzessen > shm_csr

diff shm_csr.r @ 62:aa8d37bd1930 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 05 Dec 2017 10:57:13 -0500
parents 77a7ac76c7b9
children
line wrap: on
line diff
--- a/shm_csr.r	Fri Aug 18 10:43:31 2017 -0400
+++ b/shm_csr.r	Tue Dec 05 10:57:13 2017 -0500
@@ -10,7 +10,9 @@
 empty.region.filter = args[4]
 setwd(outputdir)
 
-dat = read.table(input, header=T, sep="\t", fill=T, stringsAsFactors=F)
+#dat = read.table(input, header=T, sep="\t", fill=T, stringsAsFactors=F)
+
+dat = data.frame(fread(input, sep="\t", header=T, stringsAsFactors=F)) #fread because read.table suddenly skips certain rows...
 
 if(length(dat$Sequence.ID) == 0){
   setwd(outputdir)