 Research article
 Open Access
 Open Peer Review
 Published:
Calculating unreported confidence intervals for paired data
BMC Medical Research Methodology volume 11, Article number: 66 (2011)
Abstract
Background
Confidence intervals (or associated standard errors) facilitate assessment of the practical importance of the findings of a health study, and their incorporation into a metaanalysis. For paired design studies, these items are often not reported. Since the descriptive statistics for such studies are usually presented in the same way as for unpaired designs, direct computation of the standard error is not possible without additional information.
Methods
Elementary, wellknown relationships between standard errors and pvalues were used to develop computation schemes for paired mean difference, risk difference, risk ratio and odds ratio.
Results
Unreported confidence intervals for large sample paired binary and numeric data can be computed fairly accurately using simple methods provided the pvalue is given. In the case of paired binary data, the design based 2 × 2 table can be reconstructed as well.
Conclusions
Our results will facilitate appropriate interpretation of paired design studies, and their incorporation into metaanalyses.
Background
Not too long ago, analyses of research data comprised performing a large number of unplanned hypothesis tests, and reporting the results simply as p < 0.05 or not, often in a selective way. Described as the cult of statistical significance, this practice frequently produced flawed interpretations, and made judging their practical importance difficult. Extensive criticism has now put it out of favor [1]. Reporting guidelines and many journals currently require joint presentation of the three key statistical indices, namely, point estimate, actual pvalue, and confidence interval (CI), especially for the principal analyses of the study [2, 3]. The CI (or associated standard error) enables interpretation of the main effect measure in a more complete manner, and is needed for incorporating the study results into a metaanalysis.
Various remedies for different instances of incomplete or varied reporting of essential findings (like confidence intervals) and key data have been proposed. For independentgroups 2 × 2 tables, Pietrantonj [4] detailed and evaluated a series of methods for reconstructing the unreported data table in risk difference, risk ratio, and odds ratiobased analyses. For a factor variable presented in quantile groups, Chêne and Thompson [5] described how findings from diverse forms of analyses can be reexpressed in terms of the mean difference. Abrams et al. [6] looked at clinical trials reporting quantitative change from baseline and gave Bayesian and sensitivity analyses methods to address partial reporting. Focusing on twoperiod crossover trials, Elbourne et al. [7] presented simple formulas linking the relevant standard deviation and correlation, for both continuous and binary data. These methods can be used to estimate one when the other is given. Tang et al. [8] described and evaluated methods to find the variance of the difference between paired proportions using a specialized variance recovery method that utilizes the variances of the individual proportions.
Our paper is located within this broad field dealing with incomplete reporting of data or analysis results. We focus on pair matched designs, specifically matched and crossover clinical trials, paired cohort studies, and 1:1 matched casecontrol studies. We first note that for such studies, confidence intervalrelated reporting shortfalls persist. Mills et al. [9] scrutinized a sample of 526 randomised controlled trials. Among these, 116 or 22% were crossover trials. Of the latter, only 30% presented a CI or the standard error, 62% did not present a CI but gave enough information to enable its computation, and the remaining 8% did not report it and did not give data enabling its computation. Their paper, however, does not specify how the computability of the CI was judged. Poor reporting in crossover trials was also underscored by Elbourne et al. [7].
We assessed the current state of CI reporting in paired epidemiologic studies, albeit in a preliminary way, with a systematic survey of the PubMed database. It was searched on April 16, 2010 using the term "matched casecontrol study", then sorted by "Recently Added". Studies that classified themselves in the abstract as 1:1 matched casecontrol studies were potentially eligible. If further scrutiny revealed the actual matching not to be 1:1, the study was excluded. The first 20 studies satisfying our criteria were selected. The median number of pairs in these studies was 196 (range 7 to 42 542 pairs), and one fifth had fewer than 50 pairs. The full list of these studies is available from the authors.
Twelve studies were true casecontrol studies and eight had a pair matched design but it was not a casecontrol type. All of the 12 true casecontrol studies reported CIs, but three used inappropriate unpaired data methods. Only two of the other eight studies used paired data methods. None of them reported a CI, but seven gave sufficient data to compute it. Overall, in the twenty selflabeled pair matched casecontrol studies, correct calculation and reporting of CI was done in 9/20 (45%), and a correct CI was computed or was directly computable in 16/20 (80%).
When a paper lacks a CI or the relevant standard error, one option is to contact the authors. If that does not yield the needed data, or if the process is too timeconsuming, is there an alternative? Our paper presents simple methods, hitherto underutilized or unpublished, that can be used for this purpose.
Methods
It is helpful to note that the prevalent style of reporting descriptive statistics makes the problem of noncomputability of CIs more acute for paired designs than for independent samples designs. A comparison of two proportions usually reports the sample proportions and sample sizes, and of two means, the sample means, sizes and standard deviations, whether the design is unpaired or paired. For the former design, these quantities respectively suffice to compute the CIs for effect measures like risk difference, risk ratio, odds ratio or the difference of means. In the case of paired data, they do not. For paired binary data, we need the 2 × 2 table with the concordant and discordant pairs, and for paired continuous data, we need the mean difference and standard deviation of the differences. These entities, in part or full, are rarely reported, even when they were used to compute a pvalue via appropriate tests like the McNemar's test, paired ttest, paired ztest, or a paired exact test. If only group level proportions (or means and standard deviations) are given, additional measures, such as betweengroup correlation, are required to determine the relevant standard error [7, 8]. Such entities are hardly reported.
We deal with four common effect measures (mean difference, risk difference, risk ratio and odds ratio) for paired design studies where the CI has not been given, and the design based 2 × 2 table or the standard deviation of the differences (as the case may be) and the relevant correlation are also unknown. For each measure, we show that if the pvalue is known, the corresponding standard error and CI can be obtained from a simple computation scheme. In the case of binary data, the data table with the concordant and discordant pairs can be reconstructed as well.
Results
We first introduce the main example we use to illustrate our methods. Xie et al. [10] reported a multicenter study of critically ill surgical patients with severe sepsis (SS). Its main aim was to assess if the presence of invasive fungal infection (IFI) affected the outcome for such cases. The subjects were drawn from the surgical intensive care units of ten teaching hospitals in China. All admissions in the one year study period meeting the set criteria for sepsis were included. The data on patient characteristics, treatments, and outcomes were compiled through daily chart reviews and physician interviews.
The main study had 90 SS patients with IFI and 228 SS patients without IFI. In a substudy, 60 of the SS patients with IFI were matched, on a onetoone basis and in terms of center, sex, age and APACHE II score, with 60 SS patients without IFI. We consider the matched portion of this study, and show a part of the results in Table 1.
Paired continuous data
One variable (hospital LOS) in Table 1 is continuous; the other (mortality) is binary. For now, consider the former. The paper reports the group interquartile ranges for hospital LOS but not the standard deviations. Also, the CI for the difference in hospital LOS is not given. However, the pvalue is stated with two significant digits. Using this, we impute the needed CI as follows.
For comparing hospital LOS, the two sided p = 0.020. The associated standard normal deviate is z = 2.326. In general, let be the mean of the differences, and s _{ z } , the standard normal paired twosided ztest standard error. The mean of the differences is the difference of the means, so . From the pvalue, we get z, the corresponding twosided deviate of the standard normal distribution. Then we apply the relation
The 95% CI for the difference of the means thereby is
Suppose the paired twosided ttest was used. From the pvalue, we get t, the corresponding deviate of the t distribution with n  1 degrees of freedom. With s _{ t } as the corresponding standard error for the difference of the means, we have
The 95% CI for the difference of the means is then
where t _{ n 1 }is the 97.5th percentile of the t distribution with n  1 degrees of freedom.
Further, the standard deviation of the differences of the means is obtained by using one of the two formulas, as appropriate, given below
Applying these formulas to the hospital LOS data, we get, under the ztest, that s _{ z } = (30  20)/2.326 = 4.299 with 95% CI equal to (1.57, 18.4), and s _{ d } = 33.3. Under the ttest, we get that s _{ t } = 4.182 with 95% CI equal to (1.63, 18.4), and s _{ d } = 32.4. These 95% intervals are almost identical. Also both are quite wide, with their lower limits not far from zero days. The possibility of just a minor difference in the hospital LOS for the two groups cannot thereby be excluded.
Paired binary data
First consider paired binary data under a prospective design. The data format and the cellwise and marginal proportions are shown in Table 2. The marginal proportions are estimated by
Risk difference
Researchers usually apply two effect measures for paired prospective designs, risk difference and risk ratio. First consider the former, δ = π _{1}  π _{2}. It is estimated by
For this measure, two different variance formulas are usually used, one for null hypothesis testing and the other for confidence interval computation. These formulas are shown in the first two rows of Table 3[11–14].
Suppose the pvalue for the twosided ztest for risk difference equal to zero is known. This gives the associated standard normal deviate, z, from which we determine the null standard error as . With , , and s _{0} known, the variance formula from the second row of
Table 3 provides the three equations
Since a + b + c + d = n, we can solve for a, b, c and d. The corresponding solution scheme is in the first row of Table 4. All the numbers are rounded to the nearest integer.
We apply these formulas to the mortality data from Xie et al. Since p = 0.023, then z = 2.27. From Table 1, we find = 0.7  0.5 = 0.2. Thus, s _{0} = 0.2/2.27 = 0.088. Using these with n = 60 in the scheme of the first row of Table 4 gives a = 10, b = 8, c = 20, and d = 22. Thereby, we learn that the number of pairs in which both subjects died were 2.2 times as many as those in which both subjects remained alive.
After reconstructing the paired 2 × 2 table, we can check whether we are able to reproduce the pvalue for hospital mortality given in Table 1. Xie et al. used the McNemar test, which has test statistic equal to
to be compared with the standard normal distribution. With b = 8 and c = 20, we get T = 2.268 and p = 0.023, which is equal to the pvalue given by Xie et al.
To get the 95% CI for the risk difference, we first apply the formula in the first row of Table 3 to the reconstructed table to get . Hence the required interval is which is 0.2 ± 1.96 × 0.084 = (0.035, 0.37).
The number needed to treat (NNT) then is 1/.2 = 5, implying that for every 5 cases with SS & IFI, on average, one additional case will die as compared with those having SS & no IFI. The 95% CI for the NNT, obtained by inverting the CI for the risk difference and reversing the limits, is equal to (3, 29). This interval, like that of the risk difference, is not that precise.
Risk ratio
The risk ratio is θ = π _{2} /π _{1} which, for paired prospective data, is estimated by
The variance of the logarithm of is in the third row of Table 3. Now suppose the pvalue for the null test of the risk ratio is known. From this, we find the standard normal deviate, z.
The associated standard error then is . For this measure, the same variance formula is customarily used for null hypothesis testing and for confidence interval computation. Hence, we directly compute the 95% CI for ln θ as , and exponentiate this interval to get the 95% CI for θ.
From the variance formula, we construct four simultaneous equations in the same way as for the risk difference. The corresponding solution scheme for a, b, c and d in this case is shown in the second row of Table 4. Here also, these numbers are rounded to the nearest integer.
In the above calculations, we assumed that the published pvalue was calculated using the correct standard error for paired risk ratio. In practice, this may not occur that often, especially when the CI is unreported. If instead the pvalue was calculated using McNemar's test, we reconstruct the paired 2 × 2 table using the methods for risk difference, and then apply these numbers to get the variance estimate for the log risk ratio in Table 3. We illustrate this approach in a hypothetical situation.
Suppose we need to metaanalyse several studies with the risk ratio as the effect measure. Suppose Xie et al. [10] is one of the selected studies. The above computations then allow us to include it into this metaanalysis even though it does not report the risk ratio, and the risk ratio or its standard error are not directly computable from the paper. Using the paired data table reconstructed above (a = 10, b = 8, c = 20, and d = 22), the risk ratio is 30/42 = 0.71, and the standard error of the logrisk ratio is . We exponentiate the logarithmic interval ln (0.71) ± 1.96 × 0.15 to find the 95% CI for the risk ratio as (0.53, 0.95). This CI also indicates a low precision for the study results. The risk ratio and its standard error are now available for the metaanalysis.
Odds ratio
Now consider the odds ratio, which is mostly used in casecontrol studies. The marginal proportions in the paired data table (Table 1) are the probabilities of exposure given the outcome. Such studies usually apply a conditional form of analysis that uses only the discordant pairs whose sum is taken as fixed. The appropriate conditional probabilities for the two types of discordant pairs are π _{12}/(π _{12} + π _{21}) and π _{21}/(π _{12} + π _{21}). Further, the relevant odds ratio is ρ = π _{21}/π _{12}, which is estimated by . The (conditional) variance of the logarithm of is shown in the last row of Table 3.
For the task of reconstructing the paired 2 × 2 table, the odds ratio differs from the other two binary effect measures in an important way. Unlike for the latter, the paired data table can be computed once we know the odds ratio, the marginal proportions and the total number of pairs. This computation scheme appears in the third row of Table 4.
Consider a hypothetical casecontrol study where the odds ratio is reported , but the CI is not. Also reported are the marginal proportions ( and ) and the total number of pairs (n = 80). Using the third row of Table 4, we reconstruct the paired 2 × 2 table:
We then calculate the standard error of the log odds ratio from the last row of Table 3: . The 95% CI for the log odds ratio is ln (2.0) ± 1.96 × 0.32. After exponentiation, we find that the 95% CI for the odds ratio is (1.1, 3.7).
When only one marginal proportion is known, the appropriate pvalue is also needed. This then allows us to compute the standard error . Like the risk ratio, the same variance formula is customarily used for null hypothesis testing and for confidence interval computation, as given in the last row of Table 3. Using this, we form four simultaneous equations as done earlier. Solving these, we obtain the scheme shown in the last row of Table 4 for reconstructing the needed data table.
Suppose that in the hypothetical example above, only one of the marginal proportions was given but that p = 0.029 was reported. The standard normal deviate here is z = 2.18, and the standard error, s = 0.318, is obtained in the usual way. The 95% CI is computed as above. To reconstruct the paired 2 × 2 table, we use the last row of Table 4 to get b = (1 + 2)/0.318^{2} ≈ 30 and c = (1 + 2)/2 · 0.318^{2} ≈ 15. The quantities d = 10 and a = 25 are obtained as before.
Impact of pvalue accuracy
One important consideration is the accuracy of the stated pvalue. A small absolute change in a small p produces a large change in the z or t deviate. If the pvalue has not been stated accurately, our computation schemes can yield flawed answers. To get an initial handle on the error involved here, we performed a sensitivity analysis for the data in Table 1. For both comparisons, the pvalue was changed from 0.015 to 0.024 in increments of 0.001. For hospital LOS, the ztest scheme gave 95% CIs ranging from (1.94, 18.06) to (1.32, 18.68), and the ttest scheme gave 95% CIs ranging from (2.01, 17.99) to (1.36, 18.64). And for the same range of pvalues, the computed 95% CI for the risk difference in mortality ranged from (0.048, 0.352) to (0.035, 0.365). These changes are neither dramatic nor practically meaningful. All the intervals are close to the corresponding intervals for the observed pvalues. Our experience thus far indicates that if the pvalue is known to two significant digits, the results are sufficiently accurate, and often, one significant digit of accuracy suffices. Detailed simulation studies to resolve the concerns relating to sample size, data structure, and degree of accuracy of the pvalue are, however, called for.
Discussion
Our paper gives easy to apply computation schemes to compute confidence intervals and other entities in situations where the needed information is unavailable. For such tasks, our paper is a paireddata counterpart of Pietrantonj [4] that addresses similar problems for unpaired binary data. Note that the computational schemes we give are simpler than their unpaired variants. Further, our paper differs from previous work for paired data cited earlier in that it is not based on the knowledge of a measure of correlation, a rarely reported entity, but requires the pvalue, a commonly reported item.
The use of our computation schemes is limited by several considerations. First, for continuous data, we need to know if the paired z or t test was used to obtain the pvalue. Second, the computation scheme for the risk ratio applies only if the pvalue was obtained by the use of the risk ratio standard error. This is rarely done. pvalues for paired binary data are usually computed using a risk difference based (McNemar) test. If we then need the CI for the risk ratio, we use the risk difference scheme to reconstruct the paired data table and obtain the required standard error from this table. Third, as the paired data table can often be reconstituted with the odds ratio, the pvalue based scheme here will be rarely, if at all, used. We give it for completeness.
Other matters of concern are data structure and sample size. The data are sparse when one discordant pair is much smaller in size than the other, or when both corresponding marginal proportions are near one or zero. Suppose the given pvalue has been obtained by either the ztest or ttest for continuous data, or a test based on one of the null standard errors shown in Table 3 for binary data, and the paper notes the actual test used. In each of these cases, the computation scheme we give is valid at all sample sizes and with all data structures, whether sparse, skewed or otherwise. Validity here means that it will give the same confidence intervals as that based on the original data.
Besides these standard tests, a wide variety of other methods for computing pvalues and CIs for paired discrete and continuous data exist [15, 16]. Some of these methods have better statistical properties compared with the others. And some have complex formulas. With sparse data or small number of pairs, different methods may give appreciably different results. For such data, researchers may also use nonparametric or exact methods employing special algorithms that are not amenable to simple formulations [17]. When such methods have been used in a study with small number of pairs or sparse data, and if the CI has not been reported or the data to compute it are not available, developing a procedure to correct the deficiency is not a simple matter. Applying the computation schemes in Table 4 to such data when the pvalue has been computed using a specialized test is not advisable.
However, simulation studies show that when viewed in terms of their practical impact, the CIs for the paired data risk difference computed by several methods are fairly close to one another when the number of pairs exceeds 50, and provided the data are not too sparse [15, 16]. Based on studies of this sort, we recommend our computation schemes even when the pvalue has been computed using a nonparametric, exact, score, or some other test provided there are more than 50 pairs and the data are not too skewed or sparse. The error involved would be, we suspect, acceptable for the point of view of practical interpretation or incorporation into metaanalysis.
Conclusion
Confidence intervals allow us to judge the practical implications of a study. For a paired design, the paired data table is more informative than the marginal data summary. The standard error of difference of the means incorporates the correlation between the two paired measures. Without this information, it cannot be computed directly from the standard errors of each measure. In crossover trials, for example, it enables us to examine the presence of a treatment order effect. Reconstructing the paired data table is useful when a study analyzed in terms of the risk difference is not fully reported and incorporating it into a metaanalysis requires the risk ratio or odds ratio. That reconstruction also permits a reanalysis of the data with methods that have better statistical properties, including for performing a metaanalysis.
Our methods are valid for all types of data when the pvalue has been computed using one of the standard tests we consider, and may have reasonable accuracy even for other tests provided the number of pairs exceeds 50 and the data are not too sparse or skewed. Simulation studies to identify the types of tests, sample sizes, data structures and levels of accuracy of the pvalue under which they are acceptably accurate are, however, warranted.
Better reporting of paired data is, nonetheless, the optimal solution. When reporting a study, the point estimate, pvalue, and confidence interval for the main effect measure must be given. For paired binary designs, the full design based 2 × 2 table should be given, and for paired continuous data, the mean and standard deviation for each sample should be augmented by the standard deviation of the differences. Reporting the relevant correlation is an equivalent substitute. Improved reports will make corrective schemes like the ones we give somewhat superfluous. Until the day when such reports are the almost universal norm, however, these schemes will serve a useful purpose.
Abbreviations
 CI:

