Alternative models for the heterogeneity of mortality risks among the aged

The authors examine how sensitive the estimates of heterogeneity in the mortality risks in a population are to the choices of two types of function, "one describing the age-specific rate of increase of mortality risks for individuals and the other describing the distribution of mortality risks across individuals." U.S. data from published Medicare mortality rates for the period 1968-1978 are used to analyze total mortality among the aged. "In addition, national vital statistics data for the period 1950-1977 were used to analyze adult lung cancer mortality. For these data, the estimates of structural parameters were less sensitive to reasonable choices of the heterogeneity distribution (gamma vs. inverse Gaussian) than to reasonable choices of the hazard rate function (Gompertz vs. Weibull)."     

Bootstrap方法在死亡模型中的应用

 由于不同国家死亡率改善现象不同,世界各国所使用的死亡率模型皆不尽相同,而且不同年龄段的死亡率模型也不同。实际中,我们常常采用Gompertz模型、Makeham模型、Weibull模型等拟合高年龄段人口的死亡率,但是因高年龄段人口的死亡数据资料不够充分,较少有人以统计的观点给出模型适合性的检验过程。因此本研究提出利用Bootstrap方法检验死亡模型假设的方法,包括模型适合性的检验、参数估计、参数假设检验等。最后,本文应用中国1997－2007年65－89岁人口的粗死亡率数据,提出适合的死亡模型,然后给出利用Bootstrap方法进行死亡模型检验的全过程。     

Shape Constrained Non‐parametric Estimators of the Baseline Distribution in Cox Proportional Hazards Model

Abstract. We investigate non‐parametric estimation of a monotone baseline hazard and a decreasing baseline density within the Cox model. Two estimators of a non‐decreasing baseline hazard function are proposed. We derive the non‐parametric maximum likelihood estimator and consider a Grenander type estimator, defined as the left‐hand slope of the greatest convex minorant of the Breslow estimator. We demonstrate that the two estimators are strongly consistent and asymptotically equivalent and derive their common limit distribution at a fixed point. Both estimators of a non‐increasing baseline hazard and their asymptotic properties are obtained in a similar manner. Furthermore, we introduce a Grenander type estimator for a non‐increasing baseline density, defined as the left‐hand slope of the least concave majorant of an estimator of the baseline cumulative distribution function, derived from the Breslow estimator. We show that this estimator is strongly consistent and derive its asymptotic distribution at a fixed point.     

On finite mixture of two-component gompertz lifetime model

The Gompertz distribution has been used as a growth model, especially in epidemiological and biomedical studies. Based on Type I and II censored samples from a heterogeneous population that can be represented by a finite mixture of two-component Gompertz lifetime model, the maximum likelihood and Bayes estimates of the parameters, reliability and hazard rate functions are obtained. An approximation form due to Lindley (1980) is used in obtaining the corresponding Bayes estimates. The maximum likelihood and Bayes estimates are comparedvia a Monte Carlo simulation study.     

A Quasi‐Score Statistic for Homogeneity Testing against Covariate‐Varying Heterogeneity

In statistical modelling, it is often of interest to evaluate non‐negative quantities that capture heterogeneity in the population such as variances, mixing proportions and dispersion parameters. In instances of covariate‐dependent heterogeneity, the implied homogeneity hypotheses are nonstandard and existing inferential techniques are not applicable. In this paper, we develop a quasi‐score test statistic to evaluate homogeneity against heterogeneity that varies with a covariate profile through a regression model. We establish the limiting null distribution of the proposed test as a functional of mixtures of chi‐square processes. The methodology does not require the full distribution of the data to be entirely specified. Instead, a general estimating function for a finite dimensional component of the model, that is, of interest is assumed but other characteristics of the population are left completely unspecified. We apply the methodology to evaluate the excess zero proportion in zero‐inflated models for count data. Our numerical simulations show that the proposed test can greatly improve efficiency over tests of homogeneity that neglect covariate information under the alternative hypothesis. An empirical application to dental caries indices demonstrates the importance and practical utility of the methodology in detecting excess zeros in the data.     

