Mercurial > repos > ethevenot > profia
changeset 2:e29c563df582 draft
planemo upload for repository https://github.com/workflow4metabolomics/profia.git commit f9c59d146d83ac980eca96215d9412cf65b760a0
author | ethevenot |
---|---|
date | Thu, 29 Jun 2017 09:21:32 -0400 |
parents | bc95bcb6ead0 |
children | 5215be2bdc9e |
files | README.md profia_config.xml profia_wrapper.R runit/output/figure.pdf runit/output/information.txt runit/output/sampleMetadata.tsv runit/output/variableMetadata.tsv runit/profia_runtests.R runit/profia_tests.R |
diffstat | 9 files changed, 702 insertions(+), 561 deletions(-) [+] |
line wrap: on
line diff
--- a/README.md Wed May 03 10:39:00 2017 -0400 +++ b/README.md Thu Jun 29 09:21:32 2017 -0400 @@ -7,8 +7,8 @@ ### Description -**Version:** 3.0.4 -**Date:** 2017-05-02 +**Version:** 3.0.6 +**Date:** 2017-06-29 **Author:** Alexis Delabriere and Etienne A. Thevenot (CEA, LIST, MetaboHUB, W4M Core Development Team) **Email:** [etienne.thevenot(at)cea.fr](mailto:etienne.thevenot@cea.fr) **Citation:** Delabriere A., Hohenester U., Colsch B., Junot C., Fenaille F. and Thevenot E.A. proFIA: A data preprocessing workflow for Flow Injection Analysis coupled to High-Resolution Mass Spectrometry. *submitted*. @@ -55,6 +55,12 @@ ### News +###### CHANGES IN VERSION 3.0.6 + +NEW FEATURES + + * New (advanced) parameters available + ###### CHANGES IN VERSION 3.0.4 MINOR MODIFICATION
--- a/profia_config.xml Wed May 03 10:39:00 2017 -0400 +++ b/profia_config.xml Thu Jun 29 09:21:32 2017 -0400 @@ -1,315 +1,367 @@ -<tool id="profia" name="proFIA" version="3.0.4"> - <description>Preprocessing of FIA-HRMS data</description> - - <requirements> - <requirement type="package">r-batch</requirement> - <requirement type="package">r-FNN</requirement> - <requirement type="package">r-maxLik</requirement> - <requirement type="package">r-minpack.lm</requirement> - <requirement type="package">r-pracma</requirement> - <requirement type="package">bioconductor-xcms</requirement> - <requirement type="package">bioconductor-plasFIA</requirement> - <requirement type="package">bioconductor-proFIA</requirement> - </requirements> - - <stdio> - <exit_code range="1:" level="fatal" /> - </stdio> - - <command> - Rscript $__tool_directory__/profia_wrapper.R - - #if $inputs.input == "lib": - library $__app__.config.user_library_import_dir/$__user_email__/$inputs.library - #elif $inputs.input == "zip_file": - zipfile $inputs.zip_file - #end if - - ppmN "$ppmN" - ppmGroupN "$ppmGroupN" - fracGroupN "$fracGroupN" - kI "$kI" - - dataMatrix_out "$dataMatrix_out" - sampleMetadata_out "$sampleMetadata_out" - variableMetadata_out "$variableMetadata_out" - figure "$figure" - information "$information" - </command> - - <inputs> - <conditional name="inputs"> - <param name="input" type="select" label="Choose your input method" > - <option value="zip_file" selected="true">Zip file from your history containing your raw files</option> - <option value="lib" >Library directory name</option> - </param> - <when value="zip_file"> - <param name="zip_file" type="data" format="no_unzip.zip,zip" label="Zip file (see the details for file upload in the help section below)" /> - </when> - <when value="lib"> - <param name="library" type="text" size="40" label="Library directory name" help="The name of your directory containing all your data" > - <validator type="empty_field"/> - </param> - </when> - </conditional> - - <param name="ppmN" label="Maximum deviation between centroids during band detection (in ppm)" type="text" value = "5" help="[ppm]" /> - <param name="ppmGroupN" label="Accuracy of the mass spectrometer to be used during feature alignment (in ppm)" type="text" value = "5" help="[ppmGroup] Should be inferior or equal to the deviation parameter above." /> - <param name="fracGroupN" label=" Minimum fraction of samples in which a peak should be detected in at least one class to be kept during feature alignment" type="text" value = "0.5" help="[fracGroup]" /> - <param name="kI" label="Number of neighbour features to be used for imputation (select 0 to skip the imputation step)" type="text" value = "5" help="[k]" /> - </inputs> - - <outputs> - <data name="dataMatrix_out" label="${tool.name}_dataMatrix.tsv" format="tabular" ></data> - <data name="sampleMetadata_out" label="${tool.name}_sampleMetadata.tsv" format="tabular" ></data> - <data name="variableMetadata_out" label="${tool.name}_variableMetadata.tsv" format="tabular" ></data> - <data name="figure" label="${tool.name}_figure.pdf" format="pdf"/> - <data name="information" label="${tool.name}_information.txt" format="txt"/> - </outputs> - - <tests> - <test> - <param name="inputs|input" value="zip_file" /> - <param name="inputs|zip_file" value="input-plasFIA.zip" ftype="zip" /> - <param name="ppmN" value="2"/> - <param name="ppmGroupN" value="1"/> - <param name="fracGroupN" value="0.1"/> - <param name="kI" value="2"/> - <output name="dataMatrix_out" file="output-dataMatrix.tsv" /> - <output name="information"> - <assert_contents> - <has_text text="722 groups have been done" /> - <has_text text="3 samples x 644 variables" /> - <has_text text="78 excluded variables (near zero variance)" /> - <has_text text="2101 peaks detected" /> - </assert_contents> - </output> - </test> - </tests> - - <help> - -.. class:: infomark - -**Author** Alexis Delabriere and Etienne Thevenot (CEA, LIST, MetaboHUB Paris, etienne.thevenot@cea.fr) - ---------------------------------------------------- - -.. class:: infomark - -**Please cite** - -Delabriere A., Hohenester U., Colsch B., Junot C., Fenaille F. and Thevenot E.A. *proFIA*: A data preprocessing workflow for Flow Injection Analysis coupled to High-Resolution Mass Spectrometry. *submitted*. - ---------------------------------------------------- - -.. class:: infomark - -**R package** - -The **proFIA** package is available from the bioconductor repository `http://bioconductor.org/packages/proFIA <http://bioconductor.org/packages/proFIA>`_ - ---------------------------------------------------- - -.. class:: infomark - -**Tool updates** - -See the **NEWS** section at the bottom of this page - ---------------------------------------------------- - -========================================================== -*proFIA*: A preprocessing workflow for FIA-HRMS data -========================================================== - ------------ -Description ------------ - -**Flow Injection Analysis coupled to High-Resolution Mass Spectrometry (FIA-HRMS)** is a promising approach for **high-throughput metabolomics** (Madalinski *et al.*, 2008; Fuhrer *et al.*, 2011; Draper *et al.*, 2013). FIA- HRMS data, however, cannot be preprocessed with current software tools which rely on liquid chromatography separation, or handle low resolution data only. - -The **proFIA module is a workflow** allowing to preprocess FIA-HRMS raw data in **centroid** mode and open format (netCDF, mzData, mzXML, and mzML), and generates the table of peak intensities (**peak table**). The workflow consists in **peak detection and quantification** within individual sample files, followed by **alignment** between files in the m/z dimension, and **imputation** of the missing values in the final peak table (Delabriere *et al.*, submitted). For each ion, the graph representing the intensity as a function of time is called a **flowgram**. A flowgram can be modeled as I = kP + ME(P) + B + e, where k is the response factor (corresponding to the ionization properties of the analyte), P is the **sample peak** (normalized profile which is common for all analytes from a sample and depends on the flow injection conditions only), ME is the **matrix effect**, B is the **solvent baseline**, and e is the heteroscedastic noise. - -The generated peak table is available in the '3 table' W4M tabular format (**dataMatrix**, **sampleMetadata**, and **variableMetadata**) for downstream statistical analysis and annotation with W4M modules. - -A figure provides **diagnostics** and visualization of the preprocessed data set. - ---------------------------------------------------- - -.. class:: infomark - -**References** - -| Delabriere A., Hohenester U., Junot C. and Thevenot E.A. proFIA: A data preprocessing workflow for Flow Injection Analysis coupled to High-Resolution Mass Spectrometry. *submitted*. -| Draper J., Lloyd A., Goodacre R. and Beckmann M. (2013). Flow infusion electrospray ionisation mass spectrometry for high throughput, non-targeted metabolite fingerprinting: a review. *Metabolomics* 9, 4-29. (http://dx.doi.org/10.1007/s11306-012-0449-x) -| Fuhrer T., Dominik H., Boris B. and Zamboni N. (2011). High-throughput, accurate mass metabolome profiling of cellular extracts by flow injection-time-of-flight mass spectrometry. *Analytical Chemistry* 83, 7074-7080. (http://dx.doi.org/10.1021/ac201267k) -| Madalinski G., Godat E., Alves S., Lesage D., Genin E., Levi P., Labarre J., Tabet J., Ezan E. and Junot, C. (2008). Direct introduction of biological samples into a LTQ-orbitrap hybrid mass spectrometer as a tool for fast metabolome analysis. *Analytical Chemistry* 80, 3291-3303. (http://dx.doi.org/10.1021/ac7024915) - ---------------------------------------------------- - ------------------ -Workflow position ------------------ - -.. image:: profia_workflowPositionImage.png - :width: 600 - ------------ -Input files ------------ - -+---------------------------+------------+ -| Parameter : num + label | Format | -+===========================+============+ -| 1 : Choose your inputs | zip | -+---------------------------+------------+ - ---------------------------------------------------- - -.. class:: warningmark - -VERY IMPORTANT: Your data must be in **centroid** mode (centroidization of raw files and conversion to an open format can be achieved with the proteowizard software: http://proteowizard.sourceforge.net/). - - -You have two methods for your inputs: - | Zip file (recommended): You can put a zip file containing your inputs: myinputs.zip (containing all your conditions as sub-directories; see below). - | library folder: You must specify the name of your "library" (folder) created within your space project (for example: /projet/externe/institut/login/galaxylibrary/yourlibrary). Your library must contain all your conditions as sub-directories. - -**Steps for creating the zip file** - -**Step1: Creating your directory and hierarchize the subdirectories** - -.. class:: warningmark - -VERY IMPORTANT: If you zip your files under Windows, you must use the **7Zip** software (http://www.7-zip.org/), otherwise your zip will not be well unzipped on the platform W4M (zip corrupted bug). - -1a) Prepare a parent folder with the name of your data set (e.g., 'arabidopsis') containing your files: - | 'arabidopsis/w1.raw' - | 'arabidopsis/w2.raw' - | ... - | 'arabidopsis/m1.raw' - | 'arabidopsis/m2.raw' - | ... - | - -1b) If you have several experimental conditions resulting in distinct profiles of your samples (e.g. 'wild-type' and 'mutant' genotypes), create subfolders for your files (e.g., 'wild' and 'mutant') into your parent folder: - | 'arabidopsis/wild/w1.raw' - | 'arabidopsis/wild/w2.raw' - | ... - | 'arabidopsis/mutant/m1.raw' - | 'arabidopsis/mutant/m2.raw' - | ... - | - -**Step2: Creating a zip file** - | Zip your **parent** folder (here the 'arabidopsis' folder) containing all the subfolders and files with **7Zip**. - | - -**Step 3 : Uploading it to our Galaxy server** - | If your zip file is less than 2Gb, you get use the **Upload File** tool and the **no_unzip.zip** type to upload it. - | Otherwise if your zip file is larger than 2Gb, please refer to the HOWTO on workflow4metabolomics.org (http://application.sb-roscoff.fr/download/w4m/howto/galaxy_upload_up_2Go.pdf). - | For more informations, don't hesitate to send us an email at supportATworkflow4metabolomics.org). - | - ----------- -Parameters ----------- - -Maximum deviation between centroids during band detection; in ppm (default = 5) - | m/z tolerance of centroids corresponding to the same ion from one scan to the other. - | - -Accuracy of the mass spectrometer to be used during feature alignment; in ppm (default = 5) - | Should be inferior or equal to the deviation parameter above. - | - -Minimum fraction of samples in which a peak should be detected in at least one class to be kept during feature alignment (default = 0.5) - | Identical to the corresponding parameter in XCMS. - | - -Number of neighbour features to be used for imputation (default = 5) - | Select 0 to skip the imputation step. - | - - ------------- -Output files ------------- - -dataMatrix.tabular - | **dataMatrix** tabular separated file with the variables as rows and samples as columns. Missing values are indicated as 'NA' (i.e. when the signal was not significantly different from noise). - | - -sampleMetadata.tabular - | **sampleMetadata** tabular separated file containing the sample metadata as columns. - | - -variableMetadata.tabular - | **variableMetadata** tabular separated file containing the variable metadata as columns. The **timeShifted** flag is set to 1 when the flowgram is time shifted compared to the sample peak (probably due to liquid retention in the FI tube). The **corSampPeakMean** metric is the correlation between the feature flowgram and the sample peak (values are in [-1, 1]). A value below 0.2 suggests that the feature signal is affected by a strong matrix effect. The **meanSolvent** is the mean baseline signal in the feature flowgrams. The **signalOverSolventPvalueMean** is the mean p-value of the tests discriminating between signal and baseline solvent. - | - -figure.pdf - | Visualization and diagnostics about the preprocessed data set; **Feature quality**: Number of detected features per sample for each of the three categories: 'Well-behaved' features have a peak shape close to the sample peak (optimal FIA acquisition is achieved when the majority of the features fall into this category); 'Shifted' indicates a time shift compared to the sample peak, and probably results from retention in the FI tube; 'Significant Matrix Effect' corresponds to a correlation between the feature and the samples peaks of less than 0.2, which is usually caused by a strong matrix effect; **Sample peaks**: Visualization of the peak model for each sample; should have close shapes in case of similar FIA conditions; **m/z density**: may allow to detect a missing m/z value, and in turn, suggest that the *ppm* parameter should be modified; **PCA score plot** of the log10 intensities to detect sample outliers. - | - -information.txt - | Text file with all messages and warnings generated during the computation. - | - ---------------------------------------------------- - ---------------- -Working example ---------------- - -Figure output -============= - -.. image:: profia_workingExampleImage.png - :width: 600 - ---------------------------------------------------- - ----- -NEWS ----- - -CHANGES IN VERSION 3.0.4 -======================== - -MINOR MODIFICATION - -Details added in the documentation - -CHANGES IN VERSION 3.0.2 -======================== - -NEW FEATURE - -Parallel processing - - -CHANGES IN VERSION 3.0.0 -======================== - -NEW FEATURE - -Creation of the tool - -</help> - -<citations> - <citation type="bibtex">@Article{DelabriereSubmitted, - Title = {proFIA: A data preprocessing workflow for Flow Injection Analysis coupled to High-Resolution Mass Spectrometry}, - Author = {Delabriere, Alexis and Hohenester, Ulli and Colsch, Benoit and Junot, Christophe and Fenaille, Francois and Thevenot, Etienne A}, - Journal = {submitted}, - Year = {submitted}, - Pages = {--}, - Volume = {}, - Doi = {} - }</citation> - <citation type="doi">10.1093/bioinformatics/btu813</citation> -</citations> - -</tool> +<tool id="profia" name="proFIA" version="3.0.6"> + <description>Preprocessing of FIA-HRMS data</description> + + <requirements> + <requirement type="package">r-batch</requirement> + <requirement type="package">bioconductor-proFIA</requirement> + </requirements> + + <stdio> + <exit_code range="1:" level="fatal" /> + </stdio> + + <command> + Rscript $__tool_directory__/profia_wrapper.R + + #if $inputs.input == "lib": + library $__app__.config.user_library_import_dir/$__user_email__/$inputs.library + #elif $inputs.input == "zip_file": + zipfile $inputs.zip_file + #end if + + ppmN "$ppmN" + dmzN "$dmzN" + ppmGroupN "$ppmGroupN" + dmzGroupN "$dmzGroupN" + fracGroupN "$fracGroupN" + kI "$kI" + + #if $advCpt.opcC == "full" + bandCoverageN "$advCpt.bandCoverageN" + sizeMinN "$advCpt.sizeMinN" + scanMinI "$advCpt.scanMinI" + scanMaxI "$advCpt.scanMaxI" + #end if + + dataMatrix_out "$dataMatrix_out" + sampleMetadata_out "$sampleMetadata_out" + variableMetadata_out "$variableMetadata_out" + figure "$figure" + information "$information" + </command> + + <inputs> + <conditional name="inputs"> + <param name="input" type="select" label="Choose your input method" > + <option value="zip_file" selected="true">Zip file from your history containing your raw files</option> + <option value="lib" >Library directory name</option> + </param> + <when value="zip_file"> + <param name="zip_file" type="data" format="no_unzip.zip,zip" label="Zip file (see the details for file upload in the help section below)" /> + </when> + <when value="lib"> + <param name="library" type="text" size="40" label="Library directory name" help="The name of your directory containing all your data" > + <validator type="empty_field"/> + </param> + </when> + </conditional> + + <param name="ppmN" label="Maximum deviation between centroids during band detection (in ppm)" type="text" value = "7" help="[ppm]" /> + <param name="dmzN" label="Minimal maximum deviation between centroids during band detection (in Da)" type="text" value = "0.001" help="[dmz] shloud be at most 0.002 for high resolution" /> + <param name="ppmGroupN" label="Accuracy of the mass spectrometer to be used during feature alignment (in ppm)" type="text" value = "3" help="[ppmGroup] Should be inferior to the ppm parameter above." /> + <param name="dmzGroupN" label="Minimal accuracy of the mass spectrometer to be used during feature alignment (in Da)" type="text" value = "0.0005" help="[dmzGroup] shloud be at most 0.002 for high resolution" /> + <param name="fracGroupN" label=" Minimum fraction of samples in which a peak should be detected in at least one class to be kept during feature alignment" type="text" value = "0.5" help="[fracGroup]" /> + <param name="kI" label="Number of neighbour features to be used for imputation (select 0 to skip the imputation step)" type="text" value = "5" help="[k]" /> + + + + <conditional name="advCpt"> + <param name="opcC" type="select" label="Advanced parameters" > + <option value="default" selected="true">Use default</option> + <option value="full">Full parameter list</option> + </param> + <when value="default"/> + <when value="full"> + <param name="bandCoverageN" type="float" value="0.3" label="Minimum fraction of centroids in the estimated injection window for a band to be built" help="[bandCoverage] Must be between 0 and 1"/> + <param name="sizeMinN" type="text" value="none" label="Minimum number of consecutive centroids for a band to be built" help="[sizeMin] If set to 'none', the half of the estimated injection window will be used"/> + <param name="scanMinI" type="integer" value="1" label="First scan to be preprocessed" help="[scanMin]"/> + <param name="scanMaxI" type="text" value="none" label="Last scan to be preprocessed" help="[scanMax] Set to 'none' to preprocess up to the last acquired scan"/> + </when> + </conditional> + + + + </inputs> + + <outputs> + <data name="dataMatrix_out" label="${tool.name}_dataMatrix.tsv" format="tabular" ></data> + <data name="sampleMetadata_out" label="${tool.name}_sampleMetadata.tsv" format="tabular" ></data> + <data name="variableMetadata_out" label="${tool.name}_variableMetadata.tsv" format="tabular" ></data> + <data name="figure" label="${tool.name}_figure.pdf" format="pdf"/> + <data name="information" label="${tool.name}_information.txt" format="txt"/> + </outputs> + + <tests> + <test> + <param name="inputs|input" value="zip_file" /> + <param name="inputs|zip_file" value="input-plasFIA.zip" ftype="zip" /> + <param name="ppmN" value="2"/> + <param name="dmzN" value="0.0005"/> + <param name="ppmGroupN" value="1"/> + <param name="dmzGroupN" value="0.0005"/> + <param name="fracGroupN" value="0.1"/> + <param name="kI" value="2"/> + <output name="dataMatrix_out" file="output-dataMatrix.tsv" /> + <output name="information"> + <assert_contents> + <has_text text="722 groups have been done" /> + <has_text text="3 samples x 644 variables" /> + <has_text text="78 excluded variables (near zero variance)" /> + <has_text text="2101 peaks detected" /> + </assert_contents> + </output> + </test> + </tests> + + <help> + +.. class:: infomark + +**Author** Alexis Delabriere and Etienne Thevenot (CEA, LIST, MetaboHUB Paris, etienne.thevenot@cea.fr) + +--------------------------------------------------- + +.. class:: infomark + +**Please cite** + +Delabriere A., Hohenester U., Colsch B., Junot C., Fenaille F. and Thevenot E.A. *proFIA*: A data preprocessing workflow for Flow Injection Analysis coupled to High-Resolution Mass Spectrometry. *submitted*. + +--------------------------------------------------- + +.. class:: infomark + +**R package** + +The **proFIA** package is available from the bioconductor repository `http://bioconductor.org/packages/proFIA <http://bioconductor.org/packages/proFIA>`_ + +--------------------------------------------------- + +.. class:: infomark + +**Tool updates** + +See the **NEWS** section at the bottom of this page + +--------------------------------------------------- + +========================================================== +*proFIA*: A preprocessing workflow for FIA-HRMS data +========================================================== + +----------- +Description +----------- + +**Flow Injection Analysis coupled to High-Resolution Mass Spectrometry (FIA-HRMS)** is a promising approach for **high-throughput metabolomics** (Madalinski *et al.*, 2008; Fuhrer *et al.*, 2011; Draper *et al.*, 2013). FIA- HRMS data, however, cannot be preprocessed with current software tools which rely on liquid chromatography separation, or handle low resolution data only. + +The **proFIA module is a workflow** allowing to preprocess FIA-HRMS raw data in **centroid** mode and open format (netCDF, mzData, mzXML, and mzML), and generates the table of peak intensities (**peak table**). The workflow consists in **peak detection and quantification** within individual sample files, followed by **alignment** between files in the m/z dimension, and **imputation** of the missing values in the final peak table (Delabriere *et al.*, submitted). For each ion, the graph representing the intensity as a function of time is called a **flowgram**. A flowgram can be modeled as I = kP + ME(P) + B + e, where k is the response factor (corresponding to the ionization properties of the analyte), P is the **sample peak** (normalized profile which is common for all analytes from a sample and depends on the flow injection conditions only), ME is the **matrix effect**, B is the **solvent baseline**, and e is the heteroscedastic noise. + +The generated peak table is available in the '3 table' W4M tabular format (**dataMatrix**, **sampleMetadata**, and **variableMetadata**) for downstream statistical analysis and annotation with W4M modules. + +A figure provides **diagnostics** and visualization of the preprocessed data set. + +--------------------------------------------------- + +.. class:: infomark + +**References** + +| Delabriere A., Hohenester U., Junot C. and Thevenot E.A. proFIA: A data preprocessing workflow for Flow Injection Analysis coupled to High-Resolution Mass Spectrometry. *submitted*. +| Draper J., Lloyd A., Goodacre R. and Beckmann M. (2013). Flow infusion electrospray ionisation mass spectrometry for high throughput, non-targeted metabolite fingerprinting: a review. *Metabolomics* 9, 4-29. (http://dx.doi.org/10.1007/s11306-012-0449-x) +| Fuhrer T., Dominik H., Boris B. and Zamboni N. (2011). High-throughput, accurate mass metabolome profiling of cellular extracts by flow injection-time-of-flight mass spectrometry. *Analytical Chemistry* 83, 7074-7080. (http://dx.doi.org/10.1021/ac201267k) +| Madalinski G., Godat E., Alves S., Lesage D., Genin E., Levi P., Labarre J., Tabet J., Ezan E. and Junot, C. (2008). Direct introduction of biological samples into a LTQ-orbitrap hybrid mass spectrometer as a tool for fast metabolome analysis. *Analytical Chemistry* 80, 3291-3303. (http://dx.doi.org/10.1021/ac7024915) + +--------------------------------------------------- + +----------------- +Workflow position +----------------- + +.. image:: profia_workflowPositionImage.png + :width: 600 + +----------- +Input files +----------- + ++---------------------------+------------+ +| Parameter : num + label | Format | ++===========================+============+ +| 1 : Choose your inputs | zip | ++---------------------------+------------+ + +--------------------------------------------------- + +.. class:: warningmark + +VERY IMPORTANT: Your data must be in **centroid** mode (centroidization of raw files and conversion to an open format can be achieved with the proteowizard software: http://proteowizard.sourceforge.net/). + + +You have two methods for your inputs: + | Zip file (recommended): You can put a zip file containing your inputs: myinputs.zip (containing all your conditions as sub-directories; see below). + | library folder: You must specify the name of your "library" (folder) created within your space project (for example: /projet/externe/institut/login/galaxylibrary/yourlibrary). Your library must contain all your conditions as sub-directories. + +**Steps for creating the zip file** + +**Step1: Creating your directory and hierarchize the subdirectories** + +.. class:: warningmark + +VERY IMPORTANT: If you zip your files under Windows, you must use the **7Zip** software (http://www.7-zip.org/), otherwise your zip will not be well unzipped on the platform W4M (zip corrupted bug). + +1a) Prepare a parent folder with the name of your data set (e.g., 'arabidopsis') containing your files: + | 'arabidopsis/w1.raw' + | 'arabidopsis/w2.raw' + | ... + | 'arabidopsis/m1.raw' + | 'arabidopsis/m2.raw' + | ... + | + +1b) If you have several experimental conditions resulting in distinct profiles of your samples (e.g. 'wild-type' and 'mutant' genotypes), create subfolders for your files (e.g., 'wild' and 'mutant') into your parent folder: + | 'arabidopsis/wild/w1.raw' + | 'arabidopsis/wild/w2.raw' + | ... + | 'arabidopsis/mutant/m1.raw' + | 'arabidopsis/mutant/m2.raw' + | ... + | + +**Step2: Creating a zip file** + | Zip your **parent** folder (here the 'arabidopsis' folder) containing all the subfolders and files with **7Zip**. + | + +**Step 3 : Uploading it to our Galaxy server** + | If your zip file is less than 2Gb, you get use the **Upload File** tool and the **no_unzip.zip** type to upload it. + | Otherwise if your zip file is larger than 2Gb, please refer to the HOWTO on workflow4metabolomics.org (http://application.sb-roscoff.fr/download/w4m/howto/galaxy_upload_up_2Go.pdf). + | For more informations, don't hesitate to send us an email at supportATworkflow4metabolomics.org). + | + +---------- +Parameters +---------- + +Maximum deviation between centroids during band detection; in ppm (default = 7) + | m/z tolerance of centroids corresponding to the same ion from one scan to the other. + | + +Minimal maximum deviation between centroids during band detection; in Da (default = 0.001); to avoid bias at low mass values, the deviation is the maximum between this quantity and the deviation in ppm + | minimum m/z tolerance of centroids corresponding to the same ion from one scan to the other. + | + +Accuracy of the mass spectrometer to be used during feature alignment; in ppm (default = 3); should be less than the ppm parameter used for detection + | Should be inferior or equal to the ppm deviation parameter above. + | + +Minimal accuracy of the mass spectrometer to be used during feature alignment; in Da (default = 0.0005); to avoid bias at low mass values; the deviation is the maximum between this quantity and the deviation in ppm + | + +Minimum fraction of samples in which a peak should be detected in at least one class to be kept during feature alignment (default = 0.5) + | Identical to the corresponding parameter in XCMS. + | + +Number of neighbour features to be used for imputation (default = 5) + | Select 0 to skip the imputation step. + | + +Minimum fraction of centroids in the estimated injection window for a band to be built (default = 0.3) + | + + +Minimum number of consecutive centroids for a band to be built (default = half of the size of the estimated injection window) + | + +First scan to be preprocessed (default = 1) + | + +Last scan to be preprocessed (default = last acquisition scan) + | + +------------ +Output files +------------ + +dataMatrix.tabular + | **dataMatrix** tabular separated file with the variables as rows and samples as columns. Missing values are indicated as 'NA' (i.e. when the signal was not significantly different from noise). + | + +sampleMetadata.tabular + | **sampleMetadata** tabular separated file containing the sample metadata as columns. + | + +variableMetadata.tabular + | **variableMetadata** tabular separated file containing the variable metadata as columns. The **timeShifted** flag is set to 1 when the flowgram is time shifted compared to the sample peak (probably due to liquid retention in the FI tube). The **corSampPeakMean** metric is the correlation between the feature flowgram and the sample peak (values are in [-1, 1]). A value below 0.2 suggests that the feature signal is affected by a strong matrix effect. The **meanSolvent** is the mean baseline signal in the feature flowgrams. The **signalOverSolventPvalueMean** is the mean p-value of the tests discriminating between signal and baseline solvent. + | + +figure.pdf + | Visualization and diagnostics about the preprocessed data set; **Feature quality**: Number of detected features per sample for each of the three categories: 'Well-behaved' features have a peak shape close to the sample peak (optimal FIA acquisition is achieved when the majority of the features fall into this category); 'Shifted' indicates a time shift compared to the sample peak, and probably results from retention in the FI tube; 'Significant Matrix Effect' corresponds to a correlation between the feature and the samples peaks of less than 0.2, which is usually caused by a strong matrix effect; **Sample peaks**: Visualization of the peak model for each sample; should have close shapes in case of similar FIA conditions; **m/z density**: may allow to detect a missing m/z value, and in turn, suggest that the *ppm* parameter should be modified; **PCA score plot** of the log10 intensities to detect sample outliers. + | + +information.txt + | Text file with all messages and warnings generated during the computation. + | + +--------------------------------------------------- + +--------------- +Working example +--------------- + +Figure output +============= + +.. image:: profia_workingExampleImage.png + :width: 600 + +--------------------------------------------------- + +---- +NEWS +---- +CHANGES IN VERSION 3.0.6 +======================== + +NEW FEATURE + +dmz (and dmzGroup) parameters added for the peak detection and grouping steps; bandCoverage, sizeMin, scanMin, and scanMax added as advanced parameters for peak detection + + +CHANGES IN VERSION 3.0.4 +======================== + +MINOR MODIFICATION + +Details added in the documentation + +CHANGES IN VERSION 3.0.2 +======================== + +NEW FEATURE + +Parallel processing + + +CHANGES IN VERSION 3.0.0 +======================== + +NEW FEATURE + +Creation of the tool + +</help> + +<citations> + <citation type="bibtex">@Article{DelabriereSubmitted, + Title = {proFIA: A data preprocessing workflow for Flow Injection Analysis coupled to High-Resolution Mass Spectrometry}, + Author = {Delabriere, Alexis and Hohenester, Ulli and Colsch, Benoit and Junot, Christophe and Fenaille, Francois and Thevenot, Etienne A}, + Journal = {submitted}, + Year = {submitted}, + Pages = {--}, + Volume = {}, + Doi = {} + }</citation> + <citation type="doi">10.1093/bioinformatics/btu813</citation> +</citations> + +</tool>
--- a/profia_wrapper.R Wed May 03 10:39:00 2017 -0400 +++ b/profia_wrapper.R Thu Jun 29 09:21:32 2017 -0400 @@ -1,210 +1,262 @@ -#!/usr/bin/env Rscript - -library(batch) ## parseCommandArgs - -argVc <- unlist(parseCommandArgs(evaluate=FALSE)) - -##------------------------------ -## Initializing -##------------------------------ - - -## libraries -##---------- - -suppressMessages(library(proFIA)) - - -## constants -##---------- - -modNamC <- "proFIA" ## module name - - -## log file -##--------- - -sink(argVc["information"]) - -cat("\nStart of the '", modNamC, "' Galaxy module call: ", - format(Sys.time(), "%a %d %b %Y %X"), "\n", sep="") - - -## arguments -##---------- - - -if("zipfile" %in% names(argVc)) { - - zipfile <- argVc["zipfile"] - - ## We unzip automatically the raw files from the zip file - - if(exists("zipfile") && (zipfile!="")) { - if(!file.exists(zipfile)){ - error_message=paste("Cannot access the Zip file:", zipfile) - print(error_message) - stop(error_message) - } - - ## unzip - - suppressWarnings(unzip(zipfile, unzip="unzip")) - - ## get the directory name - - filesInZip=unzip(zipfile, list=T); - directories=unique(unlist(lapply(strsplit(filesInZip$Name,"/"), function(x) x[1]))); - directories=directories[!(directories %in% c("__MACOSX")) & file.info(directories)$isdir] - directory = "." - if (length(directories) == 1) directory = directories - - cat("files_root_directory\t",directory,"\n") - - } - -} else if ("library" %in% names(argVc)) { - - directory <- argVc["library"] - - if(!file.exists(directory)) { - - error_message=paste("Cannot access the directory:", directory,". Please check that the directory really exists.") - print(error_message) - stop(error_message) - - } - -} else { - - error_message <- "No zipfile nor input library available" - print(error_message) - stop(error_message) - -} - -##------------------------------ -## Computations -##------------------------------ - - -optWrnN <- options()$warn -options(warn = -1) - -stpI <- 1 - -cat("\n", stpI, ") Peak detection step ('proFIAset'):\n", sep = "") - -fiaset <- proFIAset(directory, - ppm = as.numeric(argVc["ppmN"]), - parallel = TRUE) - -stpI <- stpI + 1 - -cat("\n", stpI, ") Peak alignment ('group.FIA'):\n", sep = "") - -fiaset <- group.FIA(fiaset, - ppmGroup = as.numeric(argVc["ppmGroupN"]), - fracGroup = as.numeric(argVc["fracGroupN"])) - -stpI <- stpI + 1 - -cat("\n", stpI, ") Creating the peak table ('makeDataMatrix'):\n", sep = "") - -fiaset <- makeDataMatrix(fiaset, - maxo = FALSE) - -stpI <- stpI + 1 - -kI <- as.integer(argVc["kI"]) - -if(kI > 0) { - - cat("\n", stpI, ") Imputing missing values ('imputeMissingValues.WKNN_TN'):\n", sep = "") - - fiaset <- imputeMissingValues.WKNN_TN(fiaset, - k = kI) - - stpI <- stpI + 1 -} - -options(warn = optWrnN) - - -##------------------------------ -## Ending -##------------------------------ - - -## Plotting -##--------- - -cat("\n", stpI, ") Plotting ('plot'):\n", sep = "") - -pdf(argVc["figure"]) - -plot(fiaset) - -dev.off() - -stpI <- stpI + 1 - -## Printing -##--------- - -cat("\n", stpI, ") Printing ('show'):\n", sep = "") - -fiaset - -stpI <- stpI + 1 - -## Exporting -##---------- - -cat("\n", stpI, ") Exporting ('exportDataMatrix', 'exportSampleMetadata', 'exportVariableMetadata'):\n", sep = "") - -datMN <- exportDataMatrix(fiaset) -samDF <- exportSampleMetadata(fiaset) -varDF <- exportVariableMetadata(fiaset) - -if(nrow(datMN) == nrow(samDF) && ncol(datMN) == nrow(varDF)) { - datDF <- as.data.frame(t(datMN)) -} else { - datDF <- as.data.frame(datMN) -} -rownames(varDF) <- rownames(datDF) - -datDF <- cbind.data.frame(dataMatrix = rownames(datDF), - datDF) -write.table(datDF, - file = argVc["dataMatrix_out"], - quote = FALSE, - row.names = FALSE, - sep = "\t") - -samDF <- cbind.data.frame(sampleMetadata = samDF[, "sampleID"], - class = samDF[, "class"]) -write.table(samDF, - file = argVc["sampleMetadata_out"], - quote = FALSE, - row.names = FALSE, - sep = "\t") - -varDF <- cbind.data.frame(variableMetadata = rownames(varDF), - varDF) -write.table(varDF, - file = argVc["variableMetadata_out"], - quote = FALSE, - row.names = FALSE, - sep = "\t") - - -## Closing -##-------- - -cat("\nEnd of '", modNamC, "' Galaxy module call: ", - as.character(Sys.time()), "\n", sep = "") - -sink() - -rm(list = ls()) +#!/usr/bin/env Rscript + +library(batch) ## parseCommandArgs + +argVc <- unlist(parseCommandArgs(evaluate=FALSE)) + +##------------------------------ +## Initializing +##------------------------------ + + +## libraries +##---------- + +suppressMessages(library(proFIA)) + + +## constants +##---------- + +modNamC <- "proFIA" ## module name + + +## log file +##--------- + +sink(argVc["information"]) + +cat("\nStart of the '", modNamC, "' Galaxy module call: ", + format(Sys.time(), "%a %d %b %Y %X"), "\n", sep="") + + +## arguments +##---------- + + +if("zipfile" %in% names(argVc)) { + + zipfile <- argVc["zipfile"] + + ## We unzip automatically the raw files from the zip file + + if(exists("zipfile") && (zipfile!="")) { + if(!file.exists(zipfile)){ + error_message=paste("Cannot access the Zip file:", zipfile) + print(error_message) + stop(error_message) + } + + ## unzip + + suppressWarnings(unzip(zipfile, unzip="unzip")) + + ## get the directory name + + filesInZip=unzip(zipfile, list=T); + directories=unique(unlist(lapply(strsplit(filesInZip$Name,"/"), function(x) x[1]))); + directories=directories[!(directories %in% c("__MACOSX")) & file.info(directories)$isdir] + directory = "." + if (length(directories) == 1) directory = directories + + cat("files_root_directory\t",directory,"\n") + + } + +} else if ("library" %in% names(argVc)) { + + directory <- argVc["library"] + + if(!file.exists(directory)) { + + error_message=paste("Cannot access the directory:", directory,". Please check that the directory really exists.") + print(error_message) + stop(error_message) + + } + +} else { + + error_message <- "No zipfile nor input library available" + print(error_message) + stop(error_message) + +} + +##------------------------------ +## Computations +##------------------------------ + + +optWrnN <- options()$warn +options(warn = -1) + +stpI <- 1 + +cat("\n", stpI, ") Peak detection step ('proFIAset'):\n", sep = "") + +if("sizeMinN" %in% names(argVc) && argVc["sizeMinN"] != "none") { + if("scanMaxI" %in% names(argVc) && argVc["scanMaxI"] != "none") { + fiaset <- proFIAset(directory, + ppm = as.numeric(argVc["ppmN"]), + dmz = as.numeric(argVc["dmzN"]), + bandCoverage = as.numeric(ifelse("bandCoverageN" %in% names(argVc), argVc["bandCoverageN"], "0.3")), + sizeMin = as.numeric(argVc["sizeMinN"]), + scanmin = as.numeric(ifelse("scanMinI" %in% names(argVc), argVc["scanMinI"], "1")), + scanmax = as.numeric(argVc["scanMaxI"]), + parallel = FALSE) + } else { + fiaset <- proFIAset(directory, + ppm = as.numeric(argVc["ppmN"]), + dmz = as.numeric(argVc["dmzN"]), + bandCoverage = as.numeric(ifelse("bandCoverageN" %in% names(argVc), argVc["bandCoverageN"], "0.3")), + sizeMin = as.numeric(argVc["sizeMinN"]), + scanmin = as.numeric(ifelse("scanMinI" %in% names(argVc), argVc["scanMinI"], "1")), + parallel = FALSE) + } +} else { + if("scanMaxI" %in% names(argVc) && argVc["scanMaxI"] != "none") { + fiaset <- proFIAset(directory, + ppm = as.numeric(argVc["ppmN"]), + dmz = as.numeric(argVc["dmzN"]), + bandCoverage = as.numeric(ifelse("bandCoverageN" %in% names(argVc), argVc["bandCoverageN"], "0.3")), + scanmin = as.numeric(ifelse("scanMinI" %in% names(argVc), argVc["scanMinI"], "1")), + scanmax = as.numeric(argVc["scanMaxI"]), + parallel = FALSE) + } else { + fiaset <- proFIAset(directory, + ppm = as.numeric(argVc["ppmN"]), + dmz = as.numeric(argVc["dmzN"]), + bandCoverage = as.numeric(ifelse("bandCoverageN" %in% names(argVc), argVc["bandCoverageN"], "0.3")), + scanmin = as.numeric(ifelse("scanMinI" %in% names(argVc), argVc["scanMinI"], "1")), + parallel = FALSE) + } +} + +stpI <- stpI + 1 + +cat("\n", stpI, ") Peak alignment ('group.FIA'):\n", sep = "") + +fiaset <- group.FIA(fiaset, + ppmGroup = as.numeric(argVc["ppmGroupN"]), + dmz = as.numeric(argVc["dmzGroupN"]), + fracGroup = as.numeric(argVc["fracGroupN"])) + +stpI <- stpI + 1 + +cat("\n", stpI, ") Creating the peak table ('makeDataMatrix'):\n", sep = "") + +fiaset <- makeDataMatrix(fiaset, + maxo = FALSE) + +stpI <- stpI + 1 + +kI <- as.integer(argVc["kI"]) + + +###TODO add the two method for imputation. +if(kI > 0) { + + cat("\n", stpI, ") Imputing missing values ('imputeMissingValues.WKNN_TN'):\n", sep = "") + + fiaset <- imputeMissingValues.WKNN_TN(fiaset, + k = kI) + + stpI <- stpI + 1 +} + +options(warn = optWrnN) + + +##------------------------------ +## Ending +##------------------------------ + + +## Plotting +##--------- + +cat("\n", stpI, ") Plotting ('plot'):\n", sep = "") + +pdf(argVc["figure"]) + +plot(fiaset) + +dev.off() + +stpI <- stpI + 1 + +## Printing +##--------- + +cat("\n", stpI, ") Printing ('show'):\n", sep = "") + +fiaset + +stpI <- stpI + 1 + +## Exporting +##---------- + +cat("\n", stpI, ") Exporting ('exportDataMatrix', 'exportSampleMetadata', 'exportVariableMetadata'):\n", sep = "") + +datMN <- exportDataMatrix(fiaset) +samDF <- exportSampleMetadata(fiaset) +varDF <- exportVariableMetadata(fiaset) + +if(nrow(datMN) == nrow(samDF) && ncol(datMN) == nrow(varDF)) { + datDF <- as.data.frame(t(datMN)) +} else { + datDF <- as.data.frame(datMN) +} +rownames(varDF) <- rownames(datDF) + +datDF <- cbind.data.frame(dataMatrix = rownames(datDF), + datDF) +write.table(datDF, + file = argVc["dataMatrix_out"], + quote = FALSE, + row.names = FALSE, + sep = "\t") + +samDF <- cbind.data.frame(sampleMetadata = samDF[, "sampleID"], + class = samDF[, "class"]) +write.table(samDF, + file = argVc["sampleMetadata_out"], + quote = FALSE, + row.names = FALSE, + sep = "\t") + +varDF <- cbind.data.frame(variableMetadata = rownames(varDF), + varDF) +write.table(varDF, + file = argVc["variableMetadata_out"], + quote = FALSE, + row.names = FALSE, + sep = "\t") + + +## Closing +##-------- + +cat("\nEnd of '", modNamC, "' Galaxy module call: ", + as.character(Sys.time()), "\n", sep = "") + +cat("\n\n\n============================================================================") +cat("\nAdditional information about the call:\n") +cat("\n1) Parameters:\n") +print(cbind(value = argVc)) + +cat("\n2) Session Info:\n") +sessioninfo <- sessionInfo() +cat(sessioninfo$R.version$version.string,"\n") +cat("Main packages:\n") +for (pkg in names(sessioninfo$otherPkgs)) { cat(paste(pkg,packageVersion(pkg)),"\t") }; cat("\n") +cat("Other loaded packages:\n") +for (pkg in names(sessioninfo$loadedOnly)) { cat(paste(pkg,packageVersion(pkg)),"\t") }; cat("\n") + +cat("============================================================================\n") + +sink() + +rm(list = ls())
--- a/runit/output/information.txt Wed May 03 10:39:00 2017 -0400 +++ b/runit/output/information.txt Thu Jun 29 09:21:32 2017 -0400 @@ -1,5 +1,5 @@ -Start of the 'proFIA' Galaxy module call: Wed 03 May 2017 01:18:28 PM +Start of the 'proFIA' Galaxy module call: Thu 29 Jun 2017 02:34:32 PM files_root_directory plasFIA 1) Peak detection step ('proFIAset'): @@ -32,4 +32,33 @@ 7) Exporting ('exportDataMatrix', 'exportSampleMetadata', 'exportVariableMetadata'): -End of 'proFIA' Galaxy module call: 2017-05-03 13:27:09 +End of 'proFIA' Galaxy module call: 2017-06-29 14:35:49 + + + +============================================================================ +Additional information about the call: + +1) Parameters: + value +dataMatrix_out "./output/dataMatrix.tsv" +sampleMetadata_out "./output/sampleMetadata.tsv" +variableMetadata_out "./output/variableMetadata.tsv" +figure "./output/figure.pdf" +information "./output/information.txt" +zipfile "./plasfia/plasFIA.zip" +library "NULL" +ppmN "2" +dmzN "5e-04" +ppmGroupN "1" +dmzGroupN "5e-04" +fracGroupN "0.1" +kI "2" + +2) Session Info: +R version 3.3.1 (2016-06-21) +Main packages: +proFIA 1.1.10 xcms 1.50.1 Biobase 2.34.0 ProtGenerics 1.6.0 BiocGenerics 0.20.0 mzR 2.8.1 Rcpp 0.12.10 batch 1.1.4 +Other loaded packages: +RANN 2.5 quadprog 1.5.5 lattice 0.20.34 codetools 0.2.15 ropls 1.6.2 MASS 7.3.45 MassSpecWavelet 1.40.0 grid 3.3.1 plyr 1.8.4 stats4 3.3.1 pracma 1.9.9 S4Vectors 0.12.2 Matrix 1.2.8 splines 3.3.1 BiocParallel 1.8.1 RColorBrewer 1.1.2 survival 2.41.2 multtest 2.30.0 minpack.lm 1.2.1 +============================================================================
--- a/runit/output/sampleMetadata.tsv Wed May 03 10:39:00 2017 -0400 +++ b/runit/output/sampleMetadata.tsv Thu Jun 29 09:21:32 2017 -0400 @@ -1,4 +1,4 @@ -sampleMetadata samDF[, colnames(samDF) != "sampleID"] +sampleMetadata class C100a plasFIA C100b plasFIA C100c plasFIA
--- a/runit/output/variableMetadata.tsv Wed May 03 10:39:00 2017 -0400 +++ b/runit/output/variableMetadata.tsv Thu Jun 29 09:21:32 2017 -0400 @@ -13,11 +13,11 @@ M93.9950 93.9950324348781 93.994384765625 93.9957275390625 62 18 80 1 2229.79814801897 0 0 0.241612105825766 NA 22.6816537513546 0 M93.9950 M94.0000 94.000032500426 93.9993209838867 94.0006256103516 56 18 74 1 4834.00592041016 0 0 0.518934009424996 NA 12.8192023875339 0 M94.0000 M94.0078 94.0078254629064 94.0072174072266 94.008430480957 58 18 76 1 0 0 0 0.338538468735688 NA Inf 0 M94.0078 -M94.0124 94.0124465985536 94.0117416381836 94.0132141113281 64 18 82 1 2487.75439453125 0 1.41340072488116e-05 0.241464414184748 NA 16.8897639472362 0 M94.0124 +M94.0124 94.0124465985536 94.0117416381836 94.0132141113281 64 18 82 1 2487.75439453125 0 2.38033559047945e-05 0.241464414184748 NA 16.8897639472362 0 M94.0124 M94.0174 94.0173825687832 94.0166091918945 94.0180206298828 58 18 76 1 0 0 0 0.333006203441651 NA Inf 0 M94.0174 M94.0209 94.0208812366832 94.0201263427734 94.0214691162109 59 18 77 1 0 0 0 0.0445717895358037 NA Inf 0 M94.0209 M94.0262 94.0261501924491 94.025016784668 94.0277557373047 57 18 78 3 787.069051106771 0 0 0.267671736241134 0.0967164409808468 Inf 0 M94.0262 -M94.0321 94.0321010044643 94.0315017700195 94.032600402832 58 18 76 1 2696.71176147461 0 1.91379576186534e-06 0.236049991681362 NA 11.1483885383609 0 M94.0321 +M94.0321 94.0321010044643 94.0315017700195 94.032600402832 58 18 76 1 2696.71176147461 0 1.84994154861506e-06 0.236049991681362 NA 11.1483885383609 0 M94.0321 M94.0362 94.0361511230469 94.0353240966797 94.036865234375 51 29 80 1 2594.17755126953 0 0 0.396159126481773 NA 22.0100748537854 0 M94.0362 M94.0401 94.0401373142149 94.0394821166992 94.0407485961914 56 18 74 1 0 0 0 0.220869925323141 NA Inf 0 M94.0401 M94.0412 94.0412086607918 94.0401000976562 94.0417709350586 62 18 80 1 1766.9701944987 0 0 0.325103568177275 NA 35.0611157384355 0 M94.0412 @@ -35,10 +35,10 @@ M97.0761 97.0760913028895 97.0756530761719 97.0768127441406 67 18 93 3 0 0 0 0.155597889189398 0.141726241422026 Inf 1 M97.0761 M98.9189 98.9188891128755 98.9184875488281 98.9192733764648 58 18 77 3 0 0 0 0.708550417601752 0.0330146041692216 Inf 1 M98.9189 M102.0550 102.054975323093 102.054145812988 102.055130004883 85.5 18 110 2 0 0 0 -0.0763748293823842 0.12399837351998 Inf 0 M102.0550 -M102.0914 102.091388230064 102.090888977051 102.091537475586 56.5 18 77 2 2805.0125342475 0 1.77635683940025e-15 0.484732295239128 0.034646520978652 59.4814046128523 0 M102.0914 +M102.0914 102.091388230064 102.090888977051 102.091537475586 56.5 18 77 2 2805.0125342475 0 2.1094237467878e-15 0.484732295239128 0.034646520978652 59.4814046128523 0 M102.0914 M103.9195 103.919523073592 103.918594360352 103.920143127441 68 20 88 1 0 0 0 0.133331495721408 NA Inf 0 M103.9195 M103.9419 103.941940524714 103.940971374512 103.942558288574 62 20 82 1 0 0 0 0.353898796522856 NA Inf 0 M103.9419 -M104.1070 104.107025015158 104.106399536133 104.107345581055 70 18 88 2 66408.7600097656 0 0.00807211414361064 0.554019761810103 0 2.85420720740854 0 M104.1070 +M104.1070 104.107025015158 104.106399536133 104.107345581055 70 18 88 2 66408.7600097656 0 0.00808148692391319 0.554019761810103 0 2.85420720740854 0 M104.1070 M105.0699 105.069903960595 105.06950378418 105.070068359375 59 18 81 3 0 0 0 0.723271760888928 0.0189401875961994 Inf 1 M105.0699 M109.0760 109.076044828995 109.07551574707 109.076278686523 65 18 88 3 0 0 0 0.547613912169008 0.0203886595190367 Inf 1 M109.0760 M112.8957 112.895723540207 112.895500183105 112.895957946777 57 18 81 3 0 0 0 0.896542745371031 0.0124734497677377 Inf 1 M112.8957 @@ -58,21 +58,21 @@ M121.9045 121.90451741773 121.903656005859 121.905151367188 52 18 70 1 3299.53299967448 0 0 0.35979141869005 NA 25.5605654228482 0 M121.9045 M121.9157 121.915666286762 121.914970397949 121.916404724121 55 18 73 1 3720.49703543527 0 0 0.265184846004965 NA 17.5705987143516 0 M121.9157 M121.9208 121.920824323382 121.920143127441 121.921508789062 57 18 75 1 0 0 0 0.379590485496505 NA Inf 0 M121.9208 -M121.9247 121.924704031393 121.924163818359 121.92529296875 54 18 72 1 2871.18501790365 0 0.00183211477060885 0.294610695154223 NA 19.334063378455 0 M121.9247 +M121.9247 121.924704031393 121.924163818359 121.92529296875 54 18 72 1 2871.18501790365 0 0.00181624237685507 0.294610695154223 NA 19.334063378455 0 M121.9247 M121.9269 121.926917457581 121.926383972168 121.927490234375 57 18 75 1 3953.99332682292 0 0 0.411844205770503 NA 20.2418407914345 0 M121.9269 -M121.9312 121.931184823172 121.930320739746 121.931838989258 57 18 75 1 6809.88911132813 0 9.18439990726938e-10 0.32883957955863 NA 8.73740372805123 0 M121.9312 +M121.9312 121.931184823172 121.930320739746 121.931838989258 57 18 75 1 6809.88911132813 0 2.13148809713459e-09 0.32883957955863 NA 8.73740372805123 0 M121.9312 M121.9333 121.93329679049 121.932685852051 121.93376159668 62 18 80 1 4771.90551757812 0 0 0.302715322716007 NA 15.5896042025612 0 M121.9333 M121.9391 121.939097086589 121.938369750977 121.939819335938 59 18 77 1 3120.59478759766 0 0 0.492457989676642 NA 26.8791182895792 0 M121.9391 -M121.9524 121.952375531814 121.951477050781 121.95344543457 55 18 75 2 5891.23073032924 0 0.00205191314669573 0.346161763409438 0.0233719318873183 15.5795670284052 0 M121.9524 +M121.9524 121.952375531814 121.951477050781 121.95344543457 55 18 75 2 5891.23073032924 0 0.00207292995794228 0.346161763409438 0.0233719318873183 15.5795670284052 0 M121.9524 M121.9552 121.955201235684 121.954330444336 121.955688476562 55 18 73 1 3296.8867594401 0 0 0.648592985293298 NA 17.716989634183 0 M121.9552 -M121.9704 121.970371002847 121.969505310059 121.971054077148 68 20 88 1 3930.30981445312 0 9.12603326241879e-13 -0.0243177970733837 NA 8.60577780715165 0 M121.9704 -M121.9738 121.973795013428 121.972999572754 121.974403381348 62 18 80 1 3120.05792236328 0 8.47752136468038e-06 0.39289723179297 NA 39.0796586475515 0 M121.9738 +M121.9704 121.970371002847 121.969505310059 121.971054077148 68 20 88 1 3930.30981445312 0 1.00170982619829e-11 -0.0243177970733837 NA 8.60577780715165 0 M121.9704 +M121.9738 121.973795013428 121.972999572754 121.974403381348 62 18 80 1 3120.05792236328 0 8.40008728508757e-06 0.39289723179297 NA 39.0796586475515 0 M121.9738 M121.9836 121.983639053434 121.982757568359 121.984390258789 56 18 74 1 0 0 0 0.345486543396654 NA Inf 0 M121.9836 -M121.9888 121.988755337028 121.987869262695 121.989288330078 68 20 88 1 4050.15152994792 0 0.00688033038198532 0.0834597319910627 NA 12.0235260714065 0 M121.9888 -M121.9938 121.993844909668 121.99308013916 121.994384765625 60 18 78 1 4903.10241699219 0 0.000108919766856452 0.308633220363451 NA 9.26678172126117 0 M121.9938 -M121.9957 121.995664068993 121.994728088379 121.996292114258 62 18 80 1 3707.43551025391 0 0.00662612996762424 -0.136020299861705 NA 7.37780668730212 0 M121.9957 +M121.9888 121.988755337028 121.987869262695 121.989288330078 68 20 88 1 4050.15152994792 0 0.00809522593537904 0.0834597319910627 NA 12.0235260714065 0 M121.9888 +M121.9938 121.993844909668 121.99308013916 121.994384765625 60 18 78 1 4903.10241699219 0 0.000163219855365604 0.308633220363451 NA 9.26678172126117 0 M121.9938 +M121.9957 121.995664068993 121.994728088379 121.996292114258 62 18 80 1 3707.43551025391 0 0.00848493229529357 -0.136020299861705 NA 7.37780668730212 0 M121.9957 M122.0018 122.001771381923 122.000793457031 122.002372741699 58 35 93 1 0 0 0 -0.327737245344874 NA Inf 0 M122.0018 -M122.0480 122.04799041748 122.047302246094 122.048614501953 55 21 76 1 4416.31009928385 0 3.59601237676088e-12 0.362589557541678 NA 10.5154443536103 0 M122.0480 +M122.0480 122.04799041748 122.047302246094 122.048614501953 55 21 76 1 4416.31009928385 0 3.48587825271807e-12 0.362589557541678 NA 10.5154443536103 0 M122.0480 M122.0542 122.054221667832 122.053321838379 122.055038452148 59 18 80 3 1790.52235921224 0 0 0.648019126801235 0.0530872640458131 Inf 1 M122.0542 M122.0579 122.057903925578 122.057342529297 122.058380126953 50 18 68 1 0 0 0 0.284643679837114 NA Inf 0 M122.0579 M122.0747 122.074725267364 122.07405090332 122.075180053711 57 18 75 1 0 0 0 0.360307354181837 NA Inf 0 M122.0747 @@ -92,14 +92,14 @@ M132.1019 132.101878821083 132.101150512695 132.102416992188 72 18 93 3 5555.52066040039 0 0 0.906307657122881 0.0158505997084159 280.688340680345 1 M132.1019 M133.0971 133.097127339926 133.09700012207 133.097549438477 82 18 101 3 0 0 0 -0.173345431638996 0.244310413339532 Inf 1 M133.0971 M133.1052 133.10521538313 133.10481262207 133.105361938477 57 18 80 3 0 0 0 0.409149545072417 0.0816349332465311 Inf 1 M133.1052 -M134.0448 134.044755497583 134.044189453125 134.044891357422 90 26 122 2 1436.37357584635 0 1.10072731340782e-06 -0.438988624740501 0.0414663583523287 Inf 0 M134.0448 +M134.0448 134.044755497583 134.044189453125 134.044891357422 90 26 122 2 1436.37357584635 0 1.07383660419291e-06 -0.438988624740501 0.0414663583523287 Inf 0 M134.0448 M134.0964 134.096401333921 134.095794677734 134.097106933594 62 18 88 3 0 0 0 0.964456324969677 0.00799507482558049 Inf 1 M134.0964 M134.9884 134.988384731321 134.988098144531 134.988555908203 95.5 18 117 2 0 0 0 -0.531573256777832 0.0635798838633485 Inf 0 M134.9884 M135.0028 135.002830932015 135.002044677734 135.003311157227 60 18 86 3 1538.4963921441 0 0 0.921724953674066 0.0117443576005922 Inf 1 M135.0028 M135.0301 135.030094347502 135.029983520508 135.030212402344 59 25 84 1 0 0 0 0.0898480597351847 NA Inf 0 M135.0301 M135.0998 135.099750180788 135.099578857422 135.099914550781 60 18 86 3 0 0 0 0.633817132369559 0.0333361819132554 Inf 1 M135.0998 M136.0757 136.075662535082 136.075271606445 136.076187133789 71 18 94 3 0 0 0 0.931305249206762 0.0112041333734416 Inf 1 M136.0757 -M137.0573 137.057273513969 137.05712890625 137.057434082031 66 20 86 1 13551.7610909598 0 0.00944290975712225 0.0782446222599314 NA 4.09866770779716 0 M137.0573 +M137.0573 137.057273513969 137.05712890625 137.057434082031 66 20 86 1 13551.7610909598 0 0.0093071754431775 0.0782446222599314 NA 4.09866770779716 0 M137.0573 M137.0790 137.079008622603 137.078887939453 137.079147338867 58 19 77 1 0 0 0 0.534050030397736 NA Inf 0 M137.0790 M138.9065 138.906451034546 138.905853271484 138.907165527344 58 18 78 3 0 0 0 0.68790562665958 0.0309477735373366 Inf 1 M138.9065 M140.0682 140.068174501006 140.067749023438 140.068344116211 57 18 86 3 9466.70444980197 0 0 0.622536247986117 0.0298220039083829 15.9228368192193 1 M140.0682 @@ -150,7 +150,7 @@ M169.0922 169.092206682478 169.092025756836 169.092544555664 61 18 83 3 0 0 0 0.758803007884244 0.0616033213115718 Inf 1 M169.0922 M170.0424 170.042356460813 170.04216003418 170.042617797852 63 18 88 3 0 0 0 0.629777485574536 0.0718838711624807 Inf 1 M170.0424 M170.0924 170.092384338379 170.092132568359 170.09260559082 74 18 95 3 0 0 0 -0.156417733773211 0.168339791065726 Inf 1 M170.0924 -M170.0964 170.096408918792 170.096206665039 170.096618652344 59 26 85 1 25635.5879991319 0 0.00166977889507702 -0.139229768136581 NA 4.29640842461642 0 M170.0964 +M170.0964 170.096408918792 170.096206665039 170.096618652344 59 26 85 1 25635.5879991319 0 0.00170355808776779 -0.139229768136581 NA 4.29640842461642 0 M170.0964 M172.0944 172.09440111844 172.09375 172.094589233398 57 18 75 3 0 0 0 0.894796953167192 0.0230874547627367 Inf 1 M172.0944 M173.4724 173.472379181903 173.4716796875 173.472961425781 65.5 18 91 6 22129.7952564381 0 0 0.942001107207437 0.0389662769466327 17.6344446154136 1 M173.4724 M174.0386 174.038636618798 174.038330078125 174.038940429688 59 18 88 3 0 0 0 0.89512349579656 0.0170842540834756 Inf 1 M174.0386 @@ -166,7 +166,7 @@ M176.1182 176.118201870389 176.117691040039 176.118286132812 70 18 91 3 337963.152832031 0 0 0.987895793598041 0.00284698066504405 222.667788699681 1 M176.1182 M176.1280 176.12801028895 176.127426147461 176.128326416016 68 18 89 3 0 0 0 0.851530960642003 0.00779524812116477 Inf 1 M176.1280 M177.0546 177.054586839816 177.053955078125 177.054916381836 71.5 18 93 4 8010.88876994779 0 0 0.948165256999034 0.0171360945267089 273.99779209461 1 M177.0546 -M177.0869 177.086939179848 177.086395263672 177.087173461914 216 36 252 3 9160.38512369792 1 6.62292443109891e-12 NA NA 23.520555575446 1 M177.0869 +M177.0869 177.086939179848 177.086395263672 177.087173461914 216 36 252 3 9160.38512369792 1 6.77968792217598e-12 NA NA 23.520555575446 1 M177.0869 M177.0982 177.098172651854 177.097640991211 177.098434448242 78 18 99 3 0 0 0 0.494070220776238 0.0959562642562823 Inf 1 M177.0982 M177.1152 177.115224685669 177.114990234375 177.115493774414 65 18 83 3 0 0 0 0.904584032264878 0.0112768904073397 Inf 1 M177.1152 M177.1215 177.121526013167 177.12092590332 177.12174987793 68 18 88 3 0 0 0 0.981634413473625 0.0050994426900174 Inf 1 M177.1215 @@ -239,7 +239,7 @@ M224.0786 224.078640783257 224.078216552734 224.079055786133 70 18 88 3 1828.75838216146 0 0 0.923010602842713 0.0105310057686198 Inf 1 M224.0786 M224.8648 224.864822998047 224.864379882812 224.865234375 57 18 88 3 0 0 0 0.731595708083634 0.0855501852857934 Inf 1 M224.8648 M225.0345 225.034489246125 225.033798217773 225.034942626953 55 18 86 3 0 0 0 0.645015166171156 0.0345754251172273 Inf 1 M225.0345 -M225.0429 225.04286731672 225.042449951172 225.04328918457 60.5 18 80 2 33891.0204264323 0 3.79037612407984e-09 0.692319788386783 0.0108420288862525 5.61410973692305 0 M225.0429 +M225.0429 225.04286731672 225.042449951172 225.04328918457 60.5 18 80 2 33891.0204264323 0 4.06171063449534e-09 0.692319788386783 0.0108420288862525 5.61410973692305 0 M225.0429 M225.0820 225.082030785353 225.08171081543 225.08251953125 64 20 89 2 0 0 0 0.475592868108778 0.0496078488000051 Inf 0 M225.0820 M226.1049 226.104926963604 226.104446411133 226.105255126953 54 18 72 2 0 0 0 0.211976612772608 0.142168571267995 Inf 0 M226.1049 M226.9514 226.95143281168 226.951080322266 226.951782226562 54 18 73 3 112911.991391782 0 0 0.874780814193066 0.0302213467281532 6.58493759702176 1 M226.9514 @@ -295,7 +295,7 @@ M264.0556 264.055600827879 264.055023193359 264.056060791016 56 18 86 3 0 0 0 0.800486385472207 0.0805841668592475 Inf 1 M264.0556 M264.0594 264.059358799096 264.058868408203 264.059783935547 58 18 77 2 0 0 0 0.669677420477548 0.00369204662661352 Inf 0 M264.0594 M264.8525 264.852476640181 264.851959228516 264.852905273438 57 18 77 2 0 0 0 0.58943838344449 0.0111881499310966 Inf 0 M264.8525 -M265.0230 265.022975569722 265.022521972656 265.0234375 90 18 109 3 15600.1009260786 0 8.67137778914895e-07 -0.628675736411646 0.011694614163826 35.4236614886184 1 M265.0230 +M265.0230 265.022975569722 265.022521972656 265.0234375 90 18 109 3 15600.1009260786 0 8.72576002564074e-07 -0.628675736411646 0.011694614163826 35.4236614886184 1 M265.0230 M265.0529 265.052898551844 265.05224609375 265.053436279297 58 18 86 3 0 0 0 0.752584440443732 0.130323707158315 Inf 1 M265.0529 M265.2525 265.252521920028 265.251861572266 265.253051757812 63 23 88 2 6558.70487060547 0 0 0.169367880454861 0.178291019723902 6.45877462502979 0 M265.2525 M267.2682 267.268156571822 267.267639160156 267.268676757812 66 18 88 3 3654.70815429687 0 0 0.320591197514126 0.131585315749652 Inf 1 M267.2682 @@ -341,7 +341,7 @@ M297.0271 297.027119524726 297.026489257812 297.027679443359 58.5 18 79 2 0 0 0 0.317646408127947 0.0546001213907387 Inf 0 M297.0271 M298.0208 298.020847440516 298.020172119141 298.021545410156 60 18 84 3 0 0 0 0.735060344197707 0.064756142123763 Inf 1 M298.0208 M300.2896 300.289557762969 300.289001464844 300.290222167969 70 18 89 3 33603.8711066675 0 0 0.88662532357832 0.0302962881929054 123.714518118438 1 M300.2896 -M301.1409 301.140912222438 301.140167236328 301.141632080078 64.5 18 86 2 17705.1561550564 0 0.00234111973882412 0.447013765289295 0.0814584931272674 5.404218584928 0 M301.1409 +M301.1409 301.140912222438 301.140167236328 301.141632080078 64.5 18 86 2 17705.1561550564 0 0.00237001333398612 0.447013765289295 0.0814584931272674 5.404218584928 0 M301.1409 M301.2929 301.292914130471 301.292297363281 301.293609619141 70 18 92 3 9031.70613606771 0 0 0.784667192022304 0.0287785005980248 Inf 1 M301.2929 M302.2324 302.232421488702 302.231536865234 302.233001708984 64 18 87 3 0 0 0 0.842658955348086 0.0325266255988136 Inf 1 M302.2324 M302.2963 302.296271690956 302.296051025391 302.296600341797 55 18 73 1 0 0 0 0.294033201199406 NA Inf 0 M302.2963 @@ -358,7 +358,7 @@ M312.1569 312.156864149306 312.156005859375 312.157501220703 62 18 86 3 0 0 0 0.880802107443958 0.0158887984073676 Inf 1 M312.1569 M313.0882 313.088225046794 313.087615966797 313.088714599609 57 18 76 3 0 0 0 0.378550388973032 0.124230888293443 Inf 1 M313.0882 M313.2736 313.273577284112 313.273040771484 313.273986816406 61 18 80 4 36186.4700520833 0 0 0.753191656322666 0.123814657517948 20.6817266985135 1 M313.2736 -M314.2769 314.276949081186 314.276336669922 314.277587890625 66 18 86 3 8006.73060062767 0 5.31476960219379e-08 0.459632596652839 0.0803067892521975 18.3713031174221 1 M314.2769 +M314.2769 314.276949081186 314.276336669922 314.277587890625 66 18 86 3 8006.73060062767 0 5.65191864421166e-08 0.459632596652839 0.0803067892521975 18.3713031174221 1 M314.2769 M315.0673 315.067252431017 315.066131591797 315.067779541016 216 18 252 3 0 1 0 NA NA Inf 1 M315.0673 M315.0797 315.079722243089 315.079376220703 315.080078125 99 20 126 3 0 0 0 -0.68906946665802 0.0139322530260814 Inf 1 M315.0797 M316.0692 316.069240231348 316.067169189453 316.071166992188 195 36 252 6 0 1 0 NA NA Inf 1 M316.0692 @@ -433,7 +433,7 @@ M354.0300 354.029998614337 354.028900146484 354.030853271484 98 18 119 3 0 1 0 NA NA Inf 1 M354.0300 M354.0389 354.038914518536 354.038177490234 354.039855957031 81 18 101 3 0 0 0 -0.313152430438497 0.345814923787531 Inf 1 M354.0389 M354.2694 354.269373844379 354.268707275391 354.270080566406 68 18 93 3 17781.068439359 0 0 0.95211091108674 0.00298611608459795 231.419310195002 1 M354.2694 -M355.0112 355.011249735134 355.010223388672 355.012023925781 164 42 206 1 6625.73197428385 1 0.00549094935724148 NA NA 11.1275011786484 0 M355.0112 +M355.0112 355.011249735134 355.010223388672 355.012023925781 164 42 206 1 6625.73197428385 1 0.00557078474807193 NA NA 11.1275011786484 0 M355.0112 M355.0336 355.033565561584 355.032867431641 355.034210205078 70 21 107 2 0 0 0 -0.405672315398789 0.0108397116904648 Inf 0 M355.0336 M355.2715 355.271483666341 355.269317626953 355.273406982422 58 18 89 6 0 0 0 0.622783629489419 0.160540512748458 Inf 1 M355.2715 M356.2558 356.255824497768 356.254974365234 356.256561279297 58 18 78 3 0 0 0 0.480536671595142 0.100929971181428 Inf 1 M356.2558 @@ -448,7 +448,7 @@ M360.3235 360.323476776912 360.322570800781 360.324279785156 65 18 88 3 15831.1493171813 0 0 0.972567939930051 0.00758316581354671 251.688535134582 1 M360.3235 M361.3268 361.326848830532 361.326080322266 361.327697753906 67 18 88 3 0 0 0 0.876631104828427 0.0273292585552364 Inf 1 M361.3268 M362.1508 362.150843227611 362.150177001953 362.151519775391 70 18 97 3 212270.260742188 0 0 0.97582540260119 0.00400275412859612 41.8035022344611 1 M362.1508 -M362.9262 362.926198190812 362.925231933594 362.926910400391 67.5 18 87 2 10166.3879220145 0 0.00236324195672488 0.349933211676286 0.0833818333403346 9.32521639796976 0 M362.9262 +M362.9262 362.926198190812 362.925231933594 362.926910400391 67.5 18 87 2 10166.3879220145 0 0.00236019622957484 0.349933211676286 0.0833818333403346 9.32521639796976 0 M362.9262 M363.1542 363.154240396724 363.153350830078 363.154968261719 70 18 88 3 15458.3132324219 0 0 0.916885705650731 0.0286966816149146 Inf 1 M363.1542 M364.1575 364.157522201538 364.156677246094 364.158233642578 70 18 88 3 0 0 0 0.142296149146464 0.0436374058876441 Inf 1 M364.1575 M366.2314 366.231351758857 366.230529785156 366.232330322266 79.5 23 106 2 0 0 0 -0.405501621590054 0.124234688897691 Inf 0 M366.2314 @@ -458,7 +458,7 @@ M369.0006 369.000565199419 368.999572753906 369.001739501953 88 18 114 3 0 0 0 -0.654145087471826 0.0448084547448312 Inf 1 M369.0006 M369.2309 369.23093821633 369.22998046875 369.231719970703 57 18 78 2 2619.65395507812 0 0 0.421891825219025 0.041301047956437 Inf 0 M369.2309 M369.2400 369.239954148677 369.239074707031 369.2412109375 63 18 83 3 11272.796070886 0 0 0.973276109949078 0.00499418179999374 351.938542871162 1 M369.2400 -M369.3514 369.351375382831 369.350433349609 369.352294921875 65 18 86 3 7054.73975642774 0 1.92888327532614e-08 0.309239611259633 0.129624347807374 17.7434488012771 1 M369.3514 +M369.3514 369.351375382831 369.350433349609 369.352294921875 65 18 86 3 7054.73975642774 0 2.04018750782353e-08 0.309239611259633 0.129624347807374 17.7434488012771 1 M369.3514 M370.2433 370.243298124931 370.242340087891 370.244384765625 65 18 86 3 0 0 0 0.892833979125951 0.0249174597567885 Inf 1 M370.2433 M371.2381 371.238074860749 371.237152099609 371.23876953125 62 18 86 3 0 0 0 0.70781986628087 0.0836370208549056 Inf 1 M371.2381 M375.0086 375.008551890676 375.007781982422 375.009368896484 92 18 113 3 0 0 0 -0.474383073625477 0.109982135862803 Inf 1 M375.0086 @@ -511,16 +511,16 @@ M415.2114 415.211355153225 415.21044921875 415.212493896484 69 18 89 3 10266.9535039993 0 0 0.865964597278162 0.0246775136669939 113.611638510799 1 M415.2114 M415.2365 415.236474117517 415.235595703125 415.237701416016 90 18 109 3 17924.0768454843 0 1.55431223447522e-15 -0.597061058404085 0.0783990254437828 29.5702020840481 1 M415.2365 M416.2147 416.214674049449 416.213684082031 416.215576171875 65 18 91 3 0 0 0 0.75066678683018 0.143760886913443 Inf 1 M416.2147 -M416.2399 416.239894840962 416.238952636719 416.240997314453 84 18 108 2 3708.43896484375 0 0.000132672916000831 -0.580364239997042 0.059907835267169 Inf 0 M416.2399 +M416.2399 416.239894840962 416.238952636719 416.240997314453 84 18 108 2 3708.43896484375 0 0.000135386307470342 -0.580364239997042 0.059907835267169 Inf 0 M416.2399 M417.0021 417.002123564677 417.001068115234 417.002960205078 70 18 94 3 0 0 0 0.255354187797131 0.0213552418363686 Inf 1 M417.0021 M417.0302 417.030242919922 417.029388427734 417.031097412109 117 38 155 1 0 1 0 NA NA Inf 0 M417.0302 M418.9993 418.999282132662 418.998352050781 419.000274658203 71 18 91 3 0 0 0 0.0494637117593832 0.153473293881657 Inf 1 M418.9993 M421.2535 421.253461767126 421.252288818359 421.254364013672 62 18 88 3 0 0 0 0.806372839588106 0.0320789926754558 Inf 1 M421.2535 -M427.2903 427.290252086618 427.289367675781 427.291076660156 58 18 78 2 12315.2435492622 0 0.00151828375191942 -0.310957971859633 0.0814699794032974 19.6091470325956 0 M427.2903 +M427.2903 427.290252086618 427.289367675781 427.291076660156 58 18 78 2 12315.2435492622 0 0.00154461563039987 -0.310957971859633 0.0814699794032974 19.6091470325956 0 M427.2903 M429.0216 429.021555474464 429.020721435547 429.022552490234 101 19 121 3 0 0 0 -0.601328697455925 0.0566373801764399 Inf 1 M429.0216 M429.2027 429.202731892214 429.201721191406 429.203857421875 59 18 88 3 0 0 0 0.85425254420661 0.0126329682047903 Inf 1 M429.2027 M430.2061 430.20613945855 430.204833984375 430.207000732422 55 18 78 3 0 0 0 0.528782558718492 0.108466252314047 Inf 1 M430.2061 -M430.9136 430.913607377272 430.912170410156 430.914733886719 57 18 79 3 11353.0951930455 0 0.000104498070590653 0.2121107733753 0.0988316723592769 7.17762583419622 1 M430.9136 +M430.9136 430.913607377272 430.912170410156 430.914733886719 57 18 79 3 11353.0951930455 0 0.000102903302165558 0.2121107733753 0.0988316723592769 7.17762583419622 1 M430.9136 M432.2378 432.237833513532 432.236755371094 432.238800048828 61 18 87 3 0 0 0 0.915789632091554 0.0181668868762316 Inf 1 M432.2378 M432.9672 432.967202108457 432.965911865234 432.968078613281 71.5 19 109 2 0 0 0 -0.226800147328167 0.374177273074199 Inf 0 M432.9672 M432.9761 432.97612508138 432.975219726562 432.976989746094 68 18 91 3 0 0 0 0.204421438437946 0.175365733336679 Inf 1 M432.9761 @@ -537,8 +537,8 @@ M441.2530 441.252962153009 441.251800537109 441.253723144531 59 18 82 3 0 0 0 0.585250505403101 0.0380231846799636 Inf 1 M441.2530 M441.2974 441.297373644989 441.296020507812 441.298309326172 61 18 82 3 4351.5630812872 0 0 0.38713960968782 0.163991791112649 Inf 1 M441.2974 M442.9960 442.996017940251 442.994934082031 442.997039794922 134.5 18 227 2 0 1 0 NA NA Inf 0 M442.9960 -M443.2678 443.267810964374 443.2666015625 443.268768310547 89 18 111 4 49123.085530599 0 4.15239354012442e-09 -0.641449963420814 0.015365210325455 29.5195726320753 1 M443.2678 -M444.2712 444.271189026913 444.270141601562 444.272521972656 90 18 108 3 15569.8226970693 0 5.72814796593324e-05 -0.648989040526443 0.0316510789434626 16.5738420037641 1 M444.2712 +M443.2678 443.267810964374 443.2666015625 443.268768310547 89 18 111 4 49123.085530599 0 4.16431533700745e-09 -0.641449963420814 0.015365210325455 29.5195726320753 1 M443.2678 +M444.2712 444.271189026913 444.270141601562 444.272521972656 90 18 108 3 15569.8226970693 0 5.91245462555623e-05 -0.648989040526443 0.0316510789434626 16.5738420037641 1 M444.2712 M445.1766 445.176633834839 445.175506591797 445.177673339844 56 18 77 3 0 0 0 0.573261837712446 0.131487670776999 Inf 1 M445.1766 M449.2848 449.284758929715 449.283447265625 449.28564453125 60 18 80 3 3508.10217285156 0 0 0.860482946206747 0.0318459411272742 Inf 1 M449.2848 M450.2882 450.288199462891 450.287445068359 450.2890625 56 18 88 3 0 0 0 0.526622144979491 0.0267900085417479 Inf 1 M450.2882 @@ -670,7 +670,7 @@ M689.3530 689.352976769702 689.351013183594 689.354248046875 87.5 21 110 2 0 0 0 -0.388777078676843 0.0138612836936161 Inf 0 M689.3530 M689.4022 689.402190512799 689.400085449219 689.404296875 59 18 79 3 0 0 0 0.644043004138408 0.122576876367304 Inf 1 M689.4022 M695.0984 695.098381574764 695.096435546875 695.099914550781 81 19 102 3 0 0 0 -0.395403143457194 0.0560264528808235 Inf 1 M695.0984 -M703.4544 703.454356892904 703.451232910156 703.4560546875 55 18 75 3 9808.56670690466 0 2.41759945396325e-11 0.344078194120363 0.0750296818955926 10.8299646875888 1 M703.4544 +M703.4544 703.454356892904 703.451232910156 703.4560546875 55 18 75 3 9808.56670690466 0 2.46678233395414e-11 0.344078194120363 0.0750296818955926 10.8299646875888 1 M703.4544 M703.5748 703.574770078554 703.572692871094 703.576843261719 70 18 90 3 0 0 0 0.0538369468147119 0.0921709053391853 Inf 1 M703.5748 M704.3794 704.37939860026 704.377258300781 704.381286621094 62 18 81 3 0 0 0 0.37170508548455 0.116124494995131 Inf 1 M704.3794 M704.5783 704.578295255962 704.576538085938 704.579772949219 60 18 89 3 0 0 0 -0.0351561983302897 0.19672663237462 Inf 1 M704.5783
--- a/runit/profia_runtests.R Wed May 03 10:39:00 2017 -0400 +++ b/runit/profia_runtests.R Thu Jun 29 09:21:32 2017 -0400 @@ -53,7 +53,7 @@ if(.Platform$OS.type == "windows") wrapperCallC <- paste("Rscript", wrapperCallC) - + wrapperCodeN <- system(wrapperCallC) if (wrapperCodeN != 0)
--- a/runit/profia_tests.R Wed May 03 10:39:00 2017 -0400 +++ b/runit/profia_tests.R Thu Jun 29 09:21:32 2017 -0400 @@ -4,7 +4,9 @@ argLs <- list(zipfile = "./plasfia/plasFIA.zip", library = "NULL", ppmN = "2", + dmzN = "0.0005", ppmGroupN = "1", + dmzGroupN = "0.0005", fracGroupN = "0.1", kI = "2")