Nonparametric Shrinkage Estimation for Aalen's Additive Hazards Model

Aalen's nonparametric additive model in which the regression coefficients are assumed to be unspecified functions of time is a flexible alternative to Cox's proportional hazards model when the proportionality assumption is in doubt. In this paper, we incorporate a general linear hypothesis into the estimation of the time‐varying regression coefficients. We combine unrestricted least squares estimators and estimators that are restricted by the linear hypothesis and produce James‐Stein‐type shrinkage estimators of the regression coefficients. We develop the asymptotic joint distribution of such restricted and unrestricted estimators and use this to study the relative performance of the proposed estimators via their integrated asymptotic distributional risks. We conduct Monte Carlo simulations to examine the relative performance of the estimators in terms of their integrated mean square errors. We also compare the performance of the proposed estimators with a recently devised LASSO estimator as well as with ridge‐type estimators both via simulations and data on the survival of primary billiary cirhosis patients.     

Confidence Intervals for Linear Combinations of Poisson Means

Parametric confidence intervals are given for linear combinations of the means of independent Poisson variables and for their continuous versions. The performance of the intervals is assessed using simulation. A real data set is used to compare the proposed intervals with known ones. The proposed intervals are shown to be superior to known ones and comparable to exact intervals.     

Semi-polynomial copulas

In this paper, we characterise a family of bivariate copulas whose sections between the main diagonal and the border of the unit square are polynomial, generalising several families of copulas, including those with quadratic and cubic sections. We also study a measure of association and the tail dependence for this class, illustrating our results with several examples.     

Chain sampling plan for variables inspection

SUMMARY This paper extends the concept of chain sampling to variables inspection when the standard deviation of the normally distributed characteristic is known. A discussion of the shape of the known sigma single-sampling variables plan is given. The chain sampling plan for variables inspection will be useful when testing is costly or destructive.     

Tightened single-level continuous sampling plan

In this paper, a new tightening concept has been incorporated into the single-level continuous sampling plan CSP-1, such that quality degradation will warrant sampling inspection to cease beyond a certain number of sampled items, until new evidence of good quality is established. The expressions of the performance measures for this new plan, such as the operating characteristic, average outgoing quality and average fraction inspected, are derived using a Markov chain model. The advantage of the tightened CSP-1 plan is that it is possible to lower the average outgoing quality limit.     

Estimation of the number of studies with positive trends when studies with negative trends are present

A two-point estimator is proposed for the proportion of studies with positive trends among a collection of studies, some of which may demonstrate negative trends. The proposed estimator is the y-intercept of the secant line joining the points (a, F?(a)) and (b, F?(b)), where F?(p) is the empirical distribution function of p-values from one-tailed tests for positive trend derived from the individual studies. Although this estimator is negatively biased for any choice of the points 0 ≤ a < b ≤ 1, the bias is less than that of the previously proposed one-point estimator defined by setting b = 1. The bias of the two-point estimator is smallest when a and b approach the inflection point of the true distribution function, E [F?(p)]. The utility of the two-point estimator is demonstrated by using it to estimate the number of male-mouse liver carcinogens among carcinogenicity studies conducted by the National Toxicology Program.     

A new class of measurement-error models,with applications to dietary data

Measurement-error modelling occurs when one cannot observe a covariate, but instead has possibly replicated surrogate versions of this covariate measured with error. The vast majority of the literature in measurement-error modelling assumes (typically with good reason) that given the value of the true but unobserved (latent) covariate, the replicated surrogates are unbiased for latent covariate and conditionally independent. In the area of nutritional epidemiology, there is some evidence from biomarker studies that this simple conditional independence model may break down due to two causes: (a) systematic biases depending on a person's body mass index, and (b) an additional random component of bias, so that the error structure is the same as a one-way random-effects model. We investigate this problem in the context of (1) estimating distribution of usual nutrient intake, (2) estimating the correlation between a nutrient instrument and usual nutrient intake, and (3) estimating the true relative risk from an estimated relative risk using the error-prone covariate. While systematic bias due to body mass index appears to have little effect, the additional random effect in the variance structure is shown to have a potentially important effect on overall results, both on corrections for relative risk estimates and in estimating the distribution of usual nutrient intake. However, the effect of dietary measurement error on both factors is shown via examples to depend strongly on the data set being used. Indeed, one of our data sets suggests that dietary measurement error may be masking a strong risk of fat on breast cancer, while for a second data set this masking is not so clear. Until further understanding of dietary measurement is available, measurement-error corrections must be done on a study-specific basis, sensitivity analyses should be conducted, and even then results of nutritional epidemiology studies relating diet to disease risk should be interpreted cautiously.     

Quantitative convergence assessment for Markov chain Monte Carlo via cusums

Yu (1995) provides a novel convergence diagnostic for Markov chain Monte Carlo (MCMC) which provides a qualitative measure of mixing for Markov chains via a cusum path plot for univariate parameters of interest. The method is based upon the output of a single replication of an MCMC sampler and is therefore widely applicable and simple to use. One criticism of the method is that it is subjective in its interpretation, since it is based upon a graphical comparison of two cusum path plots. In this paper, we develop a quantitative measure of smoothness which we can associate with any given cusum path, and show how we can use this measure to obtain a quantitative measure of mixing. In particular, we derive the large sample distribution of this smoothness measure, so that objective inference is possible. In addition, we show how this quantitative measure may also be used to provide an estimate of the burn-in length for any given sampler. We discuss the utility of this quantitative approach, and highlight a problem which may occur if the chain is able to remain in any one state for some period of time. We provide a more general implementation of the method to overcome the problem in such cases.     

A Note on Regression-Type Estimators Using Multiple Auxiliary Information

Kiregyera (1984), Mukerjee, Rao & Vijayan (1987), and Tripathi & Ahmed (1995) considered a number of regression-type estimators where information on two auxiliary variables related to study variable is available at different levels. Mukerjee et al . (1987) suggested three estimators and computed their mean square errors, but the computations seem to be incorrect. This note corrects them, and finds their estimators are no better than that of Kiregyera (1984). The estimator suggested by Tripathi & Ahmed (1995) is the best in the sense of having the smallest mean square error.