cisAssoc {GGtools} | R Documentation

## test for variant-expression associations in cis, using VCF

### Description

test for variant-expression associations in cis, using VCF and SummarizedExperiment representations

### Usage

    cisAssoc(summex, vcf.tf, rhs = ~1,
       nperm = 3, cisradius = 50000,
       genome = "hg19", assayind = 1, lbmaf = 1e-06,
       dropUnivHet = TRUE, doEsts=FALSE)

    data(lgeu) # obtains an example SummarizedExperiment

### Arguments

`summex` | instance of `SummarizedExperiment-class`
`vcf.tf` | instance of `TabixFile`, referring to a tabix-indexed, bgzipped VCF file
`rhs` | formula ‘right hand side’ for adjustments to be made as `snp.rhs.tests` is run on each expression vector
`nperm` | number of permutations to be used for plug-in FDR computation
`cisradius` | distance in bp around each gene body to be searched for SNP association
`genome` | tag suitable for use in GenomeInfoDb structures
`assayind` | index of `assays(summex)` to use for expression data retrieval
`lbmaf` | lower bound on MAF of SNP to use
`dropUnivHet` | logical, if TRUE, will check for columns of SnpMatrix instance that possess no values other than "NA" and "A/B". See http://www.biostars.org/p/117155/#117270
`doEsts` | logical, if TRUE, will run `snp.rhs.estimates` and report beta and standard error

### Details

`snp.rhs.tests` is the workhorse for statistical modeling. VCF content is transformed to the byte-code (which allows for uncertain imputation) and used in fast testing.

### Value

a `GRanges-class` instance with mcols including chisq, permScore...

### Note

seqlevelsStyle for summex and vcf.tf content must agree

### Author(s)

VJ Carey <stvjd@channing.harvard.edu>

### Examples

     data(lgeu) # small excerpt from GEUVADIS FPKM
     tf20 = TabixFile(system.file("vcf/c20exch.vcf.gz", package="GGtools"))
     if (require(VariantAnnotation)) scanVcfHeader(tf20)
     lgeue = clipPCs(lgeu[,which(lgeu$popcode=="CEU")], 1:2)
     set.seed(4321)
     litc = cisAssoc(lgeue, tf20, nperm=2, lbmaf=.05, cisradius=50000)
     litc2 = cisAssoc(lgeue, tf20, nperm=2, lbmaf=.05, cisradius=50000, doEsts=TRUE)
     summary(litc$chisq)
     litc$pifdr = pifdr(litc$chisq, c(litc$permScore_1, litc$permScore_2))
     litc[which(litc$pifdr < .01)]

[Package GGtools version 5.16.0 Index]