Repository 'msstats'
hg clone https://toolshed.g2.bx.psu.edu/repos/galaxyp/msstats

Changeset 5:28434abe6c5c (2021-08-06)
Previous changeset 4:593839e1f2c3 (2021-02-25) Next changeset 6:b7034eff0db1 (2023-01-31)
Commit message:
"planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/msstats commit 0e253b8114e2fe6f4e33edcc5f1a4602073064c3"
modified:
msstats.xml
test-data/Comparison_plot_skyline.pdf
test-data/MSstats ProfilePlot.pdf
test-data/Profile_plot_skyline.pdf
test-data/Volcano_plot_skyline.pdf
test-data/profile_wsum_plot.pdf
added:
test-data/featurelevel_data_skyline.tabular
b
diff -r 593839e1f2c3 -r 28434abe6c5c msstats.xml
--- a/msstats.xml Thu Feb 25 08:41:37 2021 +0000
+++ b/msstats.xml Fri Aug 06 20:06:19 2021 +0000
[
b'@@ -1,7 +1,7 @@\n-<tool id="msstats" name="MSstats" version="@VERSION@.1">\n+<tool id="msstats" name="MSstats" version="@VERSION@.0">\n     <description>statistical relative protein significance analysis in DDA, SRM and DIA Mass Spectrometry</description>\n     <macros>\n-        <token name="@VERSION@">3.22.0</token>\n+        <token name="@VERSION@">4.0.0</token>\n         <xml name="useUniquePeptide">\n             <param name="useUniquePeptide" type="boolean" truevalue="TRUE" falsevalue="FALSE" checked="true" label="Remove peptides that are assigned for more than one proteins"/>\n         </xml>\n@@ -28,7 +28,6 @@\n     <command detect_errors="exit_code"><![CDATA[\n         cat \'$msstats_script\' > \'$r_script\' &&\n         Rscript \'$msstats_script\'\n-        && cat msstats*.log > \'$log\'\n     ]]></command>\n     <configfiles>\n         <configfile name="msstats_script"><![CDATA[\n@@ -66,7 +65,10 @@\n                            fewMeasurements="$input.input_options.fewMeasurements", \n                            removeMpeptides=$input.input_options.removeMpeptides,\n                            removeOxidationMpeptides=$input.input_options.removeOxidationMpeptides,\n-                           removeProtein_with1Peptide=$input.input_options.removeProtein_with1Peptide)\n+                           removeProtein_with1Peptide=$input.input_options.removeProtein_with1Peptide,\n+                           use_log_file = TRUE,\n+                           append = TRUE,\n+                           log_file_pat = "log.txt")\n \n #elif $input.input_src == \'OpenMS\'\n \n@@ -91,7 +93,10 @@\n                              useUniquePeptide=$input.input_options.useUniquePeptide, \n                              summaryforMultipleRows=$input.input_options.summaryforMultipleRows, \n                              fewMeasurements="$input.input_options.fewMeasurements", \n-                             removeProtein_with1Feature=$input.input_options.removeProtein_with1Feature)\n+                             removeProtein_with1Feature=$input.input_options.removeProtein_with1Feature,\n+                             use_log_file = TRUE,\n+                             append = TRUE,\n+                             log_file_pat = "log.txt")\n \n \n #elif $input.input_src == \'OpenSWATH\'\n@@ -114,8 +119,11 @@\n                                 useUniquePeptide=$input.input_options.useUniquePeptide,\n                                 fewMeasurements="$input.input_options.fewMeasurements",\n                                 removeProtein_with1Feature=$input.input_options.removeProtein_with1Feature,\n-                                summaryforMultipleRows=$input.input_options.summaryforMultipleRows)\n-\n+                                summaryforMultipleRows=$input.input_options.summaryforMultipleRows,\n+                                use_log_file = TRUE,\n+                                append = TRUE,\n+                                log_file_pat = "log.txt")\n+\t\t                \n #elif $input.input_src == \'Skyline\'\n \n   #if $input.skyline_input.is_of_type(\'csv\')\n@@ -142,7 +150,10 @@\n \t\t\t\tuseUniquePeptide = $input.input_options.useUniquePeptide,\n \t\t\t\tfewMeasurements="$input.input_options.fewMeasurements",\n \t\t\t\tremoveOxidationMpeptides = $input.input_options.removeOxidationMpeptides,\n-\t\t\t\tremoveProtein_with1Feature = $input.input_options.removeProtein_with1Feature)\n+\t\t\t\tremoveProtein_with1Feature = $input.input_options.removeProtein_with1Feature,\n+\t\t\t\tuse_log_file = TRUE,\n+\t\t\t\tappend = TRUE,\n+\t\t                log_file_pat = "log.txt")\n \n #end if\n \n@@ -152,7 +163,6 @@\n                           #if $dp_options.norm.normalization == \'globalStandards\'\n                           nameStandards=c($dp_options.norm.nameStandards),\n                           #end if\n-                          fillIncompleteRows=$dp_options.fillIncompleteRows,\n                           featureSubset="$dp_options.features.featureSubset",\n                           #if $dp_options.features.featureSubset == \'topN\'\n                       '..b" are censored.\n             - 0 - It assumes that all values between 0 and 1 in Intensity column are censored. If there areNAs inIntensitywith this option, NAs will be considered as random missing.\n-            - NULL - It assumes that all missing values are randomly missing.\n             - Skyline and OpenSWATH input should use '0'. MaxQuant input should use 'NA'\n-        - Cutoff value for censoring: cutoff for AFT model; only with censored intensity 'NA' or '0'; if NULL it assumes that there is no censored missing and any imputation will not be performed. In case that there are completely missing measurements in a run for a protein, any imputation will not be performed. In addition, the condition, which has no measurement at all in a protein, will be not impute.\n-\n-            - minimum value for each feature: cutoff for AFT model will be the minimum value for each feature across runs. With this option, those runs with substantial missing measurements will be biased by the cutoff value. In such case, you may remove the runs that have more than 50% missing values from the analysis.\n-            - minimum value for each run: cutoff for AFT model will be the minimum value for each run across features\n-            - smallest between minimum value of corresponding feature and minimum value of corresponding run: cutoff for AFT model will be the smallest value between minimum valueof corresponding feature and minimum value of corresponding run\n-        - Maximum quantile for deciding censored missing values: If you don\xe2\x80\x99t want to apply the threshold of noise intensity in your data, you can use maxQuantileforCensored=NULL.\n-    - Missing value imputation combination with summarization method TMP:\n-\n-        - Summarization method: TMP + censored intensity: 'NULL': It assumes that all intensities are missing at random, therefore no action with missing value imputation: No; or error with missing value imputation: Yes.\n-        - Missing value imputation: Yes + censored intensity:'NA' or '0': AFT model-based imputation using cutoff value for censoring in the AFT model\n-        - Missing value imputation: No + censored intensity:'NA' or '0': censored intensities will be replaced with the value specified  in cutoff value for censoring\n-        - Missing value imputation: No + censored intensity: NULL: no imputation\n \n - Group comparison: automatic detection of differentially abundant proteins between two conditions, conditions have to be specified with the 'comparison matrix'\n - Quantification per sample or group: choose the corresponding output option\n@@ -1115,13 +1154,13 @@\n     - MSstats log - check log file for warnings and information on the analysis steps (txt)\n     - MSstats Rscript - can be used to re-run analysis outside Galaxy or to inspect the executed code (txt)\n     - MSstats RawData - raw files combined into MSstats format (tabular)\n-    - MSstats ProcessedData - transformed, normalized, imputed intensities (tabular)\n+    - MSstats FeatureLevelData - transformed, normalized, imputed intensities (tabular)\n \n         - Intensity column:  includes original intensities values\n         - Abundance column:  contains the log2 transformed and normalized intensities and it will used for run-level summarization\n         - Censored column:  has the decision about censored missing or not, based on censored Intensity and maximum quantile for deciding censored missing values options. Abundances with TRUE value in censored column will be considered as censored missing and imputed when Missing value imputation: Yes.\n \n-    - MSstats RunlevelData - run and protein level summarized data (tabular)\n+    - MSstats ProteinLevelData - run and protein level summarized data (tabular)\n \n         - LogIntensities: log intensity summarized per run and protein, they will be used for the group comparison and summarized profile plot\n         - NumMeasuredFeature: shows how many features were used for summarization of the corresponding run and protein\n"
b
diff -r 593839e1f2c3 -r 28434abe6c5c test-data/Comparison_plot_skyline.pdf
b
Binary file test-data/Comparison_plot_skyline.pdf has changed
b
diff -r 593839e1f2c3 -r 28434abe6c5c test-data/MSstats ProfilePlot.pdf
b
Binary file test-data/MSstats ProfilePlot.pdf has changed
b
diff -r 593839e1f2c3 -r 28434abe6c5c test-data/Profile_plot_skyline.pdf
b
Binary file test-data/Profile_plot_skyline.pdf has changed
b
diff -r 593839e1f2c3 -r 28434abe6c5c test-data/Volcano_plot_skyline.pdf
b
Binary file test-data/Volcano_plot_skyline.pdf has changed
b
diff -r 593839e1f2c3 -r 28434abe6c5c test-data/featurelevel_data_skyline.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/featurelevel_data_skyline.tabular Fri Aug 06 20:06:19 2021 +0000
[
@@ -0,0 +1,46 @@
+PROTEIN PEPTIDE TRANSITION FEATURE LABEL GROUP RUN SUBJECT FRACTION originalRUN censored INTENSITY ABUNDANCE newABUNDANCE
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition1 1 1 1 121219_S_CCES_01_01_LysC_Try_1to10_Mixt_1_1raw FALSE 319070944 29.3098413927068 29.3098413927068
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition1 1 1 1 121219_S_CCES_01_01_LysC_Try_1to10_Mixt_1_1raw FALSE 114949068 27.8369589250684 27.8369589250684
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition1 1 1 1 121219_S_CCES_01_01_LysC_Try_1to10_Mixt_1_1raw FALSE 132974798 28.0471170016311 28.0471170016311
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition1 2 1 1 121219_S_CCES_01_02_LysC_Try_1to10_Mixt_1_2raw FALSE 301443168 28.0471170016311 28.0471170016311
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition1 2 1 1 121219_S_CCES_01_02_LysC_Try_1to10_Mixt_1_2raw FALSE 864180704 29.5665639874012 29.5665639874012
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition1 2 1 1 121219_S_CCES_01_02_LysC_Try_1to10_Mixt_1_2raw FALSE 99948868 26.4544931020782 26.4544931020782
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition1 3 1 1 121219_S_CCES_01_03_LysC_Try_1to10_Mixt_1_3raw FALSE 275819008 28.0471170016311 28.0471170016311
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition1 3 1 1 121219_S_CCES_01_03_LysC_Try_1to10_Mixt_1_3raw FALSE 866460512 29.6985291199933 29.6985291199933
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition1 3 1 1 121219_S_CCES_01_03_LysC_Try_1to10_Mixt_1_3raw FALSE 104015372 26.6401918724078 26.6401918724078
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition2 4 2 1 121219_S_CCES_01_04_LysC_Try_1to10_Mixt_2_1raw FALSE 282914080 28.0471170016311 28.0471170016311
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition2 4 2 1 121219_S_CCES_01_04_LysC_Try_1to10_Mixt_2_1raw FALSE 887988992 29.6972948156548 29.6972948156548
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition2 4 2 1 121219_S_CCES_01_04_LysC_Try_1to10_Mixt_2_1raw FALSE 57196342 25.7407478107147 25.7407478107147
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition2 5 2 1 121219_S_CCES_01_05_LysC_Try_1to10_Mixt_2_2raw FALSE 306392288 28.0471170016311 28.0471170016311
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition2 5 2 1 121219_S_CCES_01_05_LysC_Try_1to10_Mixt_2_2raw FALSE 562341792 28.9231842858172 28.9231842858172
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition2 5 2 1 121219_S_CCES_01_05_LysC_Try_1to10_Mixt_2_2raw FALSE 61787154 25.7371158446414 25.7371158446414
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition2 6 2 1 121219_S_CCES_01_06_LysC_Try_1to10_Mixt_2_3raw FALSE 311995872 28.0471170016311 28.0471170016311
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition2 6 2 1 121219_S_CCES_01_06_LysC_Try_1to10_Mixt_2_3raw FALSE 947982416 29.6504503596479 29.6504503596479
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition2 6 2 1 121219_S_CCES_01_06_LysC_Try_1to10_Mixt_2_3raw FALSE 57000615 25.5946394508612 25.5946394508612
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition3 7 3 1 121219_S_CCES_01_07_LysC_Try_1to10_Mixt_3_1raw FALSE 317509760 28.0471170016311 28.0471170016311
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition3 7 3 1 121219_S_CCES_01_07_LysC_Try_1to10_Mixt_3_1raw FALSE 831511536 29.4360523406791 29.4360523406791
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition3 7 3 1 121219_S_CCES_01_07_LysC_Try_1to10_Mixt_3_1raw FALSE 110326472 26.5220950579513 26.5220950579513
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition3 8 3 1 121219_S_CCES_01_08_LysC_Try_1to10_Mixt_3_2raw FALSE 253089648 28.0471170016311 28.0471170016311
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition3 8 3 1 121219_S_CCES_01_08_LysC_Try_1to10_Mixt_3_2raw FALSE 909277584 29.8921892900012 29.8921892900012
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition3 8 3 1 121219_S_CCES_01_08_LysC_Try_1to10_Mixt_3_2raw FALSE 124922268 27.0284991700692 27.0284991700692
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition3 9 3 1 121219_S_CCES_01_09_LysC_Try_1to10_Mixt_3_3raw FALSE 330089504 28.0471170016311 28.0471170016311
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition3 9 3 1 121219_S_CCES_01_09_LysC_Try_1to10_Mixt_3_3raw FALSE 735864592 29.2037000536836 29.2037000536836
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition3 9 3 1 121219_S_CCES_01_09_LysC_Try_1to10_Mixt_3_3raw FALSE 54330676 25.4440986390611 25.4440986390611
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition4 10 4 1 121219_S_CCES_01_10_LysC_Try_1to10_Mixt_4_1raw FALSE NA NA NA
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition4 10 4 1 121219_S_CCES_01_10_LysC_Try_1to10_Mixt_4_1raw FALSE 943487632 29.7700692177507 29.7700692177507
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition4 10 4 1 121219_S_CCES_01_10_LysC_Try_1to10_Mixt_4_1raw FALSE 86579608 26.3241647855115 26.3241647855115
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition4 11 4 1 121219_S_CCES_01_11_LysC_Try_1to10_Mixt_4_2raw FALSE NA NA NA
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition4 11 4 1 121219_S_CCES_01_11_LysC_Try_1to10_Mixt_4_2raw FALSE 600179488 29.4051505227937 29.4051505227937
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition4 11 4 1 121219_S_CCES_01_11_LysC_Try_1to10_Mixt_4_2raw FALSE 91340524 26.6890834804686 26.6890834804686
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition4 12 4 1 121219_S_CCES_01_12_LysC_Try_1to10_Mixt_4_3raw FALSE NA NA NA
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition4 12 4 1 121219_S_CCES_01_12_LysC_Try_1to10_Mixt_4_3raw FALSE 709016080 29.5578987363561 29.5578987363561
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition4 12 4 1 121219_S_CCES_01_12_LysC_Try_1to10_Mixt_4_3raw FALSE 87312183 26.5363352669061 26.5363352669061
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition5 13 5 1 121219_S_CCES_01_13_LysC_Try_1to10_Mixt_5_1raw FALSE NA NA NA
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition5 13 5 1 121219_S_CCES_01_13_LysC_Try_1to10_Mixt_5_1raw FALSE NA NA NA
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition5 13 5 1 121219_S_CCES_01_13_LysC_Try_1to10_Mixt_5_1raw FALSE 91599334 28.0471170016311 28.0471170016311
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition5 14 5 1 121219_S_CCES_01_14_LysC_Try_1to10_Mixt_5_2raw FALSE NA NA NA
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition5 14 5 1 121219_S_CCES_01_14_LysC_Try_1to10_Mixt_5_2raw FALSE NA NA NA
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition5 14 5 1 121219_S_CCES_01_14_LysC_Try_1to10_Mixt_5_2raw FALSE 96063704 28.0471170016311 28.0471170016311
+P32125 ADVGFLC[+57]NMLER_2 NA_NA ADVGFLC[+57]NMLER_2_NA_NA L Condition5 15 5 1 121219_S_CCES_01_15_LysC_Try_1to10_Mixt_5_3raw FALSE NA NA NA
+P32125 AFAEAMANNSFNADEK_2 NA_NA AFAEAMANNSFNADEK_2_NA_NA L Condition5 15 5 1 121219_S_CCES_01_15_LysC_Try_1to10_Mixt_5_3raw FALSE 836397712 29.6396138744808 29.6396138744808
+P32125 AGAAQTIVASQQR_2 NA_NA AGAAQTIVASQQR_2_NA_NA L Condition5 15 5 1 121219_S_CCES_01_15_LysC_Try_1to10_Mixt_5_3raw FALSE 91967458 26.4546201287815 26.4546201287815
b
diff -r 593839e1f2c3 -r 28434abe6c5c test-data/profile_wsum_plot.pdf
b
Binary file test-data/profile_wsum_plot.pdf has changed