require(metafor)
require(mada)
data(cervical)
LAG <- cervical[cervical$method==2,]
fit1 <- reitsma(LAG)
summary(fit1) # results of the bivariate meta-analysis
###
attach(LAG)
dta1 <- edta(TP,FN,TN,FP)
###
attach(dta1)
bifunnel(y,S)
Run the code above in your browser using DataLab