Combine p-values using the sum p method
sump(p)
# S3 method for sump
print(x, ...)An object of class ‘sump’ and
‘metap’, a list with entries
The transformed sum of the p-values
See details
The input vector with illegal values removed
A vector of significance values
An object of class ‘sump’
Other arguments to be passed through
Michael Dewey
Defined as (_i=1^k p_i)^kk! - k 1(_i=1^k p_i - 1)^kk! + k 2(_i=1^k p_i - 2)^kk! ... ((sum p) ^ k) / k! - (k-1)C(1) ((sum p - 1) ^ k) / k! + (k-2)C(2) ((sum p - 2) ^ k) / k! ... where there are k studies and the series continues until the numerator becomes negative edgington72ametap.
Some authors use a simpler version
(_i=1^k p_i)^kk!((sum p) ^ k) / k!
but this can be very conservative when
_i=1^k p_i > 1sum p > 1.
There seems no particular need to use this method but
it is returned as the value of conservativep
for use in checking published values.
The values of p_i should be such that 0 p_i 1 and a warning is given if that is not true. A warning is given if, possibly as a result of removing illegal values, fewer than two values remain and the return values are set to NA. A warning is given when the internal calculations are likely to have been subject to numerical error and an alternative method should be used to check the result.
The plot method for class ‘metap’ calls plotp on the valid \(p\)-values.
See also plotp
data(dat.metap)
edgington <- dat.metap$edgington
sump(edgington) # p = 0.097
Run the code above in your browser using DataLab