Mercurial > repos > drosofff > msp_sr_readmap_and_size_histograms

annotate readmap.xml @ 10:71a46afb9ce7 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 e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
author mvdbeek
date Tue, 02 Feb 2016 12:54:00 -0500
parents 44a0b0fa52a9
children edd57052f0bb
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
1 <tool id="Readmap" name="Generate readmap and histograms from alignment files" version="1.1.0">
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"/>
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
80 <param name="yrange" type="integer" size="6" value="0" label="y axis range for readmap tool" help="leave at 0 for autoscaling"/>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
81 <param name="ylabel" type="text" size="15" value="Number of reads" label="y axis label"/>
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>
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
87 <configfile name="plotCode"><![CDATA[
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
88 ## Setup R error handling to go to stderr
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
89 options( show.error.messages=F,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
90 error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
91 library(RColorBrewer)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
92 library(lattice)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
93 library(latticeExtra)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
94 library(grid)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
95 library(gridExtra)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
96
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
97 ## data frames implementation
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
98
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
99 rm=read.delim("${readmap_dataframe}", header=T, row.names=NULL)
10
71a46afb9ce7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 9
diff changeset
100 n_samples=length(unique(rm\$sample))
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
101 genes=unique(levels(rm$gene))
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
102 per_gene_readmap=lapply(genes, function(x) subset(rm, gene==x)) ####### ?
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
103 n_genes=length(per_gene_readmap)
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
104
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
105 size=read.delim("${size_distribution_dataframe}", header=T, row.names=NULL)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
106 per_gene_size=lapply(genes, function(x) subset(size, gene==x)) ###### ?
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
107
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
108 ## end of data frames implementation
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
109
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
110 ## functions
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
111
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
112 plot_readmap=function(df, ...) {
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
113 combineLimits(xyplot(count~coord|factor(sample, levels=unique(sample))+reorder(gene, count, function(x) -sum(abs(x))),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
114 data=df,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
115 type='h',
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
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)),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
117 xlab=NULL, main=NULL, ylab=NULL,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
118 as.table=T,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
119 origin = 0,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
120 horizontal=FALSE,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
121 group=polarity,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
122 col=c("red","blue"),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
123 par.strip.text = list(cex=0.7),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
124 ...))
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
125 }
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
126
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
127 plot_size_distribution= function(df, ...) {
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
128 smR.prepanel=function(x,y,...){; yscale=c(-max(abs(y)), max(abs(y)));list(ylim=yscale);}
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
129 bc= barchart(count~as.factor(size)|factor(sample, levels=unique(sample))+gene, data = df, origin = 0,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
130 horizontal=FALSE,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
131 group=polarity,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
132 stack=TRUE,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
133 col=c('red', 'blue'),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
134 cex=0.75,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
135 scales=list(y=list(tick.number=4, rot=90, relation="free", cex=0.7), x=list(cex=0.7) ),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
136 prepanel=smR.prepanel,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
137 xlab = NULL,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
138 ylab = NULL,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
139 main = NULL,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
140 as.table=TRUE,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
141 newpage = T,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
142 par.strip.text = list(cex=0.7),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
143 ...)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
144 combineLimits(bc)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
145 }
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
146
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
147 ## end of functions
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
148
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
149 ## function parameters'
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
150
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
151 par.settings.readmap=list(layout.heights=list(top.padding=0, bottom.padding=-2.5), strip.background = list(col=c("lightblue","lightgreen")) )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
152 par.settings.size=list(layout.heights=list(top.padding=-1, bottom.padding=-2.5), strip.background = list(col=c("lightblue","lightgreen")) )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
153 par.settings.combination.readmap=list(layout.heights=list(top.padding=0, bottom.padding=-3), strip.background=list(col=c("lightblue","lightgreen")) )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
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")) )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
155
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
156 ## end of function parameters'
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
157
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
158 ## GRAPHS
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
159
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
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 {
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
161 rows_per_page= 8; page_height_simple = 11.69; page_height_combi=11.69; extrarow=0 }
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
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 }
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
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 }
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
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
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
166
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
167 }
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
168
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
169 pdf(file="${readmap_PDF}", paper="special", height=page_height_simple, width=page_width)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
170 for (i in seq(1,n_genes,rows_per_page)) {
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
171 start=i
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
172 end=i+rows_per_page-1
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
173 if (end>n_genes) {end=n_genes}
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
174 if ("${yrange}" != 0) {
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
175 readmap_plot.list=lapply(per_gene_readmap[start:end], function(x) plot_readmap(x, par.settings=par.settings.readmap))
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
176 } else {
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
177 readmap_plot.list=lapply(per_gene_readmap[start:end], function(x) plot_readmap(x, ylim=c(-"{$yrange}", "{$yrange}"), par.settings=par.settings.readmap))
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
178 }
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
179 args.list=c(readmap_plot.list, list(nrow=rows_per_page, ncol=1,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
180 main=textGrob("Read Maps (nucleotide coordinates)", gp=gpar(cex=1), just="top"),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
181 left=textGrob("${ylabel}", gp=gpar(cex=1), vjust=1, rot=90)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
182 #sub=textGrob("readmap coordinates", gp=gpar(cex=.75), just="bottom")
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
183 )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
184 )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
185 do.call(grid.arrange, args.list)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
186 }
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
187 devname=dev.off()
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
188
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
189
9
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
190 pdf(file="${size_PDF}", paper="special", height=page_height_simple, width=page_width)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
191 for (i in seq(1,n_genes,rows_per_page)) {
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
192 start=i
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
193 end=i+rows_per_page-1
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
194 if (end>n_genes) {end=n_genes}
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
195 plot.list=lapply(per_gene_size[start:end], function(x) plot_size_distribution(x, par.settings=par.settings.size) )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
196 args.list=c(plot.list, list(nrow=rows_per_page, ncol=1,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
197 main=textGrob("Size distributions (in nucleotides)", gp=gpar(cex=1), just="top"),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
198 left=textGrob("${ylabel}", gp=gpar(cex=1), vjust=1, rot=90)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
199 #sub="readsize in nucleotides"
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
200 )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
201 )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
202 do.call(grid.arrange, args.list)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
203 }
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
204 devname=dev.off()
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
205
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
206 pdf(file="${combi_PDF}", paper="special", height=page_height_combi, width=page_width)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
207 for (i in seq(1,n_genes,rows_per_page/2)) {
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
208 start=i
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
209 end=i+rows_per_page/2-1
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
210 if (end>n_genes) {end=n_genes}
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
211 read_plot.list=lapply(per_gene_readmap[start:end], function(x) plot_readmap(x, par.settings=par.settings.combination.readmap))
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
212 size_plot.list=lapply(per_gene_size[start:end], function(x) plot_size_distribution(x, strip=FALSE, par.settings=par.settings.combination.size))
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
213 plot.list=rbind(read_plot.list, size_plot.list )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
214 args.list=c(plot.list, list(nrow=rows_per_page + extrarow, ncol=1,
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
215 main=textGrob("${title}", gp=gpar(cex=1), just="top"),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
216 left=textGrob("${ylabel}", gp=gpar(cex=1), vjust=1, rot=90),
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
217 sub=textGrob("${xlabel}", gp=gpar(cex=1), just="bottom")
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
218 )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
219 )
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
220 do.call(grid.arrange, args.list)
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
221 }
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
222 devname=dev.off()
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
223
44a0b0fa52a9 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
mvdbeek
parents: 8
diff changeset
224 ]]></configfile>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
225 </configfiles>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
226
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
227 <outputs>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
228 <data format="tabular" name="readmap_dataframe" label="Readmap dataframe"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
229 <data format="tabular" name="size_distribution_dataframe" label="Size distribution dataframe"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
230 <data format="pdf" name="readmap_PDF" label="Readmaps"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
231 <data format="pdf" name="size_PDF" label="Size distribution"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
232 <data format="pdf" name="combi_PDF" label="Size distribution and Readmaps"/>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
233 </outputs>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
234 <help>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
235
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
236 **What it does**
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
237
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
238 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
239 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
240 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
241
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
242
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
243 .. class:: warningmark
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
244
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
245 '''TIP''' The input data can be produced using the sRbowtie tool.
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
246
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
247 ----
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
248
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
249 '''Example'''
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
250
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
251 Query sequence::
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
252 For a SAM file as the following:
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
253
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
254 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
255
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
256 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
257
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
258 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
259
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
260 produce a plot like this:
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
261
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
262 ----
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
263
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
264 .. image:: static/images/readmap.png
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
265 :height: 800
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
266 :width: 500
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
267
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
268 </help>
1
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
269 <tests>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
270 <test>
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
271 <param name="genomeSource" value="history" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
272 <param name="ownFile" value ="transposons.fasta" ftype="fasta" />
1
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
273 <param name="series_0|input" value="sample1.srbowtie_out" ftype="tabular"/>
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
274 <param name="series_0|norm" value="1" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
275 <param name="series_1|input" value="sample2.srbowtie_out" ftype="tabular"/>
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
276 <param name="series_1|norm" value="1" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
277 <param name="series_2|input" value="sample3.srbowtie_out" ftype="tabular"/>
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
278 <param name="series_2|norm" value="1" />
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
279 <param name="minquery" value="20" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
280 <param name="maxquery" value="30" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
281 <param name="title" value="Readmaps and size distributions" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
282 <param name="xlabel" value="Coordinates/read size" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
283 <param name="ylabel" value="Number of reads" />
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
284 <param name="rows_per_page" value="8" />
1
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
285 <output name="readmap_dataframe" ftype="tabular" file="Readmap_dataframe.tab" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
286 <output name="size_distribution_dataframe" ftype="tabular" file="Size_distribution_dataframe.tab" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
287 <output name="readmap_PDF" ftype="pdf" file="Readmaps.pdf" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
288 <output name="size_PDF" ftype="pdf" file="Size_distribution.pdf" />
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
289 <output name="combi_PDF" ftype="pdf" file="Size_distribution_and_Readmaps.pdf" />
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
290 </test>
1
eee9701a7491 Uploaded
mvdbeek
parents: 0
diff changeset
291 </tests>
0
9af9983dcd02 Imported from capsule None
drosofff
parents:
diff changeset
292 </tool>