Estimation of the cumulative incidence function under multiple dependent and independent censoring mechanisms

Competing risks occur in a time-to-event analysis in which a patient can experience one of several types of events. Traditional methods for handling competing risks data presuppose one censoring process, which is assumed to be independent. In a controlled clinical trial, censoring can occur for several reasons: some independent, others dependent. We propose an estimator of the cumulative incidence function in the presence of both independent and dependent censoring mechanisms. We rely on semi-parametric theory to derive an augmented inverse probability of censoring weighted (AIPCW) estimator. We demonstrate the efficiency gained when using the AIPCW estimator compared to a non-augmented estimator via simulations. We then apply our method to evaluate the safety and efficacy of three anti-HIV regimens in a randomized trial conducted by the AIDS Clinical Trial Group, ACTG A5095.     

Modeling time series data with semi-reflective boundaries

Time series data are increasingly common in many areas of the health sciences, and in some instances, may have natural boundaries serving as performance guidelines or as thresholds associated with adverse outcomes. Such boundaries may be labeled as semi-reflective, in that the time series values have an increased chance of returning towards middle levels as the boundaries are approached, but boundaries can still be breached. In this paper we review a model that was previously proposed for such data and we investigate its statistical properties. Specifically, this model consists of a third-order auto-regressive projection component, parameterized as a constrained linear combination of linear, flat, and quadratic trends, and an error term that uses a logistic regression model for its sign. We describe and compare a previously-proposed estimation method with a modified version thereof, using computer simulations, as well as data examples from heart monitoring and from a driving simulator. We find that the two methods tend to give different results, with the modified technique having lower bias and more accurate confidence intervals than the previously-proposed method.     

Bootstrapping mean‐squared errors of robust small‐area estimators: Application to the method‐of‐payments surveys data

This paper proposes a new bootstrap procedure for mean‐squared errors of robust small‐area estimators. We formally prove the asymptotic validity of the proposed bootstrap method and examine its finite‐sample performance through Monte Carlo simulations. The results show that our procedure performs well and competes with existing ones. We also provide an application to the estimation of the total volume and value of cash, debit card, and credit card transactions in Canada as well as in its provinces and subgroups of households. In particular, we found that there is a significant average annual decline rate of 3.1% in the volume of cash transactions and that this decline is relatively higher among high‐income households living in heavily populated provinces. Our bootstrap estimator also provides indicators of quality useful in selecting the best small‐area predictor among several alternatives in practice.     

Non-nested model selection based on the quantiles and it’s application in time series

We consider the problem of model selection based on quantile analysis and with unknown parameters estimated using quantile leasts squares. We propose a model selection test for the null hypothesis that the competing models are equivalent against the alternative hypothesis that one model is closer to the true model. We follow with two applications of the proposed model selection test. The first application is in model selection for time series with non-normal innovations. The second application is in model selection in the NoVas method, short for normalizing and variance stabilizing transformation, forecast. A set of simulation results also lends strong support to the results presented in the paper.     

CUSUM Control Charts for Multivariate Poisson Distribution

A cumulative sum control chart for multivariate Poisson distribution (MP-CUSUM) is proposed. The MP-CUSUM chart is constructed based on log-likelihood ratios with in-control parameters, Θ₀, and shifts to be detected quickly, Θ₁. The average run length (ARL) values are obtained using a Markov Chain-based method. Numerical experiments show that the MP-CUSUM chart is effective in detecting parameter shifts in terms of ARL. The MP-CUSUM chart with smaller Θ₁ is more sensitive than that with greater Θ₁ to smaller shifts, but more insensitive to greater shifts. A comparison shows that the proposed MP-CUSUM chart outperforms an existing MP chart.     

Robust models in probability sampling

In the estimation of a population mean or total from a random sample, certain methods based on linear models are known to be automatically design consistent, regardless of how well the underlying model describes the population. A sufficient condition is identified for this type of robustness to model failure; the condition, which we call 'internal bias calibration', relates to the combination of a model and the method used to fit it. Included among the internally bias-calibrated models, in addition to the aforementioned linear models, are certain canonical link generalized linear models and nonparametric regressions constructed from them by a particular style of local likelihood fitting. Other models can often be made robust by using a suboptimal fitting method. Thus the class of model-based, but design consistent, analyses is enlarged to include more realistic models for certain types of survey variable such as binary indicators and counts. Particular applications discussed are the estimation of the size of a population subdomain, as arises in tax auditing for example, and the estimation of a bootstrap tail probability.     

Designing fractional two-level experiments with nested error structures

A common feature of experiments with a random blocking factor and splitplot experiments is their nested error structure. This paper proposes a general strategy to handle fractional two-level experiments with such error structures. The strategy aims to create error strata with sufficient numbers of contrasts to separate active effects from inactive effects. The strategy also details the construction of treatment generators, given the constraints of a predetermined error structure. The key elements of the strategy are illustrated with a chemical experiment that has 16 factors and 32 runs blocked according to working days, and a cheese-making experiment that has 11 factors and 128 runs, divided over milk supplies as whole plots, curds productions as subplots and sets of identically treated cheeses as sub-subplots.