Testing for a treatment effect in a heterogeneous population: a modified sign-test statistic and a leapfrog statistic

This paper proposes two non-parametric statistics that test for a treatment effect in a heterogeneous population. In the model considered, data on two examinations for both a control and a treatment group are needed to perform the test. The model allows for individual (fixed) effects that may be correlated with the choice of treatment. In addition, the model allows for an unspecified, monotonic transformation of the response variable. The techniques are illustrated by testing whether high levels of unemploymentbenefit eligibility affect the consumption patterns of unemployed American workers.     

Fitting a Curve to a Confidence Interval

Three curves are fit to a confidence interval and its estimate. These curves allow a data user to incorporate the variability of the estimation method into further calculations.     

Testing for a Unit Root in a Time Series With a Changing Mean

This study considers testing for a unit root in a time series characterized by a structural change in its mean level. My approach follows the “intervention analysis” of Box and Tiao (1975) in the sense that I consider the change as being exogenous and as occurring at a known date. Standard unit-root tests are shown to be biased toward nonrejection of the hypothesis of a unit root when the full sample is used. Since tests using split sample regressions usually have low power, I design test statistics that allow the presence of a change in the mean of the series under both the null and alternative hypotheses. The limiting distribution of the statistics is derived and tabulated under the null hypothesis of a unit root. My analysis is illustrated by considering the behavior of various univariate time series for which the unit-root hypothesis has been advanced in the literature. This study complements that of Perron (1989), which considered time series with trends.     

On uncertainty of a proportion from a stratified random sample of a small population

This article considers the uncertainty of a proportion based on a stratified random sample of a small population. Using the hypergeometric distribution, a Clopper–Pearson type upper confidence bound is presented. Another frequentist approach that uses the estimated variance of the proportion estimator is also considered as well as a Bayesian alternative. These methods are demonstrated with an illustrative example. Some aspects of planning, that is, the impact of specified strata sample sizes, on uncertainty are studied through a simulation study.     

Matching a correlation coefficient by a Gaussian copula

A Gaussian copula is widely used to define correlated random variables. To obtain a prescribed Pearson correlation coefficient of ρ_x between two random variables with given marginal distributions, the correlation coefficient ρ_z between two standard normal variables in the copula must take a specific value which satisfies an integral equation that links ρ_x to ρ_z. In a few cases, this equation has an explicit solution, but in other cases it must be solved numerically. This paper attempts to address this issue. If two continuous random variables are involved, the marginal transformation is approximated by a weighted sum of Hermite polynomials; via Mehler’s formula, a polynomial of ρ_z is derived to approximate the function relationship between ρ_x and ρ_z. If a discrete variable is involved, the marginal transformation is decomposed into piecewise continuous ones, and ρ_x is expressed as a polynomial of ρ_z by Taylor expansion. For a given ρ_x, ρ_z can be efficiently determined by solving a polynomial equation.     

On a random series in a hilbert space

Let (S_n) be partial sums of a non-degenerate sequence of Identically and independently distributed random variables taking values in a separable Hilbert space. Then for 0 ≤ β ≤ 3/2, the series converges almost nowhere. For β > 3/2 this may not be true.     

The Asymptotic Behavior of a Brownian Motion with a Drift from a Random Domain

Consider a Brownian motion with drift starting at an interior point of a random domain D in R ^d+1, d ≥ 1, let τ _D denote the first time the Brownian motion exits from D. Estimates with exact constants for the asymptotics of log P(τ _D  > T) are given for T → ∞, depending on the shape of the domain D and the order of the drift. The problem is motivated by the model in insurance and early works of Lifshits and Shi. The methods of proof are based on the calculus of variations and early works of Li, Lifshits and Shi in the drift free case.     

When is a drug not a drug?

Stephen Senn considers the strange world of bioequivalence, where studies that fail to show a difference between formulations are judged a success.     

Geometry of a corelation coefficient under a copula

Spearman's rank correlation coefficient is shown to be a measure of distance on the unit square, which characterizes the concentration of the probability density under a copula. A distance function offers insight into structuring copulas with a desired degree of association.     

On a confidence interval about a parameter estimated by a ratio of normal random variables

A confidence interval is geometrically constructed about a parameter estimated by the ratio of bivariate normal random variables. The resulting confidence interval is equivalent to that of Fieller's theorem. The geometric construction shown that such intervals are conservative. Bioassay examples are used to demonstrate the technique.     

Some Asymptotic Formulas for a Brownian Motion With a Regular Variation From a Parabolic Domain

Consider a Brownian motion with a regular variation starting at an interior point of a domain D in R^{d + 1}, d ? 1 and let τ_D denote the first time the Brownian motion exits from D. Estimates with exact constants for the asymptotics of log?P(τ_D > T) are given for T → ∞, depending on the shape of the domain D and the order of the regular variation. Furthermore, the asymptotically equivalence are obtained. The problem is motivated by the early results of Lifshits and Shi, Li in the first exit time, and Karamata in the regular variation. The methods of proof are based on their results and the calculus of variations.     

