# HG changeset patch
# User recetox
# Date 1624024399 0
# Node ID 1b9140057e5a97d251ee8ba4331a0a438f916707
# Parent  b0ba0287b76f728a6a6d5971e91b5ed266579f0e
"planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit accce68183c78552905b4446a8cab581911275ef"

diff -r b0ba0287b76f -r 1b9140057e5a aplcms_macros.xml
--- a/aplcms_macros.xml	Wed Oct 07 16:06:54 2020 +0000
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,145 +0,0 @@
-<macros>
-    <token name="@TOOL_VERSION@">6.6.6</token>
-
-    <xml name="requirements">
-        <requirements>
-            <container type="docker">recetox/aplcms:6.6.6-recetox1</container>
-        </requirements>
-    </xml>
-
-    <xml name="inputs">
-        <inputs>
-            <param name="files" type="data" format="mzdata,mzml,mzxml,netcdf" multiple="true" min="3" label="data"
-                   help="Mass spectrometry files for peak extraction." />
-            <yield />
-        </inputs>
-    </xml>
-
-    <xml name="history_db">
-        <param name="known_table" type="data" format="h5" label="known_table"
-               help="A data table containing the known metabolite ions and previously found features. The table must contain these 18 columns: chemical_formula (optional), HMDB_ID (optional), KEGG_compound_ID (optional), neutral.mass (optional), ion.type (the ion form - optional), m.z (either theoretical or mean observed m/z value of previously found features), Number_profiles_processed (the total number of processed samples to build this database), Percent_found (the percentage of historically processed samples in which the feature appeared), mz_min (minimum  observed m/z value), mz_max (maximum observed m/z value), RT_mean (mean observed retention time), RT_sd (standard deviation of observed retention time), RT_min (minimum observed retention time), RT_max (maximum observed retention time), int_mean.log. (mean observed log intensity), int_sd.log. (standard deviation of observed log intensity), int_min.log. (minimum observed log intensity), int_max.log. (maximum observed log intensity)." />
-        <section name="history_db" title="Known-Table settings">
-            <param name="match_tol_ppm" type="integer" optional="true" min="0" label="match_tol_ppm (optional)"
-                   help="The ppm tolerance to match identified features to known metabolites/features." />
-            <param name="new_feature_min_count" type="integer" value="2" min="1" label="new_feature_min_count"
-                   help="The minimum number of occurrences of a historically unseen (unknown) feature to add this feature into the database of known features." />
-        </section>
-    </xml>
-
-    <xml name="noise_filtering">
-        <section name="noise_filtering" title="Noise filtering and peak detection">
-            <param name="min_exp" type="integer" min="1" value="2"
-                   label="min_exp"
-                   help="If a feature is to be included in the final feature table, it must be present in at least this number of spectra." />
-            <param name="min_pres" type="float" value="0.5"
-                   label="min_pres"
-                   help="The minimum proportion of presence in the time period for a series of signals grouped by m/z to be considered a peak." />
-            <param name="min_run" type="float" value="12"
-                   label="min_run"
-                   help="The minimum length of elution time for a series of signals grouped by m/z to be considered a peak." />
-            <param name="mz_tol" type="float" value="1e-05"
-                   label="mz_tol"
-                   help="The m/z tolerance level for the grouping of data points. This value is expressed as the fraction of the m/z value. This value, multiplied by the m/z value, becomes the cutoff level. The recommended value is the machine's nominal accuracy level. Divide the ppm value by 1e6. For FTMS, 1e-5 is recommended." />
-            <param name="baseline_correct" type="float" value="0" label="baseline_correct"
-                   help="After grouping the observations, the highest intensity in each group is found. If the highest is lower than this value, the entire group will be deleted. The default value is NA, in which case the program uses a percentile of the height of the noise groups. If given a value, the value will be used as the threshold, and baseline.correct.noise.percentile will be ignored." />
-            <param name="baseline_correct_noise_percentile" type="float" value="0.05"
-                   label="baseline_correct_noise_percentile"
-                   help="The percentile of signal strength of those EIC that don't pass the run filter, to be used as the baseline threshold of signal strength." />
-            <param name="intensity_weighted" type="boolean" checked="false" truevalue="TRUE" falsevalue="FALSE"
-                   label="intensity_weighted"
-                   help="Whether to weight the local density by signal intensities in initial peak detection." />
-        </section>
-    </xml>
-
-    <xml name="feature_detection">
-        <section name="feature_detection" title="Feature detection">
-            <param name="shape_model" type="select" display="radio"
-                   label="shape_model"
-                   help="The mathematical model for the shape of a peak. There are two choices - bi-Gaussian and Gaussian. When the peaks are asymmetric, the bi-Gaussian is better.">
-                <option value="Gaussian">Gaussian</option>
-                <option value="bi-Gaussian" selected="true">bi-Gaussian</option>
-            </param>
-            <param name="BIC_factor" type="float" value="2.0"
-                   label="BIC_factor"
-                   help="The factor that is multiplied on the number of parameters to modify the BIC criterion. If larger than 1, models with more peaks are penalized more." />
-            <param name="peak_estim_method" type="select" display="radio"
-                   label="peak_estim_method"
-                   help="The estimation method for the bi-Gaussian peak model. Two possible values: moment and EM.">
-                <option value="moment" selected="true">Moment</option>
-                <option value="EM">EM</option>
-            </param>
-            <param name="min_bandwidth" type="float" optional="true"
-                   label="min_bandwidth (optional)"
-                   help="The minimum bandwidth to use in the kernel smoother." />
-            <param name="max_bandwidth" type="float" optional="true"
-                   label="max_bandwidth (optional)"
-                   help="The maximum bandwidth to use in the kernel smoother." />
-            <param name="sd_cut_min" type="float" value="0.01"
-                   label="sd_cut_min"
-                   help="The minimum standard deviation of a feature to be not eliminated." />
-            <param name="sd_cut_max" type="float" value="500"
-                   label="sd_cut_max"
-                   help="The maximum standard deviation of a feature to be not eliminated." />
-            <param name="sigma_ratio_lim_min" type="float" value="0.01"
-                   label="sigma_ratio_lim_min"
-                   help="The lower limit of the believed ratio range between the left-standard deviation and the right-standard deviation of the bi-Gaussian function used to fit the data." />
-            <param name="sigma_ratio_lim_max" type="float" value="100"
-                   label="sigma_ratio_lim_max"
-                   help="The upper limit of the believed ratio range between the left-standard deviation and the right-standard deviation of the bi-Gaussian function used to fit the data." />
-            <param name="component_eliminate" type="float" value="0.01"
-                   label="component_eliminate"
-                   help="In fitting mixture of bi-Gaussian (or Gaussian) model of an EIC, when a component accounts for a proportion of intensities less than this value, the component will be ignored." />
-            <param name="moment_power" type="float" value="1"
-                   label="moment_power"
-                   help="The power parameter for data transformation when fitting the bi-Gaussian or Gaussian mixture model in an EIC." />
-        </section>
-    </xml>
-
-    <xml name="peak_alignment">
-        <section name="peak_alignment" title="Peak Alignment">
-            <param name="align_chr_tol" type="float" optional="true"
-                   label="align_chr_tol (optional)"
-                   help="The retention time tolerance level for peak alignment. The default is NA, which allows the program to search for the tolerance level based on the data." />
-            <param name="align_mz_tol" type="float" optional="true"
-                   label="align_mz_tol (optional)"
-                   help="The m/z tolerance level for peak alignment. The default is NA, which allows the program to search for the tolerance level based on the data. This value is expressed as the percentage of the m/z value. This value, multiplied by the m/z value, becomes the cutoff level." />
-            <param name="max_align_mz_diff" type="float" value="0.01"
-                   label="max_align_mz_diff"
-                   help="As the m/z tolerance is expressed in relative terms (ppm), it may not be suitable when the m/z range is wide. This parameter limits the tolerance in absolute terms. It mostly influences feature matching in higher m/z range." />
-        </section>
-    </xml>
-
-    <xml name="weak_signal_recovery">
-        <section name="weak_signal_recovery" title="Weak Signal Recovery">
-            <param name="recover_mz_range" type="float" optional="true"
-                   label="recover_mz_range (optional)"
-                   help="The m/z around the feature m/z to search for observations. The default value is NA, in which case 1.5 times the m/z tolerance in the aligned object will be used." />
-            <param name="recover_chr_range" type="float" optional="true"
-                   label="recover_chr_range (optional)"
-                   help="The retention time around the feature retention time to search for observations. The default value is NA, in which case 0.5 times the retention time tolerance in the aligned object will be used." />
-            <param name="use_observed_range" type="boolean" checked="true" truevalue="TRUE" falsevalue="FALSE"
-                   label="use_observed_range"
-                   help="If the value is true, the actual range of the observed locations of the feature in all the spectra will be used." />
-            <param name="recover_min_count" type="integer" value="3"
-                   label="recover_min_count"
-                   help="The minimum number of raw data points to be considered as a true feature." />
-        </section>
-    </xml>
-
-    <token name="@GENERAL_HELP@">
-        apLCMS is a software which generates a feature table from a batch of LC/MS spectra. The m/z and retention time
-        tolerance levels are estimated from the data. A run-filter is used to detect peaks and remove noise.
-        Non-parametric statistical methods are used to find-tune peak selection and grouping. After retention time
-        correction, a feature table is generated by aligning peaks across spectra. For further information on apLCMS
-        please refer to http://web1.sph.emory.edu/apLCMS.
-    </token>
-
-    <xml name="citations">
-        <citations>
-            <citation type="doi">10.1093/bioinformatics/btp291</citation>
-            <citation type="doi">10.1186/1471-2105-11-559</citation>
-            <citation type="doi">10.1021/pr301053d</citation>
-            <citation type="doi">10.1093/bioinformatics/btu430</citation>
-        </citations>
-    </xml>
-</macros>
diff -r b0ba0287b76f -r 1b9140057e5a aplcms_unsupervised.xml
--- a/aplcms_unsupervised.xml	Wed Oct 07 16:06:54 2020 +0000
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,67 +0,0 @@
-<tool id="recetox_aplcms_unsupervised" name="apLCMS - Unsupervised" version="@TOOL_VERSION@+galaxy1">
-    <macros>
-        <import>aplcms_macros.xml</import>
-    </macros>
-
-    <expand macro="requirements" />
-
-    <command detect_errors="aggressive"><![CDATA[
-        #set file_str = str('", "').join([str($f) for $f in $files])
-
-        Rscript
-            -e 'x <- apLCMS::unsupervised(
-                    files = c("$file_str"),
-                    min_exp = $noise_filtering.min_exp,
-                    min_pres = $noise_filtering.min_pres,
-                    min_run = $noise_filtering.min_run,
-                    mz_tol = $noise_filtering.mz_tol,
-                    baseline_correct = $noise_filtering.baseline_correct,
-                    baseline_correct_noise_percentile = $noise_filtering.baseline_correct_noise_percentile,
-                    intensity_weighted = $noise_filtering.intensity_weighted,
-                    shape_model = "$feature_detection.shape_model",
-                    BIC_factor = $feature_detection.BIC_factor,
-                    peak_estim_method = "$feature_detection.peak_estim_method",
-                    min_bandwidth = $feature_detection.min_bandwidth,
-                    max_bandwidth = $feature_detection.max_bandwidth,
-                    sd_cut = c($feature_detection.sd_cut_min, $feature_detection.sd_cut_max),
-                    sigma_ratio_lim = c($feature_detection.sigma_ratio_lim_min, $feature_detection.sigma_ratio_lim_max),
-                    component_eliminate = $feature_detection.component_eliminate,
-                    moment_power = $feature_detection.moment_power,
-                    align_chr_tol = $peak_alignment.align_chr_tol,
-                    align_mz_tol = $peak_alignment.align_mz_tol,
-                    max_align_mz_diff = $peak_alignment.max_align_mz_diff,
-                    recover_mz_range = $weak_signal_recovery.recover_mz_range,
-                    recover_chr_range = $weak_signal_recovery.recover_chr_range,
-                    use_observed_range = $weak_signal_recovery.use_observed_range,
-                    recover_min_count = $weak_signal_recovery.recover_min_count
-                )'
-            -e 'apLCMS::save_peaks_to_hdf("$peaks", x)'
-    ]]></command>
-
-    <expand macro="inputs">
-        <expand macro="noise_filtering" />
-        <expand macro="feature_detection" />
-        <expand macro="peak_alignment" />
-        <expand macro="weak_signal_recovery" />
-    </expand>
-
-    <outputs>
-        <data name="peaks" format="h5" />
-    </outputs>
-
-    <tests>
-        <test>
-            <param name="files" value="mbr_test0.mzml,mbr_test1.mzml,mbr_test2.mzml" ftype="mzml"/>
-            <output name="peaks" file="peaks_unsupervised.h5" ftype="h5" compare="sim_size" delta="1000"/>
-        </test>
-    </tests>
-
-    <help>
-        This is the Unsupervised version of apLCMS which is not relying on any existing knowledge about metabolites or
-        any historically detected features. For such functionality please use the Hybrid version of apLCMS.
-
-        @GENERAL_HELP@
-    </help>
-
-    <expand macro="citations" />
-</tool>
diff -r b0ba0287b76f -r 1b9140057e5a recetox_aplcms_macros.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/recetox_aplcms_macros.xml	Fri Jun 18 13:53:19 2021 +0000
@@ -0,0 +1,159 @@
+<macros>
+    <token name="@TOOL_VERSION@">0.9.3</token>
+
+    <xml name="requirements">
+        <requirements>
+            <requirement type="set_environment">R_SCRIPT_PATH</requirement>
+            <requirement type="package" version="4.1.0">R</requirement>
+            <requirement type="package" version="4.0.1">r-arrow</requirement>
+            <requirement type="package" version="0.9.3">r-recetox-aplcms</requirement>
+        </requirements>
+    </xml>
+
+    <xml name="inputs">
+        <inputs>
+            <param name="files" type="data" format="mzdata,mzml,mzxml,netcdf" multiple="true" min="3" label="data"
+                   help="Mass spectrometry files for peak extraction." />
+            <yield />
+        </inputs>
+    </xml>
+
+    <xml name="history_db">
+        <param name="known_table" type="data" format="parquet" label="known_table"
+               help="A data table containing the known metabolite ions and previously found features. The table must contain these 18 columns: chemical_formula (optional), HMDB_ID (optional), KEGG_compound_ID (optional), neutral.mass (optional), ion.type (the ion form - optional), m.z (either theoretical or mean observed m/z value of previously found features), Number_profiles_processed (the total number of processed samples to build this database), Percent_found (the percentage of historically processed samples in which the feature appeared), mz_min (minimum  observed m/z value), mz_max (maximum observed m/z value), RT_mean (mean observed retention time), RT_sd (standard deviation of observed retention time), RT_min (minimum observed retention time), RT_max (maximum observed retention time), int_mean.log. (mean observed log intensity), int_sd.log. (standard deviation of observed log intensity), int_min.log. (minimum observed log intensity), int_max.log. (maximum observed log intensity)." />
+        <section name="history_db" title="Known-Table settings">
+            <param name="match_tol_ppm" type="integer" optional="true" min="0" label="match_tol_ppm (optional)"
+                   help="The ppm tolerance to match identified features to known metabolites/features." />
+            <param name="new_feature_min_count" type="integer" value="2" min="1" label="new_feature_min_count"
+                   help="The minimum number of occurrences of a historically unseen (unknown) feature to add this feature into the database of known features." />
+        </section>
+    </xml>
+
+    <xml name="noise_filtering">
+        <section name="noise_filtering" title="Noise filtering and peak detection">
+            <param name="min_exp" type="integer" min="1" value="2"
+                   label="min_exp"
+                   help="If a feature is to be included in the final feature table, it must be present in at least this number of spectra." />
+            <param name="min_pres" type="float" value="0.5"
+                   label="min_pres"
+                   help="The minimum proportion of presence in the time period for a series of signals grouped by m/z to be considered a peak." />
+            <param name="min_run" type="float" value="12"
+                   label="min_run"
+                   help="The minimum length of elution time for a series of signals grouped by m/z to be considered a peak." />
+            <param name="mz_tol" type="float" value="1e-05"
+                   label="mz_tol"
+                   help="The m/z tolerance level for the grouping of data points. This value is expressed as the fraction of the m/z value. This value, multiplied by the m/z value, becomes the cutoff level. The recommended value is the machine's nominal accuracy level. Divide the ppm value by 1e6. For FTMS, 1e-5 is recommended." />
+            <param name="baseline_correct" type="float" value="0" label="baseline_correct"
+                   help="After grouping the observations, the highest intensity in each group is found. If the highest is lower than this value, the entire group will be deleted. The default value is NA, in which case the program uses a percentile of the height of the noise groups. If given a value, the value will be used as the threshold, and baseline.correct.noise.percentile will be ignored." />
+            <param name="baseline_correct_noise_percentile" type="float" value="0.05"
+                   label="baseline_correct_noise_percentile"
+                   help="The percentile of signal strength of those EIC that don't pass the run filter, to be used as the baseline threshold of signal strength." />
+            <param name="intensity_weighted" type="boolean" checked="false" truevalue="TRUE" falsevalue="FALSE"
+                   label="intensity_weighted"
+                   help="Whether to weight the local density by signal intensities in initial peak detection." />
+        </section>
+    </xml>
+
+    <xml name="feature_detection">
+        <section name="feature_detection" title="Feature detection">
+            <param name="shape_model" type="select" display="radio"
+                   label="shape_model"
+                   help="The mathematical model for the shape of a peak. There are two choices - bi-Gaussian and Gaussian. When the peaks are asymmetric, the bi-Gaussian is better.">
+                <option value="Gaussian">Gaussian</option>
+                <option value="bi-Gaussian" selected="true">bi-Gaussian</option>
+            </param>
+            <param name="BIC_factor" type="float" value="2.0"
+                   label="BIC_factor"
+                   help="The factor that is multiplied on the number of parameters to modify the BIC criterion. If larger than 1, models with more peaks are penalized more." />
+            <param name="peak_estim_method" type="select" display="radio"
+                   label="peak_estim_method"
+                   help="The estimation method for the bi-Gaussian peak model. Two possible values: moment and EM.">
+                <option value="moment" selected="true">Moment</option>
+                <option value="EM">EM</option>
+            </param>
+            <param name="min_bandwidth" type="float" optional="true"
+                   label="min_bandwidth (optional)"
+                   help="The minimum bandwidth to use in the kernel smoother." />
+            <param name="max_bandwidth" type="float" optional="true"
+                   label="max_bandwidth (optional)"
+                   help="The maximum bandwidth to use in the kernel smoother." />
+            <param name="sd_cut_min" type="float" value="0.01"
+                   label="sd_cut_min"
+                   help="The minimum standard deviation of a feature to be not eliminated." />
+            <param name="sd_cut_max" type="float" value="500"
+                   label="sd_cut_max"
+                   help="The maximum standard deviation of a feature to be not eliminated." />
+            <param name="sigma_ratio_lim_min" type="float" value="0.01"
+                   label="sigma_ratio_lim_min"
+                   help="The lower limit of the believed ratio range between the left-standard deviation and the right-standard deviation of the bi-Gaussian function used to fit the data." />
+            <param name="sigma_ratio_lim_max" type="float" value="100"
+                   label="sigma_ratio_lim_max"
+                   help="The upper limit of the believed ratio range between the left-standard deviation and the right-standard deviation of the bi-Gaussian function used to fit the data." />
+            <param name="component_eliminate" type="float" value="0.01"
+                   label="component_eliminate"
+                   help="In fitting mixture of bi-Gaussian (or Gaussian) model of an EIC, when a component accounts for a proportion of intensities less than this value, the component will be ignored." />
+            <param name="moment_power" type="float" value="1"
+                   label="moment_power"
+                   help="The power parameter for data transformation when fitting the bi-Gaussian or Gaussian mixture model in an EIC." />
+        </section>
+    </xml>
+
+    <xml name="peak_alignment">
+        <section name="peak_alignment" title="Peak Alignment">
+            <param name="align_chr_tol" type="float" optional="true"
+                   label="align_chr_tol (optional)"
+                   help="The retention time tolerance level for peak alignment. The default is NA, which allows the program to search for the tolerance level based on the data." />
+            <param name="align_mz_tol" type="float" optional="true"
+                   label="align_mz_tol (optional)"
+                   help="The m/z tolerance level for peak alignment. The default is NA, which allows the program to search for the tolerance level based on the data. The tolerance is given in absolute numbers, not scaled, i.e. for 10ppm tolerance enter '1e-05'. This value, multiplied by the m/z value, becomes the cutoff level." />
+            <param name="max_align_mz_diff" type="float" value="0.01"
+                   label="max_align_mz_diff"
+                   help="As the m/z tolerance is expressed in relative terms (ppm), it may not be suitable when the m/z range is wide. This parameter limits the tolerance in absolute terms. It mostly influences feature matching in higher m/z range." />
+        </section>
+    </xml>
+
+    <xml name="weak_signal_recovery">
+        <section name="weak_signal_recovery" title="Weak Signal Recovery">
+            <param name="recover_mz_range" type="float" optional="true"
+                   label="recover_mz_range (optional)"
+                   help="The m/z around the feature m/z to search for observations. The default value is NA, in which case 1.5 times the m/z tolerance in the aligned object will be used." />
+            <param name="recover_chr_range" type="float" optional="true"
+                   label="recover_chr_range (optional)"
+                   help="The retention time around the feature retention time to search for observations. The default value is NA, in which case 0.5 times the retention time tolerance in the aligned object will be used." />
+            <param name="use_observed_range" type="boolean" checked="true" truevalue="TRUE" falsevalue="FALSE"
+                   label="use_observed_range"
+                   help="If the value is true, the actual range of the observed locations of the feature in all the spectra will be used." />
+            <param name="recover_min_count" type="integer" value="3"
+                   label="recover_min_count"
+                   help="The minimum number of raw data points to be considered as a true feature." />
+        </section>
+    </xml>
+
+    <xml name="unsupervised_outputs">
+        <data name="recovered_feature_sample_table" format="parquet" label="${tool.name} recovered_feature_sample_table on ${on_string}" />
+        <data name="aligned_feature_sample_table" format="parquet" label="${tool.name} aligned_feature_sample_table on ${on_string}" hidden="true" />
+        <collection  name="corrected_features" type="list" label="${tool.name} corrected_features on ${on_string}">
+            <discover_datasets pattern="__designation__" directory="corrected" format="parquet" />
+        </collection >
+        <collection  name="extracted_features" type="list" label="${tool.name} extracted_features on ${on_string}">
+            <discover_datasets pattern="__designation__" directory="extracted" format="parquet" />
+        </collection >
+    </xml>
+
+    <token name="@GENERAL_HELP@">
+        apLCMS is a software which generates a feature table from a batch of LC/MS spectra. The m/z and retention time
+        tolerance levels are estimated from the data. A run-filter is used to detect peaks and remove noise.
+        Non-parametric statistical methods are used to find-tune peak selection and grouping. After retention time
+        correction, a feature table is generated by aligning peaks across spectra. For further information on apLCMS
+        please refer to https://mypage.cuhk.edu.cn/academics/yutianwei/apLCMS/.
+    </token>
+
+    <xml name="citations">
+        <citations>
+            <citation type="doi">10.1093/bioinformatics/btp291</citation>
+            <citation type="doi">10.1186/1471-2105-11-559</citation>
+            <citation type="doi">10.1021/pr301053d</citation>
+            <citation type="doi">10.1093/bioinformatics/btu430</citation>
+        </citations>
+    </xml>
+</macros>
diff -r b0ba0287b76f -r 1b9140057e5a recetox_aplcms_unsupervised.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/recetox_aplcms_unsupervised.xml	Fri Jun 18 13:53:19 2021 +0000
@@ -0,0 +1,81 @@
+<tool id="recetox_aplcms_unsupervised" name="apLCMS Unsupervised" version="@TOOL_VERSION@+galaxy0">
+    <description>generate a feature table from LC/MS spectra</description>
+    <macros>
+        <import>recetox_aplcms_macros.xml</import>
+    </macros>
+
+    <expand macro="requirements" />
+
+    <command detect_errors="aggressive"><![CDATA[
+        sh ${symlink_inputs} &&
+        Rscript  -e 'source("${__tool_directory__}/main.R")' -e 'source("${run_script}")'
+    ]]></command>
+
+    <configfiles>
+        <configfile name="symlink_inputs">
+            #for $infile in $files
+                ln -s '${infile}' '${infile.element_identifier}'
+            #end for
+        </configfile>
+        <configfile name="run_script"><![CDATA[
+            #set filenames_str = str("', '").join([str($f.element_identifier) for $f in $files])
+
+            unsupervised_main(
+                sample_files = c('$filenames_str'),
+                aligned_file = '${aligned_feature_sample_table}',
+                recovered_file = '${recovered_feature_sample_table}',
+                min_exp = $noise_filtering.min_exp,
+                min_pres = $noise_filtering.min_pres,
+                min_run = $noise_filtering.min_run,
+                mz_tol = $noise_filtering.mz_tol,
+                baseline_correct = $noise_filtering.baseline_correct,
+                baseline_correct_noise_percentile = $noise_filtering.baseline_correct_noise_percentile,
+                intensity_weighted = $noise_filtering.intensity_weighted,
+                shape_model = '$feature_detection.shape_model',
+                BIC_factor = $feature_detection.BIC_factor,
+                peak_estim_method = '$feature_detection.peak_estim_method',
+                min_bandwidth = $feature_detection.min_bandwidth,
+                max_bandwidth = $feature_detection.max_bandwidth,
+                sd_cut = c($feature_detection.sd_cut_min, $feature_detection.sd_cut_max),
+                sigma_ratio_lim = c($feature_detection.sigma_ratio_lim_min, $feature_detection.sigma_ratio_lim_max),
+                component_eliminate = $feature_detection.component_eliminate,
+                moment_power = $feature_detection.moment_power,
+                align_chr_tol = $peak_alignment.align_chr_tol,
+                align_mz_tol = $peak_alignment.align_mz_tol,
+                max_align_mz_diff = $peak_alignment.max_align_mz_diff,
+                recover_mz_range = $weak_signal_recovery.recover_mz_range,
+                recover_chr_range = $weak_signal_recovery.recover_chr_range,
+                use_observed_range = $weak_signal_recovery.use_observed_range,
+                recover_min_count = $weak_signal_recovery.recover_min_count,
+                cluster = as.integer(Sys.getenv('GALAXY_SLOTS', unset = 1))
+            )
+        ]]></configfile>
+    </configfiles>
+
+    <expand macro="inputs">
+        <expand macro="noise_filtering" />
+        <expand macro="feature_detection" />
+        <expand macro="peak_alignment" />
+        <expand macro="weak_signal_recovery" />
+    </expand>
+
+    <outputs>
+        <expand macro="unsupervised_outputs" />
+    </outputs>
+
+    <tests>
+        <test>
+            <param name="files" value="mbr_test0.mzml,mbr_test1.mzml,mbr_test2.mzml" ftype="mzml" />
+            <output name="recovered_feature_sample_table" file="unsupervised_recovered_feature_sample_table.parquet" ftype="parquet" compare="sim_size" delta="1000" />
+        </test>
+    </tests>
+
+    <help>
+        This is the Unsupervised version of apLCMS which is not relying on any existing knowledge about metabolites or
+        any historically detected features. For such functionality please use the Hybrid version of apLCMS.
+
+        @GENERAL_HELP@
+    </help>
+
+    <expand macro="citations" />
+</tool>
diff -r b0ba0287b76f -r 1b9140057e5a test-data/hybrid_recovered_feature_sample_table.parquet
Binary file test-data/hybrid_recovered_feature_sample_table.parquet has changed
diff -r b0ba0287b76f -r 1b9140057e5a test-data/known_table.h5
Binary file test-data/known_table.h5 has changed
diff -r b0ba0287b76f -r 1b9140057e5a test-data/known_table.parquet
Binary file test-data/known_table.parquet has changed
diff -r b0ba0287b76f -r 1b9140057e5a test-data/peaks_hybrid.h5
Binary file test-data/peaks_hybrid.h5 has changed
diff -r b0ba0287b76f -r 1b9140057e5a test-data/peaks_unsupervised.h5
Binary file test-data/peaks_unsupervised.h5 has changed
diff -r b0ba0287b76f -r 1b9140057e5a test-data/unsupervised_recovered_feature_sample_table.parquet
Binary file test-data/unsupervised_recovered_feature_sample_table.parquet has changed