Fitting Distributions with the Polyhazard Model with Dependence

The polyhazard model with dependent causes, first introduced to fit lifetime data, generalized the traditional polyhazard model by allowing the latent causes of failure to be dependent by using copula functions. When modeling lifetime data, marginal distributions are supported on the positive reals. Dropping this restriction, the method generates a rich family of univariate distributions with asymmetries and multiple modes. We show that this new family of distributions is able to approximate other distributions proposed in the literature, such as the generalized beta-generated distributions. These distributions are fitted to three real data sets.     

Modelling Survival Events with Longitudinal Covariates Measured with Error

In survival analysis, time-dependent covariates are usually present as longitudinal data collected periodically and measured with error. The longitudinal data can be assumed to follow a linear mixed effect model and Cox regression models may be used for modelling of survival events. The hazard rate of survival times depends on the underlying time-dependent covariate measured with error, which may be described by random effects. Most existing methods proposed for such models assume a parametric distribution assumption on the random effects and specify a normally distributed error term for the linear mixed effect model. These assumptions may not be always valid in practice. In this article, we propose a new likelihood method for Cox regression models with error-contaminated time-dependent covariates. The proposed method does not require any parametric distribution assumption on random effects and random errors. Asymptotic properties for parameter estimators are provided. Simulation results show that under certain situations the proposed methods are more efficient than the existing methods.     

Failure time regression with continuous covariates measured with error

We consider failure time regression analysis with an auxiliary variable in the presence of a validation sample. We extend the nonparametric inference procedure of Zhou and Pepe to handle a continuous auxiliary or proxy covariate. We estimate the induced relative risk function with a kernel smoother and allow the selection probability of the validation set to depend on the observed covariates. We present some asymptotic properties for the kernel estimator and provide some simulation results. The method proposed is illustrated with a data set from an on-going epidemiologic study.     

Outcome-adaptive allocation with natural lead-in for three-group trials with binary outcomes

ABSTRACT

Just as Bayes extensions of the frequentist optimal allocation design have been developed for the two-group case, we provide a Bayes extension of optimal allocation in the three-group case. We use the optimal allocations derived by Jeon and Hu [Optimal adaptive designs for binary response trials with three treatments. Statist Biopharm Res. 2010;2(3):310–318] and estimate success probabilities for each treatment arm using a Bayes estimator. We also introduce a natural lead-in design that allows adaptation to begin as early in the trial as possible. Simulation studies show that the Bayesian adaptive designs simultaneously increase the power and expected number of successfully treated patients compared to the balanced design. And compared to the standard adaptive design, the natural lead-in design introduced in this study produces a higher expected number of successes whilst preserving power.     

Prediction with Additive Models—Simulation and Application with Real Data

Abstract

In this paper we study the predictor behaviour of the additive model. The prediction equation is introduced as well as the computational considerations to select the smoothing parameters through cross-validation. The additive predictor is compared with a partially linear predictor in a broad simulation study and an application to a real case, prediction of the atmospheric concentration of SO₂ in sample stations.     

Detecting clusters with increased mean using scan windows with variable radius

Applying spatiotemporal scan statistics is an effective method to detect the clustering of mean shifts in many application fields. Although several exponentially weighted moving average (EWMA) based scan statistics have been proposed, the existing methods generally require a fixed scan window size or apply the weighting technique across the temporal axis only. However, the size of shift coverage is often unavailable in practical problems. Using a mismatching scan radius may mislead the size of cluster coverage in space or delay the time to detection. This research proposed an stEWMA method by applying the weighting technique across both temporal and spatial axes with variable scan radius. The simulation analysis showed that the stEWMA method can have a significantly shorter time to detection than the likelihood ratio-based scan statistic using variable scan radius, especially when cluster coverage size is small. The application to detecting the increase of male thyroid cancer in the New Mexico state also showed the effectiveness of the proposed method.     

Three-stage group screening with unequal group sizes and with errors in observations

In this paper three-stage group screening in which group-factors contain differing number of factors is discussed. A procedure for grouping the factors in the absence of concrete prior information is described. Formulas for the expected number of runs and the expected number of incorrect decisions have also been obtained. These formulas are used to formulate criteria for optimal designs.     

Urn designs with immigration: Useful connection with continuous time stochastic processes

In this paper we present techniques that allow obtaining joint and marginal distributions of random variables arising in urn designs with immigration. The urn is embedded into a family of continuous time Markov processes. Then the generating function approach is used to obtain distributional results. A randomized version of the inverse sampling is introduced.     

Modeling Event Times with Multiple Outcomes Using the Wiener Process with Drift

Length of stay in hospital (LOS) is a widely used outcome measure in Health Services research, often acting as a surrogate for resource consumption or as a measure of efficiency. The distribution of LOS is typically highly skewed, with a few large observations. An interesting feature is the presence of multiple outcomes (e.g. healthy discharge, death in hospital, transfer to another institution). Health Services researchers are interested in modeling the dependence of LOS on covariates, often using administrative data collected for other purposes, such as calculating fees for doctors. Even after all available covariates have been included in the model, unexplained heterogeneity usually remains. In this article, we develop a parametric regression model for LOS that addresses these features. The model is based on the time, T, that a Wiener process with drift (representing an unobserved health level process) hits one of two barriers, one representing healthy discharge and the other death in hospital. Our approach to analyzing event times has many parallels with competing risks analysis (Kalbfleisch and Prentice, The Statistical Analysis of Failure Time Data, New York: John Wiley and Sons, 1980)), and can be seen as a way of formalizing a competing risks situation. The density of T is an infinite series, and we outline a proof that the density and its derivatives are absolutely and uniformly convergent, and regularity conditions are satisfied. Expressions for the expected value of T, the conditional expectation of T given outcome, and the probability of each outcome are available in terms of model parameters. The proposed regression model uses an approximation to the density formed by truncating the series, and its parameters are estimated by maximum likelihood. An extension to allow a third outcome (e.g. transfers out of hospital) is discussed, as well as a mixture model that addresses the issue of unexplained heterogeneity. The model is illustrated using administrative data.     

Ancestral Inference in Population Genetics Models with Selection (with Discussion)

A new algorithm is presented for exact simulation from the conditional distribution of the genealogical history of a sample, given the composition of the sample, for population genetics models with general diploid selection. The method applies to the usual diffusion approximation of evolution at a single locus, in a randomly mating population of constant size, for mutation models in which the distribution of the type of a mutant does not depend on the type of the progenitor allele; this includes any model with only two alleles. The new method is applied to ancestral inference for the two‐allele case, both with genic selection and heterozygote advantage and disadvantage, where one of the alleles is assumed to have resulted from a unique mutation event. The paper describes how the method could be used for inference when data are also available at neutral markers linked to the locus under selection. It also informally describes and constructs the non‐neutral Fleming–Viot measure‐valued diffusion.     

Errors-in-variables estimation with wavelets

This paper proposes a wavelet (spectral) approach to estimate the parameters of a linear regression model where the regressand and the regressors are persistent processes and contain a measurement error. We propose a wavelet filtering approach which does not require instruments and yields unbiased estimates for the intercept and the slope parameters. Our Monte Carlo results also show that the wavelet approach is particularly effective when measurement errors for the regressand and the regressor are serially correlated. With this paper, we hope to bring a fresh perspective and stimulate further theoretical research in this area.