Wiki » History » Version 1
Katie Lennard, 07/17/2020 11:32 AM
| 1 | 1 | Katie Lennard | # Wiki |
|---|---|---|---|
| 2 | |||
| 3 | # Data analysis background |
||
| 4 | |||
| 5 | The data was originally analysed by overseas collaborators (Hans mollenkopf@mpiib-berlin.mpg.de and Natalie natalie.nieuwenhuizen@uni-wuerzburg.de). Hans provided .RData objects and .RMD file summarizing his results, which were quite comprehensive (including exploratory analysis, differential expression using limma and gene module enrichment using the R package tmod). However there are very strong batch effects in this dataset that correspond to RNA extraction date in the lab. I therefore redid the analysis with combat batch correction (which was only possible by correcting the astrocyte and neuron samples separately). They were therefore analysed sepearately downstream. |
||
| 6 | |||
| 7 | # Data preprocessing |
||
| 8 | Hans performed background correction (normexp method) and normalization (quantile). The data was supplied to me as .RData objects (raw and normalized). I added ComBat batch correction. |
||
| 9 | |||
| 10 | # Downstream analysis (see attached powerpoint) |
||
| 11 | |||
| 12 | PCA, tSNE multivariate |
||
| 13 | |||
| 14 | Differential expression analysis between groups: |
||
| 15 | Astrocytes: not infected vs. H37Rv; not infected vs BC; BCG vs H37Rv |
||
| 16 | Neurons: not infected vs. H37Rv |
||
| 17 | |||
| 18 | Transcriptional module testing: using R package tmod |