Mercurial > repos > fubar > differential_count_models

changeset 128:999d4b5939bb draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author fubar
date Thu, 27 Nov 2014 19:18:10 -0500
parents 23afddb3d1e5
children a95f5496d330
files rgedgeRpaired_nocamera.xml
diffstat 1 files changed, 3 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/rgedgeRpaired_nocamera.xml	Thu Nov 27 18:39:15 2014 -0500
+++ b/rgedgeRpaired_nocamera.xml	Thu Nov 27 19:18:10 2014 -0500
@@ -538,9 +538,10 @@
     adj.p.value=p.adjust(DE\$table\$PValue, method=fdrtype),
     Dispersion=myDGEList\$tagwise.dispersion,totreads=cmrowsums
   )
+  tt = cbind(tt,URL=contigurls) # add to end so table isn't laid out strangely
+  stt = tt[order(DE\$table\$PValue),]
   print.noquote("@@ edgeR Top tags\n")
-  tt = cbind(tt,URL=contigurls) # add to end so table isn't laid out strangely
-  print.noquote(tt[1:50,])
+  print.noquote(stt[1:50,])
   deTags = rownames(uoutput[uoutput\$adj.p.value < fdrthresh,])
   nsig = length(deTags)
   print.noquote(paste('@@',nsig,'tags significant at adj p=',fdrthresh))