This function takes a list of variants and searches for them in samples from 1000 Genomes phase 3 data It then creates an LD matrix of r values (signed, and not squared) All LD values are with respect to the major alleles in the 1000G dataset. You can specify whether the allele names are displayed

ld_matrix(
variants,
with_alleles = TRUE,
pop = "EUR",
bfile = NULL,
)