39,389 research outputs found
Maximum likelihood estimation in a partially observed stratified regression model with censored data
The stratified proportional intensity model generalizes Cox's proportional
intensity model by allowing different groups of the population under study to
have distinct baseline intensity functions. In this article, we consider the
problem of estimation in this model when the variable indicating the stratum is
unobserved for some individuals in the studied sample. In this setting, we
construct nonparametric maximum likelihood estimators for the parameters of the
stratified model and we establish their consistency and asymptotic normality.
Consistent estimators for the limiting variances are also obtained
Iterated filtering methods for Markov process epidemic models
Dynamic epidemic models have proven valuable for public health decision
makers as they provide useful insights into the understanding and prevention of
infectious diseases. However, inference for these types of models can be
difficult because the disease spread is typically only partially observed e.g.
in form of reported incidences in given time periods. This chapter discusses
how to perform likelihood-based inference for partially observed Markov
epidemic models when it is relatively easy to generate samples from the Markov
transmission model while the likelihood function is intractable. The first part
of the chapter reviews the theoretical background of inference for partially
observed Markov processes (POMP) via iterated filtering. In the second part of
the chapter the performance of the method and associated practical difficulties
are illustrated on two examples. In the first example a simulated outbreak data
set consisting of the number of newly reported cases aggregated by week is
fitted to a POMP where the underlying disease transmission model is assumed to
be a simple Markovian SIR model. The second example illustrates possible model
extensions such as seasonal forcing and over-dispersion in both, the
transmission and observation model, which can be used, e.g., when analysing
routinely collected rotavirus surveillance data. Both examples are implemented
using the R-package pomp (King et al., 2016) and the code is made available
online.Comment: This manuscript is a preprint of a chapter to appear in the Handbook
of Infectious Disease Data Analysis, Held, L., Hens, N., O'Neill, P.D. and
Wallinga, J. (Eds.). Chapman \& Hall/CRC, 2018. Please use the book for
possible citations. Corrected typo in the references and modified second
exampl
Bayesian Semiparametric Multi-State Models
Multi-state models provide a unified framework for the description of the evolution of discrete phenomena in continuous time. One particular example are Markov processes which can be characterised by a set of time-constant transition intensities between the states. In this paper, we will extend such parametric approaches to semiparametric models with flexible transition intensities based on Bayesian versions of penalised splines. The transition intensities will be modelled as smooth functions of time and can further be related to parametric as well as nonparametric covariate effects. Covariates with time-varying effects and frailty terms can be included in addition. Inference will be conducted either fully Bayesian using Markov chain Monte Carlo simulation techniques or empirically Bayesian based on a mixed model representation. A counting process representation of semiparametric multi-state models provides the likelihood formula and also forms the basis for model validation via martingale residual processes. As an application, we will consider human sleep data with a discrete set of sleep states such as REM and Non-REM phases. In this case, simple parametric approaches are inappropriate since the dynamics underlying human sleep are strongly varying throughout the night and individual-specific variation has to be accounted for using covariate information and frailty terms
Bayesian Semiparametric Multi-State Models
Multi-state models provide a unified framework for the description of the evolution of discrete phenomena in continuous time. One particular example are Markov processes which can be characterised by a set of time-constant transition intensities between the states. In this paper, we will extend such parametric approaches to semiparametric models with flexible transition intensities based on Bayesian versions of penalised splines. The transition intensities will be modelled as smooth functions of time and can further be related to parametric as well as nonparametric covariate effects. Covariates with time-varying effects and frailty terms can be included in addition. Inference will be conducted either fully Bayesian using Markov chain Monte Carlo simulation techniques or empirically Bayesian based on a mixed model representation. A counting process representation of semiparametric multi-state models provides the likelihood formula and also forms the basis for model validation via martingale residual processes. As an application, we will consider human sleep data with a discrete set of sleep states such as REM and Non-REM phases. In this case, simple parametric approaches are inappropriate since the dynamics underlying human sleep are strongly varying throughout the night and individual-specific variation has to be accounted for using covariate information and frailty terms
Causal inference for continuous-time processes when covariates are observed only at discrete times
Most of the work on the structural nested model and g-estimation for causal
inference in longitudinal data assumes a discrete-time underlying data
generating process. However, in some observational studies, it is more
reasonable to assume that the data are generated from a continuous-time process
and are only observable at discrete time points. When these circumstances
arise, the sequential randomization assumption in the observed discrete-time
data, which is essential in justifying discrete-time g-estimation, may not be
reasonable. Under a deterministic model, we discuss other useful assumptions
that guarantee the consistency of discrete-time g-estimation. In more general
cases, when those assumptions are violated, we propose a controlling-the-future
method that performs at least as well as g-estimation in most scenarios and
which provides consistent estimation in some cases where g-estimation is
severely inconsistent. We apply the methods discussed in this paper to
simulated data, as well as to a data set collected following a massive flood in
Bangladesh, estimating the effect of diarrhea on children's height. Results
from different methods are compared in both simulation and the real
application.Comment: Published in at http://dx.doi.org/10.1214/10-AOS830 the Annals of
Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical
Statistics (http://www.imstat.org
- …