×

Parameter estimation for semiparametric ordinary differential equation models. (English) Zbl 07529903

Summary: We propose a new class of two-stage parameter estimation methods for semiparametric ordinary differential equation (ODE) models. In the first stage, state variables are estimated using a penalized spline approach; In the second stage, form of numerical discretization algorithms for an ODE solver is used to formulate estimating equations. Estimated state variables from the first stage are used to obtain more data points for the second stage. Asymptotic properties for the proposed estimators are established. Simulation studies show that the method performs well, especially for small sample. Real life use of the method is illustrated using Influenza specific cell-trafficking study.

MSC:

62-XX Statistics

Software:

DEDiscover; SemiPar
PDFBibTeX XMLCite
Full Text: DOI Link

References:

[1] Bansal, M.; Gatta, G. D.; di Bernardo, D., Inference of gene regulatory networks and compound mode of action from time course gene expression profiles, Bioinformatics, 22, 815-22 (2006)
[2] Bard, Y., Nonlinear parameter estimation (1974), New York: Academic Press, New York · Zbl 0345.62045
[3] Brunel, N., Parameter estimation of ODE’s via nonparametric estimators, Electronic Journal of Statistics, 2, 1242-67 (2008) · Zbl 1320.62063
[4] Cao, J.; Huang, J. Z.; Wu, H., Penalized nonlinear least squares estimation of time-varying parameters in ordinary differential equations, Journal of Computational and Graphical Statistics, 21, 1, 42-56 (2012)
[5] Chen, J.; Wu, H., Efficient local estimation for time-varying coefficients in deterministic dynamic models with applications to HIV-1 dynamics, Journal of the American Statistical Association, 103, 481, 369-84 (2008) · Zbl 1469.62365
[6] Chen, J.; Wu, H., Estimation of time-varying parameters in deterministic dynamic models with application to HIV infections, Statistica Sinica, 18, 987-1006 (2008) · Zbl 1149.62079
[7] Cheng, G.; Huang, J. Z., Bootstrap consistency for general semiparametric M-estimation, The Annals of Statistics, 38, 5, 2884-915 (2010) · Zbl 1200.62042
[8] Claeskens, G.; Krivobokova, T.; Opsomer, J., Asymptotic properties of penalized spline estimators, Biometrika, 96, 3, 529-44 (2009) · Zbl 1170.62031
[9] Craven, P.; Wahba, G., Smoothing noisy data with spline functions: estimating the correct degree of smoothing by the generalized cross-validation, Numerische Mathematik, 31, 4, 377-403 (1979) · Zbl 0377.65007
[10] D’Haeseleer, P.; Wen, X.; Fuhrman, S.; Somogyi, R., Linear modeling of MRNA expression levels during CNS development and injury, Pacific Symposium on Biocomputing, 4, 41-52 (1999)
[11] Donnet, S.; Samson, A., Estimation of parameters in incomplete data models defined by dynamical systems, Journal of Statistical Planning and Inference, 137, 9, 2815-31 (2007) · Zbl 1331.62099
[12] Efron, B., The jackknife, the bootstrap and other resampling plans (1982), Philadelphia: SIAM, Philadelphia · Zbl 0496.62036
[13] Gugushvili, S.; Klaassen, C. A. J., \(####\)-consistent parameter estimation for systems of ordinary differential equations: Bypassing numerical integration via smoothing, Bernoulli, 18, 3, 1061-98 (2012) · Zbl 1257.49033
[14] Hairer, E.; Norsett, S.; Wanner, G., Solving ordinary differential equation I: Nonstiff problems (1993), Berlin: Springer-Veralag, Berlin · Zbl 0789.65048
[15] He, X.; Xue, H.; Shi, N. Z., Sieve maximum likelihood estimation for doubly semiparametric zero-inflated Poisson models, Journal of Multivariate Analysis, 101, 9, 2026-38 (2010) · Zbl 1203.62067
[16] Huang, Y.; Liu, D.; Wu, H., Hierarchical Bayesian methods for estimation of parameters in a longitudinal HIV dynamic system, Biometrics, 62, 2, 413-23 (2006) · Zbl 1097.62128
[17] Joshi, M.; Seidel-Morgenstern, A.; Kremling, A., Exploiting the boostrap method for quantifying parameter confidence intervals in dynamic systems, Metabolic Engineering, 8, 5, 447-55 (2006)
[18] Li, L.; Brown, M. B.; Lee, K. H.; Gupta, S., Estimation and inference for a spline-enhanced population pharmacokinetic model, Biometrics, 58, 3, 601-11 (2002) · Zbl 1210.62179
[19] Li, Y.; Ruppert, D., On the asymptotics of penalized splines, Biometrika, 95, 2, 415-36 (2008) · Zbl 1437.62540
[20] Li, Z.; Osborne, M. R.; Prvan, T., Parameter estimation of ordinary differential equations, IMA Journal of Numerical Analysis, 25, 2, 264-85 (2005) · Zbl 1070.65061
[21] Liang, H.; Miao, H.; Wu, H., Estimation of constant and time-varying dynamic parameters of HIV infection in a nonlinear differential equation model, The Annals of Applied Statistics, 4, 1, 460-83 (2010) · Zbl 1189.62171
[22] Liang, H.; Wu, H., Parameter estimation for differential equation models using a framework of measurement error in regression, Journal of the American Statistical Association, 103, 484, 1570-83 (2008) · Zbl 1286.62039
[23] Ma, S.; Kosorok, M. R., Robust semiparametric M-estimation and the weighted bootstrap, Journal of Multivariate Analysis, 96, 1, 190-217 (2005) · Zbl 1073.62030
[24] Matteij, R.; Molenaar, J., Ordinary differential equations in theory and practice (2002), Philadelphia: SIAM, Philadelphia · Zbl 1016.34001
[25] Putter, H.; Heisterkamp, S.; Lange, J.; Wolf, F., A Bayesian approach to parameter estimation in HIV dynamic models, Statistics in Medicine, 21, 15, 2199-214 (2002)
[26] Qi, X.; Zhao, H., Asymptotic efficiency and finite-sample properties of the generalized profiling estimation of parameters in ordinary differential equations, The Annals of Statistics, 38, 1, 435-81 (2010) · Zbl 1181.62156
[27] Ramsay, J. O., Principal differential analysis: data reduction by differential operators, Journal of the Royal Statistical Society: Series B, 58, 495-508 (1996) · Zbl 0853.62043
[28] Ramsay, J. O.; Hooker, G.; Campbell, D.; Cao, J., Parameter estimation for differential equations: A generalized smoothing approach (with discussion), Journal of the Royal Statistical Society: Series B (Statistical Methodology), 69, 5, 741-96 (2007) · Zbl 07555374
[29] Ruppert, D.; Wand, M. P.; Carroll, R. J., Semiparametric regression (2003), Cambridge: Cambridge University Press, Cambridge · Zbl 1038.62042
[30] Schumaker, L. L., Spline functions: Basis theory (1981), New York: Wiley, New York · Zbl 0449.41004
[31] Shen, X., On methods of sieves and penalization, The Annals of Statistics, 25, 6, 2555-91 (1997) · Zbl 0895.62041
[32] Shen, X.; Wong, W. H., Convergence rate of sieve estimates, The Annals of Statistics, 22, 2, 580-615 (1994) · Zbl 0805.62008
[33] Thomaseth, K.; Kautzky-Willer, A.; Ludvik, B.; Prager, R.; Pacini, G., Integrated mathematical model to assess β-cell activity during the oral glucose test, American Journal of Physiology, 270, 3, E522-31 (1996)
[34] van der Vaart, A. W.; Wellner, J. A., Weak convergence and empirical processes (1996), New York: Springer-Verlag, New York · Zbl 0862.60002
[35] Varah, J., A spline least squares method for numerical parameter estimation in differential equations, SIAM Journal on Scientific and Statistical Computing, 3, 1, 28-46 (1982) · Zbl 0481.65050
[36] Wellner, J. A.; Zhang, Y., Two likelihood-based semiparametric estimation methods for panel count data with covariates, The Annals of Statistics, 35, 5, 2106-42 (2007) · Zbl 1126.62084
[37] Wessels, L. F.; van Someren, E. P.; Reinders, M. J., A comparison of genetic network models, Pacific Symposium on Biocomputing, 6, 508-19 (2001)
[38] Wu, H.; Xue, H.; Kumar, A., Numerical discretization-based estimation methods for ordinary differential equation models via penalized spline smoothing, Biometrics, 68, 2, 344-52 (2012) · Zbl 1274.62906
[39] Wu, H.; Kumar, A.; Miao, H.; Holden-Wiltse, J.; Mosmann, T. R.; Livingstone, A. M.; Belz, G. T.; Perelson, A. S.; Zand, M. S.; Topham, D. J., Modeling of Influenza-Specific CD8^+ T cells during the primary response indicates that the spleen is a major source of effectors, The Journal of Immunology, 187, 9, 4474-82 (2011)
[40] Xue, H.; Lam, K. F.; Li, G., Sieve maximum likelihood estimator for semiparametric regression models with current status data, Journal of the American Statistical Association, 99, 466, 346-56 (2004) · Zbl 1117.62449
[41] Xue, H.; Miao, H.; Wu, H., Sieve estimation of constant and time-varying coefficients in nonlinear ordinary differential equation models by considering both numerical error and measurement error, The Annals of Statistics, 38, 4, 2351-87 (2010) · Zbl 1203.62049
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.