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.     

Irregularly observed time series – some asymptotics and the block bootstrap

We consider time series being observed at random time points. In addition to Parzen's classical modelling by amplitude modulating sequences, we state another modelling using an integer-valued sequence as the observation times. Limiting results are presented for the sample mean and are generalized to the class of functions of smooth means. Motivated by the complicated limiting behaviour, (moving) block bootstrap possibilities are investigated. Conditional on the used modelling for the irregular spacings, one is lead to different interpretations for the block length and hence bootstrap approaches. The block length either can be interpreted as the time (resulting in an observation string of fixed length containing a random number of observations) or as the number of observations (resulting in an observation string of variable length containing a fixed number of values). Both bootstrap approaches are shown to be asymptotically valid for the sample mean. Numerical examples and an application to real-world ozone data conclude the study.     

Tables for pairwise multiple comparisons using shaffer's modifiedsequentially-rejective procedure

A sequentially rejective (SR) testing procedure introduced by Holm (1979) and modified (MSR) by Shaffer (1986) is considered for testing all pairwise mean comparisons.For such comparisons, both the SR and MSR methods require that the observed test statistics be ordered and compared, each in turn, to appropriate percentiles on Student's t distribution.For the MSR method these percentiles are based on the maximum number of true null hypotheses remaining at each stage of the sequential procedure, given prior significance at previous stages, A function is developed for determining this number from the number of means being tested and the stage of the test.For a test of all pairwise comparisons, the logical implications which follow the rejection of a null hypothesis renders the MSR procedure uniformly more powerful than the SR procedure.Tables of percentiles for comparing K means, 3 < K < 6, using the MSR method are presented.These tables use Sidak's (1967) multiplicative inequality and simplify the use of t he MSR procedure.Several modifications to the MSR are suggested as a means of further increasing the power for testing the pairwise comparisons.General use of the MSR and the corresponding function for testing other parameters besides the mean is discussed.     

MODEL-BASED VARIANCE ESTIMATION IN SURVEYS WITH STRATIFIED CLUSTERED DESIGN

A model-based method for estimating the sampling variances of estimators of (sub-)population means, proportions, quantiles, and regression parameters in surveys with stratified clustered design is described and applied to a survey of US secondary education. The method is compared with the jackknife by a simulation study. The model-based estimators of the sampling variances have much smaller mean squared errors than their jackknife counterparts. In addition, they can be improved by incorporating information about the unknown parameters (variances) from external sources. A regression-based smoothing method for estimating the sampling variances of the estimators for a large number of subpopulation means is proposed. Such smoothing may be invaluable when subpopulations are represented in the sample by only few subjects.     

Sequential testing procedures for a class of distributions representing various life-testing models

A class of probability density functions is considered, which covers several life-testing models as specific cases. Sequential probability ratio tests are developed for testing simple and composite hypotheses regarding the parameters of the probabilistic model. Expressions for the operating characteristic and the average sample number functions are derived and their behaviour is studied by means of graph-plotting.     

Comparison of ANOVA-F and ANOM tests with regard to type I error rate and test power

A Monte Carlo simulation was conducted to compare the type I error rate and test power of the analysis of means (ANOM) test to the one-way analysis of variance F-test (ANOVA-F). Simulation results showed that as long as the homogeneity of the variance assumption was satisfied, regardless of the shape of the distribution, number of group and the combination of observations, both ANOVA-F and ANOM test have displayed similar type I error rates. However, both tests have been negatively affected from the heterogeneity of the variances. This case became more obvious when the variance ratios increased. The test power values of both tests changed with respect to the effect size (Δ), variance ratio and sample size combinations. As long as the variances are homogeneous, ANOVA-F and ANOM test have similar powers except unbalanced cases. Under unbalanced conditions, the ANOVA-F was observed to be powerful than the ANOM-test. On the other hand, an increase in total number of observations caused the power values of ANOVA-F and ANOM test approach to each other. The relations between effect size (Δ) and the variance ratios affected the test power, especially when the sample sizes are not equal. As ANOVA-F has become to be superior in some of the experimental conditions being considered, ANOM is superior in the others. However, generally, when the populations with large mean have larger variances as well, ANOM test has been seen to be superior. On the other hand, when the populations with large mean have small variances, generally, ANOVA-F has observed to be superior. The situation became clearer when the number of the groups is 4 or 5.     

