Residual analysis for parametric models fit to censored data

The focus of this paper is on residual analysis for the lognormal and extreme value or Weibull models, although the proposed methods can be applied to any parametric model. Residuals developed by Barlow and Prentice (1988) for the Cox proportional hazards model are extended to the parametric model setting. Three different residuals are proposed based on this approach with two residuals measuring the impact of survival time and one measuring the impact of the covariates included in the model. In addition, a residual derived from the deviations equality presented in Efron and Johnstone (1990) and the residual proposed by Joergensen (1984) for censored data models are discussed.     

Integrated likelihoods in parametric survival models for highly clustered censored data

In studies that involve censored time-to-event data, stratification is frequently encountered due to different reasons, such as stratified sampling or model adjustment due to violation of model assumptions. Often, the main interest is not in the clustering variables, and the cluster-related parameters are treated as nuisance. When inference is about a parameter of interest in presence of many nuisance parameters, standard likelihood methods often perform very poorly and may lead to severe bias. This problem is particularly evident in models for clustered data with cluster-specific nuisance parameters, when the number of clusters is relatively high with respect to the within-cluster size. However, it is still unclear how the presence of censoring would affect this issue. We consider clustered failure time data with independent censoring, and propose frequentist inference based on an integrated likelihood. We then apply the proposed approach to a stratified Weibull model. Simulation studies show that appropriately defined integrated likelihoods provide very accurate inferential results in all circumstances, such as for highly clustered data or heavy censoring, even in extreme settings where standard likelihood procedures lead to strongly misleading results. We show that the proposed method performs generally as well as the frailty model, but it is superior when the frailty distribution is seriously misspecified. An application, which concerns treatments for a frequent disease in late-stage HIV-infected people, illustrates the proposed inferential method in Weibull regression models, and compares different inferential conclusions from alternative methods.     

Bayesian estimation of reliability using time-to-failure distribution of parametric degradation models

This paper focusses on computing the Bayesian reliability of components whose performance characteristics (degradation – fatigue and cracks) are observed during a specified period of time. Depending upon the nature of degradation data collected, we fit a monotone increasing or decreasing function for the data. Since the components are supposed to have different lifetimes, the rate of degradation is assumed to be a random variable. At a critical level of degradation, the time to failure distribution is obtained. The exponential and power degradation models are studied and exponential density function is assumed for the random variable representing the rate of degradation. The maximum likelihood estimator and Bayesian estimator of the parameter of exponential density function, predictive distribution, hierarchical Bayes approach and robustness of the posterior mean are presented. The Gibbs sampling algorithm is used to obtain the Bayesian estimates of the parameter. Illustrations are provided for the train wheel degradation data.     

Longitudinal data analysis of animal growth via multivariate dynamic models

We propose models to analyze animal growth data with the aim of estimating and predicting quantities of biological and economical interest such as the maturing rate and asymptotic weight. It is also studied the effect of environmental factors of relevant influence in the growth process. The models considered in this paper are based on an extension and specialization of the dynamic hierarchical model (Gamerman & Migon, 1993) to a non–linear growth curve setting, where some of the growth curve parameters are considered exchangeable among the units. The inference for these models are approximate conjugate analysis based on Taylor series expansions and linear Bayes procedures     

Flexible models for non-equidispersed count data: comparative performance of parametric models to deal with underdispersion

AStA Advances in Statistical Analysis - Count data as response variables are commonly modeled using Poisson regression models, which require equidispersion, i.e., equal mean and variance. However,...     

Estimating parametric semi-Markov models from panel data using phase-type approximations

Inference for semi-Markov models under panel data presents considerable computational difficulties. In general the likelihood is intractable, but a tractable likelihood with the form of a hidden Markov model can be obtained if the sojourn times in each of the states are assumed to have phase-type distributions. However, using phase-type distributions directly may be undesirable as they require estimation of parameters which may be poorly identified. In this article, an approach to fitting semi-Markov models with standard parametric sojourn distributions is developed. The method involves establishing a family of Coxian phase-type distribution approximations to the parametric distribution and merging approximations for different states to obtain an approximate semi-Markov process with a tractable likelihood. Approximations are developed for Weibull and Gamma distributions and demonstrated on data relating to post-lung-transplantation patients.     

Non parametric hypothesis tests for comparing reliability functions

In reliability and related disciplines, comparing reliability functions of two (or more) aging processes is a crucial step in the process of determining reliability and understanding an aging process. The aim of this paper is to propose a non parametric statistical methodology to compare two populations based on their mean residual life function and expected inactivity time function. We introduce some novel hypothesis testing procedures that involve both Cramér–von Mises- and Kolmogorov–Smirnov-type test statistics and their decision rules are constructed based on the asymptotic distributions of these test statistics and bootstrapping method. We study the practical behavior of the proposed testing procedures extensively through simulations. The results reveal that the proposed hypothesis testing procedures perform efficiently in identifying small and large differences. Two real-life examples are discussed to demonstrate the practical utility of the tests.     

A survey of time series analysis through parametric models

The paper is a survey of recent works on time series analysis using parametric models. The main emphasis is on linear models, in particular the ARMA model. Usual me¬thods of parameter estimation, goodness of fit tests and the choice of model order are con¬sidered. Some extensions of the methods to related problems are briefly discussed     

Non parametric regression analysis for longitudinal data with time-depending autoregressive error process

This paper considers a non parametric longitudinal model, where the within-subject correlation structure is represented by a time-depending autoregressive error process. An initial estimator without taking into account the within-subject correlation is obtained to fit the time-depending autoregressive error process. With the initial estimator, we construct a two-stage local linear estimator of the mean function. According to the asymptotic normality of the initial and two-stage estimators, it is discovered that the two-stage estimator has a smaller asymptotic variance. The simulation results show us that the two-stage estimation has some good properties. The analysis of a data set demonstrates its application.     

