Mercurial > repos > fubar > differential_count_models
diff rgedgeRpaired_nocamera.xml @ 122:6ded550dcc7d draft
Uploaded
author | fubar |
---|---|
date | Sun, 23 Nov 2014 07:18:05 -0500 |
parents | d69a656284f5 |
children | 51f998262ada |
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