Show simple item record

dc.contributor.authorXun, Xiaolei
dc.contributor.authorCao, Jiguo
dc.contributor.authorMallick, Bani
dc.contributor.authorMaity, Arnab
dc.contributor.authorCarroll, Raymond J.
dc.date.accessioned2016-02-25T13:53:34Z
dc.date.available2016-02-25T13:53:34Z
dc.date.issued2013-09
dc.identifier.citationXun X, Cao J, Mallick B, Maity A, Carroll RJ (2013) Parameter Estimation of Partial Differential Equation Models. Journal of the American Statistical Association 108: 1009–1020. Available: http://dx.doi.org/10.1080/01621459.2013.794730.
dc.identifier.issn0162-1459
dc.identifier.issn1537-274X
dc.identifier.pmid24363476
dc.identifier.doi10.1080/01621459.2013.794730
dc.identifier.urihttp://hdl.handle.net/10754/599138
dc.description.abstractPartial differential equation (PDE) models are commonly used to model complex dynamic systems in applied sciences such as biology and finance. The forms of these PDE models are usually proposed by experts based on their prior knowledge and understanding of the dynamic system. Parameters in PDE models often have interesting scientific interpretations, but their values are often unknown and need to be estimated from the measurements of the dynamic system in the presence of measurement errors. Most PDEs used in practice have no analytic solutions, and can only be solved with numerical methods. Currently, methods for estimating PDE parameters require repeatedly solving PDEs numerically under thousands of candidate parameter values, and thus the computational load is high. In this article, we propose two methods to estimate parameters in PDE models: a parameter cascading method and a Bayesian approach. In both methods, the underlying dynamic process modeled with the PDE model is represented via basis function expansion. For the parameter cascading method, we develop two nested levels of optimization to estimate the PDE parameters. For the Bayesian method, we develop a joint model for data and the PDE and develop a novel hierarchical model allowing us to employ Markov chain Monte Carlo (MCMC) techniques to make posterior inference. Simulation studies show that the Bayesian method and parameter cascading method are comparable, and both outperform other available methods in terms of estimation accuracy. The two methods are demonstrated by estimating parameters in a PDE model from long-range infrared light detection and ranging data. Supplementary materials for this article are available online. © 2013 American Statistical Association.
dc.description.sponsorshipThe research of Mallick, Carroll, and Xun was supported by grants from the National Cancer Institute (R37-CA057030) and the National Science Foundation DMS (Division of Mathematical Sciences) grant 0914951. This publication is based in part on work supported by the Award Number KUS-CI-016-04, made by King Abdullah University of Science and Technology (KAUST). Can's research is supported by a discovery grant (PIN: 328256) from the Natural Science and Engineering Research Council of Canada (NSERC). Maity's research was performed while visiting the Department of Statistics, Texas A&M University, and was partially supported by the Award Number R00ES017744 from the National Institute of Environmental Health Sciences.
dc.publisherInforma UK Limited
dc.subjectAsymptotic theory
dc.subjectBasis function expansion
dc.subjectBayesian method
dc.subjectDifferential equations
dc.subjectMeasurement error
dc.subjectParameter cascading
dc.titleParameter Estimation of Partial Differential Equation Models
dc.typeArticle
dc.identifier.journalJournal of the American Statistical Association
dc.identifier.pmcidPMC3867159
dc.contributor.institutionNovartis Pharma, Rueil-Malmaison, France
dc.contributor.institutionSimon Fraser University, Burnaby, Canada
dc.contributor.institutionTexas A and M University, College Station, United States
dc.contributor.institutionNorth Carolina State University, Raleigh, United States
kaust.grant.numberKUS-CI-016-04


This item appears in the following Collection(s)

Show simple item record