A paper titled 'Integrative network analysis of nineteen brain regions identifies molecular signatures and networks underlying selective regional vulnerability to Alzheimer?s disease' used the MSBB microarray data and reported over 6,000 DEGs across all the brain regions after taking into account all the different covariates. I've been trying to recreate their results in R but I'm getting completely different results. My model matrix looks something like this: model.matrix(~0+alzheimers+pH+Sex+Race+PMI). I was wondering if anyone had the script used by the team who published the paper or if someone had successfully managed to recreate the results? Thanks, Roberto