Example code for running PLS analysis of relative cortical thickness change and gene expression data
For details of analysis see:
Ball, Seidlitz, Beare & Seal, biorxiv 2019
Cortical remodelling in childhood is associated with genes enriched for neurodevelopmental disorders
https://www.biorxiv.org/content/10.1101/707042v1
This code performs the main PLS analysis with 3 components. Enrichment analyses can be performed using the gene component loadings with Webgestalt (http://www.webgestalt.org/) or your tool of choice.
Code was run using:
python 3.6.8
numpy 1.15.4
pandas 0.23.4
sklearn 0.20.0
Instructions:
python run_PLS_analysis.py