 Research article
 Open access
 Published:
Which is more generalizable, powerful and interpretable in metaanalyses, mean difference or standardized mean difference?
BMC Medical Research Methodology volumeĀ 14, ArticleĀ number:Ā 30 (2014)
Abstract
Background
To examine empirically whether the mean difference (MD) or the standardised mean difference (SMD) is more generalizable and statistically powerful in metaanalyses of continuous outcomes when the same unit is used.
Methods
From all the Cochrane Database (March 2013), we identified systematic reviews that combined 3 or more randomised controlled trials (RCT) using the same continuous outcome. Generalizability was assessed using the Isquared (I^{2}) and the percentage agreement. The percentage agreement was calculated by comparing the MD or SMD of each RCT with the corresponding MD or SMD from the metaanalysis of all the other RCTs. The statistical power was estimated using Zscores. Metaanalyses were conducted using both randomeffects and fixedeffect models.
Results
1068 metaanalyses were included. The I^{2} index was significantly smaller for the SMD than for the MD (Pā<ā0.0001, sign test). For continuous outcomes, the current Cochrane reviews pooled some extremely heterogeneous results. When all these or less heterogeneous subsets of the reviews were examined, the SMD always showed a greater percentage agreement than the MD. When the I^{2} index was less than 30%, the percentage agreement was 55.3% for MD and 59.8% for SMD in the randomeffects model and 53.0% and 59.8%, respectively, in the fixed effect model (both Pā<ā0.0001, sign test). Although the Zscores were larger for MD than for SMD, there were no differences in the percentage of statistical significance between MD and SMD in either model.
Conclusions
The SMD was more generalizable than the MD. The MD had a greater statistical power than the SMD but did not result in material differences.
Background
A metaanalysis aggregates indexes of effectiveness of individual trials into one pooled estimate. When the outcome of interest is a dichotomous variable, the commonly used effect sizes include the odds ratio (OR), the relative risk (RR), and the risk difference (RD). When the outcome is a continuous variable, then the effect size is commonly represented as either the mean difference (MD) or the standardised mean difference (SMD) [1].
The MD is the difference in the means of the treatment group and the control group, while the SMD is the MD divided by the standard deviation (SD), derived from either or both of the groups. Depending on how this SD is calculated, the SMD has several versions such, as Cohen's d [2], Glass's Ī [3], and Hedges' g [4].
When the outcome is measured in different units across trials, then we have no other choice but to use the SMD to combine the outcomes in the metaanalyses. On the other hand, when the outcome is measured in the same unit in every trial, theoretically, we can use either the MD or the SMD. In this latter case, there currently appears to be no unanimous agreement about which effect size is preferable, and different textbooks of metaanalyses provide differently nuanced recommendations about the selection of the appropriate effect size for continuous variables.
According to the Cochrane Handbook for Systematic Reviews of Interventions [1], the āselection of summary statistics for continuous data is principally determined by whether studies all report the outcome using the same scale when the MD can be used.ā The American Psychological Association (APA) Task Force on Statistical Inference maintains that āif the units of measurement are meaningful on a practical level (e.g. number of cigarettes smoked per day), then we usually prefer a MD to a SMDā [5]. Egger et al. writes that āthe overall treatment effect [in terms of SMD] can also be difficult to interpret as it is reported in units of standard deviation rather than in units of any of the measurement scales used in reviewā [6].
On the other hand, there are also authors who recommend the SMD along with or over the MD. The APA Publication Manual suggests that it can often be valuable to report not only the MD but also the SMD [7]. Borenstein, in his āIntroduction to MetaAnalysisā [8], wrote that if the unit is unfamiliar, the SMD serves as an easy way to judge the magnitude of the effect, thanks to the general rules of thumb described by Cohen that suggest that an SMD of 0.2 represents a āsmallā effect, an SMD of 0.5 represents a āmediumā effect, and an SMD of 0.8 represents a ālargeā effect [2]. For example, when you read that a treatment groupās mean posttreatment score on scale X was 10 points higher than that of a control group, there is no way of appreciating how much a difference this actually represents unless you are very familiar with the scale that is being used. But if the difference is expressed in terms of SMD as corresponding to an effect size of 0.5, for example, you can understand that it represents a moderate effectiveness in comparison with the control. In fact, Tian et al. noted that the SMD does not depend on the unit of measurement, and therefore the SMD has been widely used as a measure of intervention effect in many applied fields [9].
The preferability of the MD or the SMD can be examined from three aspects. First, which of the MD or the SMD is clinically more interpretable? The abovesummarised arguments made by different authors seem to concern mainly this aspect. Second, which of the MD or the SMD is more generalizable (as any summary index should have a good generalizability so that it can be applied to the next group of patients) [10]? Third, which of the MD or the SMD is statistically more powerful (as we expect metaanalyses to provide as precise an estimate of the treatment effect as possible and to be as sensitive as possible to differences among treatments)?
To the best of our knowledge, no systematic assessment focusing on the second or third aspects of the MD and the SMD has been conducted. The objective of this research was, therefore, to examine empirically which index is more generalizable and statistically powerful in metaanalyses when the same unit is used: MD or SMD?
Methods
Selection of metaanalyses
We included the following intervention metaanalyses from the Cochrane Database of Systematic Reviews (March 31, 2013) [11, 12].

