A MODEL OF SYSTEMATIC SAMPLING WITH UNEQUAL PROBABILITIES

A new method is described of drawing, without replacement, two sample units per stratum from any population. The method is developed from a consideration of the asymptotic properties of systematic sampling with unequal probabilities, as the sizes of the population units tend to zero. The essential properties of this method are very easily analysed. They also converge, over a large number of strata, to those of systematic sampling from the same strata with their population units arranged in random order. In proving this, the assumption is made that the underlying population is of the type to which it is appropriate to apply ratio estimation. The sampling method described is, however, simple enough to commend itself as an alternative to systematic sampling when the underlying population is not of this type. Consideration is given to the case where the sizes of some of the population units exceed the skip interval.     

GEOMETRIC DECAY OF THE STEADY-STATE PROBABILITIES IN A QUASI-BIRTH-AND-DEATH PROCESS WITH A COUNTABLE NUMBER OF PHASES

A sufficient condition is proved for geometric decay of the steady-state probabilities in a quasi-birth-and-death process having a countable number of phases in each level. If there is a positive number η and positive vectors x = (x _i) and y = (y _j ) satisfying some equations and inequalities, the steady-state probability π _mi decays geometrically with rate η in the sense π _mi ～ cη ^m x _i as m → ∞. As an example, the result is applied to a two-queue system with shorter queue discipline.     

ON ESTIMATING THE MEAN OF A FINITE POPULATION SAMPLED ON TWO OCCASIONS WITH VARYING PROBABILITIES

The problem considered is one of estimating the current value of a finite population mean from two samples selected with varying probabilities on the current and a previous occasion. Following the generalized least squares approach of Gurney and Daly (1965) we extend the works of Sin& (1968) and others to obtain two sampling strategies which are respectively more efficient than the ones suggested earlier by Raj (1965) and Chotai (1974).     

CONTROLLED SAMPLING WITH VARYING PROBABILITIES WITH AND WITHOUT REPLACEMENT

Goodman and Kish (1950) introduced the problem of controlled selection in the sense of decreasing the selection probability of nonpreferred combinations, such as, for example, combinations which present organisational difficulties involving additional cost, etc. The authors (Avadhani and Sukhatme (1965), Sukhatme and Avadhani (1965)) have evolved certain techniques of controlled selection which eliminate altogether some non-preferred combinations and reduce the probability of selection of the remaining combinations, if any, to the minimum possible extent without deviating from the fundamental principles of random sampling. However, it is often felt convenient in practice to draw sampling units one after another from the population rather than combinations of units. But, at present no technique by which one can select the units one after another and at the same time reduce the probability of selection of non-preferred units is available in literature.
In this paper we have suggested a solution to this problem which not only minimizes the selection probability of non-preferred units (and of samples containing predominantly large numbers of nonpreferred units) but also provides more efficient estimators than the the usual probability proportional to size (P.P.S) sampling scheme.     

THE OPTIMAL TEST FOR SELECTING THE GREATER OF TWO BINOMIAL PROBABILITIES

The efficiency of a sequential test is related to the “importance” of the trials within the test. This relationship is used to find the optimal test for selecting the greater of two binomial probabilities, p_α and p_b, namely, the stopping rule is “gambler's ruin” and the optimal discipline when p_α+p_b≤ 1 (≥ 1) is play-the-winner (loser), i.e. an α-trial which results in a success is followed by an α-trial (b-trial) whereas an α-trial which results in a failure is followed by α b-trid (α-trial) and correspondingly for b-trials.     

ON THE UNIMODALITY OF TRANSITION PROBABILITIES IN MARKOV CHAINS

Consider a discrete time Markov chain X(n) denned on {0,1,…} and let P be the transition probability matrix governing X(n). This paper shows that, if a transformed matrix of P is totally positive of order 2, then poj(n) and pio(n) are unimodal with respect to n, where pij(n) = Pr[X(n) = j |X(0) = i]. Furthermore, the modes of poj(n) and pio(n) are non-increasing in j and I, respectively, when additionally P itself is totally positive of order 2. These results are transferred to a class of semi-Markov processes via a uniformization.     

VARIANCE ESTIMATION WITH HIGHLY STRATIFIED SAMPLING DESIGNS WITH UNEQUAL PROBABILITIES

It is common practice to design a survey with a large number of strata. However, in this case the usual techniques for variance estimation can be inaccurate. This paper proposes a variance estimator for estimators of totals. The method proposed can be implemented with standard statistical packages without any specific programming, as it involves simple techniques of estimation, such as regression fitting.     

TESTING THE SIGNIFICANCE OF GENOTYPE-ENVIRONMENT INTERACTION

The AMMI (additive main effects-multiplicative interaction) model is often used to investigate interactions in two-way tables, in particular for genotype-environment interactions. Both Gollob (1968) and Mandel (1969, 1971) proposed methods for testing the significance of such interactions. These methods are compared using simulated data. Our results support Mandel's conclusions, but his method is conservative and we recommend a test proposed by Johnson & Graybill (1972).     

THE MONTE CARLO EVALUATION OF ORTHANT PROBABILITIES FOR MULTIVARIATE NORMAL DISTRIBUTIONS

The calculation of multivariate normal orthant probabilities is practically impossible when the number of variates is greater than five or six, except in very special cases. A transformation of the integral is obtained which enables quite accurate Monte Carlo estimates to be obtained for a fairly high number of dimensions, particularly if control variates are used.     

BAYESIAN ESTIMATION OF RUIN PROBABILITIES WITH A HETEROGENEOUS AND HEAVY-TAILED INSURANCE CLAIM-SIZE DISTRIBUTION

This paper describes a Bayesian approach to make inference for risk reserve processes with an unknown claim‐size distribution. A flexible model based on mixtures of Erlang distributions is proposed to approximate the special features frequently observed in insurance claim sizes, such as long tails and heterogeneity. A Bayesian density estimation approach for the claim sizes is implemented using reversible jump Markov chain Monte Carlo methods. An advantage of the considered mixture model is that it belongs to the class of phase‐type distributions, and thus explicit evaluations of the ruin probabilities are possible. Furthermore, from a statistical point of view, the parametric structure of the mixtures of the Erlang distribution offers some advantages compared with the whole over‐parametrized family of phase‐type distributions. Given the observed claim arrivals and claim sizes, we show how to estimate the ruin probabilities, as a function of the initial capital, and predictive intervals that give a measure of the uncertainty in the estimations.     

BACKWARD ESTIMATION OF STOCHASTIC PROCESSES WITH FAILURE EVENTS AS TIME ORIGINS

Stochastic processes often exhibit sudden systematic changes in pattern a short time before certain failure events. Examples include increase in medical costs before death and decrease in CD4 counts before AIDS diagnosis. To study such terminal behavior of stochastic processes, a natural and direct way is to align the processes using failure events as time origins. This paper studies backward stochastic processes counting time backward from failure events, and proposes one-sample nonparametric estimation of the mean of backward processes when follow-up is subject to left truncation and right censoring. We will discuss benefits of including prevalent cohort data to enlarge the identifiable region and large sample properties of the proposed estimator with related extensions. A SEER-Medicare linked data set is used to illustrate the proposed methodologies.     

TESTING THE LARGEST OF A SET OF CORRELATION COEFFICIENTS

A previous paper which studied the distribution of the smallest distance between N independent random points on the surface of a sphere is generalised to higher dimensions in order to study the distribution of the largest sample correlation coefficient between a set of independent normally distributed variables. Inclusion-exclusion arguments provide accurate bounds for the tail of this distribution, and by another argument more exact bounds are also found, one of which is an improvement on the result in the previous paper. Bounds are also found for the power of the test against the alternative hypothesis that one only of the population correlation coefficients is non-zero. The test is also shown to be the likelihood ratio test against the latter alternative.     

ON TRANSITION PROBABILITIES OF SKIP-FREE MARKOV CHAINS

Consider an ergodic Markov chain X(t) in continuous time with an infinitesimal matrix Q = (q_ij) defined on a finite state space {0, 1,…, N}. In this note, we prove that if X(t) is skip-free positive (negative, respectively), i.e., q_ij, = 0 for j > i+ 1 (i > j+ 1), then the transition probability p_ij(t) = Pr[X(t)=j | X(0) =i] can be represented as a linear combination of p_0N(t) (p^(m)_(N0)(t)), 0 ≤ m ≤N, where f^(m)(t) denotes the mth derivative of a function f(t) with f⁽⁰⁾(t) =f(t). If X(t) is a birth-death process, then p_ij(t) is represented as a linear combination of p_0N^(m)(t), 0 ≤m≤N - |i-j|.     

A NOTE ON OPTIMUM STRATIFICATION IN SAMPLING WITH VARYING PROBABILITIES1

Singh and Sukhatme [4] have considered the problem of optimum stratification on an auxiliary variable x when the units from the different strata are selected with probability proportional to the value of the auxiliary variable and the sample sizes for the different strata are determined by using Neyman allocation method. The present paper considers the same problem for the proportional and equal allocation methods. The rules for finding approximately optimum strata boundaries for these two allocation methods have been given. An investigation into the relative efficiency of these allocation methods with respect to the Neyman allocation has also been made. The performance of equal allocation is found to be better than that of proportional allocation and practically equivalent to the Neyman allocation.     

ON TESTING THE GOODNESS-OF-FIT OF NONLINEAR HETEROSCEDASTIC REGRESSION MODELS

For testing the adequacy of a parametric model in regression, various test statistics can be constructed on the basis of a marked empirical process of residuals. By using a discretized version of the decomposition of the corresponding Gaussian limiting process into its principal components, we obtain a test statistic with an asymptotic chi-squared distribution under the null hypothesis. We investigate the consistency of this test statistic and of the estimators needed to compute it. Numerical experiments indicate that the distributional approximations already work for small to moderate sample sizes and reveal that the test has good power properties against a variety of alternatives. The test has a simple implementation. We present an application to a real-data example for testing the adequacy of a possible heteroscedastic exponential model.     

TESTING THE EQUALITY OF PARAMETERS IN THE DISTRIBUTIONS WHEN THE RANGE DEPENDS UPON THE PARAMETER

The authors derive the null and non-null distributions of the test statistic v=y_min/y_max (where y_min= min x_ij, y_max= max x_ij, J=1,2, …, k) connected with testing the equality of scale parameters θ₁, θ₂, …θ_k in certain, class of density functions given by      

TESTING THE STABILITY OF REGRESSION COEFFICIENTS USING GENERALIZED RECURSIVE RESIDUALS

Simple but flexible methods to detect deviations from the assumption of constant coefficients in linear regression are presented. Based on recursive residuals a runs test is developed as an alternative to CUSUM- and MOSUM-techniques. Finally a simulation study gives insight into the new method.