Is it surprising that the same dataset can produce this ordination
and yet Maaslin3 would find 0 ASVs differing between these groups?
Balmoral.covariates.csv (345 Bytes)
Balmoral.features.csv (24.2 KB)
test1 ← read.csv(“Balmoral.features.csv”, header = TRUE, row.names = 1)
test2 ← read.csv(“Balmoral.covariates.csv”, header = TRUE, row.names = 1)
test2$Time_Point ← factor(test2$Time_Point)
set.seed(1024)
Balmoral.fit_out ← maaslin3(input_data = test1, input_metadata = test2,
output = ‘repro_output’,
formula = ‘~ Time_Point + Read_depth’,
normalization = ‘TSS’,
transform = ‘LOG’,
augment = TRUE,
standardize = TRUE,
max_significance = 0.1,
median_comparison_abundance = TRUE,
median_comparison_prevalence = FALSE,
max_pngs = 250,
cores = 1)
All associations had errors
or were insignificant.
