[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:
Hi Imran, There was already an answer by Michael Dewey with a useful link to MetaAnalysis Task View https://cran.r-project.org/view=MetaAnalysis Best, Gerta Am 28.09.2021 um 10:06 schrieb Awais Wahab:
Hello!
Is it possible to perform meta-analysis with available P-values in each
study as effect size, instead of HR, OR etc. ?
Thanks in advance.
Best regards,
Imran
[[alternative HTML version deleted]]
_______________________________________________ R-sig-meta-analysis mailing list R-sig-meta-analysis at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis
-- Dr. rer. nat. Gerta R?cker, Dipl.-Math. Guest Scientist Institute of Medical Biometry and Statistics, Faculty of Medicine and Medical Center - University of Freiburg Zinkmattenstr. 6a, D-79108 Freiburg, Germany Mail: ruecker at imbi.uni-freiburg.de Homepage: https://www.uniklinik-freiburg.de/imbi-en/employees.html?imbiuser=ruecker