-----Original Message-----
From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces at r-project.org] On
Behalf Of Philippe Tadger
Sent: Wednesday, 24 March, 2021 20:16
To: r-sig-meta-analysis at r-project.org
Subject: Re: [R-meta] pvalues from rma
Dear group
I would like to check how to obtain p-values form metafor::rma
I'm trying to replicate the first lines of the code from:
https://www.metafor-project.org/doku.php/tips:assembling_data_or
<https://www.metafor-project.org/doku.php/tips:assembling_data_or>
library(metafor)
dat.bcg <- dat.bcg[,c(2:7)]
#dat.bcg
dat1 <- escalc(measure="OR", ai=tpos, bi=tneg, ci=cpos, di=cneg,
data=dat.bcg)
#dat1
res1 <- rma(yi, vi, data=dat1)
#res1
#predict(res1, transf=exp, digits=2)
dat2 <- data.frame(summary(dat1))
dat2[c("yi", "ci.lb", "ci.ub")] <- data.frame(summary(dat1,
transf=exp))[c("yi", "ci.lb", "ci.ub")]
*######################The previous step need to provide pval, but is
not giving such outcome.**######################*
names(dat2)[which(names(dat2) == "yi")] <- "or"
dat2[,c("or","ci.lb","ci.ub","pval")] <-
round(dat2[,c("or","ci.lb","ci.ub","pval")], 2)
--
Kind regards
*Philippe Tadger*
Statistician (Msc.). Orthopedic Manual Therapist (Mag.).
Physical Therapist (Bsc.)
+32498774742.