Mercurial > repos > jcb-mpl > features_check_references
view features_check_references.xml @ 1:c1008b92a6c8 draft
Uploaded
author | jcb-mpl |
---|---|
date | Thu, 12 May 2022 12:15:10 +0000 |
parents | |
children |
line wrap: on
line source
<tool id="2021_05_features_check_references" name="Features -> checks references" version="0.0.1"> <description> </description> <requirements> <requirement type="package" >scilab</requirement> </requirements> <stdio> <exit_code range="1:" level="fatal" /> </stdio> <command> <![CDATA[ $__root_dir__/packages/scilab-*/bin/scilab-cli -nb -quit -f ${script_file} ]]> </command> <configfiles> <configfile name="script_file"> <![CDATA[ exec("$__tool_directory__/src/mz_9_compares_peaks.sci",-1); ... lasterror(); ... ... //load "${peaks_list}"; ... // 22avril21 peaks_list=glx_tab2div("${peaks_list}"); ... ... load "${peaks_detail}"; ... x836=glx_tab2div("${peaks_reference}"); ... diff_mz=${diff_mz}; ... k=${k}; ... x_peaks_final=mz_9_compares_peaks(peaks_list,peaks_detail,x836,diff_mz,k); ... ... div2tab(x_peaks_final,"${x_peaks_final}"); ... if ~isempty(lasterror(%f)); ... write(0,lasterror()); ... end; ]]> </configfile> </configfiles> <inputs> <param name="peaks_list" format="tabular" type="data" label="Features list" help="from function EIC->peaks:features" /> <param name="peaks_detail" format="mat" type="data" label="Detailed features" help="from function EIC->peaks:features" /> <param name="peaks_reference" format="tabular" type="data" label="Reference features" help="2columns: mz and RT, with column and line headers" /> <param name="diff_mz" value="0.0050" type="float" label="Threshold for m/z values" /> <param name="k" value="1" type="integer" label="Number of extracted features for each reference feature" help="the closest to the reference peak" /> </inputs> <outputs> <data name="x_peaks_final" format="tabular" label="Identified features" /> </outputs> <tests> <test> <param name="peaks_list" value="Peaks_list.tabular"/> <param name="peaks_detail" value="Peaks_details.mat"/> <param name="peaks_reference" value="ref_mz_rt.tab"/> <param name="diff_mz" value="0.0050"/> <param name="k" value="1"/> <output name="x_peaks_final"> <assert_contents> <has_text text="136.0616"/> <has_text text="229.135"/> </assert_contents> </output> </test> </tests> <help> **Author** Jean-Claude Boulet (INRA). --------------------------------------------------- ================================== FEATURES EXPERIMENTAL VS REFERENCE ================================== ----------- Description ----------- This function checks reference features among experimental features. ----------- Input files ----------- **Features list** The tabular file yielded by the function: EICs -> peaks:features **Detailed features** The HDF5-scilab file yielded by the function: EICs -> peaks:features **Reference features** A tabular file containing 2 columns: m/z values then retention times Row labels are the names of the reference features, e.g. chemical compounds ---------- Parameters ---------- **Threshold for m/z values** Observed features are associated to the reference features if 1) the difference in m/z values is under the threshold; 2) the reference RT falls within the range of the observed RTs. Default: 0.0050 **Number of extracted features for each reference feature** e.g. if the number is 2, the algorithm will extract the 2 experimental features closest to each reference feature. ------ Output ------ **Identified features** A tabular file with the columns: - name of the reference feature - number of the reference feature in the input file: Reference features - number of the experimental feature in the input files: Features List and Detailed features - number of identifications= the number of reference features each experimental feature has been attributed to; target = 1 - m/z value of the reference feature - m/z value of the experimental feature - difference of m/z values - RT of the reference feature - RT of the experimental feature - difference of RT - signal of the experimental feature </help> <citations> </citations> </tool>