Estimating the mean life time using right censored data

Estimation of the population mean based on right censored observations is considered. The naive sample mean will be an inconsistent and asymptotically biased estimator in this case. An estimate suggested in textbooks is to compute the area under a Kaplan–Meier curve. In this note, two more seemingly different approaches are introduced. Students’ reaction to these approaches was very positive in an introductory survival analysis course the author recently taught.     

Likelihood inference for exchangeable continuous data with covariates and varying cluster sizes; use of the Farlie–Gumbel–Morgenstern model

This article investigates the Farlie–Gumbel–Morgenstern class of models for exchangeable continuous data. We show how the model specification can account for both individual and cluster level covariates, we derive insights from comparisons with the multivariate normal distribution, and we discuss maximum likelihood inference when a sample of independent clusters of varying sizes is available. We propose a method for maximum likelihood estimation which is an alternative to direct numerical maximization of the likelihood that sometimes exhibits non-convergence problems. We describe an algorithm for generating samples from the exchangeable multivariate Farlie–Gumbel–Morgenstern distribution with any marginals, using the structural properties of the distribution. Finally, we present the results of a simulation study designed to assess the properties of the maximum likelihood estimators, and we illustrate the use of the FGM distributions with the analysis of a small data set from a developmental toxicity study.     

Small area estimation under a multivariate linear model for repeated measures data

In this article, small area estimation under a multivariate linear model for repeated measures data is considered. The proposed model aims to get a model which borrows strength both across small areas and over time. The model accounts for repeated surveys, grouped response units, and random effects variations. Estimation of model parameters is discussed within a likelihood based approach. Prediction of random effects, small area means across time points, and per group units are derived. A parametric bootstrap method is proposed for estimating the mean squared error of the predicted small area means. Results are supported by a simulation study.     

Dynamic survival models with spatial frailty

In many survival studies, covariates effects are time-varying and there is presence of spatial effects. Dynamic models can be used to cope with the variations of the effects and spatial components are introduced to handle spatial variation. This paper proposes a methodology to simultaneously introduce these components into the model. A number of specifications for the spatial components are considered. Estimation is performed via a Bayesian approach through Markov chain Monte Carlo methods. Models are compared to assess relevance of their components. Analysis of a real data set is performed, showing the relevance of both time-varying covariate effects and spatial components. Extensions to the methodology are proposed along with concluding remarks.     

A Class of Parametric Dynamic Survival Models

A class of parametric dynamic survival models are explored in which only limited parametric assumptions are made, whilst avoiding the assumption of proportional hazards. Both the log-baseline hazard and covariate effects are modelled by piecewise constant and correlated processes. The method of estimation is to use Markov chain Monte Carlo simulations Gibbs sampling with a Metropolis–Hastings step. In addition to standard right censored data sets, extensions to accommodate interval censoring and random effects are included. The model is applied to two well known and illustrative data sets, and the dynamic variability of covariate effects investigated.     

A simple variance formula for population size estimators by conditioning

This note considers the variance estimation for population size estimators based on capture–recapture experiments. Whereas a diversity of estimators of the population size has been suggested, the question of estimating the associated variances is less frequently addressed. This note points out that the technique of conditioning can be applied here successfully which also allows us to identify sources of variation: the variance due to estimation of the model parameters and the binomial variance due to sampling n units from a population of size N. It is applied to estimators typically used in capture–recapture experiments in continuous time including the estimators of Zelterman and Chao and improves upon previously used variance estimators. In addition, knowledge of the variances associated with the estimators by Zelterman and Chao allows the suggestion of a new estimator as the weighted sum of the two. The decomposition of the variance into the two sources allows also a new understanding of how resampling techniques like the Bootstrap could be used appropriately. Finally, the sample size question for capture–recapture experiments is addressed. Since the variance of population size estimators increases with the sample size, it is suggested to use relative measures such as the observed-to-hidden ratio or the completeness of identification proportion for approaching the question of sample size choice.     

Dynamic generalized extreme value modeling via particle filters

This article addresses how particle filters compare to MCMC methods for posterior density approximations of a model that allows for a dynamic state with fixed parameters and where the observation equation is nonlinear. This is a problem that was not been well studied in the specialized literature. We prove that these state and parameter estimations can be achieved via particle filter methods without the need of more expensive Forward Filtering Backward Sampling (FFBS) simulation. Estimation of a time-varying extreme value model via the generalized extreme value distribution is considered using these particle filter methods and compared to a MCMC algorithm that involves a variety of Metropolis-Hastings steps. We illustrate and compare the different methodologies with simulated data and some minimum daily stock returns occurring monthly from January 4, 1990 to December 28, 2007 using the Tokyo Stock Price Index (TOPIX).     

Robust linear mixed models for Small Area Estimation

