computeRegionSignificance: Compute significance for given regions

Description Usage Arguments Details Value Author(s) Examples

Description

For a given set of regions, region-wise pvalues and FDR is computed

Usage

1

Arguments

fit

A GenoGAM object containing the fit

regions

A GRanges object of regions of interest

smooth

Which fit should be used. The names should be equivalent to the column names of the object. Lookup with colnames(my_GenoGAM_object)

Details

For a given set of regions, region-wise pvalues are computed by applying familywise hochberg correction and taking the minimal p-value. FDR is computed by further applying Benjamini-Hochberg correction.

Value

The GRanges object from the 'region' parameter extended by two columns: pvalue and FDR

Author(s)

Georg Stricker georg.stricker@in.tum.de

Examples

1
2
3
4
5
6
## make test GenoGAM
gg <- makeTestGenoGAM()
## make region
region <- GRanges("chrXYZ", IRanges(c(2000, 4000, 6000), c(3000, 5000, 9000)))
res <- computeRegionSignificance(gg, region)
res

gstricker/fastGenoGAM documentation built on May 17, 2019, 8:56 a.m.