Closed lcolladotor closed 3 years ago
I updated a differential expression script to run on the sample files. We may want to look at the model and pca but for now this holds as an example script.
Hi @joshstolz, can you verify whether it's ok to close this issue? Thanks!
Yeah we finished that up.
From: Leonardo Collado-Torres @.***> Sent: Saturday, May 1, 2021 10:22:55 PM To: LieberInstitute/SPEAQeasy-example Cc: Josh Stolz; Mention Subject: Re: [LieberInstitute/SPEAQeasy-example] Code for differential expression analysis and EDA (#6)
Hi @joshstolzhttps://github.com/joshstolz, can you verify whether it's ok to close this issue? Thanks!
— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://github.com/LieberInstitute/SPEAQeasy-example/issues/6#issuecomment-830728652, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AOZ2YOQMNZMOLWNLJWRUNNDTLSZP7ANCNFSM4ORX77QA.
Create a directory and Rmd that uses the output of #5 as the input. Then show how to explore the data (PCA), check the ERCCs (Nick can do this part), normalize the data, perform a differential expression analysis using
limma
-voom, find enriched ontologies usingclusterProfiler
and make some heatmaps likely usingpheatmap
. This could be our own LIBD jaffelab guide for differential expression analyses. Feel free to link to Bioconductor workflows and package vignettes as well. You could also usejaffelab::agePlotter()
for example. Use the sameBiocStyle
andknitcitations
style as in #5 for the Rmd (with also the output fromsessioninfo::session_info()
).