Extracting standard errors (cont’d) : critique

In a previous post I presented a possible solution to the extraction of standard errors and hazard ratios from publications in which only rounded approximations to the risk ratios and the associated 95% confidence interval are considered.

This solution, is subject to a criticism that I will now discuss : it suffers from an internal contradiction. Specifically, while one can use two of the three pieces of data (the actual risk ratio and the two limits of the 95% CI) simultaneously, one will arrive at different and conflicting probability estimates which is clearly not what one wants!.

To show the inconsistency, we need to see that the three pieces of data H, H_{U}, H_{L} are algebraic manipulations of two random variables H, Y=exp( q_{0.975}\times se) which are interrelated via the following equations H_{U}=H \times Y, H_{L}=H / Y.

When one tries to use standard probability calculus (or rudimentary Monte Carlo) to derive the distributions of H,Y from \{H,H_{U}\} or \{H,H_{L}\} or even \{H_{L},H_{U}\} one will arrive at joint (and marginal) distributions with different support and shapes. This is illustrated in the figure below, in which yellow is the area of the joint probability density function (PDF) that assumes non-zero values, black is the area in which the PDF is zero and white is the white space that makes the figures look nice 🙂

In all three figures the blue cross marks the true H,Y used to generate the figures.

What this little exercise tells us, is that one cannot hope to arrive at an (unconditional) probability model that is internally consistent. Even though this is possible when only pairs of H, H_{U}, H_{L} are considered, the three solutions do not agree with each other. Hence, one would either need to define one’s model directly on the basis of H,Y or H,se or use the specifications discussed in the previous post in a differential fashion.
Stay tuned, more to follow 🙂

EDIT (February 9th)

The three panel image creates the (false) visual impression that the three joint PDF have the same support on the x-axis, which is definitely NOT the case. The following image corrects this impression by plotting the PDFs over the same rectangular area in the plane.


Tags: ,

2 Responses to “Extracting standard errors (cont’d) : critique”

  1. Extracting standard errors and effect estimates for meta-analysis: paging Rev Bayes | Statistical Reflections of a Medical Doctor Says:

    […] is about probability and stats as applied to solve patient problems and other stuff! « Extracting standard errors (cont’d) : critique An algebra guy’s take on the meta-analysis posts […]

  2. Extracting standard errors and treatment effects from medical journal tables (powered by R) | Statistical Reflections of a Medical Doctor Says:

    […] stated previously (here and here) the problem is to convert the reported relative risks(RR, ), 95% confidence interval () and […]

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

%d bloggers like this: