Table S1. Plasma cytokines measured at each visit by multiplex bead-based immunoassay. Table S2. Parameters of HRV (RMSSD, pNN50, HF power, and the LF/HF ratio) at baseline and follow-up. HRV parameters are compared between visits by paired analyses (Wilcoxon signed-rank test or paired t test); p values are listed. Table S3. Matrix of Spearmanâs rank correlations (R) between cytokines included in the multivariate models. p values <0.1 are included, and R values that are statistically significant (pâ<â0.05) are shown in bold type. Table S4. Associations of disease activity and flare with plasma cytokines at baseline by univariate linear regression (p values and beta regression coefficients). (DOCX 23 kb