Hi, I am running Maaslin2 in RStudio with microbial relative abundances.
However, when I run different comparisons I get no errors, but no significant results or figures are represented. I don’t know if I am doing something wrong. Please find attached the metadata and sample data. Thanks in advance.
genus_table_rel_OS_MMA.tsv (317.3 KB)
sample_metadata_MMA.tsv (2.7 KB)
input_data = system.file(“extdata”, “genus_table_rel_OS_MMA.tsv”, package=“Maaslin2”) # The abundance table file input_data
input_metadata = system.file(“extdata”, “sample_metadata_MMA.tsv”, package=“Maaslin2”) # The metadata table file input_metadata
df_input_data = read.table(file = input_data,
header = TRUE,
sep = “\t”,
row.names = 1,
stringsAsFactors = FALSE)
df_input_data[1:5, 1:5]
df_input_metadata = read.table(file = input_metadata,
header = TRUE,
sep = “\t”,
row.names = 1,
stringsAsFactors = FALSE)
df_input_metadata[1:5, ]
fit_data = Maaslin2(input_data = input_data,
input_metadata = input_metadata,
min_prevalence = 0,
normalization = “NONE”,
output = “OS_MMA_retro”,
fixed_effects = c(“OS”),
reference = c(“OS,LTS”))
fit_data2 = Maaslin2(input_data = input_data,
input_metadata = input_metadata,
normalization = “NONE”,
output = “OS_MMA_retro2”,
fixed_effects = c(“OS”),
reference = c(“OS,STS”),
min_prevalence = 0,
min_abundance = 0.00000000001)