Hierarchical models are popular in many applied statistics fields including Small Area Estimation. One well known model employed in this particular field is the Fay–Herriot model, in which unobservable parameters are assumed to be Gaussian. In Hierarchical models assumptions about unobservable quantities are difficult to check. For a special case of the Fay–Herriot model, Sinharay and Stern [2003. Posterior predictive model checking in Hierarchical models. J. Statist. Plann. Inference 111, 209–221] showed that violations of the assumptions about the random effects are difficult to detect using posterior predictive checks. In this present paper we consider two extensions of the Fay–Herriot model in which the random effects are assumed to be distributed according to either an exponential power (EP) distribution or a skewed EP distribution. We aim to explore the robustness of the Fay–Herriot model for the estimation of individual area means as well as the empirical distribution function of their ‘ensemble’. Our findings, which are based on a simulation experiment, are largely consistent with those of Sinharay and Stern as far as the efficient estimation of individual small area parameters is concerned. However, when estimating the empirical distribution function of the ‘ensemble’ of small area parameters, results are more sensitive to the failure of distributional assumptions.     

Nonparametric regression using linear combinations of basis functions

This paper discusses a Bayesian approach to nonparametric regression initially proposed by Smith and Kohn (1996. Journal of Econometrics 75: 317–344). In this approach the regression function is represented as a linear combination of basis terms. The basis terms can be univariate or multivariate functions and can include polynomials, natural splines and radial basis functions. A Bayesian hierarchical model is used such that the coefficient of each basis term can be zero with positive prior probability. The presence of basis terms in the model is determined by latent indicator variables. The posterior mean is estimated by Markov chain Monte Carlo simulation because it is computationally intractable to compute the posterior mean analytically unless a small number of basis terms is used. The present article updates the work of Smith and Kohn (1996. Journal of Econometrics 75: 317–344) to take account of work by us and others over the last three years. A careful discussion is given to all aspects of the model specification, function estimation and the use of sampling schemes. In particular, new sampling schemes are introduced to carry out the variable selection methodology.     

Model weights for model choice and averaging

A method is suggested to estimate posterior model probabilities and model averaged parameters via MCMC sampling under a Bayesian approach. The estimates use pooled output for J models (J>1) whereby all models are updated at each iteration. Posterior probabilities are based on averages of continuous weights obtained for each model at each iteration, while samples of averaged parameters are obtained from iteration specific averages that are based on these weights. Parallel sampling of models assists in deriving posterior densities for parameter contrasts between models and in assessing hypotheses regarding model averaged parameters. Four worked examples illustrate application of the approach, two involving fixed effect regression, and two involving random effects.     

Estimation in a random effects growth curve model

This paper considers estimation under the growth curve model of Potthoff and Roy (1964) with random effects. Estimation under a multivariate model is also considered. Estimation under incomplete data and estimation of random effects are also discussed. A numerical example of data on bulls is presented to illustrate these techniques.     

Multiresponse surface models with random block effects

A multiresponse experiment is one in which a number of responses are measured for each setting of a group of input (control) variables. Quite often, the experimental units are subdivided into groups — or blocks — in order to control an extraneous source of variation. This necessitates adding block effects to the hypothesized multiresponse surface model, which typically contains fixed polynomial effects. Khuri and Valeroso (1998, J. Statist. Plann. Inference, Vol. 73, pp. 7–20) discussed the analysis of such a model when the block effects are considered fixed. There are many situations, however, where these effects should more appropriately be treated as random. In the present article, we address the analysis of a multiresponse model in the latter situations. In particular, we discuss the estimation of the model's polynomial effects in two cases:

• 1.The block effects are additive in the model.
• 2.The blocks have interactive effects with the polynomial portion of the model.

Multiresponse optimization in the presence of blocks will also be discussed. A numerical example is presented for illustrative purposes.     

Conditioning of Marked Point Processes within a Bayesian Framework

Shale units with low permeability create barriers to fluid flow in a sandstone reservoir. A spatial stochastic model for the location of shale units in a reservoir is defined. The model is based on a marked point process formulation, where the marks are parameterized by random functions for the shape of a shale unit. This extends the traditional formulation in the sense that conditioning on the actual observations of the shale units is allowed in an arbitrary number of wells penetrating the reservoir. The marked point process for the shale units includes spatial interaction of units and allows a random number of units to be present. The model is defined in a Bayesian setting with prior pdfs assigned to size–shape parameters of shale units. The observations of shales in wells are associated with a likelihood function. The posterior pdf of the marked point process can only partially be developed analytically; the final solution must be determined by sampling using the Metropolis–Hastings algorithm. An example is presented, demonstrating the consequences of increasing the number of wells in which observations are made.     

Analysis of Linear Models with One Factor Having Both Fixed and Random Levels

ABSTRACT

A general theory for a case where a factor has both fixed and random effect levels is developed under one-way treatment structure model. Estimation procedures for the fixed effects and variance components are consider for the model. The testing of fixed effects is considered when the variance–covariance matrix is known and unknown. Confidence intervals for estimable functions and prediction intervals for predictable functions are constructed. The computational procedures are illustrated using data from an on-farm trial.     

