A Robust Effect Size Index

Effect size indices are useful tools in study design and reporting because they are unitless measures of association strength that do not depend on sample size. Existing effect size indices are developed for particular parametric models or population parameters. Here, we propose a robust effect size index based on M-estimators. This approach yields an index that is very generalizable because it is unitless across a wide range of models. We demonstrate that the new index is a function of Cohen’s d, R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^2$$\end{document}, and standardized log odds ratio when each of the parametric models is correctly specified. We show that existing effect size estimators are biased when the parametric models are incorrect (e.g., under unknown heteroskedasticity). We provide simple formulas to compute power and sample size and use simulations to assess the bias and standard error of the effect size estimator in finite samples. Because the new index is invariant across models, it has the potential to make communication and comprehension of effect size uniform across the behavioral sciences.


Introduction
Effect sizes are unitless indices quantifying the association strength between dependent and independent variables. These indices are critical in study design when estimates of power are desired, but the exact scale of new measurement is unknown (Cohen, 1988), and in meta-analysis, where results are compiled across studies with measurements taken on different scales or outcomes modeled differently (Chinn, 2000;Morris and DeShon, 2002). With increasing skepticism of significance testing approaches (Trafimow and Earp, 2017;Wasserstein and Lazar, 2016;Harshman et al., 2016;Wasserstein et al., 2019), effect size indices are valuable in study reporting (Fritz et al., 2012) because they are minimally affected by sample size.
Effect sizes are also important in large open source datasets because inference procedures are not designed to estimate error rates of a single dataset that is used to address many different questions across tens to hundreds of studies. While effect sizes can have similar bias to p-values when choosing among multiple hypotheses, obtaining effect size estimates for parameters specified a priori may be more useful to guide future studies than hypothesis testing because, in large datasets, p-values can be small for clinically meaningless effect sizes.
There is extensive literature in the behavioral and psychological sciences describing effect size indices and conversion formulas between different indices (see e.g., Cohen, 1988;Borenstein et al., 2009;Hedges and Olkin, 1985;Ferguson, 2009;Rosenthal, 1994;Long and Freese, 2006). Cohen 1988 defined at least eight effect size indices for different models, different types of dependent and independent variables, and provided formulas to convert between the indices. For example, Cohen's d is defined for mean differences, R 2 is used for simple linear regression, and standardized log odds ratio is used in logistic regression. Conversion formulas for some of these parametric indices are given in Table 1 and are widely used in research and software (Cohen, 1988;Borenstein et al., 2009;Lenhard and Lenhard, 2017).
Each row denotes the input argument and the column denotes the desired output value. Robust versions of classical values can be obtained by computing them as a function of S β . π 1 and π 0 denote the population proportions of each group for a two sample comparison. d is Cohen's d, f 2 β is Cohen's effect size for multiple regression, R 2 β is the partial coefficient of determination, S β is the robust index. The variables without subscripts denote the value for the full model including covariates. Conversion formulas derived by the robust index match classical formulas (Cohen, 1988;Borenstein et al., 2009;Lenhard and Lenhard, 2017) Several authors have proposed robust effect size indices based on sample quantiles (Zhang and Schoeps, 1997;Hedges and Olkin, 1984). These are robust in the sense that they do not assume a particular probability model; however, they are defined as parameters in the sense that they are a specific functional of the underlying distribution.
Despite the array of effect sizes, there are several limitations to the available indices: 1) there is no single unifying theory that links effect size indices; 2) as defined, many effect size indices do not accommodate nuisance covariates or multivariate outcomes; and 3) each index is specific to a particular population parameter. For example, Cohen's d is designed for mean differences in the absence of covariates, correlation is specific to linear regression, and existing semiparametric indices are quantile estimators. For these reasons, these classical effect size indices are not widely generalizable because their scale is dependent on the type of parameter.
In this paper, we define a new robust effect size index based on M-estimators. M-estimators are parameter estimators that can be defined as the maximizer of an estimating equation. This approach has several advantages over commonly used indices: a) The generality of M-estimators makes the index widely applicable across many types of models that satisfy mild regularity conditions, including mean and quantile estimators, so this framework serves as a canonical unifying theory to link common indices; b) the sandwich covariance estimate of M-estimators is consistent under model misspecification (MacKinnon and White, 1985;White, 1980), so the index can accommodate unknown complex relationships between second moments of multiple dependent variables and the independent variable; c) the robust effect size index is directly related to the Wald-style sandwich chi-squared statistic and is formulaically related to common indices.
Here, we describe sufficient conditions for the new effect size index to exist, describe how it relates to other indices, and show that other estimators can be biased under model misspecification.
In three examples, we show that the new index can be written as a function of Cohen's d, R 2 , and standardized log odds ratio, demonstrating that it is related to indices that were developed using intuition for specific models. In addition, we describe how to obtain a simple estimate of the index and provide functions to compute power or sample size given an effect size index and degrees of freedom of the target parameter. Finally, we use simulations to assess the bias and standard error of the proposed index estimator. An R package to estimate the index is in development; the latest release is available at https://github.com/simonvandekar/RESI.

