Two-queue polling systems with switching policy based on the queue that is not being served

We study a system of two non-identical and separate M/M/1/? queues with capacities (buffers) C₁ < ∞ and C₂ = ∞, respectively, served by a single server that alternates between the queues. The server’s switching policy is threshold-based, and, in contrast to other threshold models, is determined by the state of the queue that is not being served. That is, when neither queue is empty while the server attends Q_i (i = 1, 2), the server switches to the other queue as soon as the latter reaches its threshold. When a served queue becomes empty we consider two switching scenarios: (i) Work-Conserving, and (ii) Non-Work-Conserving. We analyze the two scenarios using Matrix Geometric methods and obtain explicitly the rate matrix R, where its entries are given in terms of the roots of the determinants of two underlying matrices. Numerical examples are presented and extreme cases are investigated.     

Optimal foldover designs

An obvious strategy for obtaining a Doptimal foldover design for p factors at two levels each in 2N runs is to fold a Doptimal main effects plan. We show that this strategy works except when N = 4t + 2 and s is even In that case there are two different classes of D-optimal main effects plans with N runs that have the same determinant. However folding them gives two different values foi the D-optimality criteiion One set of designs is D-optimal The other is not.     

Strong previsions of random elements

LetC be a class of arbitrary real random elements andP an extended real valued function onC. Two definitions of coherence forP are compared. Both definitions reduce to the classical de Finetti's one whenC includes bounded random elements only. One of the two definitions (called strong coherence) is investigated, and some criteria for checking it are provided. Moreover, conditions are given for the integral representation of a coherentP, possibly with respect to a δ-additive probability. Finally, the two definitions and the integral representation theorems are extended to the case whereC is a class of random elements taking values in a given Banach space.     

Unbiased Estimation of P(X > Y) for Exponential Populations Using Order Statistics with Application in Ranked Set Sampling

This paper addresses the problem of unbiased estimation of P[X > Y] = θ for two independent exponentially distributed random variables X and Y. We present (unique) unbiased estimator of θ based on a single pair of order statistics obtained from two independent random samples from the two populations. We also indicate how this estimator can be utilized to obtain unbiased estimators of θ when only a few selected order statistics are available from the two random samples as well as when the samples are selected by an alternative procedure known as ranked set sampling. It is proved that for ranked set samples of size two, the proposed estimator is uniformly better than the conventional non-parametric unbiased estimator and further, a modified ranked set sampling procedure provides an unbiased estimator even better than the proposed estimator.     

TESTING HYPOTHESES ABOUT THE COMMON MEAN OF TWO NORMAL POPULATIONS

Three tests are considered concerning the common mean of two normal populations: (1) an F test based on a sample from one population, (2) a test based on the addition of the F statistics from independent samples from two popultions (proposed), and (3) a test based on the maximum of the F statistics from two independent samples from two populations. A condition under which test (2) is locally more powerful than test (1) is given. As the test statistic in test (2) does not follow a standard distribution, a formula for approximating the observed significance level is provided. A simulation study is used to compare the power of these tests.     

Local comparison of empirical distributions via nonparametric regression

Given two independent samples of size n and m drawn from univariate distributions with unknown densities f and g, respectively, we are interested in identifying subintervals where the two empirical densities deviate significantly from each other. The solution is built by turning the nonparametric density comparison problem into a comparison of two regression curves. Each regression curve is created by binning the original observations into many small size bins, followed by a suitable form of root transformation to the binned data counts. Turned as a regression comparison problem, several nonparametric regression procedures for detection of sparse signals can be applied. Both multiple testing and model selection methods are explored. Furthermore, an approach for estimating larger connected regions where the two empirical densities are significantly different is also derived, based on a scale-space representation. The proposed methods are applied on simulated examples as well as real-life data from biology.     

On binomial and circular binomial distributions of orderk forl-overlapping success runs of lengthk

The number ofl-overlapping success runs of lengthk inn trials, which was introduced and studied recently, is presently reconsidered in the Bernoulli case and two exact formulas are derived for its probability distribution function in terms of multinomial and binomial coefficients respectively. A recurrence relation concerning this distribution, as well as its mean, is also obtained. Furthermore, the number ofl-overlapping success runs of lengthk inn Bernoulli trials arranged on a circle is presently considered for the first time and its probability distribution function and mean are derived. Finally, the latter distribution is related to the first, two open problems regarding limiting distributions are stated, and numerical illustrations are given in two tables. All results are new and they unify and extend several results of various authors on binomial and circular binomial distributions of orderk.     

