The Multi-Sample Block-Scalar Sphericity Test: Exact and Near-Exact Distributions for Its Likelihood Ratio Test Statistic

In this article the authors show how by adequately decomposing the null hypothesis of the multi-sample block-scalar sphericity test it is possible to obtain the likelihood ratio test statistic as well as a different look over its exact distribution. This enables the construction of well-performing near-exact approximations for the distribution of the test statistic, whose exact distribution is quite elaborate and non-manageable. The near-exact distributions obtained are manageable and perform much better than the available asymptotic distributions, even for small sample sizes, and they show a good asymptotic behavior for increasing sample sizes as well as for increasing number of variables and/or populations involved.     

Dependence bias in conventional p-charts and its correction with an approximate lot quality distribution

In this study, the variance of mean fraction nonconformance under the exact lot quality distribution with first order dependence is examined. The extent of the validity of binomial assumption with respect to the exact distribution in constructing the p-charts for dependent production processes is investigated and it is shown that neglect of dependence for small lot sizes may cause bias in the width of the area between the control limits. An approximate lot quality distribution with great versatility is proposed and is used for correcting the bias due to dependence in p-charts.     

Asymptotic approximations for the distributions of the <Emphasis Type="Italic">K</Emphasis>ø-divergence goodness-of-fit statistics

Kø-divergence’s statistic family for goodness-of-fit, under the null hypothesis, has an asymptotic chi-squared distribution; however, for small samples, the chi-squared approximation in some cases does not well agree with the exact distribution. In this paper, a closer approximation to the exact distribution is obtained by extracting the ø-dependent second order component from the distribution. Moreover, numerical results are presented for moderate sample sizes with moderate number of cells.     

Small sample distribution for a nonparametric test for trend

The exact distribution of a nonparametric test statistic for ordered alternatives, the _rank ² statistic, is computed for small sample sizes. The exact distribution is compared to an approximation.     

Tests for symmetry of two-piece normal distribution

In this paper, classical optimum tests for symmetry of two-piece normal distribution is derived. Uniformly most powerful one-sided test for the skewness parameter is obtained when the location and scale parameters are known and is compared with sequential probability ratio test. An ad-hoc test for symmetry and likelihood ratio test for symmetry for large samples, can be found in literature for this distribution. But in this paper, we derive exact likelihood ratio test for symmetry, when location parameter is known. The exact power of the test is evaluated for different sample sizes.     

On two asymptotic normal distributions for the generalized wilks lambda statistic

In this paper we.present a Normal asymptotic distribution for the logarithm of the generalized Wilks Lambda statistic based on an asymptotic distribution for the determinant of a Wishart matrix. This distribution is obtained through the combined use of Taylor expansions of random variables whose exponentials have chi-square distributions and the Lindeberg-Feller version of the Central Limit Theorem, Another asymptotic Normal distribution for the logarithm of the generalized Wilks Lambda statistic for the case when at most one of the sets has an odd number of variables is derived directly from the exact distribution. Both distributions are non-degenerate and non-singular. The first Normal distribution compares favorably with other known approximations and asymptotic distributions namely for large numbers of variables and small sample sizes, while the second Normal distribution, which has a more restricted application, compares in most cases highly favorably with other known asymptotic distributions and approximations. Finally, a method to compute approximate quantiles which lay very close and converge steadily to the exact ones is presented.     

Sample size determination for 2k-r experiments with a binomial response

This paper provides closed form expressions for the sample size for two-level factorial experiments when the response is the number of defectives. The sample sizes are obtained by approximating the two-sided test for no effect through tests for the mean of a normal distribution, and borrowing the classical sample size solution for that problem. The proposals are appraised relative to the exact sample sizes computed numerically, without appealing to any approximation to the binomial distribution, and the use of the sample size tables provided is illustrated through an example.     

Exact and Approximate Distributions of the Maximum Likelihood Estimate in a Simple Factor Analysis Model

We study a factor analysis model with two normally distributed observations and one factor. In the case when the errors have equal variance, the maximum likelihood estimate of the factor loading is given in closed form. Exact and approximate distributions of the maximum likelihood estimate are considered. The exact distribution function is given in a complex form that involves the incomplete Beta function. Approximations to the distribution function are given for the cases of large sample sizes and small error variances. The accuracy of the approximations is discussed     

