On the construction of orthogonal factorial designs of resolution 11

A new method of construction of orthogonal resolution IV designs for symmetrical and asymmetrical factorials has been presented. Many new series of orthogonal factorial designs of resolution IV can be obtained by the above general method.     

A note on incomplete block designs for symmetrical parallel line assays

Incomplete block designs estimating the relative potency free from block effects have been discussed. A method of constructing a series of such designs has been presented.     

Partial diallel grosses based on generalized right angj lar association scheme

Some new plans for partial diallel crosses (PDC) based on Generalized Right Angular (GRA) association scheme excluding parents as well as including parents have been presented. Unlike existing plans based on GRA, these plans have wider existence and provide sufficient size of PDC to have more precise conclusions. Simplified method of analysis for these plans has also been presented.     

General class of covariance structures for two or more repeated factors in longitudinal data analysis

The main difficulty in parametric analysis of longitudinal data lies in specifying covariance structure. Several covariance structures, which usually reflect one series of measurements collected over time, have been presented in the literature. However there is a lack of literature on covariance structures designed for repeated measures specified by more than one repeated factor. In this paper a new, general method of modelling covariance structure based on the Kronecker product of underlying factor specific covariance profiles is presented. The method has an attractive interpretation in terms of independent factor specific contribution to overall within subject covariance structure and can be easily adapted to standard software.     

Influence diagnostics in the Gamma regression model with adjusted deviance residuals

In this study, we develop the adjusted deviance residuals for the gamma regression model (GRM) by following Cordeiro's (2004) method. These adjusted deviance residuals under the GRM are used for influence diagnostics. A comparative analysis has been sorted out between our proposed method of the adjusted deviance residuals and an existing method for influence diagnostics. These results are illustrated by a simulation study and using a real data set. They are presented for different values of dispersion and sample sizes and indicate the significant role of the GRM inferences.     

Estimation of Weibull distribution parameters for irregular interval group failure data with unknown failure times

SUMMARY This paper presents three methods for estimating Weibull distribution parameters for the case of irregular interval group failure data with unknown failure times. The methods are based on the concepts of the piecewise linear distribution function (PLDF), an average interval failure rate (AIFR) and sequential updating of the distribution function (SUDF), and use an analytical approach similar to that of Ackoff and Sasieni for regular interval group data. Results from a large number of simulated case problems generated with specified values of Weibull distribution parameters have been presented, which clearly indicate that the SUDF method produces near-perfect parameter estimates for all types of failure pattern. The performances of the PLDF and AIFR methods have been evaluated by goodness-of-fit testing and statistical confidence limits on the shape parameter. It has been found that, while the PLDF method produces acceptable parameter estimates, the AIFR method may fail for low and high shape parameter values that represent the cases of random and wear-out types of failure. A real-life application of the proposed methods is also presented, by analyzing failures of hydrogen make-up compressor valves in a petroleum refinery.     

Adapting sample size calculations to repeated measurements in clinical trials

Many of the repeated-measures sample size calculation methods presented in the literature are not suitable when: •the different treatments are assumed to be equal on average at baseline time due to randomization, •and the experimenters are interested in a pre-specified difference to be detected after a specific time period. The method presented here has been developed for those cases where a multivariate normal distribution can reasonably be assumed. It is likelihood-based and has been designed to be flexible enough to handle repeated-measures models, including a non-linear change in time, and an arbitrary correlation structure.     

Minimum Risk Invariant Estimators of a Continuous Cumulative Distribution Function

