view minfi_qc.xml @ 76:fef50103f6b3 draft

planemo upload for repository https://github.com/kpbioteam/ewas_galaxy commit 945cecdba6528d35c6dfa866f785f182c5dd22ed
author kpbioteam
date Fri, 07 Jun 2019 08:57:26 -0400
parents 9c6fbb7d5a2a
children
line wrap: on
line source

<tool id="minfi_qc" name="Minfi QC" version="@MINFI_VERSION@">
    <description>provides a simple quality control matrix and plot</description>
    <macros>
        <import>macros.xml</import>
    </macros>
     <expand macro="requirements" />
    <command detect_errors="exit_code">
    <![CDATA[
        Rscript '$read_qc_script'
     ]]>
    </command>
    <configfiles>
        <configfile name="read_qc_script">
    <![CDATA[
require("minfi", quietly = TRUE)

MSet <- get(load('$MethylSet'))

qc <- getQC(MSet)

write.table(qc, '$qctab')

png('$qcpng')
plotQC(qc)
dev.off()
    ]]>
        </configfile>
    </configfiles>
    <inputs>
        <param type="data" name="MethylSet" format="rdata" label="MethylSet"
            help="This class holds preprocessed data for Illumina methylation microarrays, mapped to a genomic location."/>
    </inputs>
    <outputs>
        <data name="qctab" format="txt" label="Quality Control Report"/>
        <data name="qcpng" format="png" label="Quality Control Plot"/>
    </outputs>
    <tests>
        <test>
            <param name="MethylSet" value="MethylSet.rdata"/>
            <output name="qctab" file="Quality_Control_Report.txt"/>
            <output name="qcpng" file="Quality_Control_Plot.png" compare="sim_size"/>
        </test>
    </tests>
    <help><![CDATA[
Quality Control (QC) outputs plot of the log median intensity in both the methylated (M) and unmethylated (U) channels. When plotting these two medians against each other the good samples cluster together, while failed samples tend to separate and have lower median intensities.
    ]]></help>
    <expand macro="citations" />
</tool>