Testing for the partial Koziol–Green model with covariates

In this paper we consider some non-parametric goodness-of-fit statistics for testing the partial Koziol–Green regression model. In this model, the response at a given covariate value is subject to random right censoring by two independent censoring times. One of these censoring times is informative in the sense that its survival function is some power of the survival function of the response. The goodness-of-fit statistics are based on an underlying empirical process for which large sample theory is obtained.     

Optimal designs for the Michaelis–Menten model with correlated observations

In this paper we investigate the problem of designing experiments for generalized least-squares analysis in the Michaelis–Menten model. We study the structure of exact D-optimal designs in a model with an autoregressive error structure. Explicit results for locally D-optimal designs are derived for the case where two observations can be taken per subject. Additionally standardized maximin D-optimal designs are obtained in this case. The results illustrate the enormous difficulties to find exact optimal designs explicitly for nonlinear regression models with correlated observations.     

Mortality,longevity and experiments with the Lee–Carter model

The paper reviews the Lee-Carter modelling framework, illustrated with an application, and then extends the framework through the development of a wider class of generalised, parametric, non-linear models. The choice of error distribution is also generalised. These extensions permit the modelling and extrapolation of age-specific cohort effects as well as the more familiar age-specific period effects: the age-period-cohort version of the model is discussed with a worked example. The paper also provides a comparative study of simulation strategies for assessing risk in mortality rate predictions and the associated forecast estimates of life expectancy and annuity values in both period and cohort perspectives.     

r − k Class estimator in the linear regression model with correlated errors

Autocorrelation in errors and multicollinearity among the regressors are serious problems in regression analysis. The aim of this paper is to examine multicollinearity and autocorrelation problems concurrently and to compare the r ? k class estimator to the generalized least squares estimator, the principal components regression estimator and the ridge regression estimator by the scalar and matrix mean square error criteria in the linear regression model with correlated errors.     

ß-expectation tolerance region for the heteroscedastic multiple regression model with multivariate Student-t error

A ß-expectation tolerance region has been constructed for the multivariate regression model with heteroscedastic errors which follow a multivariate Student-t distribution with an unknown number of degrees of freedom. The ß-expectaion tolerance region obtained in this paper is optimal in the sense of having minimum enclosure among all such tolerance regions that guarantees that it would cover any preassigned proportions, namely, ß×100 percent of the future responses from the model.     

Bayesian–frequentist hybrid model with application to the analysis of gene copy number changes

Gene copy number (GCN) changes are common characteristics of many genetic diseases. Comparative genomic hybridization (CGH) is a new technology widely used today to screen the GCN changes in mutant cells with high resolution genome-wide. Statistical methods for analyzing such CGH data have been evolving. Existing methods are either frequentist's or full Bayesian. The former often has computational advantage, while the latter can incorporate prior information into the model, but could be misleading when one does not have sound prior information. In an attempt to take full advantages of both approaches, we develop a Bayesian-frequentist hybrid approach, in which a subset of the model parameters is inferred by the Bayesian method, while the rest parameters by the frequentist's. This new hybrid approach provides advantages over those of the Bayesian or frequentist's method used alone. This is especially the case when sound prior information is available on part of the parameters, and the sample size is relatively small. Spatial dependence and false discovery rate are also discussed, and the parameter estimation is efficient. As an illustration, we used the proposed hybrid approach to analyze a real CGH data.     

The partial Koziol–Green model with covariates

Two important models in survival analysis are that of general random censorship and the proportional hazards submodel of Koziol and Green. The difference between the two models is the way in which the lifetime variable is censored (informative versus non-informative censoring). In this paper the two viewpoints are combined into a new model which allows the lifetimes to be censored by two types of variables, one of which censors in an informative way and the other one in a non-informative way. The lifetimes and the censoring times are also allowed to depend on covariates in a very general way. The estimator for the conditional distribution of the lifetimes generalizes that of Gather and Pawlitschko (1998. Metrika 48, 189–209), who recently studied the situation without covariate information. Results obtained are the uniform strong consistency (with rate), an almost sure asymptotic representation and the weak convergence of the process.     

Optimal trajectory to the general equilibrium state of the unstable stochastic lattice gas of prey–predator model with pair-approximation

This paper focuses on the study of the optimal control problem of the unstable stochastic lattice gas of prey–predator model with pair-approximation. The Pontryagin minimum procedure is used to derive the optimal behaviour of the proposed system that makes it trajectory closer to its general equilibrium states with respect to a selected performance measure during a specified period of time T. Wide numerical simulations and illustrations are presented in graphs and tables.     

Optimum mixture designs for the log-logistic dose–response model with mixture of two similar compounds

The article studies the log-logistic class of dose–response bioassay models in the binomial set-up. The dose is identified by the potency adjusted mixing proportions of two similar compounds. Models for both absence and presence of interaction between the compounds have been considered. The aim is to investigate the D- and D_s-optimal mixture designs for the estimation of the full set of parameters or for the estimation of potency for a best guess of the parameter values. We also indicate how to find the optimal design to estimate the mixing proportions at which the probability of success attains a given value in the absence of the interaction effect.     

The Poisson Birnbaum–Saunders model with long-term survivors

