The birthday problem

Happy Birthday! If it is not your birthday today, it may well be the birthday of someone in your office, or in your class at school. Possibly, it is the birthday of two people in your office or your class. Just how big does a group have to be before its members are likely to share a birthday? Mario Cortina Borja and John Haigh explain the birthday problem.     

Number of matches and matched people in the birthday problem

The well known birthday problem asks for the probability of at least one match out of a group of n people. Also of interest are the number of matches and the number of matched people. In this paper the means and variances of the number of matches and matched people are obtained. A generalization of the use of these methods to computer storage analysis is discussed.     

Variation in caesarean delivery rates across hospitals: a Bayesian semi-parametric approach

This article presents a Bayesian semi-parametric approach for modeling the occurrence of cesarean sections using a sample of women delivering in 20 hospitals of Sardinia (Italy). A multilevel logistic regression has been fitted on the data using a Dirichlet process prior for modeling the random-effects distribution of the unobserved factors at the hospital level. Using the estimated random effects at the hospital level, a partition of the hospitals in terms of similar medical practice has been obtained that identifies different profiles of hospitals in terms of caesarean section risks. The limited number of clusters may be useful for suggesting policy implications that help to reduce the heterogeneity of caesarean delivery risks.     

Justifying an Invalid Derivation of the Probability of a Union of Events

An appealing, but invalid, derivation of the probability that at least one of n events occurs is justified, using a particular definition of subtraction of events. The probabilities that exactly m and at least m of the n events occur are derived similarly.     

A modification of the kruskal-wallis statistic for the generalized behrens-fisher problem

The problem of testing the equality of the medians of several populations is considered. Standard distribution-free procedures for this problem require that the populations have the same shape in order to maintain their nominal significance level, ever asymptotically, under the null hypothesis of equal medians , A modification of the Kruskal-Wallis test statistic is proposed which is exactly distribution-free under the usual nonparanetric asswnption that the continuous populations are identical with any shape. It is asymptotically distribution-free when the Continuous populations are asswned to be syrmnetric with equal medians.     

New simulation results for the calibration and inverse median estimation problems

This paper gives the results of a new simulation study for the familiar calibration problem and the less familiar inverse median estimation problem. The latter arises when one wishes to estimate from a linear regression analysis the value of the independent variable corresponding to a specified value of the median of the dependent variable. For example, from the results of a regression analysis between stress and time to failure, one might wish to estimate the stress at which the median time to failure is 10,000 hours. In the study, the mean square error, Pitman closeness, and probability of overestimation are compared for both the calibration problem and the inverse median estimation problem for (1) the classical estimator, (2) the inverse estimator, and (3) a modified version of an estimator proposed by Naszodi (1978) for both a small sample and a moderately large sample situation.     

基于模糊多分配的辽宁煤炭枢纽站选址研究

提出一种模糊多分配p枢纽站中位问题,其中运输成本定义为模糊变量,问题的目标函数是在给定的可信性水平下,最小化总的运输成本。对于梯形和正态运输成本,问题等价于确定的混合整数线性规划问题。在实证分析中,选取了辽宁省煤炭产业的相关面板数据,分析计算在不同可信度水平下煤炭运输枢纽站设立的数量和位置,再利用传统的优化方法(如分枝定界法)求解。经计算,这一模型和求解方法可以用来解决辽宁省煤炭运输的选址问题。     

Testing for Equivalence of Means under Heteroskedasticity by an Approximate Solution of a Partial Differential Equation of Infinite Order

Abstract. A test for two‐sided equivalence of means has been developed under the assumption of normally distributed populations with heterogeneous variances. Its rejection region is limited by functions ± h that depend on the empirical variances. h is stated implicitly by a partial differential equation, an exact solution of which would provide a test that is exactly similar at the boundary of the null hypothesis of non‐equivalence. h is approximated by Taylor series up to third powers in the reciprocal number of degrees of freedom. This suffices to obtain error probabilities of the first kind that are very close to a nominal level of α = 0 . 05 at the boundary of the null hypothesis. For more than 10 data points in each group, they range between 0.04995 and 0.05005, and are thus much more precise than those obtained by other authors.     

Some tests with specified size for the behrens-fisher problem

For a given significance level α, Welch's approximate t-test for the Behrens-Fisher Problem is modified to get a test with size α. A useful result for carrying out the Berger and Boos test is provided. Simulation results give power comparisons of several size α tests.     

A Test for the Equality of Parameters for Separate Regression Models in the Presence of Heteroskedasticity

Testing for the equality of regression coefficients across two regressions is a problem considered by analysts in a variety of fields. If the variances of the errors of the two regressions are not equal, then it is known that the standard large sample F-test used to test the equality of the coefficients is compromised by the fact that its actual size can differ substantially from the stated level of significance in small samples. This article addresses this problem and borrows from the literature on the Behrens-Fisher problem to provide some simple modifications of the large sample test which allows one to better control the probability of committing a Type I error. Empirical evidence is presented which indicates that the suggested modifications provide tests which are superior to well-known alternative tests over a wide range of the parameter space.     

An adjusted likelihood-ratio approach to the behrens-fisher problem

In many situations it is necessary to test the equality of the means of two normal populations when the variances are unknown and unequal. This paper studies the celebrated and controversial Behrens-Fisher problem via an adjusted likelihood-ratio test using the maximum likelihood estimates of the parameters under both the null and the alternative models. This procedure allows the significance level to be adjusted in accordance with the degrees of freedom to balance the risk due to the bias in using the maximum likelihood estimates and the risk due to the increase of variance. A large scale Monte Carlo investigation is carried out to show that -2 InA has an empirical chi-square distribution with fractional degrees of freedom instead of a chi-square distribution with one degree of freedom. Also Monte Carlo power curves are investigated under several different conditions to evaluate the performances of several conventional procedures with that of this procedure with respect to control over Type I errors and power.     

