pvals_lite: Computes p-values

In TomMayo/M3Ddevel: (development for M3D) Identifies differentially methylated regions across testing groups

Description

Returns p-values for each region reflecting the probability of observing the mean test-statistic of the between group comparisons among the inter-replicate comparisons. The 'lite' version works with the outcome of M3D_Wrapper_lite and the 'lite' parallel version.

Usage

pvals_lite(rrbs, CpGs, M3D_stat_lite, group1, group2, smaller = FALSE,
  comparison = "allReps", method = "empirical", closePara = 0.005)

Arguments

rrbs	An rrbs object containing methylation and coverage data as created using the BiSeq pacakge
CpGs	A GRanges object with each row being a testing region
M3D_stat_lite	A matrix containing the M3D test-statistic, as returned from the M3D_Wrapper_lite function. This has the mean of the MMDs for the comparisons between groups in the first column, and the within groups MMDs in the others, for each region in the CpGs object.
group1	The name of the first group for the comparison. This is stored in colData(rrbs)
group2	The name of the second group for the comparison. This is stored in colData(rrbs)
smaller	Determines whether the p-value is computed whether the test-statistic is greater or lesser than inter-replicate values. For our purposes, it should be set to FALSE.
comparison	Details which groups we are using to define our empirical testing distribution. The default is to use all of them, however, should the user find one group contains unusually high variability, then that group can be selected. Values are either 'allReps', 'Group1' or 'Group2'.
method	Determines which method is used to calculate p-values. 'empirical' uses the empirical distribution directly, without modelling. This is the default. 'model', fits an exponential distribution to the tail of our null distribution.
closePara	Sets a threshold for how close the exponential curve should fit the empirical distribution in the 'model' method. If the method produces errors, consider raising this parameter.

Value

Returns a list P, with 2 entries. 'FDRmean' is the Benjamini-Hochberg adjusted p-values. The unadjusted p-values are stored in 'Pmean'.

Author(s)

Tom Mayo [e-mail]

Examples

data(rrbsDemo)
data(CpGsDemo)
data(MMDlistDemo)
M3D <- MMDlistDemo$Full-MMDlistDemo$Coverage
M3D <- cbind(rowMeans(M3D[,2:5],na.rm=TRUE),M3D[,1],M3D[,6])
colnames(M3D) <- c('MeanBetween','H1-hESC1  vs  H1-hESC2','K562-1  vs  K562-2')
group1 <- unique(colData(rrbsDemo)$group)[1]
group2 <-unique(colData(rrbsDemo)$group)[2]
PDemo <- pvals(rrbsDemo, CpGsDemo, M3D,
            group1, group2, smaller=FALSE,comparison='allReps')

TomMayo/M3Ddevel documentation built on May 9, 2019, 4:53 p.m.