| Previous changeset 11:8d3ca9652a5b (2021-10-23) Next changeset 13:124f404b0fe7 (2023-02-25) |
|
Commit message:
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 6806c0677e53d52164707faeb36947987f5c500a |
|
modified:
overlapping_reads.xml signature.r signature.xml test-data/global.pdf test-data/lattice.pdf |
| b |
| diff -r 8d3ca9652a5b -r aa5e2c64dff8 overlapping_reads.xml --- a/overlapping_reads.xml Sat Oct 23 22:55:19 2021 +0000 +++ b/overlapping_reads.xml Sat Oct 22 23:49:52 2022 +0000 |
| [ |
| @@ -1,16 +1,16 @@ -<tool id="overlapping_reads" name="Get overlapping reads" version="3.3"> +<tool id="overlapping_reads" name="Get overlapping reads" version="3.4.1"> <description /> <requirements> - <requirement type="package" version="0.16.0.1=py37h45aed0b_3">pysam</requirement> - <requirement type="package" version="1.13=h8c37831_0">samtools</requirement> + <requirement type="package" version="0.18.0">pysam</requirement> </requirements> <stdio> <exit_code range="1:" level="fatal" description="Tool exception" /> </stdio> <command detect_errors="exit_code"><![CDATA[ - samtools index '$input' && + ln -f -s $input.metadata.bam_index input.bam.bai && + ln -s $input input.bam && python '$__tool_directory__'/overlapping_reads.py - --input '$input' + --input input.bam --minquery '$minquery' --maxquery '$maxquery' --mintarget '$mintarget' |
| b |
| diff -r 8d3ca9652a5b -r aa5e2c64dff8 signature.r --- a/signature.r Sat Oct 23 22:55:19 2021 +0000 +++ b/signature.r Sat Oct 22 23:49:52 2022 +0000 |
| [ |
| @@ -1,8 +1,8 @@ ## Setup R error handling to go to stderr -options(show.error.messages = F, +options(show.error.messages = FALSE, error = function() { cat(geterrmessage(), file = stderr()) - q("no", 1, F) + q("no", 1, FALSE) } ) warnings() @@ -26,12 +26,12 @@ ) parser <- OptionParser(usage = "%prog [options] file", option_list = option_list) args <- parse_args(parser) - + # data frames implementation -h_dataframe <- read.delim(args$h_dataframe, header = F) +h_dataframe <- read.delim(args$h_dataframe, header = FALSE) colnames(h_dataframe) <- c("chrom", "overlap", "sig", "z-score") h_dataframe$sig <- h_dataframe$sig * 100 # to get probs in % -z_dataframe <- read.delim(args$z_dataframe, header = F) +z_dataframe <- read.delim(args$z_dataframe, header = FALSE) colnames(z_dataframe) <- c("chrom", "overlap", "sig", "z-score") # functions @@ -77,7 +77,7 @@ pch = 19, lwd = 2, as.table = TRUE, layout = c(2, 12), - newpage = T, + newpage = TRUE, ...) plot(p) dev.off() |
| b |
| diff -r 8d3ca9652a5b -r aa5e2c64dff8 signature.xml --- a/signature.xml Sat Oct 23 22:55:19 2021 +0000 +++ b/signature.xml Sat Oct 22 23:49:52 2022 +0000 |
| [ |
| @@ -1,20 +1,20 @@ -<tool id="signature" name="Small RNA Signatures" version="3.3>"> +<tool id="signature" name="Small RNA Signatures" version="3.4.1>"> <description /> <requirements> - <requirement type="package" version="1.13=h8c37831_0">samtools</requirement> - <requirement type="package" version="1.21.3=py37h31617e3_0">numpy</requirement> - <requirement type="package" version="0.16.0.1=py37h45aed0b_3">pysam</requirement> - <requirement type="package" version="1.7.1=r40hc72bb7e_0">r-optparse</requirement> - <requirement type="package" version="0.6_29=r40hc72bb7e_1">r-latticeextra</requirement> - <requirement type="package" version="2.3=r40hc72bb7e_1003">r-gridextra</requirement> + <requirement type="package" version="1.23.4">numpy</requirement> + <requirement type="package" version="0.18.0">pysam</requirement> + <requirement type="package" version="1.7.3">r-optparse</requirement> + <requirement type="package" version="0.6_30">r-latticeextra</requirement> + <requirement type="package" version="2.3">r-gridextra</requirement> </requirements> <stdio> <exit_code range="1:" level="fatal" description="Tool exception" /> </stdio> <command detect_errors="exit_code"><![CDATA[ - samtools index '$input' && + ln -f -s $input.metadata.bam_index input.bam.bai && + ln -s $input input.bam && python '$__tool_directory__'/signature.py - --input '$input' + --input input.bam --minquery '$minquery' --maxquery '$maxquery' --mintarget '$mintarget' @@ -91,7 +91,7 @@ Compute small RNA (piRNA, siRNA, ...) signatures. This tool computes (i) the number of pairs **aligned** reads by overlap classes (in nt) and associated z-scores, -and (ii) the ping-pong signal (`Brennecke et al, 2009`_) and associated z-scores. +and (ii) the ping-pong signal (Brennecke et al. Science. 2008;322: 1387–1392) and associated z-scores. **Note** that the number of pairs of aligned reads is disctint from the number of pairs of reads when these reads can be aligned at multiple positions in the genome. The two values are equal only @@ -102,7 +102,6 @@ genome-wide level or by item (chromosomes, genes, etc.). For details on computational algorithmes for piRNA and siRNA signatures, see `Antoniewski (2014)`_. -.. _Brennecke et al, 2009: http://dx.doi.org/10.1126/science.1165171 .. _Antoniewski (2014): https://link.springer.com/protocol/10.1007%2F978-1-4939-0931-5_12 **Input** |
| b |
| diff -r 8d3ca9652a5b -r aa5e2c64dff8 test-data/global.pdf |
| b |
| Binary file test-data/global.pdf has changed |
| b |
| diff -r 8d3ca9652a5b -r aa5e2c64dff8 test-data/lattice.pdf |
| b |
| Binary file test-data/lattice.pdf has changed |