Entering edit mode
8.6 years ago
QVINTVS_FABIVS_MAXIMVS
★
2.6k
I have SHAPEIT haplotype scaffolds for a 1Mb region. I would like to determine the R2 of pairwises SNPs in a haplotype dependent manner.
Is there an option in PLINK that will produce LD statistics with haplotype in mind?
Or am I missing the point of LD and the results would be the same if the SNPs were not phased?
PLINK 1.x is not able to import phase information. LD statistics which take it into account will be a bit different.
Could you try separating the haplotypes into different plink files, then calculate it per (haplotype specific) plink file? You might have to do a lot of slicing and dicing, but you might be able to get it done this way.