Numerical discretization-based estimation methods for ordinary differential equation models via penalized spline smoothing with applications in biomedical research.

Abstract:

:Differential equations are extensively used for modeling dynamics of physical processes in many scientific fields such as engineering, physics, and biomedical sciences. Parameter estimation of differential equation models is a challenging problem because of high computational cost and high-dimensional parameter space. In this article, we propose a novel class of methods for estimating parameters in ordinary differential equation (ODE) models, which is motivated by HIV dynamics modeling. The new methods exploit the form of numerical discretization algorithms for an ODE solver to formulate estimating equations. First, a penalized-spline approach is employed to estimate the state variables and the estimated state variables are then plugged in a discretization formula of an ODE solver to obtain the ODE parameter estimates via a regression approach. We consider three different order of discretization methods, Euler's method, trapezoidal rule, and Runge-Kutta method. A higher-order numerical algorithm reduces numerical error in the approximation of the derivative, which produces a more accurate estimate, but its computational cost is higher. To balance the computational cost and estimation accuracy, we demonstrate, via simulation studies, that the trapezoidal discretization-based estimate is the best and is recommended for practical use. The asymptotic properties for the proposed numerical discretization-based estimators are established. Comparisons between the proposed methods and existing methods show a clear benefit of the proposed methods in regards to the trade-off between computational cost and estimation accuracy. We apply the proposed methods t an HIV study to further illustrate the usefulness of the proposed approaches.

journal_name

Biometrics

journal_title

Biometrics

authors

Wu H,Xue H,Kumar A

doi

10.1111/j.1541-0420.2012.01752.x

subject

Has Abstract

pub_date

2012-06-01 00:00:00

pages

344-52

issue

2

eissn

0006-341X

issn

1541-0420

journal_volume

68

pub_type