confidence interval
 IFI:

invasive fungal infection
 LOS:

length of stay
 NNT:

number needed to treat
 SS:

severe sepsis.
References
 1.
Ziliak ST, McCloskey DN: The Cult of Statistical Signficance. 2007, Ann Arbor, MI: The University of Michigan Press
 2.
Altman DG, Schulz KF, Moher D, Egger M, Davidoff F, Elbourne D, Gtzche PC, Lang T, for the CONSORT group: The revised CONSORT statement for reporting randomized trials: explanation and elaboration. Annals of Internal Medicine. 2001, 134: 663694.
 3.
Vandenbroucke JP, von Elm E, Altman DG, Gtzche PC, Mulrow CD, Pocock SJ, Poole C, Schlesselman JJ, Egger M, for the STROBE Initiative: Strengthening the reporting of observational studies in epidemiology (STROBE): explanation and elaboration. PLoS Medicine. 2007, 4 (10): e29710.1371/journal.pmed.0040297.
 4.
Pietrantonj CD: Fourfold table cell frequencies imputation in metaanalysis. Statistics in Medicine. 2006, 25: 22992322. 10.1002/sim.2287.
 5.
Chêne G, Thompson SG: Methods for summarizing the risk associations of quantitative variables in epidemiologic studies in a consistent form. American Journal of Epidemiology. 1996, 144: 610621.
 6.
