I have a list of about 250 SNPs. I want to obtain the proxy SNPs for them using the 1000 genomes. I have downloaded corresponding .vcf and .ped files from 1000 genomes. I have also used PLINK to obtain the .bim, .fam and .bed files.
How do I proceed to determine the proxy SNPs without using any software?
Thank you.
1 answer
Get LD for selected SNPs using plink --r2. Then pick any proxy SNP having highest LD score.
Why not using other software? We could, for example, use PriorityPruner which is specifically designed to pick proxies with prioritisation, i.e. when SNP1 is in LD with SNP2 and SNP3, it will prioritise based on other features of the SNP (like pvalue, frequency, annotations, etc.).
PriorityPruner is a software program which can prune a list of SNPs that are in high linkage disequilibrium (LD) with other SNPs in the list, while preferentially keeping SNPs of higher priority (e.g., the most significant SNPs in a genome-wide association study). The process of pruning SNPs based on LD is sometimes referred to as "LD clumping".
Log in to answer this question.
-
Merge only bim files with plink
written by mgdrnl •
-
Merging more than two files with plink
written by phmai1148 •
-
Reference panel data to be used for GCTA-COJO
written by Apprentice
-
How to obtain subpopulation specific allelic variants from a .vcf file using PLINK?
written by percyjackson2k20 •
-
Imputation 1000 Genome and plink file
written by mzezza •
-
Using PLINK to filter VCF files
written by Diego.Morales •
-
plink 1.9 PED to VCF, SNPs differ on input and output files
written by stolarek.ir
-
PLINK: Change affection status
written by bsmith030465
-
PLINK: merge .bed files
written by bsmith030465
-
Convert Bam To Bed/Fam/Bim
written by Anjali •
Thank you very much. I am using snpStat to calculate the LDs.