Fast computation of the exact null distribution of Spearman's ρ and Page's L statistic for samples with and without ties

We present a new algorithm for computing the exact null distribution of the Spearman rank correlation statistic ρ, which also works in the case of ties. The algorithm is based on symmetries in the representation of the probability generating function as a permanent with monomial entries. We present new critical values for sample sizes 19⩽n⩽22. Finally, we show how to derive the exact null distribution of Page's L statistic from the null distribution of ρ.     

Order statistics from discrete distributions

Several recurrence relations and identities available for single and product moments of order1 statistics in a sample size n from an arbitrary continuous distribution are extended for the discrete case,, Making use of these recurrence relations it is shown that it is sufficient to evaluate just two single moments and (n-l)/2 product moments when n is odd and two single moments and {n-2)/2 product moments when n is even, in order to evaluate the first, second and product moments of order statistics in a sample of size n drawn from an arbitrary discrete distribution, given these moments in samples of sizes n-1 and less.. A series representation for the product moments of order statistics is derived.. Besides enabling us to obtain an exact and explicit expression for the product moments of order statistics from the geometric distribution, it. makes the computation of the product moments of order statistics from other discrete distributions easy too.     

Corrected asymptotic distribution of statistics based on the multinomial law

Investigators and epidemiologists often use statistics based on the parameters of a multinomial distribution. Two main approaches have been developed to assess the inferences of these statistics. The first one uses asymptotic formulae which are valid for large sample sizes. The second one computes the exact distribution, which performs quite well for small samples. They present some limitations for sample sizes N neither large enough to satisfy the assumption of asymptotic normality nor small enough to allow us to generate the exact distribution. We analytically computed the 1/N corrections of the asymptotic distribution for any statistics based on a multinomial law. We applied these results to the kappa statistic in 2×2 and 3×3 tables. We also compared the coverage probability obtained with the asymptotic and the corrected distributions under various hypothetical configurations of sample size and theoretical proportions. With this method, the estimate of the mean and the variance were highly improved as well as the 2.5 and the 97.5 percentiles of the distribution, allowing us to go down to sample sizes around 20, for data sets not too asymmetrical. The order of the difference between the exact and the corrected values was 1/N² for the mean and 1/N³ for the variance.     

ROBUSTNESS PROPERTIES OF LOGNORMAL CONFIDENCE INTERVALS FOR LOGNORMAL AND GAMMA DISTRIBUTED DATA

ABSTRACT

The performances of six confidence intervals for estimating the arithmetic mean of a lognormal distribution are compared using simulated data. The first interval considered is based on an exact method and is recommended in U.S. EPA guidance documents for calculating upper confidence limits for contamination data. Two intervals are based on asymptotic properties due to the Central Limit Theorem, and the other three are based on transformations and maximum likelihood estimation. The effects of departures from lognormality on the performance of these intervals are also investigated. The gamma distribution is considered to represent departures from the lognormal distribution. The average width and coverage of each confidence interval is reported for varying mean, variance, and sample size. In the lognormal case, the exact interval gives good coverage, but for small sample sizes and large variances the confidence intervals are too wide. In these cases, an approximation that incorporates sampling variability of the sample variance tends to perform better. When the underlying distribution is a gamma distribution, the intervals based upon the Central Limit Theorem tend to perform better than those based upon lognormal assumptions.     

Sample size determinations when two binomial proportions are very small

Sample size determination for testing the hypothesis of equality of proportions with a specified type I and type I1 error probabilitiesis of ten based on normal approximation to the binomial distribution. When the proportionsinvolved are very small, the exact distribution of the test statistic may not follow the assumed distribution. Consequently, the sample size determined by the test statistic may not result in the sespecifiederror probabilities. In this paper the author proposes a square root formula and compares it with several existing sample size approximation methods. It is found that with small proportion (p≦.01) the squar eroot formula provides the closest approximation to the exact sample sizes which attain a specified type I and type II error probabilities. Thes quare root formula is simple inform and has the advantage that equal differencesare equally detectable.     

ON A GENERAL DISTRIBUTION THEORY FOR A CLASS OF LIKELIHOOD RATIO CRITERIA

In this paper the exact distribution of a statistic with general moment function has been derived in a gamma series and also in a beta series. It is shown that the Box expansion can be obtained from the gamma series by collecting terms of the same order.     