杂志文章
  • Ultra high-dimensional semiparametric longitudinal data analysis.

    abstract::As ultra high-dimensional longitudinal data are becoming ever more apparent in fields such as public health and bioinformatics, developing flexible methods with a sparse model is of high interest. In this setting, the dimension of the covariates can potentially grow exponentially as ...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/biom.13348

    authors: Green B,Lian H,Yu Y,Zu T

    更新日期:2020-08-04 00:00:00

  • Statistical analysis of unlabeled point sets: comparing molecules in chemoinformatics.

    abstract::We consider Bayesian methodology for comparing two or more unlabeled point sets. Application of the technique to a set of steroid molecules illustrates its potential utility involving the comparison of molecules in chemoinformatics and bioinformatics. We initially match a pair of molecules, where one molecule is regar...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.1541-0420.2006.00622.x

    authors: Dryden IL,Hirst JD,Melville JL

    更新日期:2007-03-01 00:00:00

  • Semiparametric regression of multidimensional genetic pathway data: least-squares kernel machines and linear mixed models.

    abstract::We consider a semiparametric regression model that relates a normal outcome to covariates and a genetic pathway, where the covariate effects are modeled parametrically and the pathway effect of multiple gene expressions is modeled parametrically or nonparametrically using least-squares kernel machines (LSKMs). This un...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.1541-0420.2007.00799.x

    authors: Liu D,Lin X,Ghosh D

    更新日期:2007-12-01 00:00:00

  • Small-sample inference for the comparison of means of log-normal distributions.

    abstract::We propose a likelihood-based test for comparing the means of two or more log-normal distributions, with possibly unequal variances. A modification to the likelihood ratio test is needed when sample sizes are small. The performance of the proposed procedures is compared with the F-ratio test using Monte Carlo simulati...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.0006-341X.2004.00199.x

    authors: Gill PS

    更新日期:2004-06-01 00:00:00

  • Bias in estimating association parameters for longitudinal binary responses with drop-outs.

    abstract::This paper considers the impact of bias in the estimation of the association parameters for longitudinal binary responses when there are drop-outs. A number of different estimating equation approaches are considered for the case where drop-out cannot be assumed to be a completely random process. In particular, standar...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.0006-341x.2001.00015.x

    authors: Fitzmaurice GM,Lipsitz SR,Molenberghs G,Ibrahim JG

    更新日期:2001-03-01 00:00:00

  • A century of biometrical genetics.

    abstract::We briefly review the major contribution of biometrics to genetics over the last century (population genetic models, familial correlations, segregation analysis, and gene mapping) and current areas of active research and then speculate about what problems will be tackled in the next century. ...

    journal_title:Biometrics

    pub_type: 社论,历史文章,评审

    doi:10.1111/j.0006-341x.2000.00659.x

    authors: Elston RC,Thompson EA

    更新日期:2000-09-01 00:00:00

  • Confidence interval estimation for the ratio of simple and standardized rates in cohort studies.

    abstract::Computer simulation has been used to compare four methods for calculating confidence intervals for simple rate ratios estimated from cohort studies. The method proposed by Cornfield (1956. In Proceedings of the Third Berkeley Symposium on Mathematical Statistics and Probability. Vol. IV, 135-148) for interval estimati...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:

    authors: Howe GR

    更新日期:1983-06-01 00:00:00

  • Growth curve models of repeated binary response.

    abstract::Experimental designs that include repeated measures of binary response variables over time and under different conditions are common in biology. In such settings, it is often desirable to characterize the response pattern over time. When response variables are continuous, this characterization can be made in terms of ...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:

    authors: Stanek EJ 3rd,Diehl SR

    更新日期:1988-12-01 00:00:00

  • Subsampling versus bootstrapping in resampling-based model selection for multivariable regression.

    abstract::In recent years, increasing attention has been devoted to the problem of the stability of multivariable regression models, understood as the resistance of the model to small changes in the data on which it has been fitted. Resampling techniques, mainly based on the bootstrap, have been developed to address this issue....

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/biom.12381

    authors: De Bin R,Janitza S,Sauerbrei W,Boulesteix AL

    更新日期:2016-03-01 00:00:00

  • Estimating predictors for long- or short-term survivors.

    abstract::Suppose that the response variable in a well-executed clinical or observational study to evaluate a treatment is the time to a certain event, and a set of baseline covariates or predictors was collected for each study patient. Furthermore, suppose that a significant number of study patients had nontrivial, long-term a...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.0006-341x.2003.00116.x

    authors: Tian L,Wang W,Wei LJ

    更新日期:2003-12-01 00:00:00

  • Functional multiple indicators, multiple causes measurement error models.

    abstract::Objective measures of oxygen consumption and carbon dioxide production by mammals are used to predict their energy expenditure. Since energy expenditure is not directly observable, it can be viewed as a latent construct with multiple physical indirect measures such as respiratory quotient, volumetric oxygen consumptio...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/biom.12706

    authors: Tekwe CD,Zoh RS,Bazer FW,Wu G,Carroll RJ

    更新日期:2018-03-01 00:00:00

  • A Markov model for analysing cancer markers and disease states in survival studies.

    abstract::In studies of serial cancer markers or disease states and their relation to survival, data on the marker or state are usually obtained at infrequent time points during follow-up. A Markov model is developed to assess the dependence of risk of death on marker level or disease state and inferences within this model are ...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:

    authors: Kay R

    更新日期:1986-12-01 00:00:00

  • Bayesian lasso for semiparametric structural equation models.

    abstract::There has been great interest in developing nonlinear structural equation models and associated statistical inference procedures, including estimation and model selection methods. In this paper a general semiparametric structural equation model (SSEM) is developed in which the structural equation is composed of nonpar...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.1541-0420.2012.01751.x

    authors: Guo R,Zhu H,Chow SM,Ibrahim JG

    更新日期:2012-06-01 00:00:00

  • Novel two-phase sampling designs for studying binary outcomes.

    abstract::In biomedical cohort studies for assessing the association between an outcome variable and a set of covariates, usually, some covariates can only be measured on a subgroup of study subjects. An important design question is-which subjects to select into the subgroup to increase statistical efficiency. When the outcome ...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/biom.13140

    authors: Wang L,Williams ML,Chen Y,Chen J

    更新日期:2020-03-01 00:00:00

  • Estimation of individual genetic effects from binary observations on relatives applied to a family history of respiratory illnesses and chronic lung disease of newborns.

    abstract::This paper considers methods for estimating the relationship between a binary response Y and the genetic effects responsible for a second binary trait Z. The responses Y are observed only for target individuals, and the responses Z are observed only for the relatives of these targets. The analysis consists of two part...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.0006-341x.2000.00808.x

    authors: Houwing-Duistermaat JJ,van Houwelingen HC,de Winter JP

    更新日期:2000-09-01 00:00:00

  • On logit confidence intervals for the odds ratio with small samples.

    abstract::Unless the true association is very strong, simple large-sample confidence intervals for the odds ratio based on the delta method perform well even for small samples. Such intervals include the Woolf logit interval and the related Gart interval based on adding .5 before computing the log odds ratio estimate and its st...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.0006-341x.1999.00597.x

    authors: Agresti A

    更新日期:1999-06-01 00:00:00

  • Time-varying functional regression for predicting remaining lifetime distributions from longitudinal trajectories.

    abstract::A recurring objective in longitudinal studies on aging and longevity has been the investigation of the relationship between age-at-death and current values of a longitudinal covariate trajectory that quantifies reproductive or other behavioral activity. We propose a novel technique for predicting age-at-death distribu...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.1541-0420.2005.00378.x

    authors: Müller HG,Zhang Y

    更新日期:2005-12-01 00:00:00

  • Group sequential tests for bivariate response: interim analyses of clinical trials with both efficacy and safety endpoints.

    abstract::We describe group sequential tests for a bivariate response. The tests are defined in terms of the two response components jointly, rather than through a single summary statistic. Such methods are appropriate when the two responses concern different aspects of a treatment; for example, one might wish to show that a ne...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:

    authors: Jennison C,Turnbull BW

    更新日期:1993-09-01 00:00:00

  • Some aspects of cancer epidemiology.

    abstract::Epidemiolgic studies have strongly suggested that a vast majority (80-90%) of cancers are caused by radiation, chemical and biologic agents; the remainder result from endogenous or genetic factors. Biologically, cancer is most probably the end result of a complex multistage process and therefore may be due to a sequen...

    journal_title:Biometrics

    pub_type: 杂志文章,评审

    doi:

    authors: Lilienfeld AM

    更新日期:1982-03-01 00:00:00

  • Nonparametric discrete survival function estimation with uncertain endpoints using an internal validation subsample.

    abstract::When a true survival endpoint cannot be assessed for some subjects, an alternative endpoint that measures the true endpoint with error may be collected, which often occurs when obtaining the true endpoint is too invasive or costly. We develop an estimated likelihood function for the situation where we have both uncert...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/biom.12316

    authors: Zee J,Xie SX

    更新日期:2015-09-01 00:00:00

  • Unbalanced regression analysis with residuals having a covariance structure of intra-class form.

    abstract::Let Yi be an ni X 1 vector of observations, Xi an ni X p matrix of known values, and beta an unknown p X 1 with the structure Yi = Xi beta + epsilon i, where the covariance matrix of epsilon i is of intra-class form, that is Cov (epsilon i) = sigma2[(1 - rho) Ii + rho e i e i'] where Ii is the ni X ni identity matrix ...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:

    authors: Wiorkowski JJ

    更新日期:1975-09-01 00:00:00

  • Fitting a multiplicative incidence model to age- and time-specific prevalence data.

    abstract::We discuss the assessment of age- and time-specific disease incidence using prevalence data. A method is described for conveniently fitting a discrete-time multiplicative model, subject to positivity constraints, using the EM-algorithm. Together with smoothing, it allows essentially nonparametric assessment of inciden...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:

    authors: Marschner IC

    更新日期:1996-06-01 00:00:00

  • Catch estimation with restricted randomization in the effort survey.

    abstract::One common method for estimating total catch is to multiply an estimate for CPUE, the catch per unit effort, by an estimate of total effort obtained from an independent second survey. In general, estimating total effort requires that sample times are chosen at random over the full fishing period; however, in practice,...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.0006-341x.2001.00461.x

    authors: Dauk PC,Schwarz CJ

    更新日期:2001-06-01 00:00:00

  • Sharpening bounds on principal effects with covariates.

    abstract::Estimation of treatment effects in randomized studies is often hampered by possible selection bias induced by conditioning on or adjusting for a variable measured post-randomization. One approach to obviate such selection bias is to consider inference about treatment effects within principal strata, that is, principal...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/biom.12103

    authors: Long DM,Hudgens MG

    更新日期:2013-12-01 00:00:00

  • A likelihood ratio statistic for testing goodness of fit with randomly censored data.

    abstract::A likelihood ratio statistic is proposed for testing goodness of fit with grouped data which are subject to random right censoring. It is shown that, under appropriate conditions, this statistic has an asymptotic chi-square distribution which is non-central under contiguous alternatives. A formula is given for the non...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:

    authors: Turnbull BW,Weiss L

    更新日期:1978-09-01 00:00:00

  • Use of historical marker data for assessing treatment effects in phase I/II trials when subject selection is determined by baseline marker level.

    abstract::Although the primary focus of Phase I clinical trials is to assess clinical pharmacology and possible toxicities, any information on the potential effect of treatment would be useful in helping to determine priorities between treatments for further study. We consider the scenario where data are routinely collected on ...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:

    authors: Lin HM,Hughes MD

    更新日期:1995-09-01 00:00:00

  • Improved dynamic predictions from joint models of longitudinal and survival data with time-varying effects using P-splines.

    abstract::In the field of cardio-thoracic surgery, valve function is monitored over time after surgery. The motivation for our research comes from a study which includes patients who received a human tissue valve in the aortic position. These patients are followed prospectively over time by standardized echocardiographic assess...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/biom.12814

    authors: Andrinopoulou ER,Eilers PHC,Takkenberg JJM,Rizopoulos D

    更新日期:2018-06-01 00:00:00

  • A test of homogeneity of distributions when observations are subject to measurement errors.

    abstract::When the observed data are contaminated with errors, the standard two-sample testing approaches that ignore measurement errors may produce misleading results, including a higher type-I error rate than the nominal level. To tackle this inconsistency, a nonparametric test is proposed for testing equality of two distribu...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/biom.13207

    authors: Lee D,Lahiri SN,Sinha S

    更新日期:2020-09-01 00:00:00

  • Linear mixed models with flexible distributions of random effects for longitudinal data.

    abstract::Normality of random effects is a routine assumption for the linear mixed model, but it may be unrealistic, obscuring important features of among-individual variation. We relax this assumption by approximating the random effects density by the seminonparameteric (SNP) representation of Gallant and Nychka (1987, Econome...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.0006-341x.2001.00795.x

    authors: Zhang D,Davidian M

    更新日期:2001-09-01 00:00:00

  • Model checking for ROC regression analysis.

    abstract::The receiver operating characteristic (ROC) curve is a prominent tool for characterizing the accuracy of a continuous diagnostic test. To account for factors that might influence the test accuracy, various ROC regression methods have been proposed. However, as in any regression analysis, when the assumed models do not...

    journal_title:Biometrics

    pub_type: 杂志文章

    doi:10.1111/j.1541-0420.2006.00620.x

    authors: Cai T,Zheng Y

    更新日期:2007-03-01 00:00:00