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

PCA summary plots rely on fam file over phenotype #19

Open
rkwalters opened this issue Nov 3, 2015 · 0 comments
Open

PCA summary plots rely on fam file over phenotype #19

rkwalters opened this issue Nov 3, 2015 · 0 comments

Comments

@rkwalters
Copy link
Contributor

Plots generated by pcaer (e.g. MDS plots, summary plot of PC association to phenotype) use phenotypes from the individuals FID tag (cas, con, mis) rather than from the fam file. This means that any updates to the phenotype after initial QC (e.g. for alternative outcomes, updated phenotyping, etc) that are applied to the fam file will not be reflected in the PCA output.

This is particularly problematic for the summary pdf, since the associations between the original phenotype (indicated in the FID) and the PCs will be reported rather than association with the updated/alternative phenotype, which may be misleading when using the pcaer output to identify necessary PC covariates for GWAS with the updated phenotype.

Affected code is in pca_plot_5 with the creation of the $mds_ow_file and subsequent use of the CC column in R for plotting (lines 350-600).

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

1 participant