Uniform asymptotic linearity in a regression parameter of a process based on a rank statistic

For X₁, …, X_N a random sample from a distribution F, let the process S^δ_N(t) be defined as where K²_N = σ^N_i=1(c_i ? c?)² and R x_i, + Δd, is the rank of X_i + Δd_i, among X₁ + Δd₁, …, X_N + Δd_N. The purpose of this note is to prove that, under certain regularity conditions on F and on the constants c_i and d_i, S^Δ_N (t) is asymptotically approximately a linear function of Δ, uniformly in t and in Δ, |Δ| ≤ C. The special case of two samples is considered.     

Spatio-temporal analysis of a plant disease in a non-uniform crop: a Monte Carlo approach

Identification of the type of disease pattern and spread in a field is critical in epidemiological investigations of plant diseases. For example, an aggregation pattern of infected plants suggests that, at the time of observation, the pathogen is spreading from a proximal source. Conversely, a random pattern suggests a lack of spread from a proximal source. Most of the existing methods of spatial pattern analysis work with only one variety of plant at each location and with uniform genetic disease susceptibility across the field. Pecan orchards, used in this study, and other orchard crops are usually composed of different varieties with different levels of susceptibility to disease. A new measure is suggested to characterize the spatio-temporal transmission patterns of disease; a Monte Carlo test procedure is proposed to test whether the transmission of disease is random or aggregated. In addition, we propose a mixed-transmission model, which allows us to quantify the degree of aggregation effect.     

Target detection from a classical and a Bayesian viewpoint

Two approaches have been used for designing spatial surveys to detect a target. The classical approach controls the probability of missing a target that exists; a Bayesian approach controls the probability that a target exists given that none was seen. In both cases, information about the likely size of the target can reduce sampling requirements. In this paper, previous results are summarized and then used to assess the risk that Roman remains could be present at sites scheduled for development in Greater London.     

Statistical games against a prophet-proof of a minimax conjecture

“Prophet theory” quantifies the price a statistician has to pay for his lack of information in stochastic sequences. In a recent paper, Schmitz (1991) gave a game-theoretical interpretation of this situation and he formulated in particular a minimax conjecture for the difference case. In this note we prove that conjecture and, moreover, present minimax ran domized stopping times (minimax procedures for the statistician).     

How to Approximate a Histogram by a Normal Density

Which normal density curve best approximates the sample histogram? The answer suggested here is the normal curve that minimizes the integrated squared deviation between the histogram and the normal curve. A simple computational procedure is described to produce this best-fitting normal density. A few examples are presented to illustrate that this normal curve does indeed provide a visually satisfying fit, one that is better than the traditional , s answer. Some further aspects of this procedure are investigated. In particular it is shown that there is a satisfactory answer that is independent of the bar width of the histogram. It is also noted that this graphical procedure provides highly robust estimates of the sample mean and standard deviation. We demonstrate our technique by using data including Newcomb's data of passage time of light and Fisher's iris data.     

Wavelet Estimation of a Density in a GARCH-type Model

We consider the GARCH-type model: S = σ² Z, where σ² and Z are independent random variables. The density of σ² is unknown whereas the one of Z is known. We want to estimate the density of σ² from n observations of S under some dependence assumption (the exponentially strongly mixing dependence). Adopting the wavelet methodology, we construct a nonadaptive estimator based on projections and an adaptive estimator based on the hard thresholding rule. Taking the mean integrated squared error over Besov balls, we prove that the adaptive one attains a sharp rate of convergence.     

Predicting a Future Median Life through a Power Transformation

A simple and unified prediction interval (PI) for the median of a future lifetime can be obtained through a power transformation. This interval usually possesses the correct coverage, at least asymptotically, when the transformation is known. However, when the transformation is unknown and is estimated from the data, a correction is required. A simple correction factor is derived based on large sample theory. Simulation shows that the unified PI after correction performs well. When compared with the existing frequentist PI's, it shows an equivalent or a better performance in terms of coverage probability and average length of the interval. Its nonparametric aspect and the ease of usage make it very attractive to practitioners. Real data examples are provided for illustration.     

Estimation of size and mean of a sensitive quantitative variable for a sub-group of a population

In the present paper, a new problem in survey sampling has been considered and an attempt has bcen made to find its solution. The procedures available in literature deal with the estimation of mean of a sensitive quantitative variable for the whole population. Sometimes the investigator may be interested in estimating the mean of a sensitive quantitative variable for a sub-group of the population. At the same time, he/she may also be interested in estimating the size of this sub-group. In this paper, we have obtained estimators of these parameters. Approximate bias and variance expressions are also obtained for these estimators.