This is a test version of Biostars. For the public version, visit https://www.biostars.org.
PCA to find clusters in gene expression data

Dear community,

I finally got some Gene Expression Data on some patients with cancer vs. healthy patients as control. I want to investigate if the expression data form distinct clusters. That's why I performed PCA for dimensionality reduction and plotted it with sns.scatterplot: enter image description here

You can see that patients with cancers form a cluster (blue) and the healthy patients form another (orange) cluster. But:

1.) Is it valid, to perform PCA to identify cluster or do I need to do other clustering methods, like t-SNE plots or k-means clustering? 2.) How can I show that the clusters are significantly different? Can I calculate p-values? Is it also possible to plot confidence ellipses?

I would be glad for every help!

expression clustering pca
• 1,341 views • link
updated by ATpoint ATpoint • written by Hashirama • Hashirama

This is a case where you don’t need a statistical test to find out whether the clusters are distinct.

• 0 views • link
written by 4galaxy77 4galaxy77

1 answer

All I see is two distinct populations and one outlier. Remove the outlier sample and compare orange to blue. Does that make sense?

• 0 views • link
written by ATpoint ATpoint

Thanks ATpoint for your fast and helpfull answer! But what do you actually mean by "comparing"? What are the next steps?

• 0 views • link
written by Hashirama • Hashirama

I guess you want sifferential genes? Check tools wuch as limma-voom.

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this