```r
physeq <- taxa_level(physeq, "Phylum")
env.taxa.cor <- taxa.env.correlation(physeq, grouping_column = "Country", method = "pearson", 
    pvalue.threshold = 0.05, padjust.method = "BH", adjustment = 5, num.taxa = 50, 
    select.variables = NULL)
```