1.
The outcome of the metaanalyses was a continuous outcome in the same unit.

2.
There were at least three studies contributing to that continuous outcome in question, because we need at least three studies in order to calculate the percentage agreement that we defined as one of our outcome (see below).

3.
If there were two or more outcomes that met the above criteria, we selected the outcome that contained the greatest number of studies in the review. If there were 2 or more such outcomes, we chose the outcome reported for the greatest number of patients. If the numbers of the patients were the same, we selected the outcome that appeared first in the review.
We excluded the outcomes reported in sensitivity analyses or subgroup analyses.
Search methods and data extraction
We searched the Cochrane Database of Systematic Reviews (March 31, 2013) using the term āmean differenceā in all fields to identify all metaanalyses possibly meeting the above eligibility criteria.
One author examined all the metaanalyses identified using the above search methods to decide whether the metaanalysis met our inclusion criteria and chose which outcome to focus on. In order to evaluate the reliability of the selections, we randomly selected 100 out of the original set of metaanalyses and the other authors examined the reliability of these selections. We then calculated the kappas for these selections.
Outcomes
To look at generalizability, we examined (i) the Isquared (I^{2}) index and (ii) the percentage agreement. To compare statistical power, we examined (iii) the Zscore.

(i)
I^{2} index
The I^{2} index represents the degree of heterogeneity across studies in metaanalyses [13]. It ranges from 0% to 100% and the following rough rule of thumb has been proposed [1]:
0% to 40%: might not be important
30% to 60%: may represent moderate heterogeneity
50% to 90%: may represent substantial heterogeneity
75% to 100%: considerable heterogeneity.
Because more heterogeneous results are expected to show less generalizability to similar populations, we examined the I^{2} index of each metaanalysed MD and SMD. We used the sign test to compare the I^{2} values for the MD and the SMD.

(ii)
Percentage agreement
One study was extracted from each metaanalysis. The MD and the SMD of that individual study were then compared with the metaanalytically pooled MD and the SMD of the remaining studies. Agreement was defined when the point estimate of the MD or the SMD of the individual study was included within the 95% confidence interval of the pooled MD or SMD of the remaining studies [10]. This procedure was repeated for all the studies, and the overall percentage agreement was calculated. We calculated two percentage agreement figures for each metaanalysis, using the randomeffects model [14] or the fixedeffect model [15] for the metaanalysis. In each model, we compared the percentage agreement of the MD with that of the SMD using the sign test.

