Skip to content

LDlinkR: minor issue with LDexpress - column hg38 lists hg19 positions #18

@timyers

Description

@timyers

Senf to LDlinkR maintainer email address from user ZEY, with minor edits:

I am contacting you in regards to your recently published R package LDlinkR.
......
However, I’d like to bring to your attention a minor issue that I have encountered with it that you might want to know about.

I was testing the LDexpress command as followed (requesting the SNPs to be annotated with hg38 positions in the output)
LDexpress(snps = c("rs345", "rs456"),
pop = c("CEU"),
tissue = c("LIV", "TES"),
r2d = "r2",
r2d_threshold = "0.1",
p_threshold = "0.1",
win_size = "500000",
genome_build = "grch38_high_coverage",
token = "..."
)
However, the output provided hg19 positions despite the column being named >Position_grch38_high_coverage<.

Query | RS_ID | Position_grch38_high_coverage | R2 |  
rs345 | rs756910 | chr13:32524656 | 0.473404255319149 |  
rs345 | rs756911 | chr13:32525073 | 0.473404255319149 |  
rs456 | rs17150305 | chr7:24891093 | 0.102998110552808
3 rows | 1-4 of 12 columns

The first SNV should be chr13:31950519 in hg38 build.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions