Mercurial > repos > lecorguille > xcms_plot_chromatogram
comparison xcms_plot_chromatogram.r @ 2:e6fdadaf039e draft
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f01148783819c37e474790dbd56619862960448a
author | lecorguille |
---|---|
date | Tue, 03 Apr 2018 11:37:40 -0400 |
parents | fe1f0f16d9e6 |
children | 60b156e58529 |
comparison
equal
deleted
inserted
replaced
1:86ff4497b871 | 2:e6fdadaf039e |
---|---|
26 cat("\n\n") | 26 cat("\n\n") |
27 | 27 |
28 # ----- PROCESSING INFILE ----- | 28 # ----- PROCESSING INFILE ----- |
29 cat("\tARGUMENTS PROCESSING INFO\n") | 29 cat("\tARGUMENTS PROCESSING INFO\n") |
30 | 30 |
31 #saving the specific parameters | |
32 method <- args$method; args$method <- NULL | |
33 | |
34 cat("\n\n") | 31 cat("\n\n") |
35 | 32 |
36 | 33 |
37 # ----- ARGUMENTS PROCESSING ----- | 34 # ----- ARGUMENTS PROCESSING ----- |
38 cat("\tINFILE PROCESSING INFO\n") | 35 cat("\tINFILE PROCESSING INFO\n") |
39 | 36 |
40 #image is an .RData file necessary to use xset variable given by previous tools | 37 mergeXDataReturn <- mergeXData(args) |
41 load(args$image); args$image=NULL | 38 xdata <- mergeXDataReturn$xdata |
42 if (!exists("xdata")) stop("\n\nERROR: The RData doesn't contain any object called 'xdata'. This RData should have been created by an old version of XMCS 2.*") | 39 singlefile <- mergeXDataReturn$singlefile |
43 | 40 md5sumList <- mergeXDataReturn$md5sumList |
44 # Handle infiles | 41 sampleNamesList <- mergeXDataReturn$sampleNamesList |
45 if (!exists("singlefile")) singlefile <- NULL | |
46 if (!exists("zipfile")) zipfile <- NULL | |
47 rawFilePath <- getRawfilePathFromArguments(singlefile, zipfile, args) | |
48 zipfile <- rawFilePath$zipfile | |
49 singlefile <- rawFilePath$singlefile | |
50 args <- rawFilePath$args | |
51 directory <- retrieveRawfileInTheWorkingDirectory(singlefile, zipfile) | |
52 | |
53 | 42 |
54 cat("\n\n") | 43 cat("\n\n") |
55 | 44 |
56 | 45 |
57 # ----- MAIN PROCESSING INFO ----- | 46 # ----- MAIN PROCESSING INFO ----- |
58 cat("\tMAIN PROCESSING INFO\n") | 47 cat("\tMAIN PROCESSING INFO\n") |
59 | 48 |
60 | 49 |
61 cat("\t\tDRAW GRAPHICS\n") | 50 cat("\t\tDRAW GRAPHICS\n") |
62 | 51 |
63 #@TODO: one day, use xdata instead of xset to draw the TICs and BPC or a complete other method | |
64 getPlotTICs(xdata, pdfname="TICs.pdf") | 52 getPlotTICs(xdata, pdfname="TICs.pdf") |
65 getPlotBPIs(xdata, pdfname="BPIs.pdf") | 53 getPlotBPIs(xdata, pdfname="BPIs.pdf") |
66 | 54 |
67 cat("\n\n") | 55 cat("\n\n") |
68 | 56 |