Construction of some new families of nested orthogonal arrays

Nested orthogonal arrays have been used in the design of an experimental setup consisting of two experiments, the expensive one of higher accuracy being nested in a larger and relatively less expensive one of lower accuracy. In this paper, we provide new methods for constructing two types of nested orthogonal arrays.     

An Exact Algorithm for Projection onto a Polyhedral Cone

Projection is an operation widely used in restricted statistical inferences. A polyhedral cone restriction includes many interesting problems in linear regression, order restricted inferences etc. This paper proposes an exact algorithm for the projection of a vector onto a polyhedral cone, and presents an application to second order polynomial regression subject to a non-negative, isotonic restriction.     

Strong orthogonal decompositions and non‐linear impulse response functions for infinite‐variance processes

The author proves that Wold‐type decompositions with strong orthogonal prediction innovations exist in smooth, reflexive Banach spaces of discrete time processes if and only if the projection operator generating the innovations satisfies the property of iterations. His theory includes as special cases all previous Wold‐type decompositions of discrete time processes, completely characterizes when non‐linear heavy‐tailed processes obtain a strong‐orthogonal moving average representation, and easily promotes a theory of non‐linear impulse response functions for infinite‐variance processes. The author exemplifies his theory by developing a non‐linear impulse response function for smooth transition threshold processes, and discusses how to test decomposition innovations for strong orthogonality and whether the proposed model represents the best predictor. A data set on currency exchange rates allows him to illustrate his methodology.     

Bounds for moments through a general orthogonal expansion in a pre-hilbert space-I

Bounds are obtained for the product moments of an arbitrary finite number of ordered random variables. These bounds are obtained with the help of a representation of an arbitrary function in terms of a complete orthonormal system in a pre-Hilbert space of square integrable functions defined in a k-dimensional unit cube.     

Bounds for moments through a general orthogonal expansion in a pre-hilbert space-II

In an earlier article [Canad. J. Statist., Vol, 3, No. 1, 1975, 13–34] bounds are obtained for the product moments of an arbitrary finite number of ordered random variables. These bounds are obtained with the help of a representation of an arbitrary function in terms of a complete orthonormal system in a pre-Hilbert space of square integrable functions defined in a k-dimensional unit cube. These results are extended to symmetric cases in this article.     

Evaluation of algorithms for generating Dirichlet random vectors

In this article, we describe various well-known Dirichlet generation algorithms and evaluate their performance in terms of the following criteria: (i) computer generation time, (ii) sensitivity, and (iii) goodness of fit. In addition, we examine in particular an algorithm based on transformation of beta variates and provide three useful guidelines so as to reduce its computer generation time. Simulation results show that the proposed algorithm significantly outperforms other approaches in terms of computer generation time, except in cases when all (or most) shape parameters are close to zero.     

General orthogonal designs for parameter estimation of ANOVA models under weighted sum-to-zero constraints

ABSTRACT

Traditional studies on optimal designs for ANOVA parameter estimation are based on the framework of equal probabilities of appearance for each factor's levels. However, this premise does not hold in a variety of experimental problems, and it is of theoretical and practical interest to investigate optimal designs for parameters with unequal appearing odds. In this paper, we propose a general orthogonal design via matrix image, in which all columns’ matrix images are orthogonal with each other. Our main results show that such designs have A- and E-optimalities on the estimation of ANOVA parameters which have unequal appearing odds. In addition, we develop two simple methods to construct the proposed designs. The optimality of the design is also validated by a simulation study.     

Robust algorithms for economic designing of a nonparametric control chart for abrupt shift in location

The existing statistical process control procedures typically rely on the fundamental assumption of a parametric distribution of the quality characteristic. However, when there is a lack of knowledge about the underlying distribution (as full knowledge is not available in practice), the performance of these parametric charts is very likely to be heavily degraded. Motivated by this problem, a one-sided nonparametric monitoring procedure using the single sample sign statistic is proposed for detecting a shift in the location parameter of a continuous distribution. An economic model of the control chart is developed to optimize the sample size, sampling interval, and control limits. Three data-dependent estimation approaches for the unknown parameter are evaluated and discussed. Simulation results exhibit that our proposed procedure generally performs well under a great variety of continuous distributions and hence it is recommended as an alternative scheme especially when the knowledge of the underlying distribution is imperfect. Furthermore, beneficial recommendations of estimation approach selection are provided for practical implementation of the control chart.     

Comparison of various machine learning algorithms for estimating generalized propensity score