Notation
Unless otherwise noted, capital letters denote vectors or scalars and boldface letters denote matrices; lower and upper case greek letters denote vector and matrix parameters, respectively. Let W 1 = {Y 1 , X 1 }, . . . , W n = {Y n , X n } be a sample of independent observations from W ⊂ R p with associated probability measure G and let H denote the conditional distribution of Y i given X i . Here, W i denotes a combination of a potentially multivariate outcome vector Y i with a multivariate covariate vector X i .
Let W = {W 1 , . . . , W n } denote the full dataset and θ * → (θ * ; W ) ∈ R, θ * ∈ R m be an estimating equation, where ψ is a known function. is a scalar-valued function that can be maximized to obtain the M-estimatorθ . We define the parameter θ as the maximizer of the expected value of the estimating equation under the true distribution G, and the estimatorθ isθ = arg max where α ∈ R m 0 denotes a nuisance parameter, β ∈ R m 1 is the target parameter, and m 0 +m 1 = m. We define the m × m matrices with j, kth elements which are components of the asymptotic robust covariance matrix of √ n(θ − θ).

Definition
Here, we define a robust effect size that is based on the test statistic for β 0 is a reference value used to define the index. Larger distances from β 0 represent larger effect sizes. Under the regularity conditions in the Appendix, This implies that the typical robust Wald-style statistic for the test of (4) is approximately chisquared on m 1 degrees of freedom, is the asymptotic covariance matrix ofβ, and can be derived from the covariance of (5) (Boos and Stefanski, 2013;Van der Vaart, 2000). We define the square of the effect size index as the component of the chi-squared statistic that is due to the deviation of β from the reference value: As we demonstrate in the examples below, the covariance β (θ ) serves to standardize the parameter β so that it is unitless. The regularity conditions given in Appendix are sufficient for the index to exist. The robust index, S β (θ ) := S β (θ ) 2 , is defined as the square root of S β (θ ) 2 so that the scale is proportional to that used for Cohen's d (see Example 1). This index has several advantages: it is widely applicable because it is constructed from Mestimators; it relies on a robust covariance estimate; it is directly related to the robust chi-squared statistic; it is related to classical indices, and induces several classical transformation formulas (Cohen, 1988;Borenstein et al., 2009;Lenhard and Lenhard, 2017).

An Estimator
S β (θ ) is defined in terms of parameter values and so must be estimated from data when reported in a study. Let T m 1 (θ) 2 be as defined in (6), then is consistent for S β (θ ), which follows by the consistency of the components that make up T m 1 (θ) 2 (Van der Vaart, 2000;White, 1980). We use the factor (n − m) to account for the estimation of m parameters.

Examples
In this section, we show that the robust index yields several classical effect size indices when the models are correctly specified. We demonstrate the interpretability of the effect size index through a series of examples. The following example shows that the robust index for a difference in means is proportional to Cohen's d, provided that the parametric model is correctly specified; that is, the solution to the M-estimator is equal to the MLE.
Example 1. (Difference in means) In this example, we consider a two mean model, where W i = {Y i , X i } and the conditional mean of Y i ∈ R given X i converges. That is, , and we assume the limit (9) exists. In addition, we assume P(X i = 1) = π 1 = 1 − π 0 is known and that . Note that in this example is not defined as the derivative of a log-likelihood: It defines a single parameter that is a difference in means and does not require each observation to have the same distribution. Despite this general approach, we are still able to determine the asymptotic variance of n 1/2θ , Then the robust effect size (7) is For fixed sample proportions π 0 and π 1 , when σ 2 0 = σ 2 1 , S β (θ ) is proportional to the classical index of effect size for the comparison of two means, Cohen's d (Cohen, 1988). However, S β (θ ) is more flexible: It can accommodate unequal variance among groups and accounts for the effect that unequal sample proportions has on the power of the test statistic. Thus, S is an index that accounts for all features of the study design that will affect the power to detect a difference. The robust index is proportional to the absolute value of the large sample z-statistic that does not rely on the equal variance assumption. This is what we expect in large samples when the equal variance assumption is not necessary for "correct" inference. In this example, we did not explicitly assume an identical distribution for all X i , only that the mean of the variance of Y i given X i converges in probability to a constant.
The following example derives the robust effect size for simple linear regression. This is the continuous independent variable version of Cohen's d and is related to R 2 .

