0
$\begingroup$

I have a raw counts data-set of 20,502 genes and 137 samples. I want to find out Principal Components which best explain variation between samples in different stages of tumor.

I am new to Machine learning and would like some help in selecting a criteria to reduce the number of genes before doing a PCA.

I have tried doing a DESeq2 analysis on the data-set but I am confused as to what to do next.

$\endgroup$
  • 1
    $\begingroup$ You can filter our genes that do not show variation across samples, they would not be helpful in differentiating samples anyhow. $\endgroup$ – haci Dec 5 '19 at 9:35
  • $\begingroup$ So you mean genes which do not fit the dispersion curve in DESeq2? $\endgroup$ – Pawan Verma Dec 5 '19 at 9:37
  • 1
    $\begingroup$ More like they did here: pklab.med.harvard.edu/scw2014/subpop_tutorial.html $\endgroup$ – haci Dec 5 '19 at 11:09
  • $\begingroup$ This is really helpful. Thanks a lot :) $\endgroup$ – Pawan Verma Dec 5 '19 at 11:11
  • $\begingroup$ As you found the approach useful, I will re-write this as a formal answer for future reference. $\endgroup$ – haci Dec 5 '19 at 11:17
1
$\begingroup$

You can filter our genes that do not show variation across samples, they would not differentiate samples anyhow. For the specifics on how to do so, please see the Subpopulation Analysis section of a nice single cell RNA-seq workflow from the Kharchenko Lab. I believe there won't be fundamental differences when applying this approach to bulk RNA-seq.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.