`LD_blocks.sh` is not functional
cbird808 opened this issue · 1 comments
cbird808 commented
The awk
line assumes that the second column of a *.ld
file has chom:pos info, but the second column has the major nucleotide of snp1. The fourth column (change $2
to $4
) has the snp2 chrom:pos info. Then the subsequent cut
needs to be updated to reflect the *.ld
file format