[Patchwork-commits] r156 - pkg/patchwork/R www
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Mon Jan 21 13:12:35 CET 2013
Author: sebastian_d
Date: 2013-01-21 13:12:34 +0100 (Mon, 21 Jan 2013)
New Revision: 156
Modified:
pkg/patchwork/R/patchwork.plot.r
www/pw_exec.php
Log:
uncomment line to generate normalalf
Modified: pkg/patchwork/R/patchwork.plot.r
===================================================================
--- pkg/patchwork/R/patchwork.plot.r 2013-01-04 10:48:00 UTC (rev 155)
+++ pkg/patchwork/R/patchwork.plot.r 2013-01-21 12:12:34 UTC (rev 156)
@@ -2,6 +2,11 @@
Reference=NULL,Alpha=0.0001,SD=1)
{
+ if(is.null(Normal.bam) & is.null(Reference))
+ {
+ stop("You must give either a normal bam or a reference file for patchwork to run properly! \n")
+ }
+
#Extract the name from tumor.bam
name = strsplit(tolower(Tumor.bam),".bam")
@@ -31,8 +36,8 @@
#If Normal.pileup exists, use with normal.vcf (if mpileup) or without (if pileup) to
#create normalalf in patchwork.alleledata()
- #Note that right not normalalf isnt even enabled.
- #if (!is.null(Normal.pileup)) normalalf <- patchwork.alleledata(Normal.pileup, vcf=Normal.vcf)
+ #Create normalalf
+ if (!is.null(Normal.pileup)) normalalf <- patchwork.alleledata(Normal.pileup, vcf=Normal.vcf)
alf = patchwork.alleledata(Tumor.pileup,
#normalalf=normalalf,
Modified: www/pw_exec.php
===================================================================
--- www/pw_exec.php 2013-01-04 10:48:00 UTC (rev 155)
+++ www/pw_exec.php 2013-01-21 12:12:34 UTC (rev 156)
@@ -328,7 +328,7 @@
An example run of patchwork.copynumbers(): <br />
<pre>
- patchwork.copynumbers(CNFile="path/to/prefix_copynumbers.Rdata",cn2=0.8, delta=0.28, het=0.21, hom=0.79)
+ patchwork.copynumbers(CNfile="path/to/prefix_copynumbers.Rdata",cn2=0.8, delta=0.28, het=0.21, hom=0.79)
</pre>
This will generate an additional set of plots in your working directory. <br /><br />
More information about the Patchwork-commits
mailing list