# HG changeset patch # User fubar # Date 1417133890 18000 # Node ID 999d4b5939bb551e1c49fba22a6f9a89069199af # Parent 23afddb3d1e5627727f149b569564cf8bd1b3715 Uploaded diff -r 23afddb3d1e5 -r 999d4b5939bb rgedgeRpaired_nocamera.xml --- 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))