On damage accumulation and biological aging

If biological aging is understood as some process of damage accumulation, it does not necessarily lead to increasing mortality rate. Within the framework of suggested models and relevant examples we show that even for the monotonically increasing degradation, the mortality rate can, at least, ultimately decrease. Aging properties of systems with imperfect repair are also studied. It is shown that for some models of imperfect repair the corresponding age process is monotone and stable. This means that as t→∞$t\to \infty$, degradation slows down, which results in the mortality rate deceleration and its possible convergence to a constant.     

DYNAMICAL AND STATISTICAL MODELS FOR EXPLOITED POPULATIONS†

Human activities have indirectly modified the dynamics of many populations, accelerating considerably the natural rate of species extinction and raising strong concerns about biodiversity. In many such cases, the underlying ‘natural’ dynamics of the population has been modified by human‐induced increases in mortality, even if the populations are not exploited or harvested in the strict sense. Both dynamical and statistical models are needed to investigate the consequences of human‐induced mortality on the overall dynamics of a population. This paper reviews existing approaches and the potential of recent developments to help form a conceptual and practical framework for analysing the dynamics of exploited populations. It examines both the simple case of an extra source of mortality instantaneously in time, and the theory involved when both risks compete over a continuous time scale. This basic theory is expanded to structured populations, using matrix population models, with applications to the conservation biology of long‐lived vertebrates. The type and degree of compensation expected and approaches to detect it are reviewed, and ways of handling uncertainty are discussed.     

A Bayesian approach to mixture cure models with spatial frailties for population‐based cancer relative survival data

As the treatments of cancer progress, a certain number of cancers are curable if diagnosed early. In population‐based cancer survival studies, cure is said to occur when mortality rate of the cancer patients returns to the same level as that expected for the general cancer‐free population. The estimates of cure fraction are of interest to both cancer patients and health policy makers. Mixture cure models have been widely used because the model is easy to interpret by separating the patients into two distinct groups. Usually parametric models are assumed for the latent distribution for the uncured patients. The estimation of cure fraction from the mixture cure model may be sensitive to misspecification of latent distribution. We propose a Bayesian approach to mixture cure model for population‐based cancer survival data, which can be extended to county‐level cancer survival data. Instead of modeling the latent distribution by a fixed parametric distribution, we use a finite mixture of the union of the lognormal, loglogistic, and Weibull distributions. The parameters are estimated using the Markov chain Monte Carlo method. Simulation study shows that the Bayesian method using a finite mixture latent distribution provides robust inference of parameter estimates. The proposed Bayesian method is applied to relative survival data for colon cancer patients from the Surveillance, Epidemiology, and End Results (SEER) Program to estimate the cure fractions. The Canadian Journal of Statistics 40: 40–54; 2012 © 2012 Statistical Society of Canada     

Tests for Coefficients in High‐dimensional Additive Hazard Models

We consider hypothesis testing problems for low‐dimensional coefficients in a high dimensional additive hazard model. A variance reduced partial profiling estimator (VRPPE) is proposed and its asymptotic normality is established, which enables us to test the significance of each single coefficient when the data dimension is much larger than the sample size. Based on the p‐values obtained from the proposed test statistics, we then apply a multiple testing procedure to identify significant coefficients and show that the false discovery rate can be controlled at the desired level. The proposed method is also extended to testing a low‐dimensional sub‐vector of coefficients. The finite sample performance of the proposed testing procedure is evaluated by simulation studies. We also apply it to two real data sets, with one focusing on testing low‐dimensional coefficients and the other focusing on identifying significant coefficients through the proposed multiple testing procedure.     

Determination of hazard ratio for progression‐free survival considering the tumor assessment schedule in sample size calculation

