computeRegionSignificance {GenoGAM} | R Documentation |
For a given set of regions, region-wise pvalues and FDR is computed
computeRegionSignificance(fit, regions, what = NULL)
fit |
A GenoGAM object containing the fit |
regions |
A GRanges object of regions of interest |
what |
Which fit should be used. The names should be equivalent to the column names used in the config file.
Lookup with |
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.
The GRanges object from the 'region' parameter extended by two columns: pvalue and FDR
Georg Stricker georg.stricker@in.tum.de
gg <- makeTestGenoGAM() gr <- GRanges("chrI", IRanges(1,100)) computeRegionSignificance(gg, gr)