5 Core taxa multivariate ANOVA


library(kableExtra)
library(tidyverse)
library(magrittr)

ps_prune_sample <- phyloseq::subset_samples(ps_raw, phyloseq::sample_sums(ps_raw) >15000)
ps_prune_taxa <- phyloseq::prune_taxa(phyloseq::taxa_sums(ps_prune_sample) > 0, ps_prune_sample)
ps_core <- microbiome::core(ps_prune_taxa, detection  = 0.01, prevalence = 0.01)

ps_core %>% 
betatest(group = c("nationality", "bmi_group"), distance = "bray") %>%
  set_colnames(c("Group", "Df", "SumsOfSqs", "MeanSqs", "FModel", "R2", "Pvalue")) %>% 
  kable(format = "html") %>%
  kable_styling(bootstrap_options = "basic", full_width = F, position = "float_left")  %>%
  column_spec(7, bold = TRUE, color = "black", background = "#eeeeee")
Group Df SumsOfSqs MeanSqs FModel R2 Pvalue
bmi_group 2 0.3580195 0.1790097 1.405949 0.0370904 0.196
bmi_group 73 9.2945867 0.1273231 NA 0.9629096 NA
bmi_group 75 9.6526062 NA NA 1.0000000 NA
nationality 1 1.4483886 1.4483886 13.064104 0.1500516 0.001
nationality 74 8.2042176 0.1108678 NA 0.8499484 NA
nationality 75 9.6526062 NA NA 1.0000000 NA