Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Interpretation of PC1 values among groups #92

Open
masaver opened this issue Oct 29, 2020 · 1 comment
Open

Interpretation of PC1 values among groups #92

masaver opened this issue Oct 29, 2020 · 1 comment

Comments

@masaver
Copy link

masaver commented Oct 29, 2020

Hi Gabriel,
I'm looking at a cancer transcriptomics dataset with Normal and Tumor samples, where I'm trying to quantify differences in pathways activity between Normal and Tumor tissue. I was able to run the analysis with the AESPCA_pVals() function and I'm now looking at the PC1 values for the EGFR pathway( which I undertand are a proxy for the activity of this pathway) (see picture below):

egfr

I expect/know that for the specific cancer type I'm studying the EGFR pathway is more active in Tumor than in Normal Tissue.

So my question is: How do I go from interpreting positive or negative PC1 values to to saying the X pathway is more active in a given condition ?

Thanks in advance for the help.

Best ,

-Mathias

@masaver masaver changed the title Interpretation of PC1 values amon groups Interpretation of PC1 values among groups Oct 29, 2020
@gabrielodom
Copy link
Owner

Hi Mathias,
I hate to be the bearer of bad news, but principal components are "unique up to a sign". That means that they can't be used to interpret whether or not a pathway is up-regulated or down-regulated, only that it is differentially regulated. We recommend that users target pathways for post-hoc analysis based on the p-values (or FDRs) yielded by the AESPCA_pVals() or SuperPCA_pVals() functions. In your case, it's mathematically appropriate to multiply the EGFR pathway PC by -1 and re-apply your statistical tests.

See this thread for a related question and answer: https://stats.stackexchange.com/a/30352/171235

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants