An F Approximation to the Distribution of a Linear Combination of Chi-squared Variables.

A three-parameter F approximation to the distribution of a positive linear combination of central chi-squared variables is described. It is about as easy to implement as the Satterthwaite-Welsh and Hall-Buckley-Eagleson approximations. Some reassuring properties of the F approximation are derived, and numerical results are presented. The numerical results indicate that the new approximation is superior to the Satterthwaite approximation and, for some purposes, better than the Hall-Buckley-Eagleson approximation. It is not quite as good as the Gamma-Weibull approximation due to Solomon and Stephens, but is easier to implement because iterative methods are not required.     

Minimax estimation of a lower-bounded scale parameter of a gamma distribution for scale-invariant squared-error loss

Let X have a gamma distribution with known shape parameter θr;aL and unknown scale parameter θ. Suppose it is known that θ ≥ a for some known a > 0. An admissible minimax estimator for scale-invariant squared-error loss is presented. This estimator is the pointwise limit of a sequence of Bayes estimators. Further, the class of truncated linear estimators C = {θ_ρ|θ_ρ(x) = max(a, ρ), ρ > 0} is studied. It is shown that each θ_ρ is inadmissible and that exactly one of them is minimax. Finally, it is shown that Katz's [Ann. Math. Statist., 32, 136–142 (1961)] estimator of θ is not minimax for our loss function. Some further properties of and comparisons among these estimators are also presented.     

A Parametric Model for the Interval Censored Survival Times of Acacia Mangium Plantation in a Spacing Trial

Survival times for the Acacia mangium plantation in the Segaliud Lokan Project, Sabah, East Malaysia were analysed based on 20 permanent sample plots (PSPs) established in 1988 as a spacing experiment. The PSPs were established following a complete randomized block design with five levels of spacing randomly assigned to units within four blocks at different sites. The survival times of trees in years are of interest. Since the inventories were only conducted annually, the actual survival time for each tree was not observed. Hence, the data set comprises censored survival times. Initial analysis of the survival of the Acacia mangium plantation suggested there is block by spacing interaction; a Weibull model gives a reasonable fit to the replicate survival times within each PSP; but a standard Weibull regression model is inappropriate because the shape parameter differs between PSPs. In this paper we investigate the form of the non-constant Weibull shape parameter. Parsimonious models for the Weibull survival times have been derived using maximum likelihood methods. The factor selection for the parameters is based on a backward elimination procedure. The models are compared using likelihood ratio statistics. The results suggest that both Weibull parameters depend on spacing and block.     

Restricted tests for testing independence of time to failure and cause of failure in a competing-risks model

We consider the competing-risks problem without making any assumption concerning the independence of the risks. Maximum-likelihood estimates of the cause-specific hazard rates are obtained under the condition that their ratio is monotone. We also consider the likelihood-ratio test for testing the proportionality of two cause-specific hazard rates against the alternative that the ratio of these two hazard rates is monotonic. This testing problem is equivalent to testing independence against likelihood-ratio dependence of the time to failure and the cause of failure in the competing-risks setup. We allow for random censoring on the right. The asymptotic null distribution of the test statistic is obtained and is found to be of the chi-bar-square type. The problem is extended to the case of more than two risks. A numerical example is given to illustrate the procedure.     

Counterexamples to a result of broemeling on simultaneous inferences for variance ratios of some mixed linear models

In a recent paper5 Broemeling (1978) extended his earlier work on one-sided confidence regions for the variance ratios of balanced random-effects models to the two-sided case. The extension depends on a probability Inequality which was claimed to be tru We show here that it is false, hence the proof of the main result given in Ms parer is in error W also show Lhat the ntatement of his result remains true in certain special cases.     

Utilization of treatment effect on a surrogate endpoint for planning a study to evaluate treatment effect on a final endpoint

