Mercurial > repos > shians > shrnaseq

diff hairpinTool.R @ 9:f1076bfb0ed1

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
- Fixed tool to actually make use of barcode location options - Stored variables now properly used in processHairpinReads function
author shian_su <registertonysu@gmail.com>
date Wed, 30 Jul 2014 15:13:45 +1000
parents 548802b3492f
children 8923d4ea858b
line wrap: on
line diff
--- a/hairpinTool.R	Fri May 02 17:22:24 2014 +1000
+++ b/hairpinTool.R	Wed Jul 30 15:13:45 2014 +1000
@@ -35,6 +35,7 @@
 #                            plots
 #       19.selectVals       -String specifying members selected for barcode
 #                            plots
+#    ###
 #
 # OUT:  Bar Plot of Counts Per Index
 #       Bar Plot of Counts Per Hairpin
@@ -57,7 +58,7 @@
 library(edgeR, quietly=TRUE, warn.conflicts=FALSE)
 library(limma, quietly=TRUE, warn.conflicts=FALSE)
 
-if (packageVersion("edgeR") < "3.5.23") {
+if (packageVersion("edgeR") < "3.5.27") {
   stop("Please update 'edgeR' to version >= 3.5.23 to run this tool")
 }
 
@@ -359,6 +360,7 @@
   # Use EdgeR hairpin process and capture outputs
   hpReadout <- capture.output(
   data <- processHairpinReads(fastqPath, samplePath, annoPath,
+                              barcodeStart=barStart, barcodeEnd=barEnd,
                               hairpinStart=hpStart, hairpinEnd=hpEnd, 
                               verbose=TRUE)
   )
@@ -470,6 +472,17 @@
 linkData <- rbind(linkData, newEntry)
 invisible(dev.off())
 
+# BCV Plot
+png(bcvPng, width=600, height=600)
+plotBCV(data, main="BCV Plot")
+imageData <- rbind(imageData, c("BCV Plot", "bcv.png"))
+invisible(dev.off())
+
+pdf(bcvPdf)
+plotBCV(data, main="BCV Plot")
+newEntry <- c("BCV Plot (.pdf)", "bcv.pdf")
+invisible(dev.off())
+
 if (workMode=="classic") {
   # Assess differential representation using classic exact testing methodology 
   # in edgeR