View source: R/groupedPearson.R
groupedPearson | R Documentation |
Combine p-values from grouped tests with Pearson's method. Groups are defined according to unique levels of a grouping factor.
groupedPearson(p.values, grouping, log.p = FALSE)
p.values |
A numeric vector containing p-values for individual tests. |
grouping |
A vector or factor of length equal to Alternatively, an rle object where each run corresponds to a group and specifies the entries of |
log.p |
Logical scalar indicating whether the p-values in |
Here, the joint null hypothesis for each group is that all of the individual null hypotheses are true. Pearson's method combines information from all individual nulls to determine if the joint null should be rejected. Compared to Stouffer's and Pearson's methods, Pearson's method is more sensitive to the largest individual p-value. This method is only applicable to independent tests and no weights are considered.
The representative test for each group is defined as the test with the largest p-value, as this has the greatest effect on the combined p-value. All tests for each group are considered to be influential as increasing any of them (e.g., to unity) would result in a larger combined p-value.
A list containing:
p.value
, a named numeric vector of length equal to the number of unique levels in grouping
.
This contains the Pearson combined p-value for each group, log-transformed if log.p=TRUE
.
Each entry is named according to the group.
representative
, a named integer scalar specifying the representative test for each group.
Each index refers to an entry of p.values
and is named according to its group.
influential
, a logical vector of length equal to p.values
.
Entries are TRUE
for any p-value that is deemed “influential” to the final combined p-value for its group.
Aaron Lun
Pearson K (1934). On a new method of deternining “goodness of fit.” Biometrika 26, 425-442.
parallelPearson
, for a version that operates on parallel vectors of p-values.
groupedFisher
and groupedStouffer
, for different approaches to testing a joint null of independent hypotheses.
p1 <- rbeta(100, 0.8, 1)
g <- sample(10, length(p1), replace=TRUE)
# Standard application:
out <- groupedPearson(p1, g)
str(out)
# With log p-values.
out <- groupedPearson(log(p1), g, log.p=TRUE)
str(out)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.