Linked-Tucker Factorized Individualized Regression for Paired Multivariate Categorical Outcomes
Abstract
We propose a joint individualized hurdle-ordinal regression model for paired zero-inflated ordinal outcomes with subject-specific, spatially varying, and time-varying covariate effects, motivated by the Iowa Fluoride Study (IFS). The two outcomes, dental caries and dental fluorosis, are measured repeatedly across ages at fine spatial resolution, yielding nested longitudinal data with substantial zero inflation, ordinality, and heterogeneity across individuals and locations. For each outcome, a hurdle component models disease presence, while a proportional-odds component models severity among positive observations. To parsimoniously represent the high-dimensional coefficient arrays, we introduce a linked Tucker tensor factorization. Shared subject-mode factors induce dependence between the caries and fluorosis coefficient tensors, while separate spatial factors accommodate the distinct measurement grids of tooth surfaces and tooth zones. A horseshoe prior on the core tensor elements encourages sparsity, and posterior computation is performed using the No-U-Turn Sampler in NumPyro. Population-level effect summaries are obtained by projecting individualized posterior linear predictors onto the design space, and Wasserstein barycenters aggregate these summaries across tooth locations and anatomical classes. Applied to the IFS, the model reveals spatially heterogeneous associations between early-life fluoride and dietary exposures and both outcomes. Fluoride exposure is associated with increased odds and severity of fluorosis, while soda intake consistently increases caries risk. These associations differ between presence and severity components and vary across tooth locations, ages, and subpopulations defined by prior caries status, highlighting the importance of the joint hurdle-ordinal framework for disentangling disease occurrence from disease progression in multilevel dental data.
Related
- Priors and parameterizations for Bayesian regression with categorical variables and interactions
- Minimum Volume Conformal Sets for Multivariate Regression
- Diagonalized Attention for Individualized Regression: Latent-Row Localization and Prediction
- Hybrid Variational Quantum Circuits for Multivariate Regression and High-Dimensional Data Reconstruction
- Cellwise and Casewise Robust Multivariate Regression with Inference
BibTeX
@misc{roy2026linkedtucker,
title = {{Linked-Tucker Factorized Individualized Regression for Paired Multivariate Categorical Outcomes}},
author = {Arkaprava Roy and Jeremy T. Gaskins and Steven Levy and Somnath Datta},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2605.06135},
eprint = {2605.06135},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2605.06135},
}
RIS
TY - UNPB
TI - Linked-Tucker Factorized Individualized Regression for Paired Multivariate Categorical Outcomes
AU - Arkaprava Roy
AU - Jeremy T. Gaskins
AU - Steven Levy
AU - Somnath Datta
PY - 2026
DA - 2026/09/01
PB - arXiv
DO - 10.48550/arxiv.2605.06135
UR - https://arxiv.org/abs/2605.06135
AB - We propose a joint individualized hurdle-ordinal regression model for paired zero-inflated ordinal outcomes with subject-specific, spatially varying, and time-varying covariate effects, motivated by the Iowa Fluoride Study (IFS). The two outcomes, dental caries and dental fluorosis, are measured repeatedly across ages at fine spatial resolution, yielding nested longitudinal data with substantial zero inflation, ordinality, and heterogeneity across individuals and locations. For each outcome, a hurdle component models disease presence, while a proportional-odds component models severity among positive observations. To parsimoniously represent the high-dimensional coefficient arrays, we introduce a linked Tucker tensor factorization. Shared subject-mode factors induce dependence between the caries and fluorosis coefficient tensors, while separate spatial factors accommodate the distinct measurement grids of tooth surfaces and tooth zones. A horseshoe prior on the core tensor elements encourages sparsity, and posterior computation is performed using the No-U-Turn Sampler in NumPyro. Population-level effect summaries are obtained by projecting individualized posterior linear predictors onto the design space, and Wasserstein barycenters aggregate these summaries across tooth locations and anatomical classes. Applied to the IFS, the model reveals spatially heterogeneous associations between early-life fluoride and dietary exposures and both outcomes. Fluoride exposure is associated with increased odds and severity of fluorosis, while soda intake consistently increases caries risk. These associations differ between presence and severity components and vary across tooth locations, ages, and subpopulations defined by prior caries status, highlighting the importance of the joint hurdle-ordinal framework for disentangling disease occurrence from disease progression in multilevel dental data.
ER -