I adapted a method designed by Ioannidis and Trikalinos, which compares the observed number of positive studies in a meta-analysis with the expected number, if the summary measure of effect, averaged over the individual studies, were assumed true. Excess in the observed number of positive studies, compared to the expected, is taken as evidence of publication bias. The observed number of positive studies, at a given level for statistical significance, is calculated by applying Fisher's exact test to the reported 2x2 table data of each constituent study, doubling the Fisher one-sided P-value to make a two-sided test. The corresponding expected number of positive studies was obtained by summing the statistical powers of each study. The statistical power depended on a given measure of effect which, here, was the pooled odds ratio of the meta-analysis was used. By simulating each constituent study, with the given odds ratio, and the same number of treated and non-treated as in the real study, the power of the study is estimated as the proportion of simulated studies that are positive, again by a Fisher's exact test. The simulated number of events in the treated and untreated groups was done with binomial sampling. In the untreated group, the binomial proportion was the percentage of actual events reported in the study and, in the treated group, the binomial sampling proportion was the untreated percentage multiplied by the risk ratio which was derived from the assumed common odds ratio. The statistical significance for judging a positive study may be varied and large differences between expected and observed number of positive studies around the level of 0.05 significance constitutes evidence of publication bias. The difference between the observed and expected is tested by chi-square. A chi-square test P-value for the difference below 0.05 is suggestive of publication bias, however, a less stringent level of 0.1 is often used in studies of publication bias as the number of published studies is usually small.

Author | Simon Thornley |

Date of publication | 2013-11-21 06:48:21 |

Maintainer | Simon Thornley <sithor@gmail.com> |

License | GPL-3 |

Version | 1.0 |

**BMort:** BMort

**ChisqTest_expect:** Chi-square test to test for significant difference between...

**expected_events:** calculate the number of expected events from a study included...

**plot_chase_observed_expected:** From a meta-analysis, analyse for publication bias....

**test.n:** This runs test_one function for many replicates, usually...

**test.n.treated:** Observed number of positive by testing for significance from...

**test_one:** For one constituent study and one simulation of its outcome,...

**test.one.treated:** Is one constituent study observed significant, in favour of...

PubBias

PubBias/NAMESPACE

PubBias/data

PubBias/data/BMort.rda

PubBias/R

PubBias/R/process.R
PubBias/MD5

PubBias/DESCRIPTION

PubBias/man

PubBias/man/test.one.treated.Rd
PubBias/man/ChisqTest_expect.Rd
PubBias/man/test_one.Rd
PubBias/man/test.n.treated.Rd
PubBias/man/BMort.Rd
PubBias/man/test.n.Rd
PubBias/man/plot_chase_observed_expected.Rd
PubBias/man/expected_events.Rd
Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.