General Multivariate Weibull Processes

Two multivariate stationary processes with general multivariate Weibull marginals are developed and studied. The joint distribution of the two adjacent events in the processes and the distributions of the finite sample minima as well as the geometric minima are derived. The characterization properties of these two processes are also proved.     

Bernoulli Processes with Non-Positive Correlations

We deal with sequences of dependent Bernoulli variables with non positive correlations. Some special models for 1-dependent and 2-dependent 0–1 valued variables are analyzed, namely Bernoulli variables obtained by clipping a linear combination of iid variables. Formulas describing dependence of their correlation function on the clipping parameters are derived. The lower bound for the sum of autocorrelations of Bernoulli variables obtained by clipping a Gaussian process is provided.     

A Distribution-Free Control Chart Based on Order Statistics

In this article, we introduce a new distribution-free Shewhart-type control chart that takes into account the location of a single order statistic of the test sample (such as the median) as well as the number of observations in that test sample that lie between the control limits. Exact formulae for the alarm rate, the run length distribution, and the average run length (ARL) are all derived. A key advantage of the chart is that, due to its nonparametric nature, the false alarm rate and in-control run length distribution are the same for all continuous process distributions, and so will be naturally robust. Tables are provided for the implementation of the chart for some typical ARL values and false alarm rates. The empirical study carried out reveals that the new chart is preferable from a robustness point of view in comparison to a classical Shewhart-type chart and also the nonparametric chart of Chakraborti et al. (2004 Chakraborti , S. , van der Laan , P. , van de Wiel , M. A. ( 2004 ). A class of distribution-free control charts . J. Roy. Statist. Soc. Ser. C-Appl. Statist. 53 ( 3 ): 443 – 462 .[Web of Science ®] , [Google Scholar]).     

The effects of correlation among observations on the consistency property of sample variance

It is well known that even when the sample observations are correlated and not normal the sample variance, S² converges in probability to E(S²). But the required sample size for S² to be a consistent estimator of E(S²) is an open question. Some light is shed on this question in this paper. In particular the relation between the rate of convergence and the correlation property of the observations is explored. It is shown that the retardation to the rate of convergence is not appreciable if the correlation is moderate but it can be severe for extreme correlations.     

L1 for the simple linear regression model

By modifying the direct method to solve the overdetermined linear system we are able to present an algorithm for L₁ estimation which appears to be superior computationally to any other known algorithm for the simple linear regression problem.     

Bathtub distributions: a review

Bathtub distributions are characterized by bathtub failure rate functions . These are possibly more realisitic models than the monotone failure rate models . A systematic account of such distributions is not available and this review aims to give such an account . We give some easily verifiable conditions to check the bathtub property of a distribution along with methods to construct such distributions . We also discuss some stochastic and reliablity mechanisms which lead to bathtub distributions. These include mixtures ( stochastic failure rate models ) , series system , stochastic differential equation models and so on. We also review inference on bathtub distributions. The paper concludes with a rather exhaustive list of bathtub distributions.     

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.     

Approximating M/G/1 Waiting Time Tail Probabilities

Abstract

We propose a new approximation formula for the waiting time tail probability of the M/G/1 queue with FIFO discipline and unlimited waiting space. The aim is to address the difficulty of obtaining good estimates when the tail probability has non-exponential asymptotics. We show that the waiting time tail probability can be expressed in terms of the waiting time tail probability of a notional M/G/1 queue with truncated service time distribution plus the tail probability of an extreme order statistic. The Cramér–Lundberg approximation is applied to approximate the tail probability of the notional queue. In essence, our technique extends the applicability of the Cramér–Lundberg approximation to cases where the standard Lundberg condition does not hold. We propose a simple moment-based technique for estimating the parameters of the approximation; numerical results demonstrate that our approximation can yield very good estimates over the whole range of the argument.     

Workload Process,Waiting Times,and Sojourn Times in a Discrete Time MMAP[K]/SM[K]/1/FCFS Queue

Abstract

In this paper, we study the total workload process and waiting times in a queueing system with multiple types of customers and a first-come-first-served service discipline. An M/G/1 type Markov chain, which is closely related to the total workload in the queueing system, is constructed. A method is developed for computing the steady state distribution of that Markov chain. Using that steady state distribution, the distributions of total workload, batch waiting times, and waiting times of individual types of customers are obtained. Compared to the GI/M/1 and QBD approaches for waiting times and sojourn times in discrete time queues, the dimension of the matrix blocks involved in the M/G/1 approach can be significantly smaller.     

A Retrial Queue with a Constant Retrial Rate,Server Downs and Impatient Customers

ABSTRACT

In this paper, we consider a retrial queueing system consisting of a waiting line of infinite capacity in front of a single server subject to breakdowns. A customer upon arrival may join the queue (waiting line) or go to the retrial orbit (another queue) to retry for service after a random time. Only the customer at the head of the retrial orbit is allowed to retry for service. Upon retrial, the customer enters the service if the server is idle; otherwise, it may go back to the retrial orbit or leave the system (become impatient). All the interarrival times, service times, server up times, server down times and retrial times are exponential, and all the necessary independence conditions in these variables are assumed. For this system, we provide sufficient conditions under which, for any given number of customers in the orbit, the stationary probability of the number of customers in the waiting line decays geometrically. We also provide explicitly an expression for the decay parameter.