CONFIDENCE INTERVALS FOR THE DIFFERENCE BETWEEN TWO PARTIAL AUCS

As new diagnostic tests are developed and marketed, it is very important to be able to compare the accuracy of a given two continuous‐scale diagnostic tests. An effective method to evaluate the difference between the diagnostic accuracy of two tests is to compare partial areas under the receiver operating characteristic curves (AUCs). In this paper, we review existing parametric methods. Then, we propose a new semiparametric method and a new nonparametric method to investigate the difference between two partial AUCs. For the difference between two partial AUCs under each method, we derive a normal approximation, define an empirical log‐likelihood ratio, and show that the empirical log‐likelihood ratio follows a scaled chi‐square distribution. We construct five confidence intervals for the difference based on normal approximation, bootstrap, and empirical likelihood methods. Finally, extensive simulation studies are conducted to compare the finite‐sample performances of these intervals, and a real example is used as an application of our recommended intervals. The simulation results indicate that the proposed hybrid bootstrap and empirical likelihood intervals outperform other existing intervals in most cases.     

BAYESIAN ANALYSIS FOR THE SUPERPOSITION OF TWO DEPENDENT NONHOMOGENEOUS POISSON PROCESSES

ABSTRACT

A Bayesian analysis for the superposition of two dependent nonhomogenous Poisson processes is studied by means of a bivariate Poisson distribution. This particular distribution presents a new likelihood function which takes into account the correlation between the two nonhomogenous Poisson processes. A numerical example using Markov Chain Monte Carlo method with data augmentation is considered.     

THE DISPLACED POISSON DISTRIBUTION

Consideration of the number of events in a Poisson process in excess of a threshold value r , when it is assumed that at least r events do occur, leads to the Displaced Poisson Distribution. The recurrence relationship between successive probabilities is a simple generalization of that of the Poisson distribution. Accordingly the fitting of the distribution to data is straightforward, using either a desk machine or computer.
It is a highly flexible unimodal distribution catering for both under- and over-dispersed data and generally fits Poisson type data well. Its properties have been discussed, moment and maximum likelihood estimators for the cases, r known and unknown have been given, together with several numerical examples.     

ESTIMATING THE ANGLE BETWEEN THE MEAN DIRECTIONS OF TWO SPHERICAL DISTRIBUTIONS

This paper considers the problem of calculating a confidence interval for the angular difference between the mean directions of two spherical random variables with rotationally symmetric unimodal distributions. For large sample sizes, it is shown that the asymptotic distribution of 1 – cos α, where α is the sample angular difference, is approximately exponential if the true difference is zero, and approximately normal for a ‘large’ true difference; a scaled beta approximation is determined for the general case. For small sample sizes, a bootstrap approach is recommended. The results are applied to two sets of palaeomagnetic data.     

SMOOTH TESTS FOR THE BIVARIATE POISSON DISTRIBUTION

A theorem of Rayner & Best (1989) is generalised to permit the construction of smooth tests of goodness of fit without requiring a set of orthonormal functions on the hypothesised distribution. This result is used to construct smooth tests for the bivariate Poisson distribution. The test due to Crockett (1979) is similar to a smooth test that assesses the variance structure under the bivariate Poisson model; the test due to Loukas & Kemp (1986) is related to a smooth test that seeks to detect a particular linear relationship between the variances and covariance under the bivariate Poisson model. Using focused smooth tests may be more informative than using previously suggested tests. The distribution of the Loukas & Kemp (1986) statistic is not well approximated by the x²distribution for larger correlations, and a revised statistic is suggested.     

INFERENCE ABOUT THE MEAN OF A POISSON DISTRIBUTION IN THE PRESENCE OF A NUISANCE PARAMETER

In this note we examine the problem of estimating the mean of a Poisson distribution when a nuisance parameter is present. Using a condition of Cox (1958) about ancillarity in the presence of a nuisance parameter, we justify that inference about the parameter should be carried out using the conditional distribution given the appropriate ancillary statistics. A small simulation study has been done to compare the performance of the conditional likelihood approach and the standard likelihood approach.     

