An Active Set Algorithm to Estimate Parameters in Generalized Linear Models with Ordered Predictors

Abstract

In biomedical studies, researchers are often interested in assessing the association between one or more ordinal explanatory variables and an outcome variable, at the same time adjusting for covariates of any type. The outcome variable may be continuous, binary, or represent censored survival times. In the absence of precise knowledge of the response function, using monotonicity constraints on the ordinal variables improves efficiency in estimating parameters, especially when sample sizes are small. An active set algorithm that can efficiently compute such estimators is proposed, and a characterization of the solution is provided. Having an efficient algorithm at hand is especially relevant when applying likelihood ratio tests in restricted generalized linear models, where one needs the value of the likelihood at the restricted maximizer. The algorithm is illustrated on a real life data set from oncology.

24 pages, 1 Figure, 3 Tables

References in corpus (7)

• Maximum likelihood estimation of a log-concave density and its distribution function: Basic properties and uniform consistency
• Inconsistency of bootstrap: The Grenander estimator
• Bootstrapping the Grenander estimator
• Active Set and EM Algorithms for Log-Concave Densities Based on Complete and Censored Data
• Least Squares and Shrinkage Estimation under Bimonotonicity Constraints
• Least Squares estimation of two ordered monotone regression curves
• An Algorithm to Estimate Monotone Normal Means and its Application to Identify the Minimum Effective Dose

Cited by in corpus (2)

• Sparsity with sign-coherent groups of variables via the cooperative-Lasso
• A Generic Path Algorithm for Regularized Statistical Estimation