Testing equality of cause-specific hazard rates corresponding to m competing risks among K groups

In this paper, a class of tests is developed for comparing the cause-specific hazard rates of m competing risks simultaneously in K ( 2) groups. The data available for a unit are the failure time of the unit along with the identifier of the risk claiming the failure. In practice, the failure time data are generally right censored. The tests are based on the difference between the weighted averages of the cause-specific hazard rates corresponding to each risk. No assumption regarding the dependence of the competing risks is made. It is shown that the proposed test statistic has asymptotically chi-squared distribution. The proposed test is shown to be optimal for a specific type of local alternatives. The choice of weight function is also discussed. A simulation study is carried out using multivariate Gumbel distribution to compare the optimal weight function with a proposed weight function which is to be used in practice. Also, the proposed test is applied to real data on the termination of an intrauterine device.An erratum to this article can be found at     

Optimal Weight Functions for Marginal Proportional Hazards Analysis of Clustered Failure Time Data

The choice of weights in estimating equations for multivariate survival data is considered. Specifically, we consider families of weight functions which are constant on fixed time intervals, including the special case of time-constant weights. For a fixed set of time intervals, the optimal weights are identified as the solution to a system of linear equations. The optimal weights are computed for several scenarios. It is found that for the scenarios examined, the gains in efficiency using the optimal weights are quite small relative to simpler approaches except under extreme dependence, and that a simple estimator of an exchangeable approximation to the weights also performs well.     

Some Recent Developments in Econometric Inference

Abstract

Recent results in information theory, see Soofi (1996; 2001) for a review, include derivations of optimal information processing rules, including Bayes' theorem, for learning from data based on minimizing a criterion functional, namely output information minus input information as shown in Zellner (1988; 1991; 1997; 2002). Herein, solution post data densities for parameters are obtained and studied for cases in which the input information is that in (1) a likelihood function and a prior density; (2) only a likelihood function; and (3) neither a prior nor a likelihood function but only input information in the form of post data moments of parameters, as in the Bayesian method of moments approach. Then it is shown how optimal output densities can be employed to obtain predictive densities and optimal, finite sample structural coefficient estimates using three alternative loss functions. Such optimal estimates are compared with usual estimates, e.g., maximum likelihood, two‐stage least squares, ordinary least squares, etc. Some Monte Carlo experimental results in the literature are discussed and implications for the future are provided.     

Testing recapture model MT for departure from homogeneity

This paper presents a method of testing for animal heterogeneity when the capture effort is not constant from occasion to occasion. The basic set of statistics used is a lower triangular matrix whose (I,k)th element is the proportion of those animals caught on the next (i+1)the occasion who have been caught k times until now. Both theoretical and simulation evidence is presented that the new procedure is superior to methods previously suggested. The method may breakdown for certain kinds of behavioural response.     

A New Kind of Asymptotic Quasi-Score Estimating Function

A new definition of asymptotic quasi-score sequence of estimating functions is given and studied. The relationship between asymptotic quasi-likelihood and quasi-likelihood estimates is investigated. A new practical approach for obtaining a good estimate of θ in the model y _t = f_t (θ) + m_t without any prior knowledge on the nature of E ( m ² _t |F _{t −1}) is suggested, where f_t is a predictable process and m_t is a martingale difference process. Two examples are used to show that the approach is practicable.     

Some Recent Developments in Econometric Inference

Recent results in information theory, see Soofi (1996; 2001) for a review, include derivations of optimal information processing rules, including Bayes' theorem, for learning from data based on minimizing a criterion functional, namely output information minus input information as shown in Zellner (1988; 1991; 1997; 2002). Herein, solution post data densities for parameters are obtained and studied for cases in which the input information is that in (1) a likelihood function and a prior density; (2) only a likelihood function; and (3) neither a prior nor a likelihood function but only input information in the form of post data moments of parameters, as in the Bayesian method of moments approach. Then it is shown how optimal output densities can be employed to obtain predictive densities and optimal, finite sample structural coefficient estimates using three alternative loss functions. Such optimal estimates are compared with usual estimates, e.g., maximum likelihood, two-stage least squares, ordinary least squares, etc. Some Monte Carlo experimental results in the literature are discussed and implications for the future are provided.     

