#######################################################################################
##Odds Ratio to bycatch (group)
library(meta)
library(mvmeta)
library(metafor)


dat.bcg<-read.csv("meta.csv" , header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE)
summary(dat.bcg)
names(dat.bcg)

res <- metabin(tpos,tneg,cpos,cneg, data=dat.bcg,
               sm="OR", method="I",studlab=paste(esp))
res
summary(res)
pdf(file='Figure OR3.pdf', height=120, width=120, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=10)
forest(res, xlim="s",comb.fixed=FALSE,xlab="Odds of Cacthability",text.fixed="Pooled effect estimate")
dev.off()

metainf(res, pooled="random")

#############################################################################