A single-featured EWMA-X control chart for detecting shifts in process mean and standard deviation

The combined EWMA-X chart is a commonly used tool for monitoring both large and small process shifts. However, this chart requires calculating and monitoring two statistics along with two sets of control limits. Thus, this study develops a single-featured EWMA-X (called SFEWMA-X) control chart which has the ability to simultaneously monitor both large and small process shifts using only one set of statistic and control limits. The proposed SFEWMA-X chart is further extended to monitoring the shifts in process standard deviation. A set of simulated data are used to demonstrate the proposed chart's superior performance in terms of average run length compared with that of the traditional charts. The experimental examples also show that the SFEWMA-X chart is neater and easier to visually interpret than the original EWMA-X chart.     

A note on the ratio of normal and Laplace random variables

The normal and Laplace are the two earliest known continuous distributions in statistics and the two most popular models for analyzing symmetric data. In this note, the exact distribution of the ratio | X / Y | is derived when X and Y are respectively normal and Laplace random variables distributed independently of each other. A MAPLE program is provided for computing the associated percentage points. An application of the derived distribution is provided to a discriminant problem.     

Tests of regression coefficients under ridge regression models

This paper investigates two “non-exact” t-type tests, t( k₂) and t(k₂), of the individual coefficients of a linear regression model, based on two ordinary ridge estimators. The reported results are built on a simulation study covering 84 different models. For models with large standard errors, the ridge-based t-tests have correct levels with considerable gain in powers over those of the least squares t-test, t(0). For models with small standard errors, t(k₁) is found to be liberal and is not safe to use while, t(k₂) is found to slightly exceed the nominal level in few cases. When tie two ridge tests art: not winners, the results indicate that they don't loose much against t(0).     

Matching a correlation coefficient by a Gaussian copula

A Gaussian copula is widely used to define correlated random variables. To obtain a prescribed Pearson correlation coefficient of ρ_x between two random variables with given marginal distributions, the correlation coefficient ρ_z between two standard normal variables in the copula must take a specific value which satisfies an integral equation that links ρ_x to ρ_z. In a few cases, this equation has an explicit solution, but in other cases it must be solved numerically. This paper attempts to address this issue. If two continuous random variables are involved, the marginal transformation is approximated by a weighted sum of Hermite polynomials; via Mehler’s formula, a polynomial of ρ_z is derived to approximate the function relationship between ρ_x and ρ_z. If a discrete variable is involved, the marginal transformation is decomposed into piecewise continuous ones, and ρ_x is expressed as a polynomial of ρ_z by Taylor expansion. For a given ρ_x, ρ_z can be efficiently determined by solving a polynomial equation.     

Sampling strategy for optimal classification into one of two correlated normal populations

A unit ω is to be classified into one of two correlated homoskedastic normal populations by linear discriminant function known as W classification statistic [T.W. Anderson, An asymptotic expansion of the distribution of studentized classification statistic, Ann. Statist. 1 (1973), pp. 964–972; T.W. Anderson, An Introduction to Multivariate Statistical Analysis, 2nd edn, Wiley, New York, 1984; G.J. Mclachlan, Discriminant Analysis and Statistical Pattern Recognition, John Wiley and Sons, New York, 1992]. The two populations studied here are two different states of the same population, like two different states of a disease where the population is the population of diseased patient. When a sample unit is observed in both the states (populations), the observations made on it (which form a pair) become correlated. A training sample is unbalanced when not all sample units are observed in both the states. Paired and also unbalanced samples are natural in studies related to correlated populations. S. Bandyopadhyay and S. Bandyopadhyay [Choosing better training sample for classifying an individual into one of two correlated normal populations, Calcutta Statist. Assoc. Bull. 54(215–216) (2003), pp. 167–180] studied the effect of unbalanced training sample structure on the performance of W statistics in the univariate correlated normal set-up for finding optimal sampling strategy for a better classification rate. In this study, the results are extended to the multivariate case with discussion on application in real scenario.     

Quasi-Likelihood Approach to Bioavailability and Bioequivalence Analysis