Abrams KR, Gillies CL, Lambert PC: Metaanalysis of heterogeneously reported trials assessing change from baseline. Statistics in Medicine. 2005, 24: 38233844. 10.1002/sim.2423.
 7.
Elbourne DR, Altman DG, Higgins JPT, Curtin F, Worthington HV, Vail A: Metaanalyses involving crossover trials: methodologic issues. International Journal of Epidemiology. 2002, 31: 140149. 10.1093/ije/31.1.140.
 8.
Tang ML, Ling MH, Ling L, Tian G: Confidence intervals for a difference between proportions based on paired data. Statistics in Medicine. 2010, 29: 8696.
 9.
Mills EJ, Chan AW, Wu P, Vail A, Guyatt GH, Altman DG: Design, analysis, and presentation of crossover trials. Trials. 2009, 10: 2710.1186/174562151027.
 10.
Xie GH, Fang XM, Fang Q, Wu XM, Jin YH, Wang JL, Guo QL, Gu MN, Xu QP, Wang DX, Yao SL, Yuan SY, Du ZH, Sun YB, Wang HH, Wu SJ, Cheng BL: Impact of invasive fungal infection on outcomes of severe sepsis: a multicenter matched cohort study in critically ill surgical patients. Critical Care. 2008, 12: R510.1186/cc6766.
 11.
