# HG changeset patch
# User fubar
# Date 1416745085 18000
# Node ID 6ded550dcc7d6d8b8ed5b9f3b97331dd3340ebf6
# Parent  d69a656284f526e2bfe7fc556d881a3849eab11d
Uploaded

diff -r d69a656284f5 -r 6ded550dcc7d rgedgeRpaired_nocamera.xml
--- a/rgedgeRpaired_nocamera.xml	Sun Nov 23 07:06:23 2014 -0500
+++ b/rgedgeRpaired_nocamera.xml	Sun Nov 23 07:18:05 2014 -0500
@@ -899,7 +899,7 @@
 Count_Matrix = read.table(Input,header=T,row.names=1,sep='\t') 
 snames = colnames(Count_Matrix)
 nsamples = length(snames)
-if (nsubj >  0 & nsubj != nsamples) {
+if (nsubj > 0 & nsubj != nsamples) {
 options("show.error.messages"=T)
 mess = paste('Fatal error: Supplied subject id list',paste(subjects,collapse=','),
    'has length',nsubj,'but there are',nsamples,'samples',paste(snames,collapse=','))
@@ -907,7 +907,7 @@
 quit(save="no",status=4)
 }
 if (length(subjects) != 0) {subjects = subjects[useCols]}
-Count_Matrix = Count_Matrix[,useCols] ### reorder columns
+Count_Matrix = Count_Matrix[,useCols]
 rn = rownames(Count_Matrix)
 islib = rn %in% c('librarySize','NotInBedRegions')
 LibSizes = Count_Matrix[subset(rn,islib),][1] # take first