Create a synthetic data set for testing and demo purposes, including the code to generate it. A simple approach would be to estimate the covariance matrix of the residuals of the 145 ROIs and then add the effects of age, sex, ICV from the ComBatGAM model. In addition, we could also add an effect of disease.