This function queries HaploReg web-based tool and returns results.
queryHaploreg(
query = NULL,
file = NULL,
study = NULL,
ldThresh = 0.8,
ldPop = "EUR",
epi = "vanilla",
cons = "siphy",
genetypes = "gencode",
url = Haploreg.settings[["base.url"]],
timeout = 100,
encoding = "UTF-8",
querySNP = FALSE,
fields = NULL,
verbose = FALSE
)
A data frame (table) with results similar to HaploReg uses.
Query (a vector of rsIDs).
A text file (one refSNP ID per line).
A particular study. See function getHaploRegStudyList(...)
.
Default: NULL
.
LD threshold, r2 (select NA to only show query variants). Default: 0.8.
1000G Phase 1 population for LD calculation. Can be: "AFR", "AMR", "ASN". Default: "EUR".
Source for epigenomes.
Possible values: vanilla
for ChromHMM (Core 15-state model);
imputed
for ChromHMM (25-state model using 12 imputed marks);
methyl
for H3K4me1/H3K4me3 peaks;
acetyl
for H3K27ac/H3K9ac peaks.
Default: vanilla
.
Mammalian conservation algorithm.
Possible values: gerp
for GERP,
siphy
for SiPhy-omega,
both
for both.
Default: siphy
.
Show position relative to.
Possible values: gencode
for Gencode genes;
refseq
for RefSeq genes;
both
for both.
Default: gencode
.
HaploReg url address. Default: <https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php>
A timeout
parameter for curl
.
Default: 100
sets the encoding
for correct retrieval web-page content.
Default: UTF-8
A flag indicating to return query SNPs only.
Default: FALSE
A set of fields to extract. Refer to the package vignette
for available fields. Default: All
.
Verbosing output. Default: FALSE.
library(haploR)
data <- queryHaploreg(c("rs10048158","rs4791078"))
head(data)
Run the code above in your browser using DataLab