We first normalized (mean centered and scaled) the imputed dataset. We then fit linear models for each cytokine according to the following expression: $$ Cytokine Expression = Diagnosis + Sex + Age + ...
age.M <- with(dat_M, seq.int(min(sqrt_age), max(sqrt_age), length.out = 550)) age.F <- with(dat_F, seq.int(min(sqrt_age), max(sqrt_age), length.out = 550)) par(mfrow ...