I have a bunch of independent $p$-values and now I want to combine them using the Fisher's method. However, in R, when a $p$-value is zero, the
log(p-value) becomes $-\infty$, so that the test statistic $X^2$ is $\infty$. This will give the Fisher's method $p$-value of 0.
If I check the $p$-values' histogram, they look pretty uniform, so I expect that the Fisher's method should give me a not-so-small $p$-value (should not reject the null). Is there a way to get around this issue? How about the Stouffer's method? Will that method automatically take care of the issue when $p = 0$ (or close to 0)?
Update: I tried the Stouffer's method, but still the issue of $\infty$ is not solved...