(iii)
Zscore
We also examined the Zscore or each metaanalysis in order to examine possible differences in the statistical power between the MD and the SMD using the sign test. In addition, we compared the MD or SMD values that were judged to be statistically significant at a conventional level (Pā<ā0.05) using the McNemar test.
In order to calculate the percentage agreements, the I^{2} statistics, and the Pvalue for the treatment effect between the MD and the SMD, we used the inverse variance method metaanalyses that are included in SAS, Ver. 9.3 (SAS Institute Inc.). We defined SMD as Hedgesā adjusted g to remove any upward bias that might have arisen because of small sample sizes [16].
Results
Search results
FigureĀ 1 describes the process of our search. There were 5418 reviews in the Cochrane Database of Systematic Reviews (March 31, 2013). When we searched the database using the term āmean differenceā in all the texts, 3961 reviews were selected. Of these, 1068 reviews met our eligibility criteria. The main reason why the reviews were excluded was that the outcome of the metaanalyses contained two or fewer studies. Among the 1068 metaanalyses in our dataset, only 47 (4.4%) metaanalyses reported SMD, and 1021 (95.6%) metaanalyses reported MD. In order to confirm the reliability of this identification, two authors independently examined 100 reviews that were randomly selected from the 3961 reviews. The kappa was 0.90.
The median number of the studies contributing to each outcome was 5 (interquartile range: 3 to 8). The median total of participants used to determine the outcome was 489 (interquartile range: 218 to 1146). The median absolute SMDs for the randomeffects model and for the fixedeffect model were 0.29 (interquartile range: 0.12 to 0.56) and 0.32 (interquartile range: 0.14 to 0.62), respectively. The subject areas of the 1068 reviews were as follows: internal medicine, nā=ā404; obstetrics and gynaecology (including paediatrics), nā=ā234; surgery, nā=ā204; psychiatry, nā=ā140; anaesthesiology, nā=ā33; and others, nā=ā53.
Outcomes
I^{2}index
FigureĀ 2 shows a histogram containing the I^{2} data. We found, much to our surprise, that the review authors had pooled continuous outcomes as either the MD or the SMD, even when the I^{2} statistics were extremely high. Although the median I^{2} was 47.6% for MD and 48.0% for SMD, in a matched comparison, it was significantly lower for SMD than for MD (Pā<ā0.0001, sign test).
Percentage agreement
In the randomeffects model, the median of the percentage agreement in which the point estimate of one trial was included in the 95% confidence interval of the metaanalytically pooled estimate of the remaining trials was 46.7% of the reviews if we used the SMD and 44.5% if we used the MD (Pā<ā0.0001, sign test). In the fixedeffect model, the respective figures were 27.8% and 31.9% for the SMD and the MD, respectively (Pā<ā0.0001).
The same trend was observed when we recalculated the percentage agreements by simulating a situation where the review authors limited the pooling to moderately heterogeneous (I^{2} index of both MD and SMD, <60%) or to only slightly heterogeneous (I^{2} index of both MD and SMD, <30%) instances (TableĀ 1). In each scenario, the SMD always outperformed the MD in terms of percentage agreement: the randomeffects SMD showed the highest agreement of 59.8% when we limited the metaanalyses to cases for which the I^{2} index of both the MD and the SMD was less than 30%.
Zscore
If one of the statistics was statistically more heterogeneous than another, it would be meaningless to discuss their relative detection power. To discuss the statistical power, we compared the Zscore of the MD with that of the SMD using the sign test only when the I^{2} indices of both the MD and the SMD were less than 30%, i.e. less than moderate according to the proposed interpretative guide [1].
In the randomeffects model, the median absolute Zscore for the MD was significantly higher than that for the SMD (1.59 vs. 1.61, Pā=ā0.007, sign test). In the fixedeffect model, the median Zscore for the MD was again significantly higher than that for the SMD (1.60 vs. 1.63, Pā<ā0.0001). However, when the statistical significance was compared at a conventional Pvalue of less than 0.05, no differences in the percentage of āstatistically significantā findings were observed between MD or SMD in either the randomeffects or the fixedeffect model (Pā=ā1.00 and Pā=ā0.65, respectively, McNemar test) (TableĀ 2).
Discussion
We empirically examined 1068 reviews from the Cochrane Database of Systematic Reviews. The I^{2} index, an index of heterogeneity, was significantly smaller when the metaanalysed results were expressed in terms of the SMD, rather than the MD. When each one of the included RCT was compared against the pooled results of the remaining RCTs in each metaanalysis, the SMD showed a significantly greater percentage agreement than the MD in both the randomeffects and the fixedeffect models for all degrees of heterogeneity. On the other hand, no statistically significant difference was found in terms of statistical power for identifying significant results between the MD and the SMD in either the randomeffects or the fixedeffect model.
Our research agrees with one previous study regarding heterogeneity and statistical power. This previous study examined a relatively new index of effect, known as the ratio of means, for analysing continuous outcomes in comparison with more traditional effect sizes of continuous outcomes, namely the MD and the SMD. Consequently, the study examined only the metaanalyses where the ratio of means could be calculated, and it further limited the analyses to those metaanalyses containing five or more studies and examined the randomeffects model only. However, the study also found that the Pvalue did not statistically differ between the MD and the SMD and that the SMD was less heterogeneous (defined by Pā<ā0.1 for the Q statistic) than the MD [17, 18]. When we conducted sensitivity analyses using our dataset by including only the metaanalyses that contained five or more studies, the results were also similar.
The percentage agreement figures for continuous outcomes, either in terms of MD or SMD, found in our study appeared to be very low, even when the associated I^{2} index was reasonably low. Clinicians and researchers must therefore be advised to use the most generalizable index of effectiveness, while keeping in mind that the actual degrees of expected overlap may not be as high as one would expect. We found a 5percentage point difference in the degree of expected agreement when the results were expressed as the SMD vs. when they were presented as the MD; this difference was clinically meaningful and nonnegligible.
Much to our surprise, however, the Cochrane authors often pooled their continuous outcomes as the MD and/or the SMD, even when the I^{2} statistics suggested extreme heterogeneity. While it is true that metaanalyses of continuous outcomes tend to be associated with a greater I^{2} value than those of dichotomous outcomes because of the formerās greater statistical power, the generalizability of such metaanalytic results would be highly suspected when the I^{2} values are as high as 80% or 90%.
The comparison of statistical power in the context of greater heterogeneity merits a comment. In this research, we found that the SMD was more generalizable and less heterogeneous than the MD and that there was no significant difference in statistical power between the MD and the SMD. However, strictly and logically speaking, if one statistic is statistically less generalizable and more heterogeneous than another, it would be meaningless to discuss their relative detection power. To discuss statistical power, a high percentage agreement and a low heterogeneity are essential.
As we outlined in the Introduction, whether the MD or the SMD is more clinically preferable as a summary index of metaanalyses of continuous outcomes remains controversial. When the outcome is measured in the same natural unit, such as the amount of bleeding or the number of days of hospitalisation, the MD is definitely better than the SMD from the viewpoint of interpretability. However, our studies have suggested that the SMD may be preferable than the MD from the viewpoint of generalizability. When the outcome is measured using the same patientreported outcome (PRO) measure, such as the Hamilton Rating Scale for Depression [19], the ShortForm 36 Quality of Life Questionnaire [20], or the Severity Scoring of Atopic Dermatitis [21], even though all the outcomes are measured in the same unit, the superior interpretability of the SMD is not guaranteed unless most clinicians are very familiar with that scale. And for many clinicians in most fields of medicine, such universally known and used PRO instruments are probably rare to nonexistent. In such instances, the SMD might be more interpretable than the MD for two reasons. Firstly, the SMD can be interpreted using a general rule of thumb reported by Cohen, in which an SMD of 0.2 represents a small effect, an SMD of 0.5 represents a medium effect, and an SMD of 0.8 or larger represents a large effect [2]. Second, the SMD can be directly and easily converted to a ānumber needed to treatā if the control event can be assumed [22, 23]. In all these instances of PRO instruments, the SMD appears to be more generalizable than the MD.
There is another way to increase interpretability of continuous outcomes. When the minimal important change (MIC) is known, various methods have been proposed to facilitate the interpretability of continuous outcomes [24], including conversion to MIC units and dichotomisation using the MIC threshold. Each has its own advantages and disadvantages, and a comprehensive discussion of their relative merits is out of the realm of the present study. Unfortunately, the MIC is often either not known or, if known, may not be very precise for most of the existing PRO measures.
Finally, SMD may have another important limitation. Because its value derives from the difference in the means between the treatment and control divided by their SD, if variability of the patients is artificially or accidentally reduced, SMD would be overestimated; and if its variability is increased, SMD would be underestimated [25]. However, we think that our results concerning the greater generalizability of SMD have partially resolved such concerns because, despite the variability in the SD across trials, SMD had better external applicability and can therefore be said to have been less vulnerable to over or underestimation. However, the above mentioned possibility of too small or too large SD, and correspondingly overestimated or underestimated SMD, should always be borne in mind in interpreting SMD.
Our research has two limitations. Firstly, we did not consider the influences of multiple comparisons in our analyses. However, even if we corrected the alpha level from 0.05 to 0.0042 using the Bonferroni method because we made 12 comparisons, only the difference in the absolute Zscore in the randomeffects model would lose its significance; all the other results would not change. Secondly we did not categorize the nature of the continuous outcomes. It is possible that subjective continuous outcomes may be more prone to unstable measurement and hence be more heterogeneous than objective continuous outcomes. How generalizability and power may be influenced by such differences in continuous outcomes will be an important research topic in the future.
Conclusions
When generalizability matters, SMD may be more preferred than MD as a summary measure. In order to increase interpretabillity, SMD can then be supplemented by reporting of MD or some other proposed measures [24].
References
Higgins JPT, Green S (Eds): Cochrane Handbook for Systematic Reviews of Interventions Version 5.1.0 [updated. Edited by: Higgins JPT, Green S. 2011, Available from http://www.cochranehandbook.org, March ]. The Cochrane Collaboration, 2011
Cohen J: Statistical Power Analysis for the Behavioral Sciences. 1988, Hillsdale, New Jersey: Lawrence Erlbaum Associates: Routledge
Glass GV: Primary, secondary, and metaanalysis of research. Educ Res. 1976, 5 (10): 38.
Hedges LV, Olkin I, Statistiker M, Olkin I, Olkin I: Statistical Methods for MetaAnalysis. 1985, New York: Academic Press
Wilkinson L: Statistical methods in psychology journals: guidelines and explanations. Am Psychol. 1999, 54 (8): 594
Egger M, Smith GD, Altman D: Systematic Reviews in Health Care: MetaAnalysis in Context. 2008, Wiley. com
Association AP: Publication Manual of the American Psychological Association. 2009, Washington, DC: American Psychological Association, 6
Borenstein M, Hedges LV, Higgins JP, Rothstein HR: Introduction to MetaAnalysis. 2011, Wiley. com
Tian L: Inferences on standardized mean difference: the generalized variable approach. Stat Med. 2007, 26 (5): 945953. 10.1002/sim.2589.
Furukawa TA, Guyatt GH, Griffith LE: Can we individualize the ānumber needed to treatā? An empirical study of summary effect measures in metaanalyses. Int J Epidemiol. 2002, 31 (1): 7276. 10.1093/ije/31.1.72.
JĆørgensen AW, Hilden J, GĆøtzsche PC: Cochrane reviews compared with industry supported metaanalyses and other metaanalyses of the same drugs: systematic review. BMJ. 2006, 333 (7572): 78210.1136/bmj.38973.444699.0B.
Jadad AR, Cook DJ, Jones A, Klassen TP, Tugwell P, Moher M, Moher D: Methodology and reports of systematic reviews and metaanalyses. JAMA. 1998, 280 (3): 278280. 10.1001/jama.280.3.278.
Higgins JP, Thompson SG, Deeks JJ, Altman DG: Measuring inconsistency in metaanalyses. BMJ. 2003, 327 (7414): 55710.1136/bmj.327.7414.557.
DerSimonian R, Laird N: Metaanalysis in clinical trials. Control Clin Trials. 1986, 7 (3): 177188. 10.1016/01972456(86)900462.
DeMets DL: Methods for combining randomized clinical trials: strengths and limitations. Stat Med. 1987, 6 (3): 341348. 10.1002/sim.4780060325.
Cooper H, Hedges LV, Valentine JC: Handbook of Research Synthesis and MetaAnalysis. 2009, New York: Russell Sage Foundation
Friedrich JO, Adhikari NK, Beyene J: Ratio of means for analyzing continuous outcomes in metaanalysis performed as well as mean difference methods. J Clin Epidemiol. 2011, 64 (5): 556564. 10.1016/j.jclinepi.2010.09.016.
Friedrich JO, Adhikari NK, Beyene J: Ratio of geometric means to analyze continuous outcomes in metaanalysis: comparison to mean differences and ratio of arithmetic means using empiric data and simulation. Stat Med. 2012, 31 (17): 18571886. 10.1002/sim.4501.
Hamilton M: A rating scale for depression. J Neurol Neurosurg Psychiatry. 1960, 23: 5662. 10.1136/jnnp.23.1.56.
Ware JE, Snow KK, Kosinski M: SF36 Health Survey: Manual and Interpretation Guide. 1993, London, RI: Quality Metric Incorporated
Kunz B, Oranje A, Labreze L, Stalder JF, Ring J, Taieb A: Clinical validation and guidelines for the SCORAD index: consensus report of the European Task Force on Atopic Dermatitis. Dermatology. 1997, 195 (1): 1019.
Furukawa TA: From effect size into number needed to treat. Lancet. 1999, 353 (9165): 1680
da Costa BR, Rutjes AW, Johnston BC, Reichenbach S, NĆ¼esch E, Tonia T, Gemperli A, Guyatt GH, JĆ¼ni P: Methods to convert continuous outcomes into odds ratios of treatment response and numbers needed to treat: metaepidemiological study. Int J Epidemiol. 2012, 41 (5): 14451459. 10.1093/ije/dys124.
Guyatt GH, Thorlund K, Oxman AD, Walter SD, Patrick D, Furukawa TA, Johnston BC, Karanicolas P, Akl EA, Vist G, et al: GRADE guidelines: 13 Preparing summary of findings tables and evidence profilescontinuous outcomes. J Clin Epidemiol. 2013, 66 (2): 173183. 10.1016/j.jclinepi.2012.08.001.
Thorlund K, Walter SD, Johnston BC, Furukawa TA, Guyatt GH: Pooling healthārelated quality of life outcomes in metaāanalysisāa tutorial and review of methods for enhancing interpretability. Res Synth Methods. 2011, 2 (3): 188203. 10.1002/jrsm.46.
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712288/14/30/prepub
Acknowledgements
We received no funding for this study.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The author(s) declare that they have no competing interests.
Authorsā contributions
NT and TAF conceived and designed the study. NT, YO, AT and YH extracted data. TS checked data extraction and analysed data. All authors commented on drafts of the manuscript and approved the final manuscript.
Authorsā original submitted files for images
Below are the links to the authorsā original submitted files for images.
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 credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Takeshima, N., Sozu, T., Tajika, A. et al. Which is more generalizable, powerful and interpretable in metaanalyses, mean difference or standardized mean difference?. BMC Med Res Methodol 14, 30 (2014). https://doi.org/10.1186/147122881430
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/147122881430