PCA for combined set
1
0
Entering edit mode
3.3 years ago
raalsuwaidi ▴ 100

hello all,

i am trying to plot the PCA values for our GWAS study using plink. below is the code used and results:

plink --vcf local_chr2.vcf.gz --maf 0.001 --indep-pairwise 50 5 0.2 --vcf-half-call missing --out chr2-clean plink --vcf local_chr2.vcf.gz --genome gz --out local_chr2.vcf-2 --vcf-half-call missing --extract chr2-clean.prune.in plink --vcf local_chr2.vcf.gz --extract chr2-clean.prune.in --pca var-wts --vcf-half-call missing --out chr2-clean

the chart from the R code is below

pcs = read.table( "~/Desktop/chr2-clean.eigenvec" ) plot( pcs[,3], pcs[,4], xlab = "PC1", ylab = "PC2" )

Rplot for local GWAS

then i went on to compare the PCA values against the 1000 genome. for this I merged the the two vcf files using bcftools using the following code: bcftools merge *vcf.gz -Oz -o Merged.vcf.gz

when running the plink PCA command on the merged version i got an empty PCA plot

plink --vcf Merged.vcf.gz --extract chr2-clean.prune.in --pca var-wts --vcf-half-call missing --out merged.with.1000G

PCA plot for merged vcf

any help on what i did wrong? i am so confused? is it because i used bcftools to merge? because for verification i also plotted PCA for the 1000 genome values and the diagram was fine.

plink bcftools vcf PCA • 1.4k views
ADD COMMENT
2
Entering edit mode
3.3 years ago

Have you tried using plink 2.0 for PCA? It should give you a relevant error message when it fails.

ADD COMMENT
0
Entering edit mode

The PCA on plink 1.9 didn’t fail. It just gave empty results for the combined vcf file.

ADD REPLY
1
Entering edit mode

That’s a failure.

ADD REPLY
0
Entering edit mode

Ok. I will try to run it again with version 2.0

ADD REPLY

Login before adding your answer.

Traffic: 2499 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6