Estimation in M/Er /1 queueing systems

In this paper, the maximum likelihood estimates of the parameters for the M/E_r /1 queueing model are derived when the queue size at each departure point is observed. A numerical example is generated by simulating a finite Markov chain to illustrate the methodology for estimating the parameters with variable Erlang service time distribution. The problem of hypothesis testing and simultaneous Confidence regions of the parameter is also investigated.0     

TESTING FOR INDEPEDENCE IN MULTIVARIATE EXPONENTIAL DISTRIBUTIONS2

In this paper the work of Weier & Basu (1980) is extended to a special case of the trivariate exponential distributions and to the general k-variate case. In the trivariate case several statistics are derived including one based on the likelihood ratio approach and the locally most powerful rank statistic, and power studies are carried out. The general k-variate model is derived, and testing for independence is shown to reduce to a solved problem.     

A note on bias reduction of maximum likelihood estimates for the scalar skew t distribution

The skew t distribution is a flexible parametric family to fit data, because it includes parameters that let us regulate skewness and kurtosis. A problem with this distribution is that, for moderate sample sizes, the maximum likelihood estimator of the shape parameter is infinite with positive probability. In order to try to solve this problem, Sartori (2006) has proposed using a modified score function as an estimating equation for the shape parameter. In this note we prove that the resulting modified maximum likelihood estimator is always finite, considering the degrees of freedom as known and greater than or equal to 2.     

THE EFFECT OF SAMPLE STRUCTURE ON ANALYTICAL SURVEYS1,2

The core of this paper is a dialogue between "Harry", an experienced survey statistician and "Fred", a young mathematical statistician. They have contrasting approaches to the problem of estimating a regression relationship from a stratified sample, but after one or two red herrings are dragged out, both realize that the situation is not as simple as they had supposed. The role played by the probabilities of selection is a central issue. Estimation sampling for means, totals, and ratios is also considered, and seen to be a special case of the general analytical sampling synthesis they had already agreed upon.     

Clustering gene expression time course data using mixtures of multivariate t-distributions

Clustering gene expression time course data is an important problem in bioinformatics because understanding which genes behave similarly can lead to the discovery of important biological information. Statistically, the problem of clustering time course data is a special case of the more general problem of clustering longitudinal data. In this paper, a very general and flexible model-based technique is used to cluster longitudinal data. Mixtures of multivariate t-distributions are utilized, with a linear model for the mean and a modified Cholesky-decomposed covariance structure. Constraints are placed upon the covariance structure, leading to a novel family of mixture models, including parsimonious models. In addition to model-based clustering, these models are also used for model-based classification, i.e., semi-supervised clustering. Parameters, including the component degrees of freedom, are estimated using an expectation-maximization algorithm and two different approaches to model selection are considered. The models are applied to simulated data to illustrate their efficacy; this includes a comparison with their Gaussian analogues—the use of these Gaussian analogues with a linear model for the mean is novel in itself. Our family of multivariate t mixture models is then applied to two real gene expression time course data sets and the results are discussed. We conclude with a summary, suggestions for future work, and a discussion about constraining the degrees of freedom parameter.     

Maximum penalized likelihood estimation for skew-normal and skew-t distributions

The skew-normal and the skew-t distributions are parametric families which are currently under intense investigation since they provide a more flexible formulation compared to the classical normal and t distributions by introducing a parameter which regulates their skewness. While these families enjoy attractive formal properties from the probability viewpoint, a practical problem with their usage in applications is the possibility that the maximum likelihood estimate of the parameter which regulates skewness diverges. This situation has vanishing probability for increasing sample size, but for finite samples it occurs with non-negligible probability, and its occurrence has unpleasant effects on the inferential process. Methods for overcoming this problem have been put forward both in the classical and in the Bayesian formulation, but their applicability is restricted to simple situations. We formulate a proposal based on the idea of penalized likelihood, which has connections with some of the existing methods, but it applies more generally, including the multivariate case.     

D-optimal block designs with at most six varieties

