Analyzing Supersaturated Designs by Means of an Information Based Criterion

The cost and time consumption of many industrial experimentations can be reduced using the class of supersaturated designs since this can be used for screening out the important factors from a large set of potentially active variables. A supersaturated design is a design for which there are fewer runs than effects to be estimated. Although there exists a wide study of construction methods for supersaturated designs, their analysis methods are yet in an early research stage. In this article, we propose a method for analyzing data using a correlation-based measure, named as symmetrical uncertainty. This method combines measures from the information theory field and is used as the main idea of variable selection algorithms developed in data mining. In this work, the symmetrical uncertainty is used from another viewpoint in order to determine more directly the important factors. The specific method enables us to use supersaturated designs for analyzing data of generalized linear models for a Bernoulli response. We evaluate our method by using some of the existing supersaturated designs, obtained according to methods proposed by Tang and Wu (1997 Tang , B. , Wu , C. F. J. (1997). A method for constructing supersaturated designs and its E(s ²)-optimality. Canadian Journal of Statistics 25:191–201.[Crossref], [Web of Science ®] , [Google Scholar]) as well as by Koukouvinos et al. (2008 Koukouvinos , C. , Mylona , K. , Simos , D. E. ( 2008 ). E(s ²)-optimal and minimax-optimal cyclic supersaturated designs via multi-objective simulated annealing . Journal of Statistical Planning and Inference 138 : 1639 – 1646 .[Crossref], [Web of Science ®] , [Google Scholar]). The comparison is performed by some simulating experiments and the Type I and Type II error rates are calculated. Additionally, Receiver Operating Characteristics (ROC) curves methodology is applied as an additional statistical tool for performance evaluation.     

Two-Step Residual-Based Estimation of Error Variances for Generalized Least Squares in Split-Plot Experiments

In split-plot experiments, estimation of unknown parameters by generalized least squares (GLS), as opposed to ordinary least squares (OLS), is required, owing to the existence of whole- and subplot errors. However, estimating the error variances is often necessary for GLS. Restricted maximum likelihood (REML) is an established method for estimating the error variances, and its benefits have been highlighted in many previous studies. This article proposes a new two-step residual-based approach for estimating error variances. Results of numerical simulations indicate that the proposed method performs sufficiently well to be considered as a suitable alternative to REML.     

Two Adjusted Empirical-Likelihood-Based Methods in Generalized Varying-Coefficient Partially Linear Model

An empirical likelihood method was proposed by Owen and has been extended to many semiparametric and nonparametric models with a continuous response variable. However, there has been less attention focused on the generalized regression model. This article systematically studies two adjusted empirical-likelihood-based methods in the generalized varying-coefficient partially linear models. Based on the popular profile likelihood estimation procedure, the new adjusted empirical likelihood technology for the parameter is established and the resulting statistics are shown to be asymptotically standard chi-square distributed. Further, the adjusted empirical-likelihood-based confidence regions are established, and an efficient adjusted profile empirical-likelihood-based confidence intervals/regions for any components of the parameter, which are of primary interest, is also constructed. Their asymptotic properties are also derived. Some numerical studies are carried out to illustrate the performance of the proposed inference procedures.     

Bayesian,MLE, and GMM Estimation of a Spot Rate Model

ABSTRACT

We develop Markov chain Monte Carlo algorithms for estimating the parameters of the short-term interest rate model. Using Monte Carlo experiments we compare the Bayes estimators with the maximum likelihood and generalized method of moments estimators. We estimate the model using the Japanese overnight call rate data.     

Some aspects of statistical inference on the lagrange (generalized) poisson distribution

A generalization of the Poisson distribution was defined by Consul and Jain (Ann. Math. Statist., 41, (1970)) and was obtained as a particular family of Lagrange distributions by Consul and Shenton (SIAM. J. Appl. Math., 23, (1972)). The distribution is subsequently named the generalized Poisson distribution (GPD). This GPD reduces to the Poisson distribution for ? = 0. When the data have a one-way layout structure, the asymptotically locally optimal Neyman's C(d) test is constructed and compared with the conditional test on the hypothesis H_o? = 0. Within the framework of the generalized linear models an appropriate link function is given, and the asymptotic distributions of the estimated parameters are derived.     

Analysis of accelerated life tests with weibull failure distribution via generalized linear models

