Mercurial > repos > clifinder > clifinder

changeset 18:3471a9785135 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/GReD-Clermont/CLIFinder/ commit e141f9b9038a22c92f827daa298bcc53356b4bdc"
author clifinder
date Wed, 26 Feb 2020 19:07:36 -0500
parents 3e3370387441
children 5c16e8ddff78
files script/CLIFinder_results.R
diffstat 1 files changed, 3 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/script/CLIFinder_results.R	Fri Feb 14 17:03:24 2020 -0500
+++ b/script/CLIFinder_results.R	Wed Feb 26 19:07:36 2020 -0500
@@ -110,7 +110,7 @@
   
   for (i in 0:nfastq)
   {
-    mcols(grfus)[grfus$dup>1, 11+i] <- mcols(grfus)[grfus$dup>1, 11+i] + mcols(grchim2)[position[duplicated(position[, 1]), 2], 9+i]
+    mcols(grfus)[position[duplicated(position[, 1]), 1], 11+i] <- mcols(grfus)[position[duplicated(position[, 1]), 1], 11+i] + mcols(grchim2)[position[duplicated(position[, 1]), 2], 9+i]
   }
   
   grfus2 <- grfus
@@ -123,9 +123,9 @@
   suppressWarnings(position3 <- as.data.frame(findOverlaps(gr3[gr3$dup>1], grfus2)))
   grfus$info [position3[, 2]] <- "overlap sens opposé"
   
-  min <- ddply(as.data.frame(grchim2), .(seqnames, end, strand), function(x)x[x$L1.start == min(x$L1.start), ])
+  min <- ddply(as.data.frame(grfus), .(seqnames, end, strand), function(x)x[x$L1.start == min(x$L1.start), ])
   min <- ddply(data.frame(min), .(seqnames, start, strand), function(x)x[x$L1.start == min(x$L1.start), ])
-  max <- ddply(as.data.frame(grchim2), .(seqnames, end, strand), function(x)x[x$L1.end == max(x$L1.end), ])
+  max <- ddply(as.data.frame(grfus), .(seqnames, end, strand), function(x)x[x$L1.end == max(x$L1.end), ])
   max <- ddply(data.frame(max), .(seqnames, start, strand), function(x)x[x$L1.end == max(x$L1.end), ])
   
   grfus1 <- as.data.frame(grfus)