# HG changeset patch # User stef # Date 1408547786 14400 # Node ID 19e91eb3d06bcd2b49e581823979d743b74a4502 # Parent bd81b2559c109a9d6b237f470b353ca3ac092cd4 Uploaded diff -r bd81b2559c10 -r 19e91eb3d06b QDNAseq.xml --- a/QDNAseq.xml Thu Aug 07 11:25:23 2014 -0400 +++ b/QDNAseq.xml Wed Aug 20 11:16:26 2014 -0400 @@ -1,37 +1,44 @@ - + - - + + + samtools + TOOL_DIR Quantitative copy number abberation detection - - + + QDNAseq.R $qdnaseq_cfg - + + + + .*\S - - + + + + @@ -39,11 +46,7 @@ - - - - - + @@ -55,65 +58,45 @@ - + - - - - + + + + + + + + + + + + - - - - - - - - - - - - - - - - - + + + + + - - - - - - - + - - - - - - - - - - - - - - - - + + + + + + @@ -128,109 +111,105 @@ - - - - + + + + + + + + + + + + + + + + + + + + - + + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + - + -## this file was sourced in QDNAseq R wrapper script -## in this way all galaxy params are passes to R +## Desc: this file was sourced in QDNAseq R wrapper script +## as means to pass all galaxy params to R +## ----- ## required params -as.integer( "${binSize}" ) -> binSize +## ----- +"${binSizes}" -> binSizesString "${experimentType}" -> experimentType "${jobName}" -> outputName -## extra params +## ----- +## extra main params +## ----- "${htmlFile}" -> outputHtml +"${htmlFile.id}" -> outputId +"${__new_file_path__}" -> newFilePath + "${htmlFile.files_path}" -> outputPath -"${use_own_binannotation.binannotation_file}" -> binAnnotations -as.logical( "${doCall}" ) -> doCall +as.logical( "${doSegment}" ) -> doSegment +as.logical( "${debug}" ) -> debug +## ----- +## own bin-annotations file options +## ----- +"${binannotations_source.binannotation_file}" -> binAnnotations + +## ----- ## advanced options +## ----- as.double( "${advanced.undoSD}" ) -> undoSD -as.logical( "${advanced.debug}" ) -> debug as.logical( "${advanced.blacklist}" ) -> filterBlacklistedBins as.integer( "${advanced.mappability}" ) -> mappabilityCutoff "${advanced.undo_splits}" -> undoSplits +as.logical( "${advanced.copynumbers_igv}" ) -> doOutputCopynumbersIgv -## history output params -as.logical( "${extra_history_outputs.readcounts_rds}" ) -> doOutputReadcountsRds -as.logical( "${extra_history_outputs.copynumbers_rds}" ) -> doOutputCopynumbersRds -as.logical( "${extra_history_outputs.calls_rds}" ) -> doOutputCallsRds -as.logical( "${extra_history_outputs.copynumbers_igv}" ) -> doOutputCopynumbersIgv -"${rdsReadCounts}" -> readCountsDatasetFile -"${rdsCopyNumbers}" -> copyNumbersDatasetFile -"${rdsCalls}" -> calledSegmentsDatasetFile -"${igvCopyNumbers}" -> copyNumbersIgvDatasetFile +## #for binSize in $binSizes}.split(",")# +## "${binSize}kbp_${igvCopyNumbers}" -> copyNumbersIgvDatasetFile +## #end for -## plotting params -as.integer( "${plot_options.plot_width}" ) -> PLOT_WIDTH -as.integer( "${plot_options.plot_height}" ) -> PLOT_HEIGHT -"${plot_options.exclude_chrs}" -> excludeChrsString - +## ----- +## plot options +## ----- +as.integer( "${advanced.plot_width}" ) -> PLOT_WIDTH +as.integer( "${advanced.plot_height}" ) -> PLOT_HEIGHT +"${advanced.exclude_chrs}" -> excludeChrsString + +## ----- ## input BAMs init +## ----- c() -> bamsPaths c() -> bamsNames @@ -243,27 +222,13 @@ - + - - - - extra_history_outputs['readcounts_rds'] == "TRUE" - - - - extra_history_outputs['copynumbers_rds'] == "TRUE" - - - - extra_history_outputs['calls_rds'] == "TRUE" - - - - extra_history_outputs['copynumbers_igv'] == "TRUE" - + + + @@ -295,15 +260,17 @@ .. _here: http://www.stefs.nl/wp/contact +You can **test this tool** with built-in data by selecting the option "Run with test data" and press execute. + ----- .. class:: warningmark -As there is no R 3.1.0 package yet (which is required), the **dependencies** need to be installed by hand and available to the user under which galaxy runs: R (3.1.0) and bioconductor package QDNAseq (>= 1.0.5) +As there is no R 3.1.0 galaxy-package yet (a requirement for QDNAseq), the **dependencies** need to be installed by hand and available to the user under which galaxy runs: R (3.1.0) and bioconductor package QDNAseq (>= 1.0.5). In case the default R is not 3.1.0, also the wrapper xml must be updated to include the correct path during installation of this tool. .. class:: warningmark -The input BAMs are expected to be **single end reads of 50bp length** mapped to GRCh37/hg19 genome build. Other experiment setups are currently not supported within galaxy. See the documentation of QDNAseq at bioconductor on how to deal with different setups. +The input BAMs are expected to be **single end reads of 50bp length** mapped to GRCh37/hg19 genome build. Other experiment setups are currently not tested or supported within galaxy. See the documentation of QDNAseq at bioconductor on how to deal with different setups (or keep fingers crossed ;) ) .. class:: warningmark @@ -313,13 +280,12 @@ **Citation** -For the underlying tool please cite: llari Scheinin, Daoud Sie et al. DNA copy number analysis of fresh and formalin-fixed specimens by whole-genome sequencing: improved correction of systematic biases and exclusion of problematic regions, (submitted). See also the bioconductor package_ documentation. +For the underlying QDNAseq R package please cite: llari Scheinin, Daoud Sie et al. DNA copy number analysis of fresh and formalin-fixed specimens by whole-genome sequencing: improved correction of systematic biases and exclusion of problematic regions, (submitted). See also the bioconductor package_ documentation. .. _package: http://www.bioconductor.org/packages/release/bioc/html/QDNAseq.html +.. image:: ./LGG150_copynumber_example.png :height: 240 :width: 480 -.. image:: LGG150_copynumber_example.png -.. image:: LGG150_calls_example.png