Optimal designs for the logit and probit models for binary data

We consider optimal designs for a class of symmetric models for binary data which includes the common probit and logit models. We show that for a large group of optimality criteria which includes the main ones in the literature (e.g. A-, D-, E-, F- and G-optimality) the optimal design for our class of models is a two-point design with support points symmetrically placed about the ED50 but with possibly unequal weighting. We demonstrate how one can further reduce the problem to a one-variable optimization by characterizing various of the common criteria. We also use the results to demonstrate major qualitative differences between the F - and c-optimal designs, two design criteria which have similar motivation.     

Optimal investment-consumption-insurance strategy in a continuous-time self-exciting threshold model

In this paper, we consider an optimal investment-consumption-insurance purchase problem for a wage earner. We assume that the price of the risky asset is governed by a continuous-time, finite state self-exciting threshold model. In this model, the state space of the price of the risky asset is partitioned by a set of thresholds and the parameters depend on the region which the current value of the price falls in. The wage earner’s objective is to find the optimal investment-consumption-insurance strategy that maximizes the expected discounted utilities. The optimal strategy for power utility function is derived by the martingale approach and the dynamic programming approach. Numerical examples are also provided to illustrate the effect of the thresholds.     

Design issues for population growth models

We briefly review and discuss design issues for population growth and decline models. We then use a flexible growth and decline model as an illustrative example and apply optimal design theory to find optimal sampling times for estimating model parameters, specific parameters and interesting functions of the model parameters for the model with two real applications. Robustness properties of the optimal designs are investigated when nominal values or the model is mis-specified, and also under a different optimality criterion. To facilitate use of optimal design ideas in practice, we also introduce a website for generating a variety of optimal designs for popular models from different disciplines.     

Estimation of the area under ROC curve with censored data

The area under the receiver operating characteristic curve is the most commonly used summary measure of diagnostic accuracy for a continuous-scale diagnostic test. In this paper, we develop methods to estimate the area under the curve (AUC) with censored data. Based on two different integration representations of this parameter, two nonparametric estimators are defined by the “plug in” method. Both the proposed estimators are shown to be asymptotically normal based on counting process and martingale theory. A simulation study is conducted to evaluate the performances of the proposed estimators.     

A Class of Non-parametric Tests in the Competing Risks Model for Comparing Two Samples

We consider a model when a process involving the production of elements is under inspection. The elements have possible failures due to competing risks. We assume the availability of a data set of failure times, D₁, obtained when the process is under control. Our main goal is to test if the failure rates in D₁ are equal to or less than the failure rates in another data set D₂, against "undesirable" neighbouring alternatives. A class of tests based on a two-dimensional vector statistic is obtained. Linear test statistics with weight functions giving optimal local asymptotic power are derived. Martingale techniques are used. An example illustrates the derivation of reasonable tests     

Semiparametric regression estimation for longitudinal data in models with martingale difference error's structure

An inherent characteristic of longitudinal data is the dependence among the observations within the same subject. For exhibiting dependencies among the observations within the same subject, this paper considers a semiparametric partially linear regression model for longitudinal data based on martingale difference error's structure. We establish a strong consistency for the least squares estimator of a parametric component and the estimator of a non-parametric function under some mild conditions. A simulation study shows the performance of the proposed estimator in finite samples.     

Generalized optimal design for two‐arm,randomized phase II clinical trials with endpoints from the exponential dispersion family

For two‐arm randomized phase II clinical trials, previous literature proposed an optimal design that minimizes the total sample sizes subject to multiple constraints on the standard errors of the estimated event rates and their difference. The original design is limited to trials with dichotomous endpoints. This paper extends the original approach to be applicable to phase II clinical trials with endpoints from the exponential dispersion family distributions. The proposed optimal design minimizes the total sample sizes needed to provide estimates of population means of both arms and their difference with pre‐specified precision. Its applications on data from specific distribution families are discussed under multiple design considerations. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of transformation models with right-truncated data