Example 2. (Simple linear regression) Consider the simple linear regression model
where α and β are unknown parameters, Y i ∈ R, X i ∈ R and i follows an unknown distribution with zero mean and conditional variance that can depend on where After some algebra, combining the formulas (11) and (12) gives Then (7) is The intuition of (13) is best understood by considering the homoskedastic case where E H (Y i − α − X i β) 2 = σ 2 for all i = 1, . . . , n. Then, σ 4 x /σ 2 xy β 2 = σ 2 x /σ 2 β 2 . This is similar to R 2 , except that the denominator is the variance of Y i conditional on X i instead of the marginal variance of Y i . The denominator of (13) accounts for the possible dependence between X i and Var(Y i | X i ).
In the following example, we introduce two levels of complexity by considering logistic regression with multidimensional nuisance and target parameters.
Example 3. (Logistic regression with covariates) For logistic regression, we utilize the model where Y i is a Bernoulli random variable, X i = [X i0 , X i1 ] ∈ R p−1 is a row vector, and α and β are as defined in (3). Let X = [X T 1 . . . X T n ] T ∈ R n×( p−1) and similarly define X 0 and X 1 . Let P ∈ R n×n be the matrix with If this equality does not hold, then there is under or over dispersion.
To find the robust effect size, we first need to find the covariance matrix ofβ. To simplify notation, we define the matrices for k, = 0, 1. The block matrix of J G (θ ) −1 corresponding to the parameter β is Equation (15) is the asymptotic covariance ofβ, controlling for X 0 , if model (14) is correctly specified. The robust covariance for β can be derived by finding the block matrix of J (θ ) −1 K (θ )J (θ ) −1 corresponding to β. In this general case, the asymptotic covariance matrix ofβ is If the model is correctly specified, P = Q, β (θ ) = I β (θ ) −1 , then The parameter (16) describes the effect of β controlling for the collinearity of variables of interest X 1 , with the nuisance variables, X 0 . If the collinearity is high, then the diagonal of I β (θ ) −1 will be large and the effect size will be reduced. Many suggestions have been made to compute standardized coefficients in the context of logistic regression (for a review see Menard, 2004Menard, , 2011. When m 1 = 1, the square of the robust index in this context, under correct model specification, is the square of a fully standardized coefficient and differs by a factor of √ n from the earliest proposed standardized index (Goodman, 1972). The index proposed by Goodman (1972) is simply a Wald statistic and was rightly criticized for its dependence on the sample size (Menard, 2011), despite that it correctly accounts for the fact that the variance of a binomial random variable is functions of its mean, through the use of the diagonal matrix P in the matrix I β (θ ). The robust index remediates the dependence that Goodman's standardized coefficient has on the sample size. Graphs of the robust effect size as a function of some common effect size indices (see formulas in Table 1. a Cohen's d, when π 0 = π 1 = 1/2 and σ 0 = σ 1 ; b R 2 .

Relation to Other Indices
The new index can be expressed as a function of several common effect size indices for continuous or dichotomous dependent variables when there is homoskedasticity ( Fig. 1; Table 1). The relations between effect sizes implied by the new index are equivalent to the classical conversion formulas between effect sizes (Borenstein et al., 2009;Selya et al., 2012). While the index is related to existing indices under correct model specification, the advantage of the robust index is that it is defined if the variance model is incorrectly specified. This is the case, for example, in linear regression when there is heteroskedasticity and the model assumes a single variance term for all subjects or in logistic regression when there is overdispersion. By using the formulas in Table 1, we can obtain robust versions of classical indices by writing them as functions of S β . Cohen (1988) defined ranges of meaningful effect sizes for the behavioral sciences (Table  2). These intervals can also be used to define similar regions for the robust index. These recommendations serve as a useful guide, however, ranges of meaningful effect sizes are field specific and should be based on clinical expertise and the effect an intervention could have if applied to the population of interest.

Bias of Existing Indices Under Model Misspecification
To understand the bias of the classical estimators under model misspecification, we compare the asymptotic value of the classical estimators to the effect size formulas in Table 1. Under model misspecification, the existing parametric effect size indices can be biased.
The classical estimator for R 2 converges to Taking the ratio of this value and the formula for R 2 (S) given in Table 1 gives, where variables are as defined in (12). Figure 2 plots the bias as a function of log 2 {σ 2 xy /(σ 2 x σ 2 y )}. When the variance is constant across subjects, Var(Y i | X i ) = σ 2 y , then the bias is zero. If not, then the direction of the bias of the classical estimator depends on the relationship between Var(Y i | X i ) and X i .  Percent bias for Cohen's d and R 2 . When π 1 = 1/2 or the variances are equal the classical estimator of Cohen's d is unbiased, it can be positively or negatively biased when the variances and sampling proportions are not equal. Similarly for R 2 , when Var(Y i | X i ) is constant across subjects, there is no bias (because σ 2 xy = σ 2 x σ 2 y ), but when this is not true, the classical estimator can be positively or negatively biased depending on the relationship between the variances. Variables are as defined in (12).

Determining Effect Sizes, Sample Sizes, and Power
A convenient aspect of the robust index is that it makes asymptotic power calculations easy. The formula is the same for every parameter that is a solution to an estimating equation, such as (2). For a fixed sample size and rejection threshold, power can be determined from the robust index and degrees of freedom of the chi-squared test using (6). The explicit formula for power can be written where t 1 and t 2 denote the type 1 and type 2 error rates, respectively, df denotes the degrees of freedom of the test statistic, (·; λ) denotes the cumulative distribution function of a noncentral chi-squared distribution with noncentrality parameter λ, and S β is as defined in (7). Equation (17) can be easily solved for sample size, power, error rate, or effect size, using basic statistical software with fixed values of the other variables (Fig. 3). Because the robust index is not model dependent, power curves are effectively model-free and applicable for any fixed sample size, rejection threshold, and degrees of freedom.

Simulation Analysis
We used 1,000 simulations to assess finite sample bias and standard errors of the estimator (8). Covariates of row vectors, X i , were generated from a multivariate normal distribution X i ∼ N (0, X ), where, Bias and standard error ofŜ when the data generating distribution has skew=0.63 with two nuisance covariates (m 0 = 2). S tends to be positively biased across values of S. The standard error is proportional to S and is quite large in small samples. Rhosq denotes the total squared correlation of nuisance covariates with the target variables. Rhosq does not affect the bias, standard error, or value of the effect size index because S is defined conditionally on the covariates.
with ρ 2 ∈ {0, 0.6}, m 0 ∈ {2, 5}, and m 1 ∈ {1, 3, 5}. Here, I m 0 and 1 m 0 denote the m 0 × m 0 identity matrix and a vector of ones in R m 0 , respectively. This distribution implies that the total correlation between the nuisance covariates and target covariates is equal to ρ 2 . Samples of Y i , for i = 1, . . . , n of size n ∈ {25, 50, 100, 250, 500, 1000} were generated with mean where β was determined such that S ∈ {0, 0.1, 0.25, 0.4, 0.6}. We used a gamma distribution with shape parameter a ∈ {0.5, 10} and rate equal to a/ X 2 i,m 0 +1 to generate heteroskedastic errors for Y i . For each simulation, we compute bias of the estimator (8). Only a subset of the results are reported here; however, code to run the simulations and the saved simulation results are published with this paper.
Bias and standard error of the estimator is presented for ρ 2 ∈ {0, 0.6}, for m 0 = 2 , and all values of S considered in the simulations (Fig. 4). Results demonstrate the effect size estimator is biased upwards in small samples, but the bias is close to zero for sample sizes over 500. Because the effect size is defined conditional on covariates, the existence of covariates does not affect estimation bias. The standard error of the estimator is larger in small samples and for larger values of S. When the sample size is small, n = 25, the standard error can be larger than the value of S β .

Discussion
We proposed a robust effect size index that utilizes an M-estimator framework to define an index that is generalizable across a wide range of models. The robust index provides a unifying framework for formulaically relating effect sizes across different models. The proposed index is robust to model misspecification, has an easily computable estimator, and is related to classical effect size indices. We showed that classical estimators can be asymptotically biased when the covariance model is misspecified.
The relationship between the robust index and indices based on correctly specified models (such as Cohen's d and R 2 ) is appealing because it follows intuition from other areas of robust covariance estimation. That is, when the estimating equation is proportional to the log-likelihood, then the robust index is a function of classical definitions derived from likelihood-based models. The new framework also generalizes classical indices by easily accommodating nuisance covariates and sandwich covariance estimators that are robust to heteroskedasticity. The robust index puts effect sizes for all models on the same scale so that asymptotically accurate power analyses can be performed for model parameters using a single framework.
One important feature of the proposed index is that it is defined conditional on covariates. While the effect size lies on a standardized scale that is related directly to the power of the test, the inclusion of covariates affects the interpretation of the index because it is defined conditional on the covariates. For this reason, careful consideration of the target parameter is necessary for accurate interpretation and comparison across studies that present the robust index. Marginal estimators (without conditioning on covariates) should be considered if the investigator is interested in the general effect across a given population.
Several limitations may inspire future research topics: Like p-values, estimates of effect size indices can be subject to bias by data dredging. Also, the estimator can be biased in small samples because the index is based on asymptotic results. Thus, methods for bias adjustment or low mean squared error estimators could be considered to adjust the effects of data dredging or small sample sizes. Here, we considered an M-estimator framework, but a semiparametric or robust likelihood framework may have useful properties as well (Royall and Tsou, 2003;Blume et al., 2007). We believe this index serves as a first step in constructing a class of general robust effect size estimators that can make communication of effect sizes uniform across models in the behavioral sciences.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.