1  library(tidyverse)
 2  
 3  enname <- function(x) {
 4    names(x) <- x
 5    x
 6  }
 7  
 8  sim <- read_csv("analysis/features_normalized/mrxcat_simulation/features.csv") %>%
 9    separate(file, into=c("noise", "seed"), sep="_") %>%
10    mutate(patient="sim")
11  
12  ranks_sim <- read_csv("analysis/tables/mrxcat_simulation/ranks.csv") %>% select(-filter, -feature_class, -feature) %>% mutate(patient="sim")
13  
14  acdc <- dir("analysis/features_normalized/ACDC/", full.names = TRUE) %>% map_df(read_csv) %>%
15    separate(file, into=c("patient", "noise", "seed"), sep="_")
16  
17  ranks_acdc <- dir("analysis/tables/ACDC", full.names = TRUE) %>%
18    enname() %>%
19    map_df(read_csv, .id="patient") %>%
20    mutate(patient = str_remove(patient, ".*ranks_"), patient=str_remove(patient, ".csv"))
21  
22  bae <- dir("analysis/features_normalized/subject/", full.names = TRUE) %>% map_df(read_csv) %>%
23    separate(file, into=c("patient", "noise", "seed"), sep="_")
24  
25  ranks_bae <- dir("analysis/tables/subject", full.names = TRUE) %>%
26    enname() %>%
27    map_df(read_csv, .id="patient") %>%
28    mutate(patient = str_remove(patient, ".*ranks_"), patient=str_remove(patient, ".csv"))
29  
30  combined <- bind_rows(sim, acdc, bae)
31  ranks_combined <- bind_rows(ranks_sim, select(ranks_acdc, colnames(ranks_sim)), select(ranks_bae, colnames(ranks_sim)))
32  
33  combined %>% write_tsv("analysis/combined/features.tsv")
34  ranks_combined %>% write_tsv("analysis/combined/ranks.tsv")