This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Prune SNPs in LD with a list of SNPs

Hi,

I have a list of SNPs. I want to find the SNPs in a genotype (Plink format) data set, which are in LD (> 0.2) with the list of SNPs. Any help on how to achieve this is appreciated.

Thanks.

snp plink ld
• 3,511 views • link
written by aritra90 aritra90

I couldn't figure out the right command to enter a SNP list and a genotype file and obtain another snp list which are in high LD. For now I've been using indep-pairwise to remove the SNPs in high LD beforehand and then do downstream analyses. However, I wanted to remove SNPs which are in LD with a particular subset of SNPs (provided in a file such as what --extract/--exclude takes in as input)

• 1 views • link
written by aritra90 aritra90

Not sure but you may try the --ld-snp-list flag, which accepts a file-listing of input SNPs for the purpose of identifying LD.

• 1 views • link
written by Kevin Blighe Kevin Blighe

0 answers

No answers yet.

Log in to answer this question.

More posts like this