The problem of nonparametric minimum risk invariant estimation has engaged a good deal of attention in the literature and minimum risk invariant estimators (MRIE's) have been constructed for some special statistical models. We present a new and simple method of obtaining the MRIE's of a continuous cumulative distribution function (cdf) under a general invariant loss function. All the MRIE's, which are known from the literature, can be constructed by the method presented in the article, in particular, under the weighted quadratic, LINEX and entropy loss functions. This method enables also to construct the MRIE's in nonparametric statistical models which have not been considered until now. In particular, considering a family of nonparametric precautionary loss functions, a new class of MRIE's of the cdf has been found. We also give some general remarks on obtaining the MRIE's and a review concerning minimaxity and admissibility of MRIE's.     

Sample size calculation for a proportional hazards mixture cure model with nonbinary covariates

Sample size calculation is a critical issue in clinical trials because a small sample size leads to a biased inference and a large sample size increases the cost. With the development of advanced medical technology, some patients can be cured of certain chronic diseases, and the proportional hazards mixture cure model has been developed to handle survival data with potential cure information. Given the needs of survival trials with potential cure proportions, a corresponding sample size formula based on the log-rank test statistic for binary covariates has been proposed by Wang et al. [25]. However, a sample size formula based on continuous variables has not been developed. Herein, we presented sample size and power calculations for the mixture cure model with continuous variables based on the log-rank method and further modified it by Ewell's method. The proposed approaches were evaluated using simulation studies for synthetic data from exponential and Weibull distributions. A program for calculating necessary sample size for continuous covariates in a mixture cure model was implemented in R.     

Balanced block designs

In the present paper an attempt has been made to characterize and unify the three different concepts of balancing in incomplete block designs, namely (i) variance balance, (ii) efficiency balance and (iii) pairwise balance. Simple characterizations of variance balance and efficiency balance have been given using the P matrix. A method of constructing efficiency balanced (EB) and variance balanced designs has also been presented.     

An iterative Monte Carlo method for nonconjugate Bayesian analysis

The Gibbs sampler has been proposed as a general method for Bayesian calculation in Gelfand and Smith (1990). However, the predominance of experience to date resides in applications assuming conjugacy where implementation is reasonably straightforward. This paper describes a tailored approximate rejection method approach for implementation of the Gibbs sampler when nonconjugate structure is present. Several challenging applications are presented for illustration.     

Bayesian projection of the acquired immune deficiency syndrome epidemic

Short-term projections of the acquired immune deficiency syndrome (AIDS) epidemic in England and Wales have been regularly updated since the publication of the Cox report in 1988. The key approach for those updates has been the back-calculation method, which has been informally adapted to acknowledge various sources of uncertainty as well as to incorporate increasingly available information on the spread of the human immunodeficiency virus (HIV) in the population. We propose a Bayesian formulation of the back-calculation method which allows a formal treatment of uncertainty and the inclusion of extra information, within a single coherent composite model. Estimation of the variably dimensioned model is carried out by using reversible-jump Markov chain Monte Carlo methods. Application of the model to data for homosexual and bisexual males in England and Wales is presented, and the role of the various sources of information and model assumptions is appraised. Our results show a massive peak in HIV infections around 1983 and suggest that the incidence of AIDS has now reached a plateau, although there is still substantial uncertainty about the future.     

Influence diagnostic analysis in the possibly heteroskedastic linear model with exact restrictions

The local influence method has proven to be a useful and powerful tool for detecting influential observations on the estimation of model parameters. This method has been widely applied in different studies related to econometric and statistical modelling. We propose a methodology based on the Lagrange multiplier method with a linear penalty function to assess local influence in the possibly heteroskedastic linear regression model with exact restrictions. The restricted maximum likelihood estimators and information matrices are presented for the postulated model. Several perturbation schemes for the local influence method are investigated to identify potentially influential observations. Three real-world examples are included to illustrate and validate our methodology.     

Non-existence of an optimum estimator in a class of ratio estimators

A number of estimators formulated in the field of the ratio method of estimation has been presented. A class of estimators encompassing these estimators is constructed. It is noted that an optimum estimator does not exist uniformly in this class. The “Optimum” so obtained reduces to the usual regression estimator.     

Unbiased estimation following selection for bernoulli populations

Although a large number of selection procedures have been published in the statistics literature, the selection approach has received only limited use in applications. One drawback to the use of such procedures has been the lack of parameter estimates, which prevents quantitative comparisons among the treatments. To partially address this criticism, we present a general method for constructing unbiased estimators of the success probabilities after the termination of a sequential experiment involving two or more Bernoulli populations. Some theoretical properties are presented, and examples are provided for several different selection procedures.     

Understanding Canonical Correlation through the General Linear Model and Principal Components

Canonical correlation has been little used and little understood, even by otherwise sophisticated analysts. An alternative approach to canonical correlation, based on a general linear multivariate model, is presented. Properties of principal component analysis are used to help explain the method. Standard computational methods for full rank canonical correlation, techniques for canonical correlation on component scores, and canonical correlation with less than full rank are discussed. They are seen to be essentially equivalent when the model equation for canonical correlation on component scores is presented. The two approaches to less than full rank situations are equivalent in some senses, but quite different in usefulness, depending on the application. An example dataset is analyzed in detail to help demonstrate the conclusions.     

Bayesian Estimation of Bivariate Exponential Distributions Based on Linex and Quadratic Loss Functions: A Survival Approach with Censored Samples

In this article, four bivariate exponential (BVE) distributions with subject to right censoring samples are presented. Bayesian estimates of the parameters of BVE are obtained through Linex and quadratic loss functions. Gamma prior distribution has been suggested to reforming the posterior function. The estimations and standard errors of parameters have also been obtained through simulation method. Markov chain Monte Carlo (MCMC) method is employed for the case of Block-Buse bivariate distribution because there was no closed form for estimator criteria. Simulation studies have been conducted to show that the computation parts can be implemented easily and comparing the estimated values due to two methods and with the true values as well.     

Distance-based outlier detection for high dimension,low sample size data

Despite the popularity of high dimension, low sample size data analysis, there has not been enough attention to the sample integrity issue, in particular, a possibility of outliers in the data. A new outlier detection procedure for data with much larger dimensionality than the sample size is presented. The proposed method is motivated by asymptotic properties of high-dimensional distance measures. Empirical studies suggest that high-dimensional outlier detection is more likely to suffer from a swamping effect rather than a masking effect, thus yields more false positives than false negatives. We compare the proposed approaches with existing methods using simulated data from various population settings. A real data example is presented with a consideration on the implication of found outliers.     

Random balanced resampling: A new method for estimating variance components in unbalanced designs

Variance components in factorial designs with balanced data are commonly estimated by equating mean squares to expected mean squares. For unbalanced data, the usual extensions of this approach are the Henderson methods, which require formulas that are rather involved. Alternatively, maximum likelihood estimation based on normality has been proposed. Although the algorithm for maximum likelihood is computationally complex, programs exist in some statistical packages. This article introduces a simpler method, that of creating a balanced data set by resampling from the original one. Revised formulas for expected mean squares are presented for the two-way case; they are easily generalized to larger factorial designs. The results of a number of simulation studies indicate that, in certain types of designs, the proposed method has performance advantages over both the Henderson Method I and maximum likelihood estimators.     

Taguchi methods: linear graphs of high resolution

An important reason behind the success of the Taguchi methodology in qual- ity assurance has been the use of statistical methods, presented in a way that is accessible to the nonexpert user. Among the tools used to simplify the sta- tistical design of experiments has been the linear graph, apparently introduced by Taguchi. However, he did not consider the resolution of the corresponding designs (the higher the resolution, the more accurate the conclusions). For example, it will be shown that half of the linear graphs given by Taguchi for the L₁₆(2¹⁵) orthogonal array correspond to designs of resolution III, when designs of resolution IV are available (with the same lines in the linear graphs but with different assignments to the columns of the orthogonal array). A nontraditional but very straightforward method is presented for obtaining the alias chains and the linear graphs corresponding to an orthogonal array. The procedure can be easily understood and employed by nonstatisticians to find an experimental design of the highest possible resolution. The design can be used to obtain products or processes that are robust to variation.