The number of decision sets involving n means

A decision set is the result of the multiple comparison of n means. The means are placed in order in a diagram and a line is drawn under those that are not significantly different. In this paper, we show how the number of decision sets involving n means can be calculated.     

A modified subset selection formulation with special reference to one-way and two-way layout experiments

The usual formulation of subset selection due to Gupta (1956) requires a minimum guaranteed probability of a correct selection. The modified formulation of the present paper includes an additional requirement that the expected number of the nonbest populations be bounded above by a specified constant when the best and the next best populations are ‘sufficiently’ apart. A class of procedures is defined and the determination of the minimum sample size required is discussed. The specific problems discussed for normal populations include selection in terms of means and variances, and selection in terms of treatment effects in a two-way layout.     

Estimating the population median using only a central subset of observations

In estimating the population median, it is common to encounter estimators which are linear combinations of a small number of central observations. Sample medians, sample quasi medians, trimmed means, jackknifed (and delete‐d jackknifed) medians and jackknifed quasi medians are all familiar examples. The objective of this paper is to show that within this class the quasi medians turn out to have the best asymptotic mean squared error.     

Comparing Nonparametric Tests of Equality of Means for Randomized Block Designs

A number of nonparametric tests are compared empirically for a randomized block layout. We assess tests appropriate for when the data are not consistent with normality or when outliers invalidate traditional analysis of variance (ANOVA) tests. The objective is to assess, within this setting, tests that use ranks within blocks, the rank transform procedure that ranks the complete sample and continuous analogs of the Cochran–Mantel–Haenszel tests. The usual linear model is assumed, and our primary foci are tests of equality of means and component tests that assess linear and quadratic trends in the means. These tests include the traditional Page and Friedman tests. We conclude that the rank transform tests have competitive power and warrant greater use than is currently apparent.     

A comparison of selection procedures for selecting the best normal population under heterogeneity of variance

This study examines the comparative probabilities of making a correct selection when using the means procedure (M), the medians procedure (D) and the rank-sum procedure (S) to correctly select the normal population with the largest mean under heterogeneity of variance. The comparison is conducted by using Monte-Carlo simulation techniques for 3, 4, and 5 normal populations under the condition that equal sample sizes are taken from each population. The population means and standard deviations are assumed to be equally-spaced. Two types of heterogeneity of variance are considered: (1) associating larger means with larger variances, and (2) associating larger means with smaller variances.     

Power and Sample Size Determination When the Alternative Hypotheses are Given in Quantiles

The power of normal-theory tests about means depends on a noncentrality parameter which is a function of the unknown parameter σ. In order to calculate power and to solve sample-size problems based on power, differences between hypothesized and alternative values of the means are frequently selected as a multiple of σ, a choice which eliminates σ from the noncentrality parameter and permits a solution. Perhaps a more natural (but equivalent) way to express alternatives is to give one or more means as the quantile of order p (say Q_p ) of a distribution with another mean. As we will demonstrate, this kind of alternative also eliminates σ from the problem.     

Post-hoc selection of covariates in randomized experiments

Monte Carlo methods are used to compere a number of adaptive strategies for deciding which of several covariates to incorporate into the analysis of a randomized experiment.Sixteen selection strategies in three categories are considered: 1)select covariates correlated with the response, 2)select covariates with means differing across groups, and 3)select covariates with means differing across groups that are also correlated with the response. The criteria examined are the type I error rate of the test for equality of adjusted group means and the variance of the estimated treatment effect. These strategies can result in either inflated or deflated type I errors, depending on the method and the population parameters. The adaptive methods in the first category some times yieldpoint estimates of the treatment effect more precise than estimators derive dusing either all or none of the covariates.     

Nonparametric bootstrap of sample means of positive-definite matrices with an application to diffusion-tensor-imaging data analysis

