scripts/collect_n_alignments.R

In richelbilderbeek/wiritttea: WIRITTTE analysis

# Collects the number of species trees used in all .RDa files in one file
library(wiritttea)
options(warn = 2)

args <- commandArgs(trailingOnly = TRUE)
if (length(args) == 0) {
  stop("Supply a source folder as a first argument, e.g. '~/wirittte_data/stub'")
}
if (length(args) == 1) {
  stop("Supply a target filename as a second argument, e.g. '~/n_alignments_stub.csv'")
}
if (length(args) != 2) {
  stop("Supply two parameters: a source folder and a target filename, ",
    "e.g. '~/wirittte_data/stub ~/parameters_stub.csv'")
}

path_data <- args[1]
n_alignments_filename <- args[2]

args <- commandArgs(trailingOnly = TRUE)
if (length(args) > 0) path_data <- args[1]
if (length(args) > 1) n_alignments_filename <- args[2]

print(paste("path_data:", path_data))
print(paste("n_alignments_filename:", n_alignments_filename))

print("Collecting .RDa files")
my_filenames <- list.files(path_data, pattern = "*.RDa", full.names = TRUE)

print("Collecting # alignments")
df_n_alignments <- wiritttea::collect_files_n_alignments(filenames = my_filenames)

print("Saving # alignments")
write.csv(df_n_alignments, n_alignments_filename)