REVIEW 2 major objections 1 minor 23 references
Direct and Indirect Discrimination in Generalized Linear Models
T0 review · 2 major / 1 minor · reviewed 2026-06-26 · grok-4.3
Pith's one-line read A moment decomposition of GLM predictions isolates four linear disparity channels plus two curvature terms from the inverse link.
desk verdict The paper supplies explicit moment decompositions for GLM fairness that add curvature coupling and amplification terms, but the whole thing rests on an untested extension of the Wasserstein criterion past the linear-Gaussian case. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The leading decomposition D1(f) of the empirical output-scale criterion U2(f), which extends the reduced Wasserstein two-moment criterion by preserving linear channels and adding inverse-link curvature terms.
What would settle it
A numerical check in which the sum of the six decomposed components deviates materially from the observed difference in average GLM predictions between two groups on the same data.
Extended reading notes
Core claim
In the exact linear-Gaussian benchmark the Wasserstein barycentric criterion reduces to a two-moment criterion and decomposes into direct mean, indirect mean, interaction and structural components. For GLMs the leading decomposition D1(f) preserves the four linear channels and adds two curvature components induced by the inverse link: curvature coupling and curvature amplification. Explicit formulas are derived for logistic, Poisson and Tweedie specifications.
Load-bearing premise
The moment-based reduction of the Wasserstein barycentric criterion remains a sufficient diagnostic for group disparities once the nonlinear inverse link is present.
Editorial extensions
If this is right
- Disparities in GLM predictions can be attributed separately to explicit group effects, proxy covariates, covariance structure and nonlinear link effects.
- Explicit formulas allow direct computation of curvature coupling and amplification for logistic, Poisson and Tweedie models.
- The decomposition supplies a tractable actuarial diagnostic on real datasets such as medical-expenditure surveys.
Reading between the lines
- The same moment decomposition could be applied to other GLM families or to models with multiple inverse links to test whether curvature terms scale with variance.
- Comparing the leading D1(f) term against full distributional parity measures on the same data would quantify how much information is lost by the moment reduction.
- The framework suggests auditing fitted models by recomputing the decomposition after removing each sensitive or proxy variable in turn.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a moment-based decomposition framework for group disparities in fitted GLM predictions. In the linear-Gaussian benchmark the Wasserstein barycentric criterion reduces to a two-moment criterion that decomposes into direct mean, indirect mean, interaction, and structural components. For GLMs the authors distinguish the empirical output-scale criterion U_2(f), a within-group proxy ilde U_2(f), and a leading decomposition D_1(f) that retains the four linear channels while adding two curvature components (coupling and amplification) induced by the inverse link; explicit formulas are supplied for logistic, Poisson, and Tweedie families and the diagnostic is illustrated on medical-expenditure survey data. The framework is presented as a tractable actuarial tool rather than a full distributional-parity characterization or legal test.
Significance. If the explicit formulas are correct and D_1(f) remains informative, the work supplies actuaries with closed-form expressions that separate direct, indirect, structural, and nonlinear-link contributions to prediction disparities in common GLM specifications. The provision of exact formulas for three families and the explicit acknowledgment of scope limitations are concrete strengths that could support practical adoption.
major comments (2)
- [Abstract / derivation of D_1(f)] Abstract and derivation of D_1(f): the central claim that the leading decomposition preserves the four linear channels and adds explicit curvature terms rests on an unverified reduction from the Wasserstein barycentric criterion once the nonlinear inverse link is present; the manuscript supplies no numerical checks, simulation studies, or comparisons against the full distributional criterion to confirm that omitted higher-order link effects do not dominate in typical GLM regimes.
- [Illustration section] Illustration section: the medical-expenditure example is described only at a high level; without reported values of U_2(f), ilde U_2(f), and the separate curvature components it is impossible to assess whether the added curvature terms materially alter the linear-channel decomposition in a real GLM fit.
minor comments (1)
- [Methods] Notation for U_2(f) and ilde U_2(f) is introduced in the abstract but the precise definitions and their relationship to the Wasserstein criterion should be restated at the beginning of the methods section for clarity.
Simulated Author's Rebuttal
We thank the referee for the constructive comments. We address each major point below and indicate the revisions we will make.
read point-by-point responses
-
Referee: [Abstract / derivation of D_1(f)] Abstract and derivation of D_1(f): the central claim that the leading decomposition preserves the four linear channels and adds explicit curvature terms rests on an unverified reduction from the Wasserstein barycentric criterion once the nonlinear inverse link is present; the manuscript supplies no numerical checks, simulation studies, or comparisons against the full distributional criterion to confirm that omitted higher-order link effects do not dominate in typical GLM regimes.
Authors: The derivation of D_1(f) is obtained via a first-order Taylor expansion of the inverse link around the group means; this expansion exactly retains the four linear channels and isolates the two leading curvature terms by construction. The manuscript positions D_1(f) explicitly as a leading-term diagnostic rather than an exact match to the full Wasserstein criterion outside the linear-Gaussian case. We agree that a numerical check would strengthen the claim and will add a short simulation study comparing D_1(f) to the full distributional criterion under representative GLM regimes in the revision. revision: yes
-
Referee: [Illustration section] Illustration section: the medical-expenditure example is described only at a high level; without reported values of U_2(f), ilde U_2(f), and the separate curvature components it is impossible to assess whether the added curvature terms materially alter the linear-channel decomposition in a real GLM fit.
Authors: We agree that the illustration section would be more informative with explicit numerical values. The revised manuscript will report the computed values of U_2(f), ilde U_2(f), and the individual components of D_1(f) for the medical-expenditure data, together with a short discussion of the relative size of the curvature contributions. revision: yes
Circularity Check
No circularity: decomposition constructed from GLM structure and Wasserstein benchmark
full rationale
The paper starts from the Wasserstein barycentric criterion (reducing to moments only in linear-Gaussian case), then defines empirical output-scale U2(f), within-group proxy, and leading D1(f) by adding explicit curvature-coupling and curvature-amplification terms induced by the inverse link. Closed-form expressions for logistic, Poisson and Tweedie follow directly from standard GLM mean-variance and link properties rather than re-expressing any fitted parameter as a 'prediction' of itself. No self-citation chain, uniqueness theorem, or ansatz-smuggling is invoked; the abstract explicitly disclaims full distributional parity outside the benchmark, confirming the construction does not collapse to its inputs by definition.
Assumptions & free parameters
assumptions (1)
- domain assumption Wasserstein barycentric criterion measures distributional demographic-parity violation
Cite this review
Pith. "Pith review of Direct and Indirect Discrimination in Generalized Linear Models." pith.science (2026). https://pith.science/paper/V7GFCY2P
@misc{pith2026260623116,
author = {Pith},
title = {Pith review of: Direct and Indirect Discrimination in Generalized Linear Models},
year = {2026},
howpublished = {\url{https://pith.science/paper/V7GFCY2P}},
note = {Machine review of arXiv:2606.23116}
}
abstract
Generalized linear models are central to actuarial modelling of binary risk, claim frequency, utilization, and cost-related outcomes. Yet fairness diagnostics often rely on linear-model intuitions, although GLM predictions are obtained by transporting a latent score through a nonlinear inverse link. We develop a moment-based decomposition framework for diagnosing group disparities in fitted GLM predictions. In an exact linear-Gaussian benchmark, the Wasserstein barycentric criterion for distributional demographic-parity violation reduces to a two-moment criterion and decomposes into direct mean, indirect mean, interaction, and structural components. For GLMs, we distinguish the empirical output-scale criterion $U_2(f)$, a within-group proxy $\widetilde U_2(f)$, and a leading decomposition $D_1(f)$. This leading term preserves the four linear channels and adds two curvature components induced by the inverse link: curvature coupling and curvature amplification. We derive explicit formulas for logistic, Poisson, and Tweedie specifications and illustrate the diagnostic on medical-expenditure survey data. The framework is not a legal test of discrimination, nor a full characterization of distributional parity outside the linear-Gaussian case. It is a tractable actuarial diagnostic for identifying whether fitted prediction disparities arise from explicit sensitive effects, proxy-mediated covariate profiles, covariance-structure differences, or nonlinear link effects.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
Fairness and risk: An ethical argument for a group fairness definition insurers can use
Baumann, J., Loi, M., 2023. Fairness and risk: An ethical argument for a group fairness definition insurers can use. Philosophy & Technology 36, 45
2023
-
[2]
Wage discrimination: Reduced form and structural estimates
Blinder, A.S., 1973. Wage discrimination: Reduced form and structural estimates. The Journal of Human Resources 8, 436–455
1973
-
[3]
Insurance, Biases, Discrimination and Fairness
Charpentier, A., 2024. Insurance, Biases, Discrimination and Fairness. Springer Actuarial, Springer
2024
-
[4]
A fair price to pay: Exploiting causal graphs for fairness in insurance
Côté, O., Côté, M.P., Charpentier, A., 2025. A fair price to pay: Exploiting causal graphs for fairness in insurance. Journal of Risk and Insurance doi:10.1111/jori.12503
-
[5]
Generalized linear models for insurance data
De Jong, P., Heller, G.Z., 2008. Generalized linear models for insurance data. Cambridge University Press
2008
-
[6]
Effective Statistical Learning Methods for Actuaries I: GLMs and Extensions
Denuit, M., Hainaut, D., Trufin, J., 2019. Effective Statistical Learning Methods for Actuaries I: GLMs and Extensions. Springer
2019
-
[7]
Generalized Linear Models with Examples in R
Dunn, P.K., Smyth, G.K., 2018. Generalized Linear Models with Examples in R. 2 ed., Springer. 60
2018
-
[8]
An extension of the blinder-oaxaca decomposition technique to logit and probit models
Fairlie, R.W., 2005. An extension of the blinder-oaxaca decomposition technique to logit and probit models. Journal of economic and social measurement 30, 305–316
2005
Show all 23 references
-
[9]
Decomposition methods in eco- nomics, in: Ashenfelter, O., Card, D
Fortin, N., Lemieux, T., Firpo, S., 2011. Decomposition methods in eco- nomics, in: Ashenfelter, O., Card, D. (Eds.), Handbook of Labor Economics. Elsevier. volume 4A, pp. 1–102
2011
-
[10]
Regression modeling with actuarial and financial appli- cations
Frees, E.W., 2009. Regression modeling with actuarial and financial appli- cations. Cambridge University Press
2009
-
[11]
Predictive modeling applica- tions in actuarial science
Frees, E.W., Derrig, R.A., Meyers, G., 2014. Predictive modeling applica- tions in actuarial science. volume 1. Cambridge University Press
2014
-
[12]
Predictive modeling applica- tions in actuarial science
Frees, E.W., Derrig, R.A., Meyers, G., 2016. Predictive modeling applica- tions in actuarial science. volume 2. Cambridge University Press
2016
-
[13]
Demographic parity constrained minimax optimal regression under linear model, in: Advances in Neural Information Processing Systems, pp
Fukuchi, K., Sakuma, J., 2023. Demographic parity constrained minimax optimal regression under linear model, in: Advances in Neural Information Processing Systems, pp. 8653–8689
2023
-
[14]
Generalized Additive Models
Hastie, T.J., Tibshirani, R.J., 1990. Generalized Additive Models. Chapman and Hall
1990
-
[15]
The blinder–oaxaca decomposition for linear regression models
Jann, B., 2008. The blinder–oaxaca decomposition for linear regression models. The Stata Journal 8, 453–479
2008
-
[16]
Modern Actuarial Risk Theory: Using R
Kaas, R., Goovaerts, M., Dhaene, J., Denuit, M., 2008. Modern Actuarial Risk Theory: Using R. Springer
2008
-
[17]
Components of a difference between two rates
Kitagawa, E.M., 1955. Components of a difference between two rates. Journal of the American Statistical Association 50, 1168–1194
1955
-
[18]
Discrimination-free insurance pricing
Lindholm, M., Richman, R., Tsanakas, A., Wüthrich, M.V., 2022. Discrimination-free insurance pricing. ASTIN Bulletin 52, 55–89. doi:10. 1017/S0515036121000234. 61
2022
-
[19]
Generalized Linear Models
McCullagh, P., Nelder, J.A., 1989. Generalized Linear Models. 2 ed., Chapman and Hall
1989
-
[20]
Male-female wage differentials in urban labor markets
Oaxaca, R., 1973. Male-female wage differentials in urban labor markets. International Economic Review 14, 693–709
1973
-
[21]
Non-life insurance pricing with generalized linear models
Ohlsson, E., Johansson, B., 2010. Non-life insurance pricing with generalized linear models. volume 174. Springer
2010
-
[22]
Decomposing direct and indirect biases in linear models under demographic parity constraint, in: Proceedings of the AAAI Conference on Artificial Intelligence
Tierny, B., Charpentier, A., Hu, F., 2026. Decomposing direct and indirect biases in linear models under demographic parity constraint, in: Proceedings of the AAAI Conference on Artificial Intelligence
2026
-
[23]
Generalized Additive Models: An Introduction with R
Wood, S.N., 2017. Generalized Additive Models: An Introduction with R. 2 ed., Chapman and Hall/CRC. Appendix A. Empirical Design and Reproducibility This appendix describes the empirical design used in Sections 2.5, 4.4, 5.5, and 6.4. The numerical results can be reproduced by...
2017
Reviewed June 26, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.