Network analysis
We used the qgraph package in R (v3.6.1) to perform network analysis. A network is a set of nodes and a set of edges, where each node represents either an immune parameter or an operational taxonomic unit (OTU) from microbiome analysis whereas the edges represent associations amongst them. Pearson correlation coefficients were used to quantify the strength of associations between combinations of immune parameters or OTUs.