This paper presents nonparametric two-sample bootstrap tests for means of random symmetric positive-definite (SPD) matrices according to two different metrics: the Frobenius (or Euclidean) metric, inherited from the embedding of the set of SPD metrics in the Euclidean set of symmetric matrices, and the canonical metric, which is defined without an embedding and suggests an intrinsic analysis. A fast algorithm is used to compute the bootstrap intrinsic means in the case of the latter. The methods are illustrated in a simulation study and applied to a two-group comparison of means of diffusion tensors (DTs) obtained from a single voxel of registered DT images of children in a dyslexia study.     

New approximate confidence intervals for the difference between two Poisson means and comparison

The problem of estimating the difference between two Poisson means is considered. A new moment confidence interval (CI), and a fiducial CI for the difference between the means are proposed. The moment CI is simple to compute, and it specializes to the classical Wald CI when the sample sizes are equal. Numerical studies indicate that the moment CI offers improvement over the Wald CI when the sample sizes are different. Exact properties of the CIs based on the moment, fiducial and hybrid methods are evaluated numerically. Our numerical study indicates that the hybrid and fiducial CIs are in general comparable, and the moment CI seems to be the best when the expected total counts from both distributions are two or more. The interval estimation procedures are illustrated using two examples.     

Some inference results on pr(x < y) in the bivariate exponential model

This paper provides three different estimators for Pr(X < Y) when X and Y have a bivariate exponential distribution. The asymptotic variances of the three estimators are also derived. A test for the equality of the means of X and Y and confidence limits for the difference of the two means are presented. Our results are directly applicable in a reliability context with underlying bivariate exponential distribution.     

A fiducial approach to multiple comparisons

Comparing treatment means from populations that follow independent normal distributions is a common statistical problem. Many frequentist solutions exist to test for significant differences amongst the treatment means. A different approach would be to determine how likely it is that particular means are grouped as equal. We developed a fiducial framework for this situation. Our method provides fiducial probabilities that any number of means are equal based on the data and the assumed normal distributions. This methodology was developed when there is constant and non-constant variance across populations. Simulations suggest that our method selects the correct grouping of means at a relatively high rate for small sample sizes and asymptotic calculations demonstrate good properties. Additionally, we have demonstrated the flexibility in the methods ability to calculate the fiducial probability for any number of equal means. This was done by analyzing a simulated data set and a data set measuring the nitrogen levels of red clover plants that were inoculated with different treatments.     

Analysis of Means Using Ranks

A nonparametric procedure, called the analysis of means using ranks (ANOMR), is proposed for testing the equality of several population means. The ANOMR procedure may be used graphically in the form of a Shewhart control chart and so has the advantage of pinpointing which population mean, if any, is significantly different from the others. Exact and asymptotic critical values are given for the implementation of ANOMR. Results from a Monte Carlo power study are presented which indicate that for light-tailed distributions such as the uniform and the normal, ANOMR is only slightly less powerful than the parametric competitive procedures based on analysis of variance and analysis of means. For heavy-tailed distributions such as the Cauchy, ANOMR is shown to provide greater power than the parametric procedures. The results also indicate that for both light and heavy-tailed distributions the use of the ANOMR test instead of the Kruskal-Wallis test leads to only a small loss of power for a range of alternatives.     

Simple and effective number-of-bins circumference selectors for a histogram

Two very effective data-based procedures which are simple and fast to compute are proposed for selecting the number of bins in a histogram. The idea is to choose the number of bins that minimizes the circumference (or a bootstrap estimate of the expected circumference) of the frequency histogram. Contrary to most rules derived in the literature, our method is therefore not dependent on precise asymptotic analyses. It is shown by means of an extensive Monte-Carlo study that our selectors perform well in comparison with recently suggested selectors in the literature, for a wide range of density functions and sample sizes. The behaviour of one of the proposed rules is also illustrated on real data sets.     

Estimation and hypothesis in linear models-A reparameterization approach to the cell mlans modll

A general approach to analysis of fixed effects models through a reparameterization of the cell means model is outlined. After applying a Q-op3rator, distribution theory for the model is developed. Methods of estimation and hypothesis testing are given and the exact form of each hypothesis is shown. The implementation of these results in a currently available computer package is discussed.