As I was reading through the meta-analysis posts in order to correct various typos, the forgotten non-probability me woke up and raised the following question:

**What if** one were to treat the reported RR (), 95% confidence interval () and p-value () as the true values of the non-reported quantities, in essence ignoring the round-off error?

Could this lead to a (?simpler) solution bypassing the need for Monte Carlo? What this solution would look like and how it differs (implementationally) from the Bayesian one ? More importantly how does it hold up against the Bayesian solution?

The solution involves solving algebraically for the various quantities; if more than one estimate can be obtained (e.g. the log-relative risk can be obtained from the confidence interval as well as the reported RR) then one averages over these estimates. So the solution looks something like this:

- First we obtain an estimate for from the reported RR:
- Secondly we transform the confidence interval to estimates for :

- Averaging over yields
- A second estimate for the is obtained from the p-value using the quantile function of the normal distribution:
- Averaging over yields

Of note, the second step is used in the Monte Carlo algorithm when generating samples for the quantities from the confidence interval, so that the algorithm solution is living in the Bayesian one. Furthermore, the algebraic solution utilizes statistical concepts (the mean) to arbitrate between alternative estimates of the same quantities.

As the algebraic solution does indeed look simple enough to be programmed in Excel š one decide to use this rather than the complicated Monte Carlo one. So in a subsequent post we will pit these two solutions in a “death match” to decide whether the extra complexity of the Bayesian solution is worth the extra programming effort.

### Like this:

Like Loading...

*Related*

Tags: Bayesian, Meta-Analysis, statistics

This entry was posted on April 29, 2013 at 5:46 am and is filed under Bayesian, meta-analysis. You can follow any responses to this entry through the RSS 2.0 feed.
You can leave a response, or trackback from your own site.

November 10, 2013 at 10:05 pm |

[…] in order to recover the output of the statistical software. The simplest approach is one based on algebraic manipulationsĀ and disregards the finite precision of the medical journal table, treating the reported figures as […]