```r
co_occr <- co_occurence_network(physeq, grouping_column = "Country", rhos = 0.35, 
    method = "cor", qval_threshold = 0.05, select.condition = "V", scale.vertex.size = 4, 
    scale.edge.width = 15, plotNetwork = T, plotBetweennessEeigenvalue = F)
```