Read in GWAS summary data with filters on datasets (if multiple datasets per file) and/or chromosome/position, rsids or pvalues. Chooses most optimal choice for the detected operating system. Typically chrompos searches are the fastest. On Windows, rsid or pvalue filters from a file will be slow.
query_gwas(
vcf,
chrompos = NULL,
rsid = NULL,
pval = NULL,
id = NULL,
rsidx = NULL,
pvali = NULL,
build = "GRCh37",
os = Sys.info()[["sysname"]],
proxies = "no",
bfile = NULL,
dbfile = NULL,
tag_kb = 5000,
tag_nsnp = 5000,
tag_r2 = 0.6,
threads = 1
)Path or URL to GWAS-VCF file or VCF object e.g. output from VariantAnnotation::readVcf() or create_vcf()
Either vector of chromosome and position ranges e.g. "1:1000" or "1:1000-2000", or data frame with columns chrom, start, end.
Vector of rsids
P-value threshold (NOT -log10)
If multiple GWAS datasets in the vcf file, the name (sample ID) from which to perform the filter
Path to rsidx index file
Path to pval index file
="GRCh37" Build of vcffile
The operating system. Default is as detected. Determines the method used to perform query
="no" If SNPs are absent then look for proxies (yes) or not (no). Can also mask all target SNPs and only return proxies (only), for testing purposes. Currently only possible if querying rsid.
=path to plink bed/bim/fam ld reference panel
=path to sqlite tag snp database
=5000 Proxy parameter
=5000 Proxy parameter
=0.6 Proxy parameter
=1 NUmber of threads
vcf object