Applying skovgaard's modified directed likelihood statistic to mixed linear models

The introduction of software to calculate maximum likelihood estimates for mixed linear models has made likelihood estimation a practical alternative to methods based on sums of squares. Likelihood based tests and confidence intervals, however, may be misleading in problems with small sample sizes. This paper discusses an adjusted version of the directed log-likelihood statistic for mixed models that is highly accurate for testing one parameter hypotheses. Indroduced by Skovgaard (1996, Journal of the Bernoulli Society,2,145-165), we show in mixed models that the statistic has a simple conpact from that may be obtained from standard software. Simulation studies indicate that this statistic is more accurate than many of the specialized procedure that have been advocated.     

Assessing the bias of maximum likelihood estimates of contaminated garch models

It is well known that Gaussian maximum likelihood estimates of time series models are not robust. In this paper we prove this is also the case for the Generalized Autoregressive Conditional Heteroscedastic (GARCH) models. By expressing the Gaussian maximum likelihood estimates as Ψ estimates and by assuming the existence of a contaminated process, we prove they possess zero breakdown point and unbounded influence curves. By simulating GARCH processes under several proportions of contaminations we assess how much biased the maximum likelihood estimates may become and compare these results to a robust alternative. The t-student maximum likelihood estimates of GARCH models are also considered.     

Exact inference and prediction for k-sample exponential case under type-ii censoring

The exact inference and prediction intervals for the K-sample exponential scale parameter under doubly Type-II censored samples are derived using an algorithm of Huffer and Lin [Huffer, F.W. and Lin, C.T., 2001, Computing the joint distribution of general linear combinations of spacings or exponen-tial variates. Statistica Sinica, 11, 1141–1157.]. This approach provides a simple way to determine the exact percentage points of the pivotal quantity based on the best linear unbiased estimator in order to develop exact inference for the scale parameter as well as to construct exact prediction intervals for failure times unobserved in the ith sample. Similarly, exact prediction intervals for failure times of units from a future sample can also be easily obtained.     

On a redundant repairable system with switching failure: Bayesian approach

System characteristics of a redundant repairable system are studied from a Bayesian viewpoint with different types of priors assumed for the unknown parameters. The system consists of two primary units, one standby unit, and one repair facility which is activated when switching to standby fails. Times to failure and times to repair of the operating units are assumed to follow exponential distributions. When time to failure and time to repair have uncertain parameters, a Bayesian approach is adopted to evaluate system characteristics. Monte Carlo simulation is used to derive the posterior distribution for the mean time to system failure and steady-state availability. Some numerical experiments are performed to illustrate the results derived in this paper.     

Algorithm for determining whether various two-level fractional factorial split-plot row–column designs are non-isomorphic

The confounding and aliasing scheme for fractional factorial split-plot designs with the units within each wholeplot arranged in rows and columns is described and illustrated. Isomorphism for this design type is described, together with a procedure which considers extensions of the concepts of wordlength patterns and letter patterns that can be used to test isomorphism between designs. Using in part this isomorphism testing procedure, a construction algorithm that may be used to obtain a complete set of such non-isomorphic two-level designs is described. Software based on this construction algorithm was used to obtain a complete set of non-isomorphic designs for up to five wholeplot factors, five subplot factors and up to 64 runs, which is presented as a table of designs. To aid the experimenter in distinguishing between competing designs, the estimation capacity sequence for each design is presented.     

Estimating monotone convex functions via sequential shape modification

We propose a sequential method to estimate monotone convex functions that consists of: (i) monotone regression via solving a constrained least square (LS) problem and (ii) convexification of the monotone regression estimate via solving a uniform approximation problem with associated constraints. We show that this method is faster than the constrained LS method. The ratio of computation time increases as data size increases. Moreover, we show that, under an appropriate smoothness condition, the uniform convergence rate achieved by the proposed method is nearly comparable to the best achievable rate for a non-parametric estimate which ignores the shape constraint. Simulation studies show that our method is comparable to the constrained LS method in estimation error. We illustrate our method by analysing ground water level data of wells in Korea.     

A Monte Carlo comparison of some variance estimators of the Horvitz–Thompson estimator

At the design and estimation stage of a survey, large survey organization often uses auxiliary information. This article discusses various procedures for improving variance estimation of the Horvitz–Thompson estimator of a finite population total with the aid of auxiliary information. To study the design-based properties of the proposed variance estimators relative to the standard one, a small scale Monte Carlo study is performed.     

A simple method for estimating parameters of the location–scale distribution family

The purpose of this paper is to estimate the parameters of the location–scale distribution family. As a special case, the method is used for estimating the parameters of the normal distribution and Cauchy distribution. For the Cauchy distribution, neither the moment estimation method nor the maximum likelihood estimation method works properly for estimating the parameters. The quantiles for obtaining confidence intervals and point estimates for the parameters of the two-parameter Cauchy distribution are given in the paper. It is shown that the estimators obtained in this paper are unbiased with respect to the median and possess some optimal properties.     

Correlation estimation in Downton's bivariate exponential distribution using incomplete samples

Suppose (X, Y) has a Downton's bivariate exponential distribution with correlation ρ. For a random sample of size n from (X, Y), let X _r:n be the rth X-order statistic and Y _[r:n] be its concomitant. We investigate estimators of ρ when all the parameters are unknown and the available data is an incomplete bivariate sample made up of (i) all the Y-values and the ranks of associated X-values, i.e. (i, Y _[i:n]), 1≤i≤n, and (ii) a Type II right-censored bivariate sample consisting of (X _i:n , Y _[i:n]), 1≤i≤r<n. In both setups, we use simulation to examine the bias and mean square errors of several estimators of ρ and obtain their estimated relative efficiencies. The preferred estimator under (i) is a function of the sample correlation of (Y _i:n , Y _[i:n]) values, and under (ii), a method of moments estimator involving the regression function is preferred.     

Modified BLUEs and BLIEs of the location and scale parameters and the population mean using ranked set sampling

Ranked set sampling (RSS) is a sampling procedure that can be used to improve the cost efficiency of selecting sample units of an experiment or a study. In this paper, RSS is considered for estimating the location and scale parameters a and b>0, as well as the population mean from the family F((x?a)/b). Modified best linear unbiased estimators (BLUEs) and best linear invariant estimators (BLIEs) are considered. Numerical computations with different location-scale distributions and different sample sizes are conducted to assess the efficiency of the suggested estimators. It is found that the modified BLIEs are uniformly higher than that of BLUEs for all distributions considered in this study. The modified BLUE and BLIE are more efficient when the underlying distribution is symmetric.