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