@@ -117,19 +117,19 @@ library(vegan)
117117** Load species abundance, taxonomical annotation and metadata file**
118118
119119``` {r}
120- TRAINING_DIR <- "/workspaces/dsp_metagenomics_training/"
120+ training_dir <- "/workspaces/dsp_metagenomics_training/"
121121```
122122
123123``` {r}
124124# Abundance table
125- abundance <- readRDS(file = paste0(TRAINING_DIR , "data/MetaphlanAbundance_Species.rds"))
125+ abundance <- readRDS(file = paste0(training_dir , "data/MetaphlanAbundance_Species.rds"))
126126rownames(abundance) <- gsub("_SRR_db1.metaphlan", "", rownames(abundance))
127127
128128# Taxonomical annotation
129- annotation <- readRDS(file = paste0(TRAINING_DIR , "data/MetaphlanAnnotations_Species.rds"))
129+ annotation <- readRDS(file = paste0(training_dir , "data/MetaphlanAnnotations_Species.rds"))
130130
131131# Metadata
132- metadata <- read.table(file = paste0(TRAINING_DIR , "data/metadata.tsv"), header = TRUE,
132+ metadata <- read.table(file = paste0(training_dir , "data/metadata.tsv"), header = TRUE,
133133 sep = "\t",
134134 quote = "",
135135 row.names = NULL)
@@ -141,7 +141,7 @@ abundance <- abundance[rownames(metadata), ]
141141
142142# Creating a directory for the PCoA related results
143143``` {r}
144- results_dir <- paste0(TRAINING_DIR , "results/report/02_Beta-diversity/")
144+ results_dir <- paste0(training_dir , "results/report/02_Beta-diversity/")
145145
146146dir.create(results_dir, recursive = TRUE, showWarnings = FALSE)
147147```
0 commit comments