| Sign In to gain access to subscriptions and/or personal tools. |
Odds Ratio, Delta, ETS Classification, and Standardization Measures of DIF Magnitude for Binary Logistic RegressionIndiana University
Merck & Co., Inc
Regenstrief Institute, Inc. and Indiana University
Previous methodological and applied studies that used binary logistic regression (LR) for detection of differential item functioning (DIF) in dichotomously scored items either did not report an effect size or did not employ several useful measures of DIF magnitude derived from the LR model. Equations are provided for these effect size indices. Using two large data sets, the authors demonstrate the usefulness of these effect sizes for judging practical importance: the LR adjusted odds ratio and its conversions to the delta metric, the Educational Testing Service (ETS) classification system, and the p metric; the LR model-based standardization indices, using various weights for averaging stratum-specific differences in fitted probabilities; and a p metric classification system. Pros and cons of these effect sizes are discussed. Recommendations are offered. These LR effect sizes will be valuable to practitioners, particularly for preventing flagging of statistically significant but practically unimportant DIF in large samples.
Key Words: Keywords: differential item functioning logistic regression effect sizes In differential item functioning (DIF) analyses, groups are compared on item performance after adjusting for overall performance on the measured trait (Holland & Wainer, 1993). Since Swaminathan and Rogers (1990) applied the binary logistic regression (LR) procedure to the detection of DIF in dichotomous test items, the LR method has become increasingly popular for this purpose. However, Swaminathan and colleagues focused on hypothesis testing (Narayanan & Swaminathan, 1996; Rogers & Swaminathan, 1993; Swaminathan & Rogers, 1990). It is important to incorporate an effect size into flagging rules, especially in large samples, because high power can yield significance for practically unimportant effect sizes (e.g., Kirk, 1996). Several methodological and applied studies investigating binary LR for DIF have flagged items for DIF based only on statistical significance (Clauser, Nungester, Mazor, & Ripkey, 1996; Huang & Dunbar, 1998; Kwak, Davenport, & Davison, 1998; Marshall, Mungas, Weldon, Reed, & Haan, 1997; Mazor, Kanjee, & Clauser, 1995; Whitmore & Schumacker, 1999; Woodard, Auchus, Godsall, & Green, 1998). Previous attempts to report effect sizes for binary LR have included using the LR Wald chi-square value (Huang & Dunbar, 1998), reporting raw or standardized LR coefficients on the log odds scale (Borsboom, Mellenbergh, & Heerden, 2002; Clauser & Mazor, 1998; Millsap & Everson, 1993; Swanson, Clauser, Case, Nungester, & Featherman, 2002), presenting R2-like measures (Swanson et al., 2002; Zumbo, 1999), calculating the partial gamma (Groenvold, Bjorner, Klee, & Kreiner, 1995), listing eta-squared (Whitmore & Schumacker, 1999), adopting a chance-corrected proportion of correct classification (Hess, Olejnik, & Huberty, 2001), and plotting fitted probabilities or fitted logits (Schmitt, Holland, & Dorans, 1993). These attempts contributed to DIF literature. However, none of these works focused on several intuitive effect sizes that can easily be derived from binary LR: the adjusted odds ratio, delta statistic, Educational Testing Service (ETS) classification system, adjusted odds ratio reported on the p metric, and model-based standardization indices of conditional differences in proportions. We found only one DIF study that reported odds ratios for binary LR (Volk, Cantor, Steinbauer, & Cass, 1997). The purposes of this article are to (a) provide and explain the equations for obtaining these useful effect sizes for the LR procedure, (b) demonstrate the application of these effect sizes, and (c) present the pros and cons of these effect sizes and offer guidance in how to use them. We focus here on effect sizes for uniform DIF. We are investigating LR effect sizes for nonuniform DIF. Although a strength of LR is powerful detection of nonuniform DIF, corresponding effect sizes requires more research because the choice of weights for averaging stratum-specific measures is especially critical when interactions are present (e.g., Mosteller & Tukey, 1977).
The Logistic Regression (LR) Procedure for DIF Detection In the binary LR model, the probability of endorsing a dichotomously scored item is
and the log odds (or logit) of endorsing the item is modeled as
where ln is the natural logarithm, x is a measure of overall proficiency (usually total score), g is a dummy variable representing group membership (traditionally, 1 = reference group, 0 = focal group), xg is the interaction term between total score and group membership, and
Effect Sizes for LR Based on the Conditional-Log-Odds-Ratio Definition of DIF
Odds ratios range from 0 to
Another option for an effect size is to transform
It is apparent that LR-D-DIF is a simple linear rescaling of the regression coefficient, In addition, one can calculate the ETS classification system (Dorans & Holland, 1993):
Notice that assigning Categories A and B entails using LR to test Ho :
It is also possible for reporting purposes to convert a conditional log-odds-ratio-based index to the metric of differences in item-proportion-correct called the p metric. We use the formula that Dorans and Holland (1993) used to convert
where,
The term Pr
Effect Sizes for LR Based on the Conditional-Difference-in-Proportions Definition of DIF
where PfmLR and PrmLR are predicted from the LR model.
This index is reminiscent of item response theory (IRT) model-based standardization (Wainer, 1993) except instead of integrating over
The weights (wm) for averaging conditional differences in proportions in the STD procedure have traditionally been based on intuitive rationale. In DIF studies, wm is often chosen to be the number of focal group examinees at each stratum (Nfm) (Dorans & Kulick, 1986). Other plausible weights include (Dorans & Holland, 1993; Mosteller & Tukey, 1977) (a) the number of reference group examinees at each stratum (Nrm), (b) the number of the total examinees at each stratum (Ntm), or (c) the relative frequency of some real or hypothetical standard group. One could also use Cochrans (1954) statistically driven weights (Dorans & Holland, 1993):
Another option, available in the STDIF software (Robin, 2001; Zenisky, Hambleton, & Robin, 2003), is the equal weight (wm = 1), which yields an unweighted average.
We performed a gender (1 = male, 0 = female) DIF analysis in two data sets. The first data set was the Supplement on Aging (SOA) to the 1984 National Health Interview Survey (U.S. Department of Health and Human Services, 1997). This study was designed to assess the future needs of the elderly in the United States. Participants were 55 and older (n = 12, 943). We analyzed 23 dichotomous functional status items. Each item measured whether participants reported a problem (1 = yes, 0 = no) performing an activity. The second data set was from the Established Populations for Epidemiologic Studies of the Elderly (EPESE). Persons (age 65) were interviewed to identify predictors of mortality (Taylor, Wallace, Ostfeld, & Blazer, 1998). We analyzed the 20 dichotomous items of the Center for Epidemiologic Studies Depression Scale (CES-D; Radloff, 1977) obtained at baseline on 3,401 participants from the Duke site. The CES-D is a widely used self-report measure of depressive symptomatology for the general population. Each item was scored for presence (1) or absence (0) of a depressive symptom.
Statistical Methods
Result of Using LR Statistical Test Alone to Detect Uniform DIF
Effect Sizes for LR Based on the Conditional-Log-Odds-Ratio Definition of DIF We will interpret the effect sizes in Table 1 from left to right, beginning with the LR odds ratio ( LR). By sorting on ascending LR-D-DIF (equivalently, descending LR), items were conveniently grouped by direction and magnitude of DIF. Thus, for the 11 statistically significant items, the 6 items at the top were more greatly endorsed by men and the 5 items at the bottom were more greatly endorsed by women after adjusting for total score. The LR for these 11 items varied in strength from 1.55 to 2.94 (men displayed greater functional problems after adjustment) and from 0.74 to 0.21 (women displayed greater functional problems after adjustment) (Table 1). For example, the LR model estimated that the odds that men reported having a problem with lifting and carrying 25 pounds was about one fifth (0.21) times the odds that women reported this problem, adjusted for overall functional status. The LR-estimated odds of having a problem with using the telephone was almost 3 (2.94) times greater for men compared to women, controlling for overall functional status (Table 1). However, the odds of having a problem with walking was only 1 times greater for men. Therefore, LR indicates that not all 11 statistically significant items exhibited equally important DIF magnitude.
Likewise, LR-D-DIF (
Effect Sizes for LR Based on the Conditional-Difference-in-Proportions Definition of DIF
Abbreviated Results for EPESE Data
Sensitivity of Results Results were very similar after deleting examinees at the floor and ceiling. We computed Cochrans (1954) test criterion by specifying wm = cm in LR-STD-P-DIF and by using LR-predicted proportions in the standard error; the observed significance was extremely similar to the LR Wald observed significance for all items in both data sets (differed at most by .008). This is not surprising given that Cochran (1954) derived these weights for a test criterion that would be powerful for detecting an alternative hypothesis of a constant difference on either the logit or probit scale. Thus, in LR, although Cochran weights are an option when computing STD-LR-P-DIF, the Cochran test might be an unnecessary adjunct to the LR Wald test.
Choosing an Effect Size: Pros and Cons The effect sizes can be contrasted on a number of dimensions. First, as for ease of interpretation, indices reported on the delta and p metric are symmetrical around their null value of zero, which facilitates interpreting DIF in opposite directions. However, those experienced with interpreting odds ratios may find LR easier to interpret than LR-D-DIF, which is on the ETS-preferred delta metric. For data conforming to the two-parameter logistic (2PL) IRT model, one advantage of LR-D-DIF is that the MH-D-DIF parameter ( 2PL) can be written as a linear rescaling of the difference between b parameters (Roussos, Schnipke, & Pashley, 1999):4 2PL = 4a(bR–bF ). The MH parameter also shares the advantage of being related, although nonlinearly, to IRT b-DIF (Roussos et al., 1999). The p metric is probably the most universally understood metric, conveniently connected to total and true score metrics. Practitioners should choose an effect size that they and their readership can easily interpret.
Second, practitioners should choose an effect size whose metric for defining departures from the null hypothesis supplies the most valid definition of DIF for their purpose. Specifically, relative to conditional odds ratios (
Third, in terms of fundamental connections to the LR model,
Fourth, as far as ease of programming, standard software for LR automatically provides Fifth, the purpose of weights in standardization is not only to standardize according to the distribution of interest but also to yield smaller weight to sparse strata that provide less precise information. Using equal weights could be dangerous if one or more strata are sparse. (None of the values for total score were sparse for the present data sets.) If one chooses an outside (real or hypothetical) standard distribution, one must be careful to not combine large weights with ill-determined differences in proportions (Mosteller & Tukey, 1977).
Recommendations on How to Use the Effect Sizes
Second, practitioners must decide what values of the effect size represent negligible, moderate, and large magnitudes for the intended purpose. For example, ETS uses thresholds of 1.0 and 1.5 on the absolute value of the delta metric, which are equivalent to odds ratios greater than 1.53 (or less than 0.65) and greater than 1.89 (or less than 0.53), respectively. Users of STD procedures often use .05 and .10 thresholds on the absolute value of the p metric. In the medical sciences, the Third, one can take steps to facilitate interpretations. One could calculate the reciprocal of odds ratios less than one. By calculating 1/.74 = 1.35, one can readily see that .74 for Item 21 in Table 1 is not as strong as 1.55 for Items 5 and 16. One can use graphs (e.g., scatter, line, and bar), which help discern relative distances between DIF magnitudes. In addition, sorting items by direction and magnitude of DIF in tables, as we did here, aids interpretation. Fourth, these effect sizes can be used to facilitate comparisons of DIF procedures. One could compare MH, LR, SIBTEST, STD, and IRT procedures on the p metric (using LR-P-DIF or LR-STD-P-DIF for LR). Likewise, one could compare procedures on the odds ratio or delta metric, where STD-P-DIF and the latent-true-score adjusted difference in proportions of SIBTEST are converted using a formula similar to Equation 22 in Dorans and Holland (1993).
Limitations
Conclusions
1 The original proposal was to use the two df simultaneous test of uniform and nonuniform differential item functioning (DIF); however, when only uniform DIF is present, including the interaction term in the test may decrease power (Swaminathan & Rogers, 1990).
2 We considered the ML subscript (maximum likelihood estimation); however, the LR subscript in Equation 3 reminds practitioners that the odds ratio was estimated by assuming a logistic regression (LR) model.
3 Jodoin and Gierl (2001) suggested that R2-like indices are preferable to effect sizes based on
4 In this formula (i.e., Equation 16 in Roussos, Schnipke, & Pashley, 1999), item discrimination (a) for the two-parameter logistic (2PL) item response theory (IRT) model varies over items and the MH delta-DIF (MH-D-DIF) parameter is conditional on theta, whereas in Equation 13 in Donoghue, Holland, and Thayer (1993), a is constant across items because the MH-D-DIF parameter is conditional on observed total score where the corresponding IRT model is the Rasch model.
PATRICK O. MONAHAN is assistant professor, Division of Biostatistics, Department of Medicine, School of Medicine, Indiana University, 410 West 10th Street Suite 3000, Indianapolis, IN 46202;pmonahan{at}iupui.edu. His area of interest is measurement and statistics applied to the behavioral and social sciences.
COLLEEN A. McHORNEY, PhD, is director of outcomes research at Merck & Co., Inc., WP39-166, 770 Sumneytown Pike, West Point, PA 19486-0004. Her areas of expertise relate to the measurement and evaluation of patient-reported outcomes, including health status, quality of life, patient satisfaction, and patient preferences.
TIMOTHY E. STUMP is statistician, Regenstrief Institute, Inc. and the Indiana University Center for Aging Research;tstump{at}regenstrief.org. His area of interest is measurement and statistics in the medical sciences.
ANTHONY J. PERKINS is a statistical consultant for the Regenstrief Institute, Inc. and the Indiana University Center for Aging Research;tperkins348{at}sbcglobal.net. His area of interest is item bias in quality of life instruments.
This research was supported by NIA Grant R01 AG022067, NCI Grant R03 CA 113099-01, and the Mary Margaret Walther Program for Cancer Care Research. Suggestions by the editor and two anonymous reviewers led to improved presentation. Manuscript received July 15, 2004. Accepted for publication August 2, 2005.
Borsboom, D, Mellenbergh, GJ, & Heerden, Jv. (2002). Different kinds of DIF: A distinction between absolute and relative forms of measurement invariance and bias. Applied Psychological Measurement, 26, 433-450. 2 tests. Biometrics, 10, 417-451.[Medline]
[Order article via Infotrieve]Donoghue, JR, Holland, PW, & Thayer, DT. In Holland, PW, & Wainer, H (Eds.). (1993). A Monte Carlo study of factors that affect the Mantel-Haenszel and standardization measures of differential item functioning. Differential item functioning (p. 137-166). Hillsdale, NJ: Lawrence Erlbaum.Dorans, NJ, & Holland, PW. In Holland, PW, & Wainer, H (Eds.). (1993). DIF detection and description: Mantel-Haenszel and standardization. Differential item functioning (p. 35-66). Hillsdale, NJ: Lawrence Erlbaum.Dorans, NJ, & Kulick, E. (1986). Demonstrating the utility of the standardization approach to assessing unexpected differential item performance on the Scholastic Aptitude Test. Journal of Educational Measurement, 23, 355-368.[CrossRef][Web of Science]Groenvold, M, Bjorner, JB, Klee, MC, & Kreiner, S. (1995). Test for item bias in a quality of life questionnaire. Journal of Clinical Epidemiology, 48, 805-816.[CrossRef][Web of Science][Medline]
[Order article via Infotrieve]Hess, B, Olejnik, S, & Huberty, CJ. (2001, April). The efficacy of two improvement-over-chance effect size measures for two-group univariate comparisons under variance heterogeneity and nonnormality, Paper presented at the annual meeting of the American Educational Research Association, Seattle, WAHolland, PW, & Thayer, DT. In Wainer, H, & Braun, HI (Eds.). (1988). Differential item performance and the Mantel-Haenszel procedure. Test validity (p. 129-145). Hillsdale, NJ: Lawrence Erlbaum.Holland, PW, & Wainer, H (Eds.). (1993). Differential item functioning. Hillsdale, NJ: Lawrence Erlbaum.Hosmer, DW, & Lemeshow, S. (2000). Applied logistic regression (2.). New York: John Wiley.Huang, C-Y, & Dunbar, SB. (1998, April). Factors influencing the reliability of DIF detection methods, Paper presented at the annual meeting of the American Educational Research Association, San Diego, CA.Jodoin, MG, & Gierl, MJ. (2001). Evaluating Type I error and power rates using an effect size measure with the logistic regression procedure for DIF detection. Applied Measurement in Education, 14, 329-349.[CrossRef][Web of Science]Kirk, RE. (1996). Practical significance: A concept whose time has come. Educational and Psychological Measurement, 56, 746-759.[Abstract]Kwak, N, Davenport, EC, Jr, & Davison, ML. (1998, April). A comparative study of observed score approaches and purification procedures for detecting differential item functioning, Paper presented at the annual meeting of the National Council on Measurement in Education, San Diego, CA.Marshall, SC, Mungas, D, Weldon, M, Reed, B, & Haan, M. (1997). Differential item functioning in the Mini-Mental State Examination in English- and Spanish-speaking older adults. Psychology and Aging, 12, 718-725.[CrossRef][Web of Science][Medline]
[Order article via Infotrieve]Mazor, KM, Kanjee, A, & Clauser, BE. (1995). Using logistic regression and the Mantel-Haenszel with multiple ability estimates to detect differential item functioning. Journal of Educational Measurement, 32, 131-144.[CrossRef][Web of Science]Millsap, RE, & Everson, HT. (1993). Methodological review: Statistical approaches for assessing measurement bias. Applied Psychological Measurement, 17, 297-334.
Journal of Educational and Behavioral Statistics, Vol. 32, No. 1,
92-109 (2007) This article has been cited by other articles:
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



0 is the intercept (
j) are estimated on the log odds scale. The exponential of 
. Values of
LR further from 1.0 represent greater DIF magnitude. An odds ratio and its reciprocal are equivalent in strength but not symmetrical in distance from the null value of 1.0 (e.g., 4.0 and 0.25).
MH): 
1.0. Practitioners can perform the latter test in LR by testing Ho: |
LR of 1.0 equals 

is the predicted proportion of examinees endorsing the item in the reference group based on 
, averaging occurs over total scores. Historically, absolute values between .05 and .10 are inspected to ensure that no possible DIF is overlooked, and absolute values above .10 are considered more unusual and should be examined (
65) were interviewed to identify predictors of mortality (
times greater for men. Therefore,
MH parameter also shares the advantage of being related, although nonlinearly, to IRT b-DIF (
2 tests. Biometrics, 10, 417-451.






