plink計算兩個SNP位點的連鎖不平衡值（LD）

PLINK提供了“--ld”的參數計算兩個SNP位點的連鎖不平衡值。

命令如下：

plink --file file --ld rs123 rs134 --out rs123_rs134

 

生成如下數據：

--ld rs123 rs134:

R-sq = 0.0313386       D' = 1

Haplotype Frequency Expectation under LE
--------- --------- --------------------
TG 0 0.022549
CG 0.171717 0.149168
TA 0.131313 0.108764
CA 0.696970 0.719518

In phase alleles are TA/CG