I adapted a method designed by Ioannidis and Trikalinos, which compares the observed number of positive studies in a metaanalysis 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 onesided Pvalue to make a twosided 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 metaanalysis was used. By simulating each constituent study, with the given odds ratio, and the same number of treated and nontreated 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 chisquare. A chisquare test Pvalue 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.
Package details 


Author  Simon Thornley 
Date of publication  20131121 06:48:21 
Maintainer  Simon Thornley <sithor@gmail.com> 
License  GPL3 
Version  1.0 
Package repository  View on CRAN 
Installation 
Install the latest version of this package by entering the following in R:

Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.