[R-meta] pvalues from rma

[Philippe Tadger]

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.

	[[alternative HTML version deleted]]