Bayesian analysis of reliability in multicomponent stress-strength models

This paper provides a Bayesian treatment of the problem of inference about the reliability of a multicomponent stress-strength system which functions if s or more of k identical components simultaneously operate. All stresses and strengths are assumed to be independent, exponentially distributed random variables. Exact and approximate asymptotic posterior distributions for the reliability are derived, and the results are illustrated by a numerical example.     

Estimation of zero-inflated parameter-driven models via data cloning

ABSTRACT

In this paper, we propose the use of the Data Cloning (DC) approach to estimate parameter-driven zero-inflated Poisson and Negative Binomial models for time series of counts. The data cloning algorithm obtains the familiar maximum likelihood estimators and their standard errors via a fully Bayesian estimation. This provides some computational ease as well as inferential tools such as confidence intervals and diagnostic methods which, otherwise, are not readily available for parameter-driven models. To illustrate the performance of the proposed method, we use Monte Carlo Simulations and real data on asthma-related emergency department visits in the Canadian province of Ontario.     

Composable models for online Bayesian analysis of streaming data

Data is rapidly increasing in volume and velocity and the Internet of Things (IoT) is one important source of this data. The IoT is a collection of connected devices (things) which are constantly recording data from their surroundings using on-board sensors. These devices can record and stream data to the cloud at a very high rate, leading to high storage and analysis costs. In order to ameliorate these costs, the data is modelled as a stream and analysed online to learn about the underlying process, perform interpolation and smoothing and make forecasts and predictions. Conventional state space modelling tools assume the observations occur on a fixed regular time grid. However, many sensors change their sampling frequency, sometimes adaptively, or get interrupted and re-started out of sync with the previous sampling grid, or just generate event data at irregular times. It is therefore desirable to model the system as a partially and irregularly observed Markov process which evolves in continuous time. Both the process and the observation model are potentially non-linear. Particle filters therefore represent the simplest approach to online analysis. A functional Scala library of composable continuous time Markov process models has been developed in order to model the wide variety of data captured in the IoT.     

Markov chain models for multivariate repeated binary data analysis

Repeated categorical outcomes frequently occur in clinical trials. Muenz and Rubinstein (1985) presented Markov chain models to analyze binary repeated data in a breast cancer study. We extend their method to the setting when more than one repeated outcome variable is of interest. In a randomized clinical trial of breast cancer, we investigate the dependency of toxicities on predictor variables and the relationship among multiple toxic effects.     

Ignorability conditions for frequentist non parametric analysis of conditional distributions with incomplete data

Rubin (1976 Rubin, D.B. (1976). Inference and missing data. Biometrika 63(3):581–592.[Crossref], [Web of Science ®] , [Google Scholar]) derived general conditions under which inferences that ignore missing data are valid. These conditions are sufficient but not generally necessary, and therefore may be relaxed in some special cases. We consider here the case of frequentist estimation of a conditional cdf subject to missing outcomes. We partition a set of data into outcome, conditioning, and latent variables, all of which potentially affect the probability of a missing response. We describe sufficient conditions under which a complete-case estimate of the conditional cdf of the outcome given the conditioning variable is unbiased. We use simulations on a renal transplant data set (Dienemann et al.) to illustrate the implications of these results.     

Goodness-of-fit tests for parametric models in censored regression

The authors propose a goodness-of-fit test for parametric regression models when the response variable is right-censored. Their test compares an estimation of the error distribution based on parametric residuals to another estimation relying on nonparametric residuals. They call on a bootstrap mechanism in order to approximate the critical values of tests based on Kolmogorov-Smirnov and Cramér-von Mises type statistics. They also present the results of Monte Carlo simulations and use data from a study about quasars to illustrate their work.     

An adaptive-to-model test for partially parametric single-index models

Residual marked empirical process-based tests are commonly used in regression models. However, they suffer from data sparseness in high-dimensional space when there are many covariates. This paper has three purposes. First, we suggest a partial dimension reduction adaptive-to-model testing procedure that can be omnibus against general global alternative models although it fully use the dimension reduction structure under the null hypothesis. This feature is because that the procedure can automatically adapt to the null and alternative models, and thus greatly overcomes the dimensionality problem. Second, to achieve the above goal, we propose a ridge-type eigenvalue ratio estimate to automatically determine the number of linear combinations of the covariates under the null and alternative hypotheses. Third, a Monte-Carlo approximation to the sampling null distribution is suggested. Unlike existing bootstrap approximation methods, this gives an approximation as close to the sampling null distribution as possible by fully utilising the dimension reduction model structure under the null model. Simulation studies and real data analysis are then conducted to illustrate the performance of the new test and compare it with existing tests.     

A lack-of-fit test for parametric zero-inflated Poisson models

Count data often contain many zeros. In parametric regression analysis of zero-inflated count data, the effect of a covariate of interest is typically modelled via a linear predictor. This approach imposes a restrictive, and potentially questionable, functional form on the relation between the independent and dependent variables. To address the noted restrictions, a flexible parametric procedure is employed to model the covariate effect as a linear combination of fixed-knot cubic basis splines or B-splines. The semiparametric zero-inflated Poisson regression model is fitted by maximizing the likelihood function through an expectation–maximization algorithm. The smooth estimate of the functional form of the covariate effect can enhance modelling flexibility. Within this modelling framework, a log-likelihood ratio test is used to assess the adequacy of the covariate function. Simulation results show that the proposed test has excellent power in detecting the lack of fit of a linear predictor. A real-life data set is used to illustrate the practicality of the methodology.