To design a phase III study with a final endpoint and calculate the required sample size for the desired probability of success, we need a good estimate of the treatment effect on the endpoint. It is prudent to fully utilize all available information including the historical and phase II information of the treatment as well as external data of the other treatments. It is not uncommon that a phase II study may use a surrogate endpoint as the primary endpoint and has no or limited data for the final endpoint. On the other hand, external information from the other studies for the other treatments on the surrogate and final endpoints may be available to establish a relationship between the treatment effects on the two endpoints. Through this relationship, making full use of the surrogate information may enhance the estimate of the treatment effect on the final endpoint. In this research, we propose a bivariate Bayesian analysis approach to comprehensively deal with the problem. A dynamic borrowing approach is considered to regulate the amount of historical data and surrogate information borrowing based on the level of consistency. A much simpler frequentist method is also discussed. Simulations are conducted to compare the performances of different approaches. An example is used to illustrate the applications of the methods.     

Additional approximations to the distribution of a structural coefficient estimator for the case of two included endogenous variables

Four new approximations t o the exact distribution of the two-stage l e a s t squares estimator of astructuralcoefficient for

the case of two included endogeneous variables are introduced and compared with the others in the literatur e . Two of the new approximations are based on the Pearson distribution and are found to be adequate throughout the parameter space. A normal approximation using exact moments and an approximation based on the saddlepoint method (Holly and Phillips,1979) are found to be

poor for a wide range of parameter values.     

Prediction of Times to Failure of Censored Items for a Simple Step-Stress Model with Regular and Progressive Type I Censoring from the Exponential Distribution

The problem of predicting times to failure of units from the Exponential Distribution which are censored under a simple step-stress model is considered in this article. We discuss two types of censoring—regular and progressive Type I—and two kinds of predictors—the maximum likelihood predictors (MLP) and the conditional median predictors (CMP) for each type of censoring. Numerical examples are used to illustrate the prediction methods. Using simulation studies, mean squared prediction error (MSPE) and prediction intervals are generated for these examples. MLP and the CMP are then compared with respect to MSPE and the prediction interval.     

Sample size needed to get given ratio of endpoints for confidence interval of standard deviation in a normal distribution

Confidence interval (CI) for a standard deviation in a normal distribution, based on pivotal quantity with a Chi-square distribution, is considered. As a measure of CI quality, the ratio of its endpoints is taken. There are given formulas for sample sizes so that this ratio does not exceed a fixed value. Both equally tailed and minimum ratio of endpoint CIs are considered.     

Recurrence relations for single and product moments of order statistics from a generalized half logistic distribution with applications to inference

In this paper, we derive several recurrence relations satisfied by the single and product moments of order statistics from a generalized half logistic distribution. These generalize the corresponding results for the half logistic distribution established by Balakrishnan (1985). The relations established in this paper will enable one to compute the single and product moments of all order statistics for all sample sizes in a simple recursive manner; this may be done for any choice of the shape parameter k. These moments can then be used to determine the best linear unbiased estimators of location and scale parameters from complete as well as Type-I1 censored samples.     

Application of a predictive distribution formula to Bayesian computation for incomplete data models

We consider exact and approximate Bayesian computation in the presence of latent variables or missing data. Specifically we explore the application of a posterior predictive distribution formula derived in Sweeting And Kharroubi (2003), which is a particular form of Laplace approximation, both as an importance function and a proposal distribution. We show that this formula provides a stable importance function for use within poor man’s data augmentation schemes and that it can also be used as a proposal distribution within a Metropolis-Hastings algorithm for models that are not analytically tractable. We illustrate both uses in the case of a censored regression model and a normal hierarchical model, with both normal and Student t distributed random effects. Although the predictive distribution formula is motivated by regular asymptotic theory, it is not necessary that the likelihood has a closed form or that it possesses a local maximum.     

A time-varying measurement error model for age of onset of a psychiatric diagnosis: applied to first depressive episode diagnosed by the Composite International Diagnostic Interview (CIDI)