The outcomes AUCT (area-under-curve from time zero to time t) of n individuals randomized to one of two groups TR or RT, where the group name denotes the order in which the subjects receive a test formulation (T) or a reference formulation (R), are used to assess average bioequivalence for the two formulations. The classical method is the mixed model, for example, proc mixed or proc glm with random statement in SAS can be used to analyze this type of data. This is equivalent to the marginal likelihood approach in a normal–normal model. There are some limitations for this approach. It is not appropriate if the random effect is not normally distributed. In this article, we introduce a hierarchical quasi-likelihood approach. Instead of assuming the random effect is normal, we make assumptions only about the mean and the variance function of the random effect. Our method is flexible to model the random effect. Since we can estimate the random effect for each individual, we can check the adequacy of the distribution assumption about the random effect. This method can also be used to handle high-dimensional crossover data. Simulation studies are conducted to check the finite sample performance of the method under various conditions and two real data examples are used for illustration.     

Admissibility intervals for linear correlation coefficients

The correlation coefficient is widely used to quantify the degree of association between two quantitative variables. By resorting to the geometric representation of the linear correlation coefficient, it is possible to calculate the upper and lower bounds of the correlation coefficient between two variables x ₁,x ₂ when the correlation coefficients with a third variable x ₃ are available. Implications in observational studies, where x ₃ could be a proxy of a target variable x ₂, whose direct measurement is too expensive or impractical, are discussed.     

THE DIFFERENCE BETWEEN TWO POISSON EXPECTATIONS1

A number of statistics have been suggested for testing the difference between two Poisson expectations. On the basis of size and power calculations a new statistic, v= (2x+ 3/4)1/2 - (2y + 3/4)1/2, where x and y are the two observed Poisson values, is recommended for testing such differences. This statistic is shown to be similar in performance to the familiar u=(x-y)f(x+y)1/2, except in the tails of the distribution where it is superior.     

Sample size calculation for an agreement study

It is often necessary to compare two measurement methods in medicine and other experimental sciences. This problem covers a broad range of data. Many authors have explored ways of assessing the agreement of two sets of measurements. However, there has been relatively little attention to the problem of determining sample size for designing an agreement study. In this paper, a method using the interval approach for concordance is proposed to calculate sample size in conducting an agreement study. The philosophy behind this is that the concordance is satisfied when no more than the pre‐specified k discordances are found for a reasonable large sample size n since it is much easier to define a discordance pair. The goal here is to find such a reasonable large sample size n. The sample size calculation is based on two rates: the discordance rate and tolerance probability, which in turn can be used to quantify an agreement study. The proposed approach is demonstrated through a real data set. Copyright © 2009 John Wiley & Sons, Ltd.     

Generalized run tests for statistical process control

ABSTRACT

In a sequence of elements, a run is defined as a maximal subsequence of like elements. The number of runs or the length of the longest run has been widely used to test the randomness of an ordered sequence. Based on two different sampling methods and two types of test statistics used, run tests can be classified into one of four cases. Numerous researchers have derived the probability distributions in many different ways, treating each case separately. In the paper, we propose a unified approach which is based on recurrence arguments of two mutually exclusive sub-sequences. We also consider the sequence of nominal data that has more than two classes. Thus, the traditional run tests for a binary sequence are special cases of our generalized run tests. We finally show that the generalized run tests can be applied to many quality management areas, such as testing changes in process variation, developing non-parametric multivariate control charts, and comparing the shapes and locations of more than two process distributions.     

Tests on the Equality of Two Multiple Correlation Coefficients

Asymptotic tests are suggested for testing the equality of two multiple correlation coefficients calculated from a single sample from a multivariate normal distribution. An F test is possible only when the two dependent variables coincide and one set of independent variables is a subset of the second set. Tests are compared by simulation for situations in which the F test is inapplicable. Special attention is paid to cases in which asymptotic normality of the test statistics does not hold.     

Reliability analysis of systems with components having two dependent subcomponents

In this article, a system that consists of n independent components each having two dependent subcomponents (A_i, B_i), i = 1, …, n is considered. The system is assumed to compose of components that have two correlated subcomponents (A_i, B_i), and functions iff both systems of subcomponents A₁, A₂, …, A_n and B₁, B₂, …, B_n work under certain structural rules. The expressions for reliability and mean time to failure of such systems are obtained. A sufficient condition to compare two systems of bivariate components in terms of stochastic ordering is also presented.