In this paper, we conducted a simulation study to evaluate the performance of four algorithms: multinomial logistic regression (MLR), bagging (BAG), random forest (RF), and gradient boosting (GB), for estimating generalized propensity score (GPS). Similar to the propensity score (PS), the ultimate goal of using GPS is to estimate unbiased average treatment effects (ATEs) in observational studies. We used the GPS estimates computed from these four algorithms with the generalized doubly robust (GDR) estimator to estimate ATEs in observational studies. We evaluated these ATE estimates in terms of bias and mean squared error (MSE). Simulation results show that overall, the GB algorithm produced the best ATE estimates based on these evaluation criteria. Thus, we recommend using the GB algorithm for estimating GPS in practice.     

Roundoff error for floating point representation of real data

In the present paper we define aprobabilistic model for the study ofthe relative roundoff error in the floating point representation ofreal data. In particular, starting from the assumption that real data is uniformly distributed in the interval (0,1), each floating point is assigned aprobability value corresponding to the proportion of real data assigned this floating point value. For each real number x～ U(0,1), we compute the relative roundoff error in the floating point representation and we calculate i t s expected value and variance.     

Confidence interval estimation for the population coefficient of variation using ranked set sampling: a simulation study

In this paper, an evaluation of the performance of several confidence interval estimators of the population coefficient of variation (τ) using ranked set sampling compared to simple random sampling is performed. Two performance measures are used to assess the confidence intervals for τ, namely: width and coverage probabilities. Simulated data were generated from normal, log-normal, skew normal, Gamma, and Weibull distributions with specified population parameters so that the same values of τ are obtained for each distribution, with sample sizes n=15, 20, 25, 50, 100. A real data example representing birth weight of 189 newborns is used for illustration and performance comparison.     

On an application of the diffusion approximation of the generalized birth and death process for stochastic population projection

Assuming that both birth and death rates are density and time dependent, a diffusion approximation of the generalized birth and death process has been considered in this paper to obtain a suitable stochastic population model describing the population size and its moments. A simple method of estimating the parameters of the model Is discussed. The predictions of the expected size of the population, and the variance are made and compared with the corresponding census figures as well as with another deterministic projection series made for the corresponding period.     

Crossover designs based on type I orthogonal arrays for a self and simple mixed carryover effects model with correlated errors

We investigate the performance of crossover designs based on type I orthogonal arrays for a self and simple mixed carryover effects model in the presence of correlated errors. Assuming that between-subject errors are independent while within-subject errors behave according to the stationary first-order autoregressive and moving average processes, analytical optimality results for 3-period designs are established and, as an illustration, numerical details for a number of 4-period cases are tabulated.     

An asymptotic approximation for the permanent of a doubly stochastic matrix

A determinantal approximation is obtained for the permanent of a doubly stochastic matrix. For moderate-deviation matrix sequences, the asymptotic relative error is of order O(n^?1).     

Multivariate analysis of variance for functional data

Functional data are being observed frequently in many scientific fields, and therefore most of the standard statistical methods are being adapted for functional data. The multivariate analysis of variance problem for functional data is considered. It seems to be of practical interest similarly as the one-way analysis of variance for such data. For the MANOVA problem for multivariate functional data, we propose permutation tests based on a basis function representation and tests based on random projections. Their performance is examined in comprehensive simulation studies, which provide an idea of the size control and power of the tests and identify differences between them. The simulation experiments are based on artificial data and real labeled multivariate time series data found in the literature. The results suggest that the studied testing procedures can detect small differences between vectors of curves even with small sample sizes. Illustrative real data examples of the use of the proposed testing procedures in practice are also presented.     

Relative efficienctes of sampling plans for selecting a small number of units from a rectangular region

We examine the efficiency of several sampling plans for use in certain agricultural, ecological and environmental studies. One concern for such studies is that plots that arephysically close might be more similar than distant plots. We considered sampling plansthat are designed to generate samples that represent the entire population while avoidingthe selection of units that provide essentially redundant information. All plans havethe property that they avoid the simultaneous selection of units that are, in some sense,neighboring units. By means of a simulation study, the efficiency of these plans iscompared to simple random Aampling Factors that influence the relative efficiencies areexamined. This is done for a number of different populations, representing variouspossible patterns for a response variable.     

Posterior moments of scalar functions of a covariance matrix

Given multivariate normal data and a certain spherically invariant prior distribution on the covariance matrix, it is desired to estimate the moments of the posterior marginal distributions of some scalar functions of the covariance matrix by importance sampling. To this end a family of distributions is defined on the group of orthogonal matrices and a procedure is proposed for selecting one of these distributions for use as a weighting distribution in the importance sampling process. In an example estimates are calculated for the posterior mean and variance of each element in the covariance matrix expressed in the original coordinates, for the posterior mean of each element in the correlation matrix expressed in the original coordinates, and for the posterior mean of each element in the covariance matrix expressed in the coordinates of the principal variables.