A substantial degree of uncertainty exists surrounding the reconstruction of events based on memory recall. This form of measurement error affects the performance of structured interviews such as the Composite International Diagnostic Interview (CIDI), an important tool to assess mental health in the community. Measurement error probably explains the discrepancy in estimates between longitudinal studies with repeated assessments (the gold-standard), yielding approximately constant rates of depression, versus cross-sectional studies which often find increasing rates closer in time to the interview. Repeated assessments of current status (or recent history) are more reliable than reconstruction of a person's psychiatric history based on a single interview. In this paper, we demonstrate a method of estimating a time-varying measurement error distribution in the age of onset of an initial depressive episode, as diagnosed by the CIDI, based on an assumption regarding age-specific incidence rates. High-dimensional non-parametric estimation is achieved by the EM-algorithm with smoothing. The method is applied to data from a Norwegian mental health survey in 2000. The measurement error distribution changes dramatically from 1980 to 2000, with increasing variance and greater bias further away in time from the interview. Some influence of the measurement error on already published results is found.     

Revisit of a randomized response model for estimating a rare sensitive attribute under probability proportional to size sampling using Poisson probability distribution

The present article deals with the estimation of mean number of individuals possess a rare sensitive attribute using Poisson probability distribution, when the population consists of clusters. Unbiased estimation procedures for the mean number of individuals have been suggested and their properties are discussed when the parameter of a rare non-sensitive unrelated attribute is assumed to be known as well as unknown. The suggested estimation procedure is further discussed for situation of stratified cluster population. Empirical studies are carried out to show the dominance of proposed method and resultant estimators over a well-known contemporary estimator.     

Estimation of a rare sensitive attribute for two-stage randomized response model in probability proportional to size sampling using Poisson probability distribution

This paper addresses the estimation of the mean number of individuals in the population who possess a rare sensitive attribute using Poisson distribution for the situations of (i) clustered population and (ii) stratified population with clusters are strata units. Properties of the proposed estimation procedures have been discussed when the proportion of a rare unrelated non-sensitive attribute is assumed to be known as well as unknown. Empirical studies are carried out to support the theoretical results which showed dominance over Lee et al. [Estimation of a rare sensitive attribute in probability proportional to size measures using Poisson distribution. Statistics (Ber). 2014;48(3):685–709] estimation procedures.     

EM algorithm-based likelihood estimation for a generalized Gompertz regression model in presence of survival data with long-term survivors: an application to uterine cervical cancer data

In this paper we develop a regression model for survival data in the presence of long-term survivors based on the generalized Gompertz distribution introduced by El-Gohary et al. [The generalized Gompertz distribution. Appl Math Model. 2013;37:13–24] in a defective version. This model includes as special case the Gompertz cure rate model proposed by Gieser et al. [Modelling cure rates using the Gompertz model with covariate information. Stat Med. 1998;17:831–839]. Next, an expectation maximization algorithm is then developed for determining the maximum likelihood estimates (MLEs) of the parameters of the model. In addition, we discuss the construction of confidence intervals for the parameters using the asymptotic distributions of the MLEs and the parametric bootstrap method, and assess their performance through a Monte Carlo simulation study. Finally, the proposed methodology was applied to a database on uterine cervical cancer.     

Recurrence relations for single and product moments of progressively Type-II censored order statistics from a generalized half-logistic distribution with application to inference

In this paper, we establish several recurrence relations for the single and product moments of progressively Type-II right-censored order statistics from a generalized half-logistic distribution. The use of these relations in a systematic recursive manner enables the computation of all the means, variances, and covariances of progressively Type-II right-censored order statistics from the generalized half-logistic distribution for all sample sizes n, effective sample sizes m, and all progressive censoring schemes (R ₁, …, R _m ). The results established here generalize the corresponding results for the usual order statistics due to Balakrishnan and Sandhu [Recurrence relations for single and product moments of order statistics from a generalized half-logistic distribution with applications to inference, J. Stat. Comput. Simul. 52 (1995), pp. 385–398.]. The moments so determined are then utilized to derive the best linear unbiased estimators of the scale and location–scale parameters of the generalized half-logistic distribution. The best linear unbiased predictors of censored failure times are discussed briefly. Finally, a numerical example is presented to illustrate the inferential method developed here.