I was looking for: fisher.comb <- function (pvalues) { df=length(pvalues) ch2=(-2*sum(log(pvalues))) return pchisq(ch2, df=df, lower.tail=FALSE) } http://en.wikipedia.org/wiki/Fisher%27s_method for the second part: the calculation of the p-value from the result of fischer combination test . Kind regards Knut