Optimal Composite Markers for Time‐Dependent Receiver Operating Characteristic Curves with Censored Survival Data

Abstract. To increase the predictive abilities of several plasma biomarkers on the coronary artery disease (CAD)‐related vital statuses over time, our research interest mainly focuses on seeking combinations of these biomarkers with the highest time‐dependent receiver operating characteristic curves. An extended generalized linear model (EGLM) with time‐varying coefficients and an unknown bivariate link function is used to characterize the conditional distribution of time to CAD‐related death. Based on censored survival data, two non‐parametric procedures are proposed to estimate the optimal composite markers, linear predictors in the EGLM model. Estimation methods for the classification accuracies of the optimal composite markers are also proposed. In the article we establish theoretical results of the estimators and examine the corresponding finite‐sample properties through a series of simulations with different sample sizes, censoring rates and censoring mechanisms. Our optimization procedures and estimators are further shown to be useful through an application to a prospective cohort study of patients undergoing angiography.     

Analysis of Poisson varying-coefficient models with autoregression

In the regression analysis of time series of event counts, it is of interest to account for serial dependence that is likely to be present among such data as well as a nonlinear interaction between the expected event counts and predictors as a function of some underlying variables. We thus develop a Poisson autoregressive varying-coefficient model, which introduces autocorrelation through a latent process and allows regression coefficients to nonparametrically vary as a function of the underlying variables. The nonparametric functions for varying regression coefficients are estimated with data-driven basis selection, thereby avoiding overfitting and adapting to curvature variation. An efficient posterior sampling scheme is devised to analyse the proposed model. The proposed methodology is illustrated using simulated data and daily homicide data in Cali, Colombia.     

Generalized linear models and transformations for unreplicated factorial experiments in the presence of dispersion effects

For the analysis of replicated designs, many different methods have been suggested. These allow for the estimation of functional dependencies between mean and variance as well as possible dispersion effects within the same model framework. However, in the situation of unreplicated designs, most methods known so far rely on the assumption of constant variances, or a functional relationship between mean and variance as the only source of heteroscedasticity. In this paper, we propose two methods for dealing with unreplicated data, when dispersion effects might also be of importance. One of these is an extension of the Box–Cox-method [Box, G.E.P., Cox, D.R., 1964. An analysis of transformations. Journal of the Royal Statistical Society B 26, 211–252], the other is based on double generalized linear models. Both these methods turn out to yield approximately equivalent results in the case of comparable assumptions, whereas the double generalized linear model is the more general one and allows further extensions. If this class of models is assumed, consistency, asymptotic efficiency and normality of the resulting estimates are shown.     

QUALITY IMPROVEMENT FROM DISASSEMBLY-REASSEMBLY EXPERIMENTS

The goal of achieving high quality products has led to an emphasis on reducing variation in performance characteristics. It may often happen that one of the product's components is responsible for much of the observed variation. This research is stimulated by the problem of detecting a component that impairs quality by systematically inflating the variance in a product that is assembled from “interchangeable components.” We consider the class of “disassembly-reassembly” experiments, in which components are swapped among assemblies. The specific units used in the experiment are sampled from a large population of units, so it is natural to measure the influence of each factor by its variance component. We present the model for these experiments as a special case of the mixed linear model, compare several estimators for the variance components and consider the problem of testing hypotheses to identify troublesome components.     

Interblock information in multidimensional block designs

In many experimental situations, d-way heterogeneity among experimental units may be controlled through use of multiple blocking criteria. In some cases it is reasonable to regard some or all of the block effects as random. Then the model is mixed and observations within blocks are correlated. Very general estimators of treatment effects and their dispersion matrix with recovery of interblock information are provided. They apply to designs with d > 1 blocking criteria that may be crossed, nested, or a combination thereof. These general results may be specialized to provide analyses of new classes of MBD's or used directly for numerical analyses of designs in the general class, perhaps through use as the basis for very general computer programs. Estimation of variance components is discussed, and an example is provided to illustrate adaptation of the general results.     

On the estimation of homogeneous population size from a complex dual-record system

ABSTRACT

A dual-record system (DRS) (equivalently two sample capture–recapture experiments) model, with time and behavioural response variation, has attracted much attention specifically in the domain of official statistics and epidemiology, as the assumption of list independence often fails. The relevant model suffers from parameter identifiability problem, and suitable Bayesian methodologies could be helpful. In this article, we formulate population size estimation in DRS as a missing data problem and two empirical Bayes approaches are proposed along with the discussion of an existing Bayes treatment. Some features and associated posterior convergence for these methods are mentioned. Investigation through an extensive simulation study finds that our proposed approaches compare favourably with the existing Bayes approach for this complex model depending upon the availability of directional nature of underlying behavioural response effect. A real-data example is given to illustrate these methods.