Mercurial > repos > drosofff > msp_sr_readmap_and_size_histograms

annotate readmap.xml @ 18:893560ece89f draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit 9237338d798251fb2667280d597746e852f3ffcc-dirty
author mvdbeek
date Wed, 03 Feb 2016 11:34:35 -0500
parents 70f4385534f9
children 99f478136171
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
18
893560ece89f planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit 9237338d798251fb2667280d597746e852f3ffcc-dirty
mvdbeek
parents: 16
diff changeset
1 <tool id="Readmap" name="Generate readmap and histograms from alignment files" version="1.1.3">
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
2 <description>from sRbowtie aligment</description>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
3 <requirements>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
4 <requirement type="package" version="0.12.7">bowtie</requirement>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
5 <requirement type="package" version="0.7.7">pysam</requirement>
6
0f7cf7fb3465 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit fe40dec87779c1fcfbd03330e653aa886f4a2cda
drosofff
parents: 5
diff changeset
6 <requirement type="package" version="3.1.2">R</requirement>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
7 <requirement type="package" version="2.14">biocbasics</requirement>
1
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
8 <requirement type="package" version="1.9">numpy</requirement>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
9 </requirements>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
10 <command interpreter="python">
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
11 readmap.py
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
12 #if $refGenomeSource.genomeSource == "history":
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
13 --reference_fasta ## sys.argv[2]
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
14 $refGenomeSource.ownFile ## index source
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
15 #else:
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
16 #silent reference= filter( lambda x: str( x[0] ) == str( $refGenomeSource.series[0].input.dbkey ), $__app__.tool_data_tables[ 'bowtie_indexes' ].get_fields() )[0][-1]
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
17 --reference_bowtie_index
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
18 $reference
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
19 #end if
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
20 --rcode
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
21 $plotCode
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
22 --output_readmap
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
23 $readmap_dataframe
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
24 --output_size_distribution
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
25 $size_distribution_dataframe
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
26 --minquery
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
27 $minquery
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
28 --maxquery
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
29 $maxquery
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
30 --input
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
31 #for $i in $refGenomeSource.series
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
32 $i.input
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
33 #end for
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
34 --ext
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
35 #for $i in $refGenomeSource.series
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
36 $i.input.ext
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
37 #end for
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
38 --label
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
39 #for $i in $refGenomeSource.series
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
40 "$i.input.name"
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
41 #end for
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
42 --normalization_factor
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
43 #for $i in $refGenomeSource.series
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
44 $i.norm
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
45 #end for
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
46 #if $gff:
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
47 --gff
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
48 $gff
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
49 #end if
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
50
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
51 </command>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
52 <inputs>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
53 <conditional name="refGenomeSource">
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
54 <param name="genomeSource" type="select" label="Will you select a reference genome from your history or use a built-in index?" help="Built-ins were indexed using default options">
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
55 <option value="indexed">Use a built-in index</option>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
56 <option value="history">Use one from the history</option>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
57 </param>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
58 <when value="indexed">
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
59 <repeat name="series" title="Add alignment files">
4
dffa22efc6a8 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 3
diff changeset
60 <param name="input" type="data" label="Select multiple alignments to parse" format="tabular,sam,bam">
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
61 <validator type="dataset_metadata_in_data_table" table_name="bowtie_indexes" metadata_name="dbkey" metadata_column="0" message="database not set for this bowtie output. Select the database(=genome used for matching) manually, or select a reference fasta from your history."/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
62 </param>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
63 <param name="norm" type="float" value="1" label="Indicate a normalization factor to compare multiple aligments"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
64 </repeat>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
65 </when>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
66 <when value="history">
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
67 <param name="ownFile" type="data" format="fasta" label="Select a fasta file, that served as the reference index for the alignments" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
68 <repeat name="series" title="Add alignment files">
4
dffa22efc6a8 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 3
diff changeset
69 <param name="input" type="data" label="Select multiple alignments to parse" format="tabular,sam,bam"/>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
70 <param name="norm" type="float" value="1" label="Indicate a normalization factor to compare multiple aligments"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
71 </repeat>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
72 </when>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
73 </conditional>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
74 <param name="gff" type="data" format="gff3" optional="true" label="Optional: select a GFF to investigate regions of interest" help="GFF must match genome build"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
75 <!-- <validator type="dataset_metadata_in_data_table" table_name="bowtie_indexes" metadata_name="dbkey" metadata_column="0" message="GFF database and alignment file databse do not match!"/> -->
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
76 <param name="minquery" type="integer" size="3" value="18" label="Min size of query small RNAs" help="'18' = 18 nucleotides"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
77 <param name="maxquery" type="integer" size="3" value="28" label="Max size of query small RNAs" help="'28' = 28 nucleotides"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
78 <param name="title" type="text" size="15" value= "Readmaps and size distributions" label="Main Titles"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
79 <param name="xlabel" type="text" size="15" value="Coordinates/read size" label="x axis label"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
80 <param name="ylabel" type="text" size="15" value="Number of reads" label="y axis label"/>
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
81 <param name="yrange" type="integer" size="3" value="0" label="y axis range for readmaps. 0 means auto-scaling."/>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
82 <param name="rows_per_page" type="text" size="9" value="8" label="How many items to display per page?">
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
83 <validator type="in_range" min="6" max="20" message="Select between 6 and 20 rows, as the readability will suffer otherwise."/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
84 </param>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
85 </inputs>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
86 <configfiles>
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
87 <configfile name="plotCode">
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
88 ## Setup R error handling to go to stderr
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
89 options( show.error.messages=F,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
90 error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
91 library(RColorBrewer)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
92 library(lattice)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
93 library(latticeExtra)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
94 library(grid)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
95 library(gridExtra)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
96
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
97 ## data frames implementation
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
98
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
99 rm=read.delim("${readmap_dataframe}", header=T, row.names=NULL)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
100 n_samples=length(unique(rm\$sample))
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
101 genes=unique(levels(rm\$gene))
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
102 per_gene_readmap=lapply(genes, function(x) subset(rm, gene==x)) ####### ?
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
103 n_genes=length(per_gene_readmap)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
104
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
105 size=read.delim("${size_distribution_dataframe}", header=T, row.names=NULL)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
106 per_gene_size=lapply(genes, function(x) subset(size, gene==x)) ###### ?
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
107
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
108 ## end of data frames implementation
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
109
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
110 ## functions
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
111
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
112 plot_readmap=function(df, ...) {
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
113 combineLimits(xyplot(count~coord|factor(sample, levels=unique(sample))+reorder(gene, count, function(x) -sum(abs(x))),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
114 data=df,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
115 type='h',
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
116 scales= list(relation="free", x=list(rot=0, cex=0.7, axs="i", tck=0.5), y=list(tick.number=4, rot=90, cex=0.7)),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
117 xlab=NULL, main=NULL, ylab=NULL,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
118 as.table=T,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
119 origin = 0,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
120 horizontal=FALSE,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
121 group=polarity,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
122 col=c("red","blue"),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
123 par.strip.text = list(cex=0.7),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
124 ...))
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
125 }
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
126
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
127 plot_size_distribution= function(df, ...) {
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
128 smR.prepanel=function(x,y,...){; yscale=c(-max(abs(y)), max(abs(y)));list(ylim=yscale);}
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
129 bc= barchart(count~as.factor(size)|factor(sample, levels=unique(sample))+gene, data = df, origin = 0,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
130 horizontal=FALSE,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
131 group=polarity,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
132 stack=TRUE,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
133 col=c('red', 'blue'),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
134 cex=0.75,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
135 scales=list(y=list(tick.number=4, rot=90, relation="free", cex=0.7), x=list(cex=0.7) ),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
136 prepanel=smR.prepanel,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
137 xlab = NULL,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
138 ylab = NULL,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
139 main = NULL,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
140 as.table=TRUE,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
141 newpage = T,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
142 par.strip.text = list(cex=0.7),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
143 ...)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
144 combineLimits(bc)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
145 }
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
146
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
147 ## end of functions
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
148
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
149 ## function parameters'
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
150
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
151 par.settings.readmap=list(layout.heights=list(top.padding=0, bottom.padding=-2.5), strip.background = list(col=c("lightblue","lightgreen")) )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
152 par.settings.size=list(layout.heights=list(top.padding=-1, bottom.padding=-2.5), strip.background = list(col=c("lightblue","lightgreen")) )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
153 par.settings.combination.readmap=list(layout.heights=list(top.padding=0, bottom.padding=-3), strip.background=list(col=c("lightblue","lightgreen")) )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
154 par.settings.combination.size=list(layout.heights=list(top.padding=-2, bottom.padding=-0.5), strip.background=list(col=c("lightblue", "lightgreen")) )
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
155
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
156 ## end of function parameters'
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
157
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
158 ## GRAPHS
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
159
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
160 if (n_genes > 7) {page_height_simple = 11.69; page_height_combi=11.69; rows_per_page=${rows_per_page}; extrarow=0 } else {
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
161 rows_per_page= 8; page_height_simple = 11.69; page_height_combi=11.69; extrarow=0 }
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
162 ## rows_per_page= 8; page_height_simple = 11.69/7*n_genes; page_height_combi=11.69/9*(n_genes*2); extrarow=0 }
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
163 ## rows_per_page= n_genes; page_height_simple = 11.69/n_genes/4; page_height_combi=11.69/(n_genes*2); extrarow=1 }
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
164 if (n_samples > 4) {page_width = 8.2677*n_samples/4} else {page_width = 8.2677*n_samples/3} # to test
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
165
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
166 pdf(file="${readmap_PDF}", paper="special", height=page_height_simple, width=page_width)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
167 for (i in seq(1,n_genes,rows_per_page)) {
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
168 start=i
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
169 end=i+rows_per_page-1
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
170 if (end>n_genes) {end=n_genes}
15
16d0811f7a09 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 14
diff changeset
171 if (${yrange} == 0) { readmap_plot.list=lapply(per_gene_readmap[start:end], function(x) plot_readmap(x, par.settings=par.settings.readmap)) } else {
16d0811f7a09 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 14
diff changeset
172 readmap_plot.list=lapply(per_gene_readmap[start:end], function(x) plot_readmap(x, ylim=c(-${yrange}, ${yrange}) , par.settings=par.settings.readmap)) }
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
173 args.list=c(readmap_plot.list, list(nrow=rows_per_page, ncol=1,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
174 main=textGrob("Read Maps (nucleotide coordinates)", gp=gpar(cex=1), just="top"),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
175 left=textGrob("${ylabel}", gp=gpar(cex=1), vjust=1, rot=90)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
176 #sub=textGrob("readmap coordinates", gp=gpar(cex=.75), just="bottom")
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
177 )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
178 )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
179 do.call(grid.arrange, args.list)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
180 }
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
181 devname=dev.off()
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
182
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
183
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
184 pdf(file="${size_PDF}", paper="special", height=page_height_simple, width=page_width)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
185 for (i in seq(1,n_genes,rows_per_page)) {
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
186 start=i
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
187 end=i+rows_per_page-1
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
188 if (end>n_genes) {end=n_genes}
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
189 plot.list=lapply(per_gene_size[start:end], function(x) plot_size_distribution(x, par.settings=par.settings.size) )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
190 args.list=c(plot.list, list(nrow=rows_per_page, ncol=1,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
191 main=textGrob("Size distributions (in nucleotides)", gp=gpar(cex=1), just="top"),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
192 left=textGrob("${ylabel}", gp=gpar(cex=1), vjust=1, rot=90)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
193 #sub="readsize in nucleotides"
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
194 )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
195 )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
196 do.call(grid.arrange, args.list)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
197 }
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
198 devname=dev.off()
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
199
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
200 pdf(file="${combi_PDF}", paper="special", height=page_height_combi, width=page_width)
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
201 for (i in seq(1,n_genes,rows_per_page/2)) {
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
202 start=i
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
203 end=i+rows_per_page/2-1
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
204 if (end>n_genes) {end=n_genes}
15
16d0811f7a09 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 14
diff changeset
205 if (${yrange} == 0) {readmap_plot.list=lapply(per_gene_readmap[start:end], function(x) plot_readmap(x, par.settings=par.settings.readmap)) } else {
13
355940295e76 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 12
diff changeset
206 readmap_plot.list=lapply(per_gene_readmap[start:end], function(x) plot_readmap(x, ylim=c(-${yrange}, ${yrange}), par.settings=par.settings.readmap)) }
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
207 size_plot.list=lapply(per_gene_size[start:end], function(x) plot_size_distribution(x, strip=FALSE, par.settings=par.settings.combination.size))
14
f3e338e433ef planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 13
diff changeset
208 plot.list=rbind(readmap_plot.list, size_plot.list )
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
209 args.list=c(plot.list, list(nrow=rows_per_page + extrarow, ncol=1,
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
210 main=textGrob("${title}", gp=gpar(cex=1), just="top"),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
211 left=textGrob("${ylabel}", gp=gpar(cex=1), vjust=1, rot=90),
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
212 sub=textGrob("${xlabel}", gp=gpar(cex=1), just="bottom")
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
213 )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
214 )
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
215 do.call(grid.arrange, args.list)
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
216 }
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
217 devname=dev.off()
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
218
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
219
12
2fc0d4756048 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 11
diff changeset
220 </configfile>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
221 </configfiles>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
222
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
223 <outputs>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
224 <data format="tabular" name="readmap_dataframe" label="Readmap dataframe"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
225 <data format="tabular" name="size_distribution_dataframe" label="Size distribution dataframe"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
226 <data format="pdf" name="readmap_PDF" label="Readmaps"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
227 <data format="pdf" name="size_PDF" label="Size distribution"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
228 <data format="pdf" name="combi_PDF" label="Size distribution and Readmaps"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
229 </outputs>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
230 <help>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
231
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
232 **What it does**
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
233
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
234 Takes one or more alignment files (BAM, SAM or tabular bowtie output) as input and produces a "Readmap",
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
235 where by default for each "chromosome" the position of the read is recorded on the x-axis, and the y-axis indicates
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
236 the number of reads per position. Reads that map in sense are on the top, reads that map antisense are on the bottom.
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
237
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
238
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
239 .. class:: warningmark
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
240
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
241 '''TIP''' The input data can be produced using the sRbowtie tool.
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
242
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
243 ----
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
244
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
245 '''Example'''
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
246
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
247 Query sequence::
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
248 For a SAM file as the following:
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
249
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
250 5 16 2L_79 24393 255 17M * 0 0 CCTTCATCTTTTTTTTT IIIIIIIIIIIIIIIII XA:i:0 MD:Z:17 NM:i:0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
251
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
252 11 0 2R_1 12675 255 21M * 0 0 AAAAAAAACGCGTCCTTGTGC IIIIIIIIIIIIIIIIIIIII XA:i:0 MD:Z:21 NM:i:0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
253
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
254 2 16 2L_5 669 255 23M * 0 0 TGTTGCTGCATTTCTTTTTTTTT IIIIIIIIIIIIIIIIIIIIIII XA:i:0 MD:Z:23 NM:i:0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
255
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
256 produce a plot like this:
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
257
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
258 ----
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
259
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
260 .. image:: static/images/readmap.png
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
261 :height: 800
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
262 :width: 500
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
263
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
264 </help>
1
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
265 <tests>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
266 <test>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
267 <param name="genomeSource" value="history" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
268 <param name="ownFile" value ="transposons.fasta" ftype="fasta" />
1
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
269 <param name="series_0|input" value="sample1.srbowtie_out" ftype="tabular"/>
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
270 <param name="series_0|norm" value="1" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
271 <param name="series_1|input" value="sample2.srbowtie_out" ftype="tabular"/>
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
272 <param name="series_1|norm" value="1" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
273 <param name="series_2|input" value="sample3.srbowtie_out" ftype="tabular"/>
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
274 <param name="series_2|norm" value="1" />
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
275 <param name="minquery" value="20" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
276 <param name="maxquery" value="30" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
277 <param name="title" value="Readmaps and size distributions" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
278 <param name="xlabel" value="Coordinates/read size" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
279 <param name="ylabel" value="Number of reads" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
280 <param name="rows_per_page" value="8" />
1
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
281 <output name="readmap_dataframe" ftype="tabular" file="Readmap_dataframe.tab" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
282 <output name="size_distribution_dataframe" ftype="tabular" file="Size_distribution_dataframe.tab" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
283 <output name="readmap_PDF" ftype="pdf" file="Readmaps.pdf" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
284 <output name="size_PDF" ftype="pdf" file="Size_distribution.pdf" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
285 <output name="combi_PDF" ftype="pdf" file="Size_distribution_and_Readmaps.pdf" />
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
286 </test>
1
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
287 </tests>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
288 </tool>