Self-modeling ordinal regression with time invariant covariates

Abstract

In a prostate cancer study, the severity of genito-urinary (bladder) toxicity is assessed for patients who were given different doses of radiation. The ordinal responses (severity of side effects) are recorded longitudinally along with the cancer stage of a patient. Differences among the patients due to time-invariant covariates are captured by the parameters. To build up a suitable framework for an analysis of such data, we propose the use of self-modeling ordinal longitudinal model where the conditional cumulative probabilities for a category of an outcome have a relation with shape-invariant model. Since patients suffering from a common disease usually exhibit a similar pattern, it is natural to build up a nonlinear model that is shape invariant. The model is essentially semi-parametric where the population time curve is modeled with penalized regression spline. Monte Carlo expectation maximization technique is used to estimate the parameters of the model. A simulation study is also carried out to justify the methodology used.

Keywords

Longitudinal data self-modeling nonparametric penalized spline

Get full access to this article

View all access options for this article.

References

Lawton

Sylvestre

Maggio

. Self modeling non-linear regression. Technometrics 1972; 14: 513–532.

Kneip

Gasser

. Convergence and consistency results for self-modeling nonlinear regression. Ann Stat 1988; 16: 82–112.

Altman

Villarreal

. Self-modelling regression for longitudinal data with time-invariant covariates. Can J Stat 2004; 32: 251–268.

Ruppert

Wand

Carroll

. Semiparametric regression, Cambridge: Cambridge University Press, 2003.

Ruppert D and Carroll RJ. Penalized regression splines. Technical Report, Operations Research, Cornell University, 1997.

Sheiner

Beal

. Evaluation of methods for estimating population pharmacokinetic parameters. I. Michaelis–Menten model: Routine clinical pharmacokinetic data. J Pharm Biopharm 1980; 8: 553–571.

Lindstrom

Bates

. Nonlinear mixed effects models for repeated measures data. Biometrics 1990; 46: 673–687.

Wolfinger

. Laplacian’s approximation for nonlinear mixed models. Biometrika 1993; 80: 791–795.

Vonesh

. A note on Laplacians approximation in nonlinear mixed effects models. Biometrika 1996; 83: 447–452.

10.

Davidian

Giltinan

. Nonlinear models for repeated measures data, New York: Chapman & Hall, 1995.

11.

Pinheiro

Bates

. Approximations to the log-likelihood function in the nonlinear mixed-effects model. J Comput Graph Statist 1995; 4: 12–35.

12.

Lindstrom

Bates

. Newton-Raphson and EM algorithms for linear mixed-effects models for repeated-measures data. J Am Stat Assoc 1988; 83: 1014–1022.

13.

Wei

GCG

Tanner

. Calculating the content and the boundary of the highest posterior density region via data augmentation. Biometrika 1990; 77: 649–652.

14.

McCullagh

. Regression models for ordinal data (with discussion). J Roy Stat Soc B 1980; 42: 109–142.

15.

Kneip

Engel

. Model estimation in nonlinear regression under shape invariance. Ann Stat 1995; 23: 551–570.

16.

Wahba

. A comparison of GCV and GML for choosing the smoothing parameter in the generalized spline smoothing problem. Ann Stat 1985; 13: 1378–1402.

17.

Robinson GK. That BLUP is a good thing: The estimation of random eects. Statist Sci 1991; 6: 15–51.

18.

McCulloch CE and Searle SR. Generalized, Linear, and Mixed Models. New York: John Wiley & Sons, 2001.

19.

Wang Y. Smoothing spline models with correlated random errors. J Am Stat Assoc 1998; 93: 341–348.

20.

Brumback BA and Rice JA. Smoothing spline models for the analysis of nested and crossed samples of curves. J Am Stat Assoc 1998; 93: 961–976.

Self-modeling ordinal regression with time invariant covariates – An application to prostate cancer data

Abstract

Keywords

Get full access to this article

References