In many applications, statistical data are frequently observed subject to a retrospective sampling criterion resulting in right-truncated data. In this article, a general class of semiparametric transformation models that include proportional hazards model and proportional odds model as special cases is studied for the analysis of right-truncated data. We proposed two estimators for regression coefficients. The first estimator is based on martingale estimating equations. The second estimator is based on the conditional likelihood function given the truncation times. The asymptotic properties of both estimators are derived. The finite sample performance is examined through a simulation study.     

Limit theorems for dependent Bernoulli variables with statistical inference

We consider a class of dependent Bernoulli variables where the conditional success probability is a linear combination of the last few trials and the original success probability. We obtain its limit theorems including the strong law of large numbers, weak invariance principle, and law of the iterated logarithm. We also derive some statistical inference results which make the model applicable. Simulation results are exhibited as well to show that with small sample size the convergence rate is satisfying and the proposed estimators behave well.     

Zero-inflated proportion data models applied to a biological control assay

Biological control of pests is an important branch of entomology, providing environmentally friendly forms of crop protection. Bioassays are used to find the optimal conditions for the production of parasites and strategies for application in the field. In some of these assays, proportions are measured and, often, these data have an inflated number of zeros. In this work, six models will be applied to data sets obtained from biological control assays for Diatraea saccharalis , a common pest in sugar cane production. A natural choice for modelling proportion data is the binomial model. The second model will be an overdispersed version of the binomial model, estimated by a quasi-likelihood method. This model was initially built to model overdispersion generated by individual variability in the probability of success. When interest is only in the positive proportion data, a model can be based on the truncated binomial distribution and in its overdispersed version. The last two models include the zero proportions and are based on a finite mixture model with the binomial distribution or its overdispersed version for the positive data. Here, we will present the models, discuss their estimation and compare the results.     

Standardized Martingale Residuals Applied to Grouped Left Truncated Observations of Dementia Cases

The use of martingale residuals have been proposed for modelchecking and also to get a non-parametric estimate of the effectof an explanatory variable. We apply this approach to an epidemiologicalproblem which presents two characteristics: the data are lefttruncated due to delayed entry in the cohort; the data are groupedinto geographical units (parishes). This grouping suggests anatural way of smoothing the graph of residuals which is to computethe sum of the residuals for each parish. It is also naturalto present a graph with standardized residuals. We derive thevariances of the estimated residuals for left truncated datawhich allows computing the standardized residuals. This methodis applied to the study of dementia in a cohort of old people,and to the possible effect of the concentration of aluminum andsilica in drinking water on the risk of developing dementia.     

A weak convergence criterion for constructing changes of measure

Based on a weak convergence argument, we provide a necessary and sufficient condition that guarantees that a nonnegative local martingale is indeed a martingale. Typically, conditions of this sort are expressed in terms of integrability conditions (such as the well-known Novikov condition). The weak convergence approach that we propose allows to replace integrability conditions by a suitable tightness condition. We then provide several applications of this approach ranging from simplified proofs of classical results to characterizations of processes conditioned on first passage time events and changes of measures for jump processes.     

Testing the difference between two independent regression models

ABSTRACT

In some situations, for example, in biology or psychology studies, we wish to determine whether the linear relationship between response variable and predictor variables differs in two populations. The analysis of the covariance (ANCOVA) or, equivalently, the partial F-test approaches are the commonly used methods. In this study, the asymptotic distribution for the difference between two independent regression coefficients was established. The proposed method was used to derive the asymptotic confidence set for the difference between coefficients and hypothesis testing for the equality of the two regression models. Then a simulation study was conducted to compare the proposed method with the partial F method. The performance of the new method was comparable with that of the partial F method.     

On almost sure convergence for sums of stochastic sequence

In this work we discuss almost sure convergence for sums of arbitrarily dependent stochastic sequence under different conditions of Chung’s type. Our approach is based on the stopping time technique and the theorem of convergence for martingale difference sequence. Meanwhile, the results here include some relevant classical conclusions.