A new ranked set sampling protocol for the signed rank test

In this paper, we construct a new ranked set sampling protocol that maximizes the Pitman asymptotic efficiency of the signed rank test. The new sampling design is a function of the set size and independent order statistics. If the set size is odd and the underlying distribution is symmetric and unimodal, then the new sampling protocol quantifies only the middle observation. On the other hand, if the set size is even, the new sampling design quantifies the two middle observations. This data collection procedure for use in the signed rank test outperforms the data collection procedure in the standard ranked set sample. We show that the exact null distribution of the signed rank statistic W_RSS⁺ based on a data set generated by the new ranked set sample design for odd set sizes is the same as the null distribution of the simple random sample signed rank statistic W_SRS⁺ based on the same number of measured observations. For even set sizes, the exact null distribution of W_RSS⁺ is simulated.     

An approximation to the exact distribution of the wilcoxon-mann-whitney rank sum test statistic

An approximation to the exact distribution of the Wilcoxon rank sum test (Mann-Whitney U-test) and the Siegel-Tukey test based on a linear combination of the two-sample t-test applied to ranks and the normal approximation is compared with the usual normal approximation. The normal approximation results in a conservative test in the tails while the linear combination of the test statistics provides a test that has a very high percentage of agreement with tables of the exact distribution. Sample sizes 3≤m, n≤50 were considered.     

Asymptotic distributions of the sphericity test in a complex multivariate normal distribution

In this paper, asymptotic expansions of the null and non-null distributions of the sphericity test criterion in the case of a complex multivariate normal distribution are obtained for the first time in terms of beta distributions. In the null case, it is found that the accuracy of the approximation by taking the first term alone in the asymptotic series is sufficient for practical purposes. In fact for p - 2. the asymptotic expansion reduces to the first term which is also the exact distribution in this case. Applications of the results to the area of inferences on multivariate time series are also given.     

Some Properties of the Exact and Score Methods for Binomial Proportion and Sample Size Calculation

In this article, we point out some interesting relations between the exact test and the score test for a binomial proportion p. Based on the properties of the tests, we propose some approximate as well as exact methods of computing sample sizes required for the tests to attain a specified power. Sample sizes required for the tests are tabulated for various values of p to attain a power of 0.80 at level 0.05. We also propose approximate and exact methods of computing sample sizes needed to construct confidence intervals with a given precision. Using the proposed exact methods, sample sizes required to construct 95% confidence intervals with various precisions are tabulated for p = .05(.05).5. The approximate methods for computing sample sizes for score confidence intervals are very satisfactory and the results coincide with those of the exact methods for many cases.     

A note on spearman's footrule

We obtain exact tables of the null distribution of Spearman's footrule for sample sizes n = 4(1)40 by using a certain Markov chain property, and we investigate the adequacy of approximations to the distribution.     

Maximizing generalized linear mixed model likelihoods with an automated Monte Carlo EM algorithm

Two new implementations of the EM algorithm are proposed for maximum likelihood fitting of generalized linear mixed models. Both methods use random (independent and identically distributed) sampling to construct Monte Carlo approximations at the E-step. One approach involves generating random samples from the exact conditional distribution of the random effects (given the data) by rejection sampling, using the marginal distribution as a candidate. The second method uses a multivariate t importance sampling approximation. In many applications the two methods are complementary. Rejection sampling is more efficient when sample sizes are small, whereas importance sampling is better with larger sample sizes. Monte Carlo approximation using random samples allows the Monte Carlo error at each iteration to be assessed by using standard central limit theory combined with Taylor series methods. Specifically, we construct a sandwich variance estimate for the maximizer at each approximate E-step. This suggests a rule for automatically increasing the Monte Carlo sample size after iterations in which the true EM step is swamped by Monte Carlo error. In contrast, techniques for assessing Monte Carlo error have not been developed for use with alternative implementations of Monte Carlo EM algorithms utilizing Markov chain Monte Carlo E-step approximations. Three different data sets, including the infamous salamander data of McCullagh and Nelder, are used to illustrate the techniques and to compare them with the alternatives. The results show that the methods proposed can be considerably more efficient than those based on Markov chain Monte Carlo algorithms. However, the methods proposed may break down when the intractable integrals in the likelihood function are of high dimension.