Skip to content
Prev 3308 / 5632 Next

[R-meta] Fwd:

Hello!
I need some help please.
I have 5-6 one sided P-values from different studies and need to do their
meta-analysis.

I am using these (Example) codes in R as below, but they are not working,

*> library(metap)*

*pvals <- c(0.1, 0.1, 0.9, 0.9, 0.9, 0.9)  *

*> istwo <- c(TRUE, FALSE, TRUE, FALSE, TRUE, FALSE)*

*> toinvert <- c(FALSE, TRUE, FALSE, FALSE, TRUE, TRUE)*

*> two2one(pvals, two = istwo, invert = toinvert)*

*> print(validity)*

*> plotp(validity)*

*> schweder(validity)*


*OR this below,*



*Fisher's method, T = -2*log(p-value1)-2*log(p-value2)-...-2*log(p-valuem)
where m is the number of p-values. Under the null hypothesis, this will be
distributed as chi-squared with m degrees of freedom, so you would then do
T.pval = 1-pchisq(T,m)*


Are these codes correct ?  I am not sure where exactly to put my p-values.
It will be great if someone can help me with the correct codes.


I will be much grateful for your kind support.

best regards.

On Tue, Sep 28, 2021 at 11:27 AM Dr. Gerta R?cker <
ruecker at imbi.uni-freiburg.de> wrote: