I am a freshmen with plink. I downloaded the 1000 genome phase 3 data and try to get all the snp pairs with ld score > 0.8 in a matrix format. I use the following command:
plink --vcf ALL.chr1.phase3_shapeit2_mvncall_integrated_v5a.20130502.genotypes.vcf.gz --keep-fam EUR.id --geno 0.05 --hwe 1e-6 --make-bed --out tmp
but when I tied to filter the r2 value and get matrix output with this command:
plink --bfile tmp --r2 square --ld-window-r2 0.8 --out EUR_chr1
it showed that "Error: --ld-window-r2 flag cannot be used with the --r2 matrix output modifiers."
if I didn't use "square" argument, I only got all the snp pairs with ld score > 0.8 in a table format, like "EUR_chr1.ld". Is there any way to change "EUR_chr1.ld" into matrix with plink?
If I could get the matrix directly, that would be better.
thanks for any helps.