Interval Estimation for the Difference of Two Independent Variances

Analytical methods for interval estimation of differences between variances have not been described. A simple analytical method is given for interval estimation of the difference between variances of two independent samples. It is shown, using simulations, that confidence intervals generated with this method have close to nominal coverage even when sample sizes are small and unequal and observations are highly skewed and leptokurtic, provided the difference in variances is not very large. The method is also adapted for testing the hypothesis of no difference between variances. The test is robust but slightly less powerful than Bonett's test with small samples.     

Estimation for the three-parameter gamma distribution based on progressively censored data

Some work has been done in the past on the estimation for the three-parameter gamma distribution based on complete and censored samples. In this paper, we develop estimation methods based on progressively Type-II censored samples from a three-parameter gamma distribution. In particular, we develop some iterative methods for the determination of the maximum likelihood estimates (MLEs) of all three parameters. It is shown that the proposed iterative scheme converges to the MLEs. In this context, we propose another method of estimation which is based on missing information principle and moment estimators. Simple alternatives to the above two methods are also suggested. The proposed estimation methods are then illustrated with a numerical example. We also consider the interval estimation based on large-sample theory and examine the actual coverage probabilities of these confidence intervals in case of small samples using a Monte Carlo simulation study.     

Transformed Logit Confidence Intervals for Small Populations in Single Capture–Recapture Estimation

The good performance of logit confidence intervals for the odds ratio with small samples is well known. This is true unless the actual odds ratio is very large. In single capture–recapture estimation the odds ratio is equal to 1 because of the assumption of independence of the samples. Consequently, a transformation of the logit confidence intervals for the odds ratio is proposed in order to estimate the size of a closed population under single capture–recapture estimation. It is found that the transformed logit interval, after adding .5 to each observed count before computation, has actual coverage probabilities near to the nominal level even for small populations and even for capture probabilities near to 0 or 1, which is not guaranteed for the other capture–recapture confidence intervals proposed in statistical literature. Thus, given that the .5 transformed logit interval is very simple to compute and has a good performance, it is appropriate to be implemented by most users of the single capture–recapture method.     

Hierarchical Bayes estimation in small area estimation using cross-sectional and time-series data

Bayesian methods have been extensively used in small area estimation. A linear model incorporating autocorrelated random effects and sampling errors was previously proposed in small area estimation using both cross-sectional and time-series data in the Bayesian paradigm. There are, however, many situations that we have time-related counts or proportions in small area estimation; for example, monthly dataset on the number of incidence in small areas. This article considers hierarchical Bayes generalized linear models for a unified analysis of both discrete and continuous data with incorporating cross-sectional and time-series data. The performance of the proposed approach is evaluated through several simulation studies and also by a real dataset.     

Robust estimation of the exponential mean parameter for small samples: Complete and censored data

This paper studies methods for simple estimation of the exponential mean parameter in small samples in the presence of outliers. Existing estimation methods are discussed. Adaptation of these methods to allow for Type I censoring is investigated. New robust procedures are proposed. A series of simulation experiments Indicate trimming provides significant protection against outliers while the premium is usually small when trimming uncontarninated samples. A linearly weighted mean is recommended for uncontarninated samples, both censored and complete. In larger samples, (n - 20), the proposed Huber-type estimator performs quite well in all situations of censoring and contarnination     

Confidence limits for stress–strength reliability involving Weibull models

The problem of interval estimation of the stress–strength reliability involving two independent Weibull distributions is considered. An interval estimation procedure based on the generalized variable (GV) approach is given when the shape parameters are unknown and arbitrary. The coverage probabilities of the GV approach are evaluated by Monte Carlo simulation. Simulation studies show that the proposed generalized variable approach is very satisfactory even for small samples. For the case of equal shape parameter, it is shown that the generalized confidence limits are exact. Some available asymptotic methods for the case of equal shape parameter are described and their coverage probabilities are evaluated using Monte Carlo simulation. Simulation studies indicate that no asymptotic approach based on the likelihood method is satisfactory even for large samples. Applicability of the GV approach for censored samples is also discussed. The results are illustrated using an example.     

Inference for the Birnbaum–Saunders Lifetime Regression Model with Applications

The Birnbaum–Saunders distribution is a positively skewed distribution that is frequently used for analyzing lifetime data. Regression analysis is widely used in this context when some covariates are involved in the life-test. In this article, we discuss the maximum likelihood estimation of the model parameters and associated inference. We discuss the likelihood-ratio tests for some hypotheses of interest as well as some interval estimation methods. A Monte Carlo simulation study is then carried out to examine the performance of the proposed estimators and the interval estimation methods. Finally, some numerical data analyses are done for illustrating all the inferential methods developed here.     

