Publication bias

par(mfrow=c(1,2))
funnel(result, atransf=exp, xlab="Relative Risk")
funnel(result, level=c(90, 95, 99), shade=c("white","gray","darkgray"),
refline=0, atransf=exp, xlab="Relative Risk", at=log(c(.25, .5, 1, 2, 4)), yaxis="seinv")

베거 검정

ranktest(result)

## 
## Rank Correlation Test for Funnel Plot Asymmetry
## 
## Kendall's tau = 0.0256, p = 0.9524

에거 검정

regtest(result)

## 
## Regression Test for Funnel Plot Asymmetry
## 
## model:     mixed-effects meta-regression model
## predictor: standard error
## 
## test for funnel plot asymmetry: z = -0.7148, p = 0.4747

H0 : Publication bias가 없다.

H1 : Publication bias가 있다.

->Publication bias가 없다.

forest(result,atransf=exp)