ESTIMATION OF THE PARAMETER OF A POISSON DISTRIBUTION USING A LINEX LOSS FUNCTION

This paper considers estimation of the parameter of a Poisson distribution using Varian's (1975) asymmetric LINEX loss function L (δ) = b{exp(aδ) - aδ - 1}, where δ is the estimation error and b > 0, a 0. It is shown that for a < 0, the sample mean X¯ is admissible whereas for a > 0, X¯ is dominated by c*X¯, where c*= (n/a)log(1+a/n). Practical implications of this result are indicated. More general results, concerning the admissibility of estimators of the form cX¯+ d are also presented.     

POISSON VERSUS BINOMIAL: APPOINTMENT OF JUDGES TO THE U.S. SUPREME COURT

The problem of discriminating between the Poisson and binomial models is discussed in the context of a detailed statistical analysis of the number of appointments of the U.S. Supreme Court justices from 1789 to 2004. Various new and existing tests are examined. The analysis shows that both simple Poisson and simple binomial models are equally appropriate for describing the data. No firm statistical evidence in favour of an exponential Poisson regression model was found. Two attendant results were obtained by simulation: firstly, that the likelihood ratio test is the most powerful of those considered when testing for the Poisson versus binomial and, secondly, that the classical variance test with an upper-tail critical region is biased.     

SOME PROPERTIES OF A RANDOM LINEAR DIFFERENCE EQUATION1

Some simple conditions are given for the absolute continuity of the limiting distribution of a random linear difference equation. These results are applied to the super-critical Bellman-Harris branching process with immigration. When the coefficients of the difference equation are non-negative and there is no limiting distribution, it is shown that the asymptotic behaviour of the solutions is the same as that of the partial sums of a divergent random power series. A number of limit theorems are given for the latter situation.     

THE RELATIONSHIPS BETWEEN SKEWNESS AND KURTOSIS

Theoretical considerations of kurtosis, whether of partial orderings of distributions with respect to kurtosis or of measures of kurtosis, have tended to focus only on symmetric distributions. With reference to historical points and recent work on skewness and kurtosis, this paper defines anti-skewness and uses it as a tool to discuss the concept of kurtosis in asymmetric univariate distributions. The discussion indicates that while kurtosis is best considered as a property of symmetrised versions of distributions, symmetrisation does not simply remove skewness. Skewness, anti-skewness and kurtosis are all inter-related aspects of shape. The Tukey g and h family and the Johnson Su family are considered as examples.     

THE CHOICE BETWEEN SETS OF REGRESSORS

This paper examines the choice of critical values for testing both non-sequential and nested sequential sets of constraints in the standard linear regression model. Modest increases in (e.g.) t-ratio critical values relative to their one-off values are often sufficient to maintain proper size. A Bayesian decision-theoretic approach, highlighted by the Schwarz (1978) criterion, provides a framework for deriving consistency and asymptotic local power properties of both forms of testing (data mining) algorithms.     

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.     

A TEST FOR THE NONPARAMETRIC TWO SAMPLE SCALE PROBLEM

We propose a distribution-free test for the nonparametric two sample scale problem. Unlike the other tests for this problem, we do not assume that the two distribution functions have a common median. We assume that they have a common quantile of order a (not necessarily 1/2). The test statistic is a modification of the Sukhatme statistic for the scale problem and the Wilcoxon-Mann-Whitney statistic for stochastic dominance. It is shown that the new test is uniformly more efficient (in the Pitman sense) than the Sukhatme test and has very good efficiency when compared to the Mood test.     

ON THE USE OF X2 AS A TEST OF INDEPENDENCE IN CONTINGENCY TABLES WITH SMALL CELL EXPECTATIONS

When an r×c contingency table has many cells having very small expectations, the usual χ² approximation to the upper tail of the Pearson χ² goodness-of-fit statistic becomes very conservative. The alternatives considered in this paper are to use either a lognormal approximation, or to scale the usual χ² approximation. The study involves thousands of tables with various sample sizes, and with tables whose sizes range from 2×2 through 2×10×10. Subject to certain restrictions the new scaled χ² approximations are recommended for use with tables having an average cell expectation as small as 0·5.