[R] weighted inverse chi-square method for combining p-values
pburns at pburns.seanet.com
Wed Jun 25 12:53:30 CEST 2008
Are you sure you want to use Fisher's method?
page 13 shows an example where Fisher's method
is not attractive, but Stouffer's method gives reasonable
results. Weighting in Stouffer's method is straightforward.
patrick at burns-stat.com
+44 (0)20 8525 0696
(home of S Poetry and "A Guide for the Unwilling S User")
Andreas Cederfeldt wrote:
> This is more of a general question than a pure R one, but I hope that is OK.
> I want to combine one-tailed independent p-values using the weighted version
> of fisher's inverse chi-square method. The unweighted version is pretty
> straightforward to implement. If x is a vector with p-values, then I guess
> that this will do for the unweighted version:
> statistic <- -2*sum(log(x))
> comb.p <- 1-pchisq(statistic,2*length(x))
> Has anyone written some R code for the weighted version?
> Anders Cederfeldt
> [[alternative HTML version deleted]]
> R-help at r-project.org mailing list
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help