Mercurial > repos > recetox > recetox_aplcms_recover_weaker_signals

annotate utils.R @ 9:a9bb2ccc53de draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
author recetox
date Wed, 11 Oct 2023 11:18:27 +0000
parents 99d118321d5f
children e0536ff73f36
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
1 library(recetox.aplcms)
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
2
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
3 get_env_sample_name <- function() {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
4 sample_name <- Sys.getenv("SAMPLE_NAME", unset = NA)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
5 if (nchar(sample_name) == 0) {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
6 sample_name <- NA
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
7 }
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
8 if (is.na(sample_name)) {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
9 message("The mzML file does not contain run ID.")
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
10 }
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
11 return(sample_name)
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
12 }
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
13
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
14 save_sample_name <- function(df, sample_name) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
15 attr(df, "sample_name") <- sample_name
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
16 return(df)
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
17 }
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
18
8
99d118321d5f planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 0873c036d734782d2647b0eea84cd1a0658f5e58
recetox
parents: 6
diff changeset
19 restore_sample_name <- function(df) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
20 return(df$sample_id[1])
8
99d118321d5f planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 0873c036d734782d2647b0eea84cd1a0658f5e58
recetox
parents: 6
diff changeset
21 }
99d118321d5f planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 0873c036d734782d2647b0eea84cd1a0658f5e58
recetox
parents: 6
diff changeset
22
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
23 load_sample_name <- function(df) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
24 sample_name <- attr(df, "sample_name")
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
25 if (is.null(sample_name)) {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
26 return(NA)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
27 } else {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
28 return(sample_name)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
29 }
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
30 }
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
31
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
32 save_data_as_parquet_file <- function(data, filename) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
33 arrow::write_parquet(data, filename)
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
34 }
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
35
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
36 load_data_from_parquet_file <- function(filename) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
37 return(arrow::read_parquet(filename))
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
38 }
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
39
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
40 load_parquet_collection <- function(files) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
41 features <- lapply(files, arrow::read_parquet)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
42 features <- lapply(features, tibble::as_tibble)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
43 return(features)
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
44 }
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
45
4
755ff0afaaff planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 724ecb1b81ebd8a67488b8a9397177b2ff0357db
recetox
parents: 3
diff changeset
46 save_parquet_collection <- function(feature_tables, sample_names, subdir) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
47 dir.create(subdir)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
48 for (i in seq_len(length(feature_tables))) {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
49 filename <- file.path(subdir, paste0(sample_names[i], ".parquet"))
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
50 feature_table <- as.data.frame(feature_tables[[i]])
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
51 feature_table <- save_sample_name(feature_table, sample_names[i])
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
52 arrow::write_parquet(feature_table, filename)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
53 }
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
54 }
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
55
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
56 sort_by_sample_name <- function(tables, sample_names) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
57 return(tables[order(sample_names)])
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
58 }
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
59
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
60 save_tolerances <- function(table, tol_file) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
61 mz_tolerance <- c(table$mz_tol_relative)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
62 rt_tolerance <- c(table$rt_tol_relative)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
63 arrow::write_parquet(data.frame(mz_tolerance, rt_tolerance), tol_file)
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
64 }
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
65
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
66 save_aligned_features <- function(aligned_features, metadata_file, rt_file, intensity_file) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
67 save_data_as_parquet_file(aligned_features$metadata, metadata_file)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
68 save_data_as_parquet_file(aligned_features$rt, rt_file)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
69 save_data_as_parquet_file(aligned_features$intensity, intensity_file)
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
70 }
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
71
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
72 select_table_with_sample_name <- function(tables, sample_name) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
73 sample_names <- lapply(tables, load_sample_name)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
74 index <- which(sample_names == sample_name)
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
75 if (length(index) > 0) {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
76 return(tables[[index]])
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
77 } else {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
78 stop(sprintf(
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
79 "Mismatch - sample name '%s' not present in %s",
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
80 sample_name, paste(sample_names, collapse = ", ")
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
81 ))
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
82 }
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
83 }
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
84
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
85 select_adjusted <- function(recovered_features) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
86 return(recovered_features$adjusted_features)
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
87 }
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
88
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
89 known_table_columns <- function() {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
90 c(
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
91 "chemical_formula", "HMDB_ID", "KEGG_compound_ID", "mass", "ion.type",
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
92 "m.z", "Number_profiles_processed", "Percent_found", "mz_min", "mz_max",
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
93 "RT_mean", "RT_sd", "RT_min", "RT_max", "int_mean(log)", "int_sd(log)",
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
94 "int_min(log)", "int_max(log)"
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
95 )
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
96 }
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
97
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
98 save_known_table <- function(table, filename) {
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
99 columns <- known_table_columns()
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
100 arrow::write_parquet(table$known_table[columns], filename)
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
101 }
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
102
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
103 read_known_table <- function(filename) {
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
104 arrow::read_parquet(filename, col_select = known_table_columns())
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
105 }
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
106
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
107 save_pairing <- function(table, filename) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
108 df <- table$pairing %>%
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
109 as_tibble() %>%
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
110 setNames(c("new", "old"))
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
111 arrow::write_parquet(df, filename)
0
067a308223e3 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 19de0924a65bc65cbbf7c1fc17e9b5348305f95c
recetox
parents:
diff changeset
112 }
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
113
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
114 join_tables_to_list <- function(metadata, rt_table, intensity_table) {
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
115 features <- new("list")
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
116 features$metadata <- metadata
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
117 features$intensity <- intensity_table
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
118 features$rt <- rt_table
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
119 return(features)
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
120 }
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
121
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
122 validate_sample_names <- function(sample_names) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
123 if ((any(is.na(sample_names))) || (length(unique(sample_names)) != length(sample_names))) {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
124 stop(sprintf(
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
125 "Sample names absent or not unique - provided sample names: %s",
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
126 paste(sample_names, collapse = ", ")
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
127 ))
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
128 }
2
472dc85ce7c5 planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 506df2aef355b3791567283e1a175914f06b405a
recetox
parents: 0
diff changeset
129 }
3
c69a12bfc2fb planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 4254c6f6f05915ceab5af5d7040ac2c44a8364b4
recetox
parents: 2
diff changeset
130
c69a12bfc2fb planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 4254c6f6f05915ceab5af5d7040ac2c44a8364b4
recetox
parents: 2
diff changeset
131 determine_sigma_ratios <- function(sigma_ratio_lim_min = NA, sigma_ratio_lim_max = NA) {
9
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
132 if (is.na(sigma_ratio_lim_min)) {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
133 sigma_ratio_lim_min <- 0
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
134 }
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
135 if (is.na(sigma_ratio_lim_max)) {
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
136 sigma_ratio_lim_max <- Inf
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
137 }
a9bb2ccc53de planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 97249a1af94ac5c387e1ede274dec5364f71cde9
recetox
parents: 8
diff changeset
138 return(c(sigma_ratio_lim_min, sigma_ratio_lim_max))
3
c69a12bfc2fb planemo upload for repository https://github.com/RECETOX/galaxytools/tree/master/tools/recetox_aplcms commit 4254c6f6f05915ceab5af5d7040ac2c44a8364b4
recetox
parents: 2
diff changeset
139 }