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 = "http://archive.broadinstitute.org/mammals/haploreg/haploreg.php",
timeout = 100, encoding = "UTF-8", querySNP = FALSE, fields = NULL,
verbose = FALSE)
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: <http://archive.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.
A data frame (table) with results similar to HaploReg uses.
# NOT RUN {
library(haploR)
data <- queryHaploreg(c("rs10048158","rs4791078"))
head(data)
# }
Run the code above in your browser using DataLab