A Note on the Uniformity Assumption in the Birthday Problem

If 23 people are selected at random from some large population whose birthdays are uniformly distributed throughout the year, then it is well-known that the chances are better than even that at least two of them will have the same birthday. In this note it is shown that this is true for any distribution of birthdays.     

Large-scale multiple testing under dependence

Summary.  The paper considers the problem of multiple testing under dependence in a compound decision theoretic framework. The observed data are assumed to be generated from an underlying two-state hidden Markov model. We propose oracle and asymptotically optimal data-driven procedures that aim to minimize the false non-discovery rate FNR subject to a constraint on the false discovery rate FDR. It is shown that the performance of a multiple-testing procedure can be substantially improved by adaptively exploiting the dependence structure among hypotheses, and hence conventional FDR procedures that ignore this structural information are inefficient. Both theoretical properties and numerical performances of the procedures proposed are investigated. It is shown that the procedures proposed control FDR at the desired level, enjoy certain optimality properties and are especially powerful in identifying clustered non-null cases. The new procedure is applied to an influenza-like illness surveillance study for detecting the timing of epidemic periods.     

Statistical inference on difference or ratio of means from heteroscedastic normal populations

In this paper, we propose a two-stage sampling method (Stein, 1945) for making inferences on the means of two normal populations with heteroscedastic variances which is usually referred to as the Behrens–Fisher problem (6 and 7). The two-stage method can be used to control Type-II error and the length of confidence interval conditional on a specified Type-I error or confidence level. Simulation results show that the two-stage method improves the powers of the tests and its performance exhibits superior robustness properties over the existing single-stage methods. Moreover, our method can be easily extended to compare multiple populations.     

Testing the Equality of Covariance Operators in Functional Samples

Abstract. We propose a non‐parametric test for the equality of the covariance structures in two functional samples. The test statistic has a chi‐square asymptotic distribution with a known number of degrees of freedom, which depends on the level of dimension reduction needed to represent the data. Detailed analysis of the asymptotic properties is developed. Finite sample perfo‐rmance is examined by a simulation study and an application to egg‐laying curves of fruit flies.     

Locally optimal invariant tests for change in level

The classical change point problem is considered, from the invariance point of view. Locally optimal invariant tests are derived for the change in level, when the initial level and the common variance are assumed to be unknown. The tests derived by Chernoff and Zacks (1964) and Gardner (1969), for the change in level, when variance is known, are shown to be locally optimal invariant tests.     

The probability generating function of empty cell variable in a randomized occupancy problem

Let M_o denote the number of empty cells when n distinguishable balls are distributed independently and at random in ra cells such that each ball stays with probability p in its cell, and falls through with probability 1-p. We find the probability generating function of M_o by solving a partial differential equation satisfied by a suitable generating function. The corresponding function for the classical case p = 1 is well-known, but obtained by different methods.     

Modifications of the Bonferroni-Holm procedure for a multi-way ANOVA

We present step-wise test procedures based on the Bonferroni-Holm principle for multi-way ANOVA-type models. It is shown for two plausible modifications that the multiple level α is preserved. These theoretical results are supplemented by a simulation study, in a two-way ANOVA setting, to compare the multiple procedures with respect to their simultaneous power and the relative frequency of correctly rejected false hypotheses. Financial support of the Deutsche Forschungsgemeinschaft is gratefully acknowledged.     

Anna Karenina and the two envelopes problem

The Anna Karenina principle is named after the opening sentence in the eponymous novel: Happy families are all alike; every unhappy family is unhappy in its own way. The two envelopes problem (TEP) is a much-studied paradox in probability theory, mathematical economics, logic and philosophy. Time and again a new analysis is published in which an author claims finally to explain what actually goes wrong in this paradox. Each author (the present author included) emphasises what is new in their approach and concludes that earlier approaches did not get to the root of the matter. We observe that though a logical argument is only correct if every step is correct, an apparently logical argument which goes astray can be thought of as going astray at different places. This leads to a comparison between the literature on TEP and a successful movie franchise: it generates a succession of sequels, and even prequels, each with a different director who approaches the same basic premise in a personal way. We survey resolutions in the literature with a view to synthesis, correct common errors, and give a new theorem on order properties of an exchangeable pair of random variables, at the heart of most TEP variants and interpretations. A theorem on asymptotic independence between the amount in your envelope and the question whether it is smaller or larger shows that the pathological situation of improper priors or infinite expectation values has consequences as we merely approach such a situation.     

One-way ANOVA with Unequal Variances

We consider the one-way ANOVA problem of testing the equality of several normal means when the variances are not assumed to be equal. This is a generalization of the Behrens-Fisher problem, but even in this special case there is no exact test and the actual size of any test depends on the values of the nuisance parameters. Therefore, controlling the actual size of the test is of main concern. In this article, we first consider a test using the concept of generalized p-value. Extensive simulation studies show that the actual size of this test does not exceed the nominal level, for practically all values of the nuisance parameters, but the test is not too conservative either, in the sense that the actual size of the test can be very close to the nominal level for some values of the nuisance parameters. We then use this test to propose a simple F-test, which has similar properties but avoids the computations associated with generalized p-values. Because of its simplicity, both conceptually as well as computationally, this F-test may be more useful in practice, since one-way ANOVA is widely used by practitioners who may not be familiar with the generalized p-value and its computational aspects.