The authors report a low number of overlapping DE gene features between the sgACC and the amygdala, but at the same time an absence of strong evidence for an interaction effect between phenotype and brain region. To investigate this pattern further, the authors should correlate fold changes (between cases and controls) across all gene features between the two brain regions.
[x] Create scatterplots of DE t-stats between sACC and Amyd at the gene, exon, exon-exon and transcript levels.
Likely requires using /dcl01/lieber/ajaffe/lab/zandiHyde_bipolar_rnaseq/case_control/bipolarControl_deStats_byRegion_qSVAjoint_withAnnotation.rda.
Likely requires using
/dcl01/lieber/ajaffe/lab/zandiHyde_bipolar_rnaseq/case_control/bipolarControl_deStats_byRegion_qSVAjoint_withAnnotation.rda
.Could use
ggplot2
and facetting, similar to say https://github.com/LieberInstitute/brainseq_phase2/blob/master/twas/pdf/pgc2_twas_z.pdf or base R like in https://github.com/LieberInstitute/brainseq_phase2/tree/master/check_noQsva#with-vs-without-qsva.