Efficient industrial experiments for reliability analysis of manufactured goods may consist in subjecting the units to higher stress levels than those of the usual working conditions. This results in the so called "accelerated life tests" where, for each pre-fixed stress level, the experiment ends after the failure of a certain pre-fixed proportion of units or a certain test time is reached. The aim of this paper is to determine estimates of the mean lifetime of the units under usual working conditions from censored failure data obtained under stress conditions. This problem is approached through generalized linear modelling and related inferential techniques, considering a Weibull failure distribution and a log-linear stress-response relationship. The general framework considered has as particular cases, the Inverse Power Law model, the Eyring model, the Arrhenius model and the generalized Eyring model. In order to illustrate the proposed methodology, a numerical example is provided.     

Efficiency and Validity Analyses of Two-Stage Estimation Procedures and Derived Testing Procedures in Quantitative Linear Models with AR(1) Errors

Abstract

In a quantitative linear model with errors following a stationary Gaussian, first-order autoregressive or AR(1) process, Generalized Least Squares (GLS) on raw data and Ordinary Least Squares (OLS) on prewhitened data are efficient methods of estimation of the slope parameters when the autocorrelation parameter of the error AR(1) process, ρ, is known. In practice, ρ is generally unknown. In the so-called two-stage estimation procedures, ρ is then estimated first before using the estimate of ρ to transform the data and estimate the slope parameters by OLS on the transformed data. Different estimators of ρ have been considered in previous studies. In this article, we study nine two-stage estimation procedures for their efficiency in estimating the slope parameters. Six of them (i.e., three noniterative, three iterative) are based on three estimators of ρ that have been considered previously. Two more (i.e., one noniterative, one iterative) are based on a new estimator of ρ that we propose: it is provided by the sample autocorrelation coefficient of the OLS residuals at lag 1, denoted r(1). Lastly, REstricted Maximum Likelihood (REML) represents a different type of two-stage estimation procedure whose efficiency has not been compared to the others yet. We also study the validity of the testing procedures derived from GLS and the nine two-stage estimation procedures. Efficiency and validity are analyzed in a Monte Carlo study. Three types of explanatory variable x in a simple quantitative linear model with AR(1) errors are considered in the time domain: Case 1, x is fixed; Case 2, x is purely random; and Case 3, x follows an AR(1) process with the same autocorrelation parameter value as the error AR(1) process. In a preliminary step, the number of inadmissible estimates and the efficiency of the different estimators of ρ are compared empirically, whereas their approximate expected value in finite samples and their asymptotic variance are derived theoretically. Thereafter, the efficiency of the estimation procedures and the validity of the derived testing procedures are discussed in terms of the sample size and the magnitude and sign of ρ. The noniterative two-stage estimation procedure based on the new estimator of ρ is shown to be more efficient for moderate values of ρ at small sample sizes. With the exception of small sample sizes, REML and its derived F-test perform the best overall. The asymptotic equivalence of two-stage estimation procedures, besides REML, is observed empirically. Differences related to the nature, fixed or random (uncorrelated or autocorrelated), of the explanatory variable are also discussed.     

Bayes Prediction for a Heteroscedastic Regression Superpopulation Model Using Balanced Loss Function

We consider Prais–Houthakker heteroscedastic normal regression model having variance of the dependent variable same as square of its expectation. Bayes predictors for the regression coefficient and the mean of a finite population are derived using Zellner's balanced loss function. Bayes predictive expected losses are obtained and compared with those of classical predictors and Bayes predictors under squared error loss function to examine their loss robustness.     

Regression Properties for Asymmetric Generalized Scale Mixtures of Multivariate Gaussian Variables

Analytical properties of regression and the variance–covariance matrix of asymmetric generalized scale mixture of multivariate Gaussian variables are presented. The analysis includes an in-depth analytical investigation of the first two conditional moments of the mixing variable. Exact computable expressions for the prediction and the conditional variance are presented for the generalized hyperbolic distribution using the inversion theorem for Fourier transforms. An application to financial log returns is demonstrated via the classical Euler approximation. The methodology is illustrated by analyzing the regression of intraday log returns for CISCO against the corresponding data from S&P 500.     

Accent on Teaching Materials: Notices and Reviews of Important Teaching Materials

This article examines the probabilities of outcomes from rolling dice with the dimension 1 × 1 × r for various values of r. Experiments were conducted by school students and university students. The results of the experiments are given and the probabilities examined using a generalized linear model. Notes are also made about the value of the experiment in teaching the groups of students.