Altman DG, Machin D, Bryant TN, Gardner MJ, (eds): Statistics with confidence. 2000, BMJ Books, 2
 12.
Fleiss JL, Levin B, Paik MC: Statistical Methods for Rates and Proportions. 2003, Hoboken, NJ: John Wiley & Sons, Inc, 3
 13.
Væth M, Poulsen C: Comments on a commentary: statistical evaluation of split mouth carries trials. Community Dentistry and Oral Epidemiology. 1998, 26: 8083. 10.1111/j.16000528.1998.tb01930.x.
 14.
Rothman KJ, Greenland S: Modern Epidemiology. 1998, Philadelphia, PA: Lippincott Williams & Wilkins, 2
 15.
May WL, Johnson WD: Confidence intervals for difference in correlated binary proportions. Statistics in Medicine. 1997, 16: 21272136. 10.1002/(SICI)10970258(19970930)16:18<2127::AIDSIM633>3.0.CO;2W.
 16.
Newcombe RG: Improved confidence intervals for the difference between binomial proportions based on paired data. Statistics in Medicine. 1998, 17 (22): 26352650. 10.1002/(SICI)10970258(19981130)17:22<2635::AIDSIM954>3.0.CO;2C.
 17.
Hirji KF: Exact Analysis of Discrete Data. 2006, Boca Raton, FL: Chapman and Hall/CRC Press
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712288/11/66/prepub
Acknowledgements
We thank the three referees for helpful comments on the initial draft.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
KFH conceived the study and wrote the initial and final drafts. MWF performed the literature survey, did additional computations, gave comments and made additions to the final draft. Both authors read and approved the final manuscript.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Received
Accepted
Published
DOI
Keywords
 Risk Ratio
 Invasive Fungal Infection
 Computation Scheme
 Risk Difference
 Variance Formula