4
$\begingroup$

I'm certain there's a unique solution to this, and I think I've worked it out before but now it has me pulling my hair out:

Given the margins of a 2x2 contingency table such as the prevalences of a binary exposure and a binary outcome, as well as their odds ratio, how can you calculate the four cell proportions (i.e. the reverse of the usual problem of calculating an odds ratio)? e.g.

                outcome     no outcome
exposed         a           b               Pr(E)
not exposed     c           d               Pr(Ê)=1-Pr(E)
                Pr(O)       Pr(Ô)=1-Pr(O)   100%

given the equations:

a + b = Pr(E)
b + c = Pr(Ê) = 1 - Pr(E)
a + c = Pr(O)
b + d = Pr(Ô) = 1 - Pr(O)
(ad) / (bc) = OR

and the constraints that all cell values and odds ratio are positive, and that a + b + c + d sum to 100% (i.e. everything is expressed as as proportion, although you can expand the equations to include an arbitrary known population count if you wish) express each cell value in terms of Pr(E), Pr(O) and OR.

$\endgroup$
1
5
$\begingroup$

Write $\rho$ for the odds ratio, $\beta=\Pr(E)$, $\gamma=\Pr(O)$. Four independent equations are

$$\cases{a+b=\beta \\ a+c=\gamma \\ a+b+c+d=1 \\ ad = \rho bc.}$$

Adding the first two shows

$$b+c = \beta + \gamma - 2a. \tag{1}$$

Multiplying the first two equations and using $(1)$ yields

$$bc = \beta\gamma - (\beta+\gamma)a + a^2.\tag{2}$$

Multiplying the third equation by $a$, using the fourth to re-express $ad$ and plugging in $(1)$ and $(2)$ gives

$$a^2 + a(\beta+\gamma-2a) + \rho(\beta\gamma - (\beta+\gamma)a + a^2) = a.$$

In a more standard form, this a zero of the quadratic

$$(\rho-1)a^2 + [(\beta+\gamma)(1-\rho)-1]a + \rho\beta\gamma.\tag{3}$$

Provided $\beta, \gamma,\rho$ are consistent with some $2\times 2$ table, there will be at least one real zero of $(3)$, easily found using any quadratic formula. For either zero, solutions for the remaining entries are readily found from the first three of the original equations as

$$\cases{b=\beta-a\\c=\gamma-a\\d = 1+a-\beta-\gamma.}$$


There will be at most one valid solution for $a$, determined by the non-negativity of the coefficients. Here is an R implementation of the solution as the function f along with a test using randomly generated tables. The test outputs the random table, its reconstructed value from f, and a measure of the difference between them. By wrapping the test in replicate, I have run it 10,000 times. The final output gives the largest difference found: up to floating point error it equals zero, demonstrating the correctness of this approach.

f <- function(beta, gamma, rho, eps=1e-15) {
  a <- rho-1
  b <- (beta+gamma)*(1-rho)-1
  c_ <- rho*beta*gamma

  if (abs(a) < eps) {
    z <- -c_ / b
  } else {
    d <- b^2 - 4*a*c_
    if (d < eps*eps) s <- 0 else s <- c(-1,1)
    z <- (-b + s*sqrt(max(0, d))) / (2*a)
  }
  y <- vapply(z, function(a) zapsmall(matrix(c(a, gamma-a, beta-a, 1+a-beta-gamma), 2, 2)),
              matrix(0.0, 2, 2))
  i <- apply(y, 3, function(u) all(u >= 0))
  return(y[,,i])
}
set.seed(17)
i<-0
sim <- replicate(1e4, {
  while(TRUE) {
    x <- matrix(round(rexp(4), 2), 2, 2)
    if(all(rowSums(x) > 0) && all(colSums(x) > 0) && x[1,2]*x[2,1] > 0) break
  }
  x <- x / sum(x)
  beta <- rowSums(x)[1]
  gamma <- colSums(x)[1]
  rho <- x[1,1]*x[2,2] / (x[1,2]*x[2,1])

  y<- f(beta, gamma, rho)
  delta <- try(zapsmall(c(1, sqrt(crossprod(as.vector(x-y)))))[2])
  if ("try-error" %in% class(delta)) cat("Error processing ", x, "\n")
  delta
})
max(sim)
$\endgroup$
3
  • $\begingroup$ Thank you! Realising that it could be solved with a quadratic formula was exactly what I needed. $\endgroup$
    – JayDee
    Oct 24 '16 at 9:42
  • $\begingroup$ To summarise the above for those less technically minded, the steps are: $\endgroup$
    – JayDee
    Oct 25 '16 at 8:38
  • $\begingroup$ 1: Use simultaneous equations to solve for any one cell value (combining OR formula with sum formulae) 2: The solution will be some horrendous non-linear term, but you can express it as a quadratic equation, i.e. Ax^2 + Bx + C = 0, where x represents one of your cell values. Plug this into the quadratic formula: x = (-B - sqrt(B^2 - 4AC)) / 2A to solve for x. Note the sign before the sqrt is a minus, rather than the usual plus/minus. This keeps the value of x positive. 3: plug the cell value represented by x back into your simple linear sum equations to solve the rest of the cell values. $\endgroup$
    – JayDee
    Oct 25 '16 at 8:48

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.