Progression‐free survival is recognized as an important endpoint in oncology clinical trials. In clinical trials aimed at new drug development, the target population often comprises patients that are refractory to standard therapy with a tumor that shows rapid progression. This situation would increase the bias of the hazard ratio calculated for progression‐free survival, resulting in decreased power for such patients. Therefore, new measures are needed to prevent decreasing the power in advance when estimating the sample size. Here, I propose a novel calculation procedure to assume the hazard ratio for progression‐free survival using the Cox proportional hazards model, which can be applied in sample size calculation. The hazard ratios derived by the proposed procedure were almost identical to those obtained by simulation. The hazard ratio calculated by the proposed procedure is applicable to sample size calculation and coincides with the nominal power. Methods that compensate for the lack of power due to biases in the hazard ratio are also discussed from a practical point of view.     

Model‐Based Non‐parametric Variance Estimation for Systematic Sampling

Abstract. Systematic sampling is frequently used in surveys, because of its ease of implementation and its design efficiency. An important drawback of systematic sampling, however, is that no direct estimator of the design variance is available. We describe a new estimator of the model‐based expectation of the design variance, under a non‐parametric model for the population. The non‐parametric model is sufficiently flexible that it can be expected to hold at least approximately in many situations with continuous auxiliary variables observed at the population level. We prove the model consistency of the estimator for both the anticipated variance and the design variance under a non‐parametric model with a univariate covariate. The broad applicability of the approach is demonstrated on a dataset from a forestry survey.     

Modelling survival data using mixtures of frailties

Frailty models are often used to describe the extra heterogeneity in survival data by introducing an individual random, unobserved effect. The frailty term is usually assumed to act multiplicatively on a baseline hazard function common to all individuals. In order to apply the frailty model, a specific frailty distribution has to be assumed. If at least one of the latent variables is continuous, the frailty must follow a continuous distribution. In this paper, a finite mixture of continuous frailty distributions is used in order to describe situations in which one (or more) of the latent variables separates the population in study into two (or more) subpopulations. Closure properties of the unobserved quantity are given along with the maximum-likelihood estimates under the most common choices of frailty distributions. The model is illustrated on a set of lifetime data.     

The McDonald Gompertz distribution: Properties and applications

This article introduces a five-parameter lifetime model called the McDonald Gompertz (McG) distribution to extend the Gompertz, generalized Gompertz, generalized exponential, beta Gompertz, and Kumaraswamy Gompertz distributions among several other models. The hazard function of new distribution can be increasing, decreasing, upside-down bathtub, and bathtub shaped. We obtain several properties of the McG distribution including moments, entropies, quantile, and generating functions. We provide the density function of the order statistics and their moments. The parameter estimation is based on the usual maximum likelihood approach. We also provide the observed information matrix and discuss inferences issues. The flexibility and usefulness of the new distribution are illustrated by means of application to two real datasets.     

Exploring Parameter Relations for Multi‐Stage Models in Stage‐Wise Constant and Time Dependent Hazard Rates

Single cohort stage‐frequency data are considered when assessing the stage reached by individuals through destructive sampling. For this type of data, when all hazard rates are assumed constant and equal, Laplace transform methods have been applied in the past to estimate the parameters in each stage‐duration distribution and the overall hazard rates. If hazard rates are not all equal, estimating stage‐duration parameters using Laplace transform methods becomes complex. In this paper, two new models are proposed to estimate stage‐dependent maturation parameters using Laplace transform methods where non‐trivial hazard rates apply. The first model encompasses hazard rates that are constant within each stage but vary between stages. The second model encompasses time‐dependent hazard rates within stages. Moreover, this paper introduces a method for estimating the hazard rate in each stage for the stage‐wise constant hazard rates model. This work presents methods that could be used in specific types of laboratory studies, but the main motivation is to explore the relationships between stage maturation parameters that, in future work, could be exploited in applying Bayesian approaches. The application of the methodology in each model is evaluated using simulated data in order to illustrate the structure of these models.     

A Proportional Hazard Marshall–Olkin Extended Family of Distributions and its Application to Gompertz Distribution