For given positive integers v, b, and k (all of them ≥2) a block design is a k × b array of the variety labels 1,…,v with blocks as columns. For the usual one-way heterogeneity model in standard form the problem is studied of finding a D-optimal block design for estimating the variety contrasts, when no balanced block design (BBD) exists. The paper presents solutions to this problem for v≤6. The results on D-optimality are derived from a graph-theoretic context. Block designs can be considered as multigraphs, and a block design is D-optimal iff its multigraph has greatest complexity (=number of spanning trees).     

A discrete/continuous choice model on a nonconvex budget set

Decreasing block rate pricing is a nonlinear price system often used for public utility services. Residential gas services in Japan and the United Kingdom are provided under this price schedule. The discrete/continuous choice approach is used to analyze the demand under decreasing block rate pricing. However, the nonlinearity problem, which has not been examined in previous studies, arises because a consumer’s budget set (a set of affordable consumption amounts) is nonconvex, and hence, the resulting model includes highly nonlinear functions. To address this problem, we propose a feasible, efficient method of demand estimation on the nonconvex budget. The advantages of our method are as follows: (i) the construction of an Markov chain Monte Carlo algorithm with an efficient blanket based on the Hermite–Hadamard integral inequality and the power-mean inequality, (ii) the explicit consideration of the (highly nonlinear) separability condition, which often makes numerical likelihood maximization difficult, and (iii) the introduction of normal disturbance into the discrete/continuous choice model on the nonconvex budget set. The proposed method is applied to estimate the Japanese residential gas demand function and evaluate the effect of price schedule changes as a policy experiment.     

Computing the conditional stationary distribution in Markov chains of level-dependent M/G/1-type

This paper considers the computation of the conditional stationary distribution in Markov chains of level-dependent M/G/1-type, given that the level is not greater than a predefined threshold. This problem has been studied recently and a computational algorithm is proposed under the assumption that matrices representing downward jumps are nonsingular. We first show that this assumption can be eliminated in a general setting of Markov chains of level-dependent G/G/1-type. Next we develop a computational algorithm for the conditional stationary distribution in Markov chains of level-dependent M/G/1-type, by modifying the above-mentioned algorithm slightly. In principle, our algorithm is applicable to any Markov chain of level-dependent M/G/1-type, if the Markov chain is irreducible and positive-recurrent. Furthermore, as an input to the algorithm, we can set an error bound for the computed conditional distribution, which is a notable feature of our algorithm. Some numerical examples are also provided.     

ON THE ESTIMATION OF TOTAL WEIGHT IN SINGULAR SPRING BALANCE WEIGHING DESIGNS UNDER THE CO VARIANCE MATRIX OF ERRORS s2G

The problem of the estimation of the linear combination of weights, c′w, in a singular spring balance weighing design when the error structure takes the form E(ee′) =s?²G has been studied. A lower bound for the variance of the estimated linear combination of weights is obtained and a necessary and sufficient condition for this lower bound to be attained is given. The general results are applied to the case of the total of the weights. For a specified form for G, some optimum spring balance weighing designs for the estimated total weight are found.     

A new alternative to the standard F test for clustered data

The data collection process and the inherent population structure are the main causes for clustered data. The observations in a given cluster are correlated, and the magnitude of such correlation is often measured by the intra-cluster correlation coefficient. The intra-cluster correlation can lead to an inflated size of the standard F test in a linear model. In this paper, we propose a solution to this problem. Unlike previous adjustments, our method does not require estimation of the intra-class correlation, which is problematic especially when the number of clusters is small. Our simulation results show that the new method outperforms the existing methods.     

An L-statistic approach to a test of exponentiality against IFR alternatives

In this note we consider the problem of testing exponentiality against IFR alternatives. A measure of deviation from exponentiality is developed and a test statistic constructed on the basis of this measure. It is shown that the test statistic is an L-statistic. The asymptotic as well as the exact distribution of the test statistic is obtained and the test is shown to be consistent.     

Nonparametric model checking for k-out-of-n systems

