I have a list of z-scores (obtained via effect size/std.err of effect size). I want to combine all the z-scores via Stouffer's method but without weights. Till now I have seen "Stouffer.test" implemented in metaseq package but it requires weights along with the z-score. I have tried "sumz" method implemented in metap package which according to my understanding can work without weights as follows:
`sumz(p, weights = NULL)` -> where p is the vector of values; in my case z-scores
My understanding of Stouffer's method is
sum of all the z-scores/sqrt of total number of samples
when I compared the results of the "sumz" method as described above, it is not the same as the above formula when I calculated it in excel using the following
My question is : Is there any R function which works like above (sum of all the z-scores/sqrt of total number of samples) so that I can cross check the results? Or I have to do it manually?