In this paper, we propose a cure rate survival model by assuming that the number of competing causes of the event of interest follows the Poisson distribution and the time to event has the Birnbaum–Saunders (BS) distribution. We define the Poisson BS distribution and provide two useful representations for its density function which facilitate to obtain some mathematical properties. Two closed-form expressions for the moments of the new distribution are given. We estimate the parameters of the model with cure rate using maximum likelihood. For different parameter settings, sample sizes and censoring percentages, several simulations are performed. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and present some ways to perform a global influence study. We analyse a real data set from the medical area.     

Optimality of the quasi-score estimator in a mean–variance model with applications to measurement error models

We consider a regression of y$y$ on x$x$ given by a pair of mean and variance functions with a parameter vector θ$\theta$ to be estimated that also appears in the distribution of the regressor variable x$x$. The estimation of θ$\theta$ is based on an extended quasi-score (QS) function. We show that the QS estimator is optimal within a wide class of estimators based on linear-in-y$y$ unbiased estimating functions. Of special interest is the case where the distribution of x$x$ depends only on a subvector α$\alpha$ of θ$\theta$, which may be considered a nuisance parameter. In general, α$\alpha$ must be estimated simultaneously together with the rest of θ$\theta$, but there are cases where α$\alpha$ can be pre-estimated. A major application of this model is the classical measurement error model, where the corrected score (CS) estimator is an alternative to the QS estimator. We derive conditions under which the QS estimator is strictly more efficient than the CS estimator.     

On the dependent competing risks using Marshall–Olkin bivariate Weibull model: Parameter estimation with different methods

Competing risks models are of great importance in reliability and survival analysis. They are often assumed to have independent causes of failure in literature, which may be unreasonable. In this article, dependent causes of failure are considered by using the Marshall–Olkin bivariate Weibull distribution. After deriving some useful results for the model, we use ML, fiducial inference, and Bayesian methods to estimate the unknown model parameters with a parameter transformation. Simulation studies are carried out to assess the performances of the three methods. Compared with the maximum likelihood method, the fiducial and Bayesian methods could provide better parameter estimation.     

A two-stage estimation in the Clayton–Oakes model with marginal linear transformation models for multivariate failure time data

This paper considers the analysis of multivariate survival data where the marginal distributions are specified by semiparametric transformation models, a general class including the Cox model and the proportional odds model as special cases. First, consideration is given to the situation where the joint distribution of all failure times within the same cluster is specified by the Clayton–Oakes model (Clayton, Biometrika 65:141–151, l978; Oakes, J R Stat Soc B 44:412–422, 1982). A two-stage estimation procedure is adopted by first estimating the marginal parameters under the independence working assumption, and then the association parameter is estimated from the maximization of the full likelihood function with the estimators of the marginal parameters plugged in. The asymptotic properties of all estimators in the semiparametric model are derived. For the second situation, the third and higher order dependency structures are left unspecified, and interest focuses on the pairwise correlation between any two failure times. Thus, the pairwise association estimate can be obtained in the second stage by maximizing the pairwise likelihood function. Large sample properties for the pairwise association are also derived. Simulation studies show that the proposed approach is appropriate for practical use. To illustrate, a subset of the data from the Diabetic Retinopathy Study is used.     

A new cure rate model based on the Yule–Simon distribution with application to a melanoma data set

In this paper, a new survival cure rate model is introduced considering the Yule–Simon distribution [12 H.A. Simon, On a class of skew distribution functions, Biometrika 42 (1955), pp. 425–440.[Crossref], [Web of Science ®] , [Google Scholar]] to model the number of concurrent causes. We study some properties of this distribution and the model arising when the distribution of the competing causes is the Weibull model. We call this distribution the Weibull–Yule–Simon distribution. Maximum likelihood estimation is conducted for model parameters. A small scale simulation study is conducted indicating satisfactory parameter recovery by the estimation approach. Results are applied to a real data set (melanoma) illustrating the fact that the model proposed can outperform traditional alternative models in terms of model fitting.     

Influence diagnostics for Grubbs’s model with asymmetric heavy-tailed distributions

Grubbs’s model (Grubbs, Encycl Stat Sci 3:42–549, 1983) is used for comparing several measuring devices, and it is common to assume that the random terms have a normal (or symmetric) distribution. In this paper, we discuss the extension of this model to the class of scale mixtures of skew-normal distributions. Our results provide a useful generalization of the symmetric Grubbs’s model (Osorio et al., Comput Stat Data Anal, 53:1249–1263, 2009) and the asymmetric skew-normal model (Montenegro et al., Stat Pap 51:701–715, 2010). We discuss the EM algorithm for parameter estimation and the local influence method (Cook, J Royal Stat Soc Ser B, 48:133–169, 1986) for assessing the robustness of these parameter estimates under some usual perturbation schemes. The results and methods developed in this paper are illustrated with a numerical example.     

Approximate estimations in capture–recapture model with heteroscedastic measurement errors

In a capture–recapture experiment, the number of measurements for individual covariates usually equals the number of captures. This creates a heteroscedastic measurement error problem and the usual surrogate condition does not hold in the context of a measurement error model. This study adopts a small measurement error assumption to approximate the conventional estimating functions and the population size estimator. This study also investigates the biases of the resulting estimators. In addition, modifications for two common approximation methods, regression calibration and simulation extrapolation, to accommodate heteroscedastic measurement error are also discussed. These estimation methods are examined through simulations and illustrated by analysing a capture–recapture data set.