It is an important problem in reliability analysis to decide whether for a given k-out-of-n system the static or the sequential k-out-of-n model is appropriate. Often components are redundantly added to a system to protect against failure of the system. If the failure of any component of the system induces a higher rate of failure of the remaining components due to increased load, the sequential k-out-of-n model is appropriate. The increase of the failure rate of the remaining components after a failure of some component implies that the effects of the component redundancy are diminished. On the other hand, if all the components have the same failure distribution and whenever a failure occurs, the remaining components are not affected, the static k-out-of-n model is adequate. In this paper, we consider nonparametric hypothesis tests to make a decision between these two models. We analyze test statistics based on the profile score process as well as test statistics based on a multivariate intensity ratio and derive their asymptotic distribution. Finally, we compare the different test statistics.     

Computing the probability of hash table/urn overflow

We analyze the probability of a random distribution of n balls into m urns of size b resulting in no overflows. This solves the computational problem associated with a classical combinatorial extreme-value distribution. The problem arose during the analysis of a technique, called perfect hashing, for organizing data in computer files. The results and techniques presented can be used to solve several problems in the analysis of hashing techniques     

Moments of truncated normal/independent distributions

In this work we have considered the problem of finding the moments of a doubly truncated member of the class of normal/independent distributions. We obtained a general result and then use it to derive the moments in the case of doubly truncated versions of Pearson type VII distribution, slash distribution, contaminated normal distribution, double exponential distribution and variance gamma distribution. We also give an application of some actuarial data.     

A LINEAR MODEL WITH ERRORS LACKING A VARIANCE1

The problem discussed is that of estimating β= (β¹, …, β_k) in the model Y=βX +ε when X has a specified multivariate distribution and the error ε does not necessarily have a finite second moment, for example, ε symmetric stable. We construct a moment estimator based on the empirical characteristic function and establish asymptotic unbiassedness and normality. Most of the paper is concerned with the case when X is normal. Forms of the suggested estimator are given in (2.5), (4.6) and (5.5).     

Quality/Value Relationship for Imperfect Information in the Umbrella Problem

The so-called umbrella problem (or cost-loss ratio situation), in which an individual must decide whether to take an umbrella in the face of uncertainty concerning whether it will rain today, is sometimes used as a textbook example of decision making under uncertainty. This problem is extended to provide a simple demonstration of the way in which the economic value of imperfect information changes as its quality increases.     

The Slippage Problem for the Circular Normal Distribution

The slippage problem occurs when an unspecified observation in a given random sample is from a distribution other than that for all the remaining observations. This paper studies the problem in terms of the 'slip' in the mean direction of a circular normal distribution. The slippage problem is first treated as a multiple decision problem with a prior which is invariant under the permutations of the hypotheses. The probabilities of accepting the various hypotheses for the Bayes rule with respect to this prior are explicitly obtained. The likelihood ratio tests for this slippage problem, for the cases when the mean directions are both known and unknown, are shown to be easily computable. The tests are illustrated through two well-known datasets. The performances of a range of tests are compared using extensive simulation.     

E-optimal design in irregular BIBD settings

When the necessary conditions for a BIBD are satisfied, but no BIBD exists, there is no simple answer for the optimal design problem. This paper identifies the E-optimal information matrices for any such irregular BIBD setting when the number of treatments is no larger than 100. A- and D-optimal designs are typically not E-optimal. An E-optimal design for 15 treatments in 21 blocks of size 5 is found.     

Testing for Heteroscedasticity and/or Autocorrelation in Longitudinal Mixed Effect Nonlinear Models with AR(1) Errors

The effect of influental observation son the parameter estimates of ordinary least squares regression models has received considerable a t t e n t i o n fn the last decade. However, very little attention has been given to the problem of influential observation sinthea naysis of variace . The purpose of this paper is to show by way of examples that in fluential observations can alter the conclusions of tests of hypotheses in the analysis of variance . Regression diagno stics for identifying both extreme points and out liers can be used toreveal potential data and design problems.