In this paper, we have presented a proportional hazard version of the Marshall–Olkin extended family of distributions. This family of distributions has been compared in terms of stochastic orderings with the Marshall-Olkin extended family of distributions. Considering the Gompertz distribution as the baseline, the monotonicity of the resulting failure rate is shown to be either increasing or bathtub, even though the Gompertz distribution has an increasing failure rate. The maximum likelihood estimation of the parameters has been studied and a data set, involving the serum–reversal times, has been analyzed and it has been shown that the model presented in this paper fit better than the Gompertz or even the Mrashall–Olkin Gompertz distribution. The extension presented in this paper can be used in other family of distributions as well.     

The Effects of Variations in the Components of Change on the Size and Structure of The Australian Population: 1911-66

The main aim of the study was to estimate separately the effects of the variations in the components of population change and the effects of their mutual interactions on the size and age structure and other characteristics of the Australian population during 1911-66. The method proposed, called here the factorial projection method, is to project the population over the period under consideration, under different assumptions according to a 2×2×2 factorial design, utilizing the observed variations in fertility, mortality and migration. It was found that, apart from the effects of variations in the components occurring individually, the two factor interaction effects– the interaction effects of the changes occurring simultaneously in two of the components (viz. fertility and mortality or mortality and migration or fertility and migration)–on the population size, were considerable. The contribution of migration to the working age group was greater than its effect on the total size, and mortality improvement played a more prominent role in increasing the size of the old age population, while the reduction in fertility caused a higher proportion in the old age group and thus increased the mean age of the population both in the case of males and females. The two factor interaction effects were not negligible in changing the size of the school age population, population in the working age group and also in the old age group. But the interaction of migration and the decline in fertility was the only interaction effect which was important in changing the percentage age distribution and the mean age of the male and female population.     

Dynamic mixing probability measures of mixtures

ABSTRACT

Lifetime of heterogeneous population can be modeled as mixture of a family of lifetime distributions according to a mixing probability measure. With the help of dynamic mixing measure, the hazard rate of the mixture can also be expressed as the mixture of the hazard rates of the lifetime distributions. Various local stochastic orderings are defined in this article. Applying these local stochastic orderings, we can explore the behavior of the dynamic mixing measures locally and then compare the hazard rates of two heterogeneous populations in both the local and global ways.     

Why a time effect often has a limited impact on capture‐recapture estimates in closed populations

The author is concerned with log‐linear estimators of the size N of a population in a capture‐recapture experiment featuring heterogeneity in the individual capture probabilities and a time effect. He also considers models where the first capture influences the probability of subsequent captures. He derives several results from a new inequality associated with a dispersive ordering for discrete random variables. He shows that in a log‐linear model with inter‐individual heterogeneity, the estimator N is an increasing function of the heterogeneity parameter. He also shows that the inclusion of a time effect in the capture probabilities decreases N in models without heterogeneity. He further argues that a model featuring heterogeneity can accommodate a time effect through a small change in the heterogeneity parameter. He demonstrates these results using an inequality for the estimators of the heterogeneity parameters and illustrates them in a Monte Carlo experiment     

Least squares estimation of varying‐coefficient hazard regression with application to breast cancer dose‐intensity data

To enhance modeling flexibility, the authors propose a nonparametric hazard regression model, for which the ordinary and weighted least squares estimation and inference procedures are studied. The proposed model does not assume any parametric specifications on the covariate effects, which is suitable for exploring the nonlinear interactions between covariates, time and some exposure variable. The authors propose the local ordinary and weighted least squares estimators for the varying‐coefficient functions and establish the corresponding asymptotic normality properties. Simulation studies are conducted to empirically examine the finite‐sample performance of the new methods, and a real data example from a recent breast cancer study is used as an illustration. The Canadian Journal of Statistics 37: 659–674; 2009 © 2009 Statistical Society of Canada     

Minimax A‐ and D‐optimal integer‐valued wavelet designs for estimation

The author discusses integer‐valued designs for wavelet estimation of nonparametric response curves in the possible presence of heteroscedastic noise using a modified wavelet version of the Gasser‐Müller kernel estimator or weighted least squares estimation. The designs are constructed using a minimax treatment and the simulated annealing algorithm. The author presents designs for three case studies in experiments for investigating Gompertz's theory on mortality rates, nitrite utilization in bush beans and the impact of crash helmets in motorcycle accidents.