Notes on estimation of the proportion ratio in the presence of carryover effects under the AB/BA crossover trial

The authors focus discussion on estimation of the proportion ratio (PR) in the presence of residual effects under the AB/BA design. Under a random effects multiplicative risk model, we develop three point estimators and three interval estimators accounting for residual effects for the PR. Using Monte Carlo simulations, we compare the performance of these point (or interval) estimators with point (or interval) estimators assuming no residual effects with respect to the bias and mean-squared-error (or the coverage probability and average length). The authors use the data taken from an AB/BA trial comparing two new inhalation devices to illustrate the use of these estimators.     

A comparison of two varying coefficient meta-analysis methods for an average risk difference

Two interval estimation methods for a general linear function of binomial proportions have been proposed. One method [Zou GY, Huang W, Zhang X. A note on confidence interval estimation for a linear function of binomial proportions. Comput Statist Data Anal. 2009;53:1080–1085] combines Wilson interval estimates of individual proportions, and the other method [Price RM, Bonett DG. An improved confidence interval for a linear function of binomial proportions. Comput Statist Data Anal. 2004;45:449–456] uses an adjusted Wald interval. Both methods are appropriate in varying coefficient meta-analysis models where the risk differences are allowed to vary across studies. The two methods were compared in a simulation study under realistic meta-analysis conditions and the adjusted Wald method was found to have the best performance characteristics.     

Bayes/frequentist compromise decision rules for Gaussian sampling

Bayesian methods have the potential to confer substantial advantages over frequentist when the assumed prior is approximately correct, but otherwise can perform poorly. Therefore, estimators and other inferences that strike a compromise between Bayes and frequentist optimality are attractive. To evaluate potential trade-offs, we study Bayes vs. frequentist risk under Gaussian sampling for families of point estimators and interval estimators. Bayes/frequentist compromises for interval estimation are more challenging than for point estimation, since performance involves an interplay between coverage and length. Each family allows ‘purchasing’ improved frequentist performance by allowing a small increase in Bayes risk over the Bayes rule. Any degree of increase can be specified, thus enabling greater or lesser trade-offs between Bayes and frequentist risk.     

Estimation for the three-parameter inverse Gaussian distribution under progressive Type-II censoring

Inverse Gaussian distribution has been used widely as a model in analysing lifetime data. In this regard, estimation of parameters of two-parameter (IG2) and three-parameter inverse Gaussian (IG3) distributions based on complete and censored samples has been discussed in the literature. In this paper, we develop estimation methods based on progressively Type-II censored samples from IG3 distribution. In particular, we use the EM-algorithm, as well as some other numerical methods for determining the maximum-likelihood estimates (MLEs) of the parameters. The asymptotic variances and covariances of the MLEs from the EM-algorithm are derived by using the missing information principle. We also consider some simplified alternative estimators. The inferential methods developed are then illustrated with some numerical examples. We also discuss the interval estimation of the parameters based on the large-sample theory and examine the true coverage probabilities of these confidence intervals in case of small samples by means of Monte Carlo simulations.     

Estimation of the Reliability of a Stress-Strength System from Power Lindley Distributions

In this paper, we are interested in the estimation of the reliability parameter R = P(X > Y) where X, a component strength, and Y, a component stress, are independent power Lindley random variables. The point and interval estimation of R, based on maximum likelihood, nonparametric and parametric bootstrap methods, are developed. The performance of the point estimate and confidence interval of R under the considered estimation methods is studied through extensive simulation. A numerical example, based on a real data, is presented to illustrate the proposed procedure.     

Zone control charts with estimated parameters for detecting prespecified changes in linear profiles

ABSTRACT

In profile monitoring, control charts are proposed to detect unanticipated changes, and it is usually assumed that the in-control parameters are known. However, due to the characteristics of a system or process, the prespecified changes would appear in the process. Moreover, in most applications, the in-control parameters are usually unknown. To overcome these issues, we develop the zone control charts with estimated parameters to detect small shifts of these prespecified changes. The effects of estimation error have been investigated on the performance of the proposed charts. To account for the practitioner-to-practitioner variability, the expected average run length (ARL) and the standard deviation of the average run length (SDARL) is used as the performance metrics. Our results show that the estimation error results in the significant variation in the ARL distribution. Furthermore, in order to adequately reduce the variability, more phase I samples are required in terms of the SDARL metric than that in terms of the expected ARL metric. In addition, more observations on each sampled profile are suggested to improve the charts' performance, especially for small phase I sample sizes. Finally, an illustrative example is given to show the performance of the proposed zone control charts.     

Prediction intervals for general balanced linear random models

The main interest of prediction intervals lies in the results of a future sample from a previously sampled population. In this article, we develop procedures for the prediction intervals which contain all of a fixed number of future observations for general balanced linear random models. Two methods based on the concept of a generalized pivotal quantity (GPQ) and one based on ANOVA estimators are presented. A simulation study using the balanced one-way random model is conducted to evaluate the proposed methods. It is shown that one of the two GPQ-based and the ANOVA-based methods are computationally more efficient and they also successfully maintain the simulated coverage probabilities close to the nominal confidence level. Hence, they are recommended for practical use. In addition, one example is given to illustrate the applicability of the recommended methods.     

Robustness of confidence intervals for scale parameters based on m-estimators

The robustness of confidence intervals for a scale parameter based on M-esimators is studied, especially in small size samples. The coverage probablity is used as measure of robustness. A theorem for a lower bound of the minimum coverage probability of M-estimators is presented and it is applied in order to examine the behavior of the standard deviation and the median absolute deviation, as interval estimators. This bound can confirm the robustness of any other scale M-estimator in interval estimation. The idea of stretching is used to formulate the family of distributions that are considered as underlying. Critical values for the confidence interval are computed where it is needed, that is for the median absolute deviation in the Normal, Uniform and Cauchy distribution and for the standard deviation in the Uniform and Cauchy distribution. Simulation results have been achieved for the estimation of the coverage probabilities and the critical values.     

The Performance of Two Data-Generation Processes for Data with Specified Marginal Treatment Odds Ratios

Monte Carlo simulation methods are increasingly being used to evaluate the property of statistical estimators in a variety of settings. The utility of these methods depends upon the existence of an appropriate data-generating process. Observational studies are increasingly being used to estimate the effects of exposures and interventions on outcomes. Conventional regression models allow for the estimation of conditional or adjusted estimates of treatment effects. There is an increasing interest in statistical methods for estimating marginal or average treatment effects. However, in many settings, conditional treatment effects can differ from marginal treatment effects. Therefore, existing data-generating processes for conditional treatment effects are of little use in assessing the performance of methods for estimating marginal treatment effects. In the current study, we describe and evaluate the performance of two different data-generation processes for generating data with a specified marginal odds ratio. The first process is based upon computing Taylor Series expansions of the probabilities of success for treated and untreated subjects. The expansions are then integrated over the distribution of the random variables to determine the marginal probabilities of success for treated and untreated subjects. The second process is based upon an iterative process of evaluating marginal odds ratios using Monte Carlo integration. The second method was found to be computationally simpler and to have superior performance compared to the first method.     

一种简单有效的二项分布比例参数似然比置信区间的求法

对二项分布比例参数p的似然比置信区间,提出一种简便求解方法。在平均覆盖率、平均区间长度及区间长度的95%置信区间准则下与WScore、Plus4、Jeffreys置信区间进行模拟比较。试验表明,在二项分布b(n,p)的参数n≥20且p∈(0.1,0.9)时,该方法获取的似然比置信区间性能优良。当点估计p值不是接近于0或1且n≥20时,推荐使用本方法获取p的置信区间。     

Notes on interval estimation of the proportion ratio under a three-treatment three-period crossover trial

When comparing two experimental treatments with a placebo, we focus our attention on interval estimation of the proportion ratio (PR) of patient responses under a three-period crossover design. We propose a random effects exponential multiplicative risk model and derive asymptotic interval estimators in closed form for the PR between treatments and placebo. Using Monte Carlo simulations, we compare the performance of these interval estimators in a variety of situations. We use the data comparing two different doses of an analgesic with placebo for the relief of primary dysmenorrhea to illustrate the use of these interval estimators and the difference in estimates of the PR and odds ratio (OR) when the underlying relief rates are not small.     

Semiparametric Inference for ROC Curves with Censoring

Abstract.  Comparison of two samples can sometimes be conducted on the basis of analysis of receiver operating characteristic (ROC) curves. A variety of methods of point estimation and confidence intervals for ROC curves have been proposed and well studied. We develop smoothed empirical likelihood-based confidence intervals for ROC curves when the samples are censored and generated from semiparametric models. The resulting empirical log-likelihood function is shown to be asymptotically chi-squared. Simulation studies illustrate that the proposed empirical likelihood confidence interval is advantageous over the normal approximation-based confidence interval. A real data set is analysed using the proposed method.