Mercurial > repos > fubar > differential_count_models

changeset 122:6ded550dcc7d draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author fubar
date Sun, 23 Nov 2014 07:18:05 -0500
parents d69a656284f5
children 51f998262ada
files rgedgeRpaired_nocamera.xml
diffstat 1 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- 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