ldByGene {ldblock}R Documentation

obtain LD statistics in region specified by a gene model

Description

Obtain LD statistics in region specified by a gene model.

Usage

ldByGene(sym = "MMP24", 
  vcf = system.file("vcf/c20exch.vcf.gz", package = "gQTLstats"), 
  flank = 1000, 
  vcfSLS = "NCBI", 
  genomeSLS = "hg19", 
  stats = "D.prime", depth = 10)

Arguments

sym

A standard gene symbol for use with genemodel

vcf

Path to a tabix-indexed VCF file

flank

number of basepairs to flank gene model for search

vcfSLS

seqlevelsStyle (SLS) token for VCF; will be imposed on gene model

genomeSLS

character tag for genome, to be used with readVcf

stats

passed to ld

depth

passed to ld

Value

sparse matrix representation of selected LD statistic, as returned by ld

Examples

if (.Platform$OS.type == "unix") {
 ld1 = ldByGene(depth=150)
 image(ld1[1:200,1:200], col.reg=heat.colors(120), colorkey=TRUE,
   main="SNPs in MMP24 (chr20)") 
 }

[Package ldblock version 1.12.1 Index]