Estimating the Incidence of Hidden,Untreated Opiate Use

Illegal drug use is a hidden phenomenon, making it extremely difficult to obtain accurate estimates of untreated use. Yet it is this very estimate of incidence that is required by policy makers and service providers who are responsible for allocating increasingly scarce resources in times of worldwide recession. Using Irish data on more than 6,000 cases of clients presenting for first treatment, a Gamma incubation period distribution describing progression to first opiate treatment is fitted. An analytical solution of the back calculation method and a linear Volterra integral equation are used to project back and estimate the hidden, untreated population of opiate users.     

Minutes of Committee of Presidents of Statistical Societies

The general approach to generating random variates through transformations with multiple roots is discussed. Multinomial probabilities are determined for the selection of the different roots. An application of the general result yields a new and simple technique for the generation of variates from the inverse Gaussian distribution.     

Evaluating and improving risk formulas for allocating limited budgets to expensive risk-reduction opportunities

Simple risk formulas, such as risk = probability × impact, or risk = exposure × probability × consequence, or risk = threat × vulnerability × consequence, are built into many commercial risk management software products deployed in public and private organizations. These formulas, which we call risk indices, together with risk matrices, “heat maps,” and other displays based on them, are widely used in applications such as enterprise risk management (ERM), terrorism risk analysis, and occupational safety. But, how well do they serve to guide allocation of limited risk management resources? This article evaluates and compares different risk indices under simplifying conditions favorable to their use (statistically independent, uniformly distributed values of their components; and noninteracting risk‐reduction opportunities). Compared to an optimal (nonindex) approach, simple indices produce inferior resource allocations that for a given cost may reduce risk by as little as 60% of what the optimal decisions would provide, at least in our simple simulations. This article suggests a better risk reduction per unit cost index that achieves 98–100% of the maximum possible risk reduction on these problems for all budget levels except the smallest, which allow very few risks to be addressed. Substantial gains in risk reduction achieved for resources spent can be obtained on our test problems by using this improved index instead of simpler ones that focus only on relative sizes of risk (or of components of risk) in informing risk management priorities and allocating limited risk management resources. This work suggests the need for risk management tools to explicitly consider costs in prioritization activities, particularly in situations where budget restrictions make careful allocation of resources essential for achieving close‐to‐maximum risk‐reduction benefits.     

Study of redundancy of vector variables in canonical correlations

Fujikoshi (1982) obtained the necessary and sufficient conditions for the increased number of variables in the two sets of vectors not affecting the original nonzero canonical correlations and used these to obtain the likelihood ratio test procedure. He assumed a nonsingular covariance matrix due to random variables. Here, we study the same problem when the covariance matrix is singular and establish some further results. In this study, we note that the unit canonical correlations have to be separated in some of the situations. These results are valid for complex random vector variables and in some situations, the test for redundancy is given for complex random variables.     

A solution to the multivariate behrens-fisher problem

Some new algebra on pattern and transition matrices is used to determine the degrees of freedom and the parameter matrix, if the distribution of a linear sum of Wishart matrices is approximated by a single Wishart distribution. This approximation is then used to find a solution to the multivariate Behrens-Fisher problem similar to the Welch (1947) solution in the univariate case.     

On the distributions of the inner distance functions in the multivariate normal populations with unequal covariance matrices

The sample distance functions between an observation and a population were deduced by the likelihood procedures for discrimination problem in the case of several normal populations with unequal covariance matrices(1986). The present paper gives the exact MGFs of the distance functions for the case that the observation and the sample come from the same population and the limiting distributions of the distance functions by using the MCFs.     

A characterization of the independence - distribution - preserving covariance structure for the multivariate maximum squared - radii statistic

Necessary and sufficient conditions on the observation covariance structure and on the set of linear transformations are given for which the distribution of the multivariate maximum squared - radii statistic for detecting a single multivariate outlier is invariant from the distribution assuming the usual independence covariance structure. Thus, we extend the work of Baksalary and Puntanen (1990), who have given necessary and sufficient conditions for an independence-distribution-preserving covariance structure for Grubbs' statistic for detecting a univariate outlier. We also extend the work of Marco, Young, and Turner (1987) and Pavur and Young (1991), who have given sufficient conditions for an independence-distribution-preserving dependency structure for the multivariate squared - radii statistic.     

A Note on Sufficient Conditions for Valid Unmodified t Testing in Correlation Analysis with Autocorrelated and Heteroscedastic Sample Data

Traditionally, sphericity (i.e., independence and homoscedasticity for raw data) is put forward as the condition to be satisfied by the variance–covariance matrix of at least one of the two observation vectors analyzed for correlation, for the unmodified t test of significance to be valid under the Gaussian and constant population mean assumptions. In this article, the author proves that the sphericity condition is too strong and a weaker (i.e., more general) sufficient condition for valid unmodified t testing in correlation analysis is circularity (i.e., independence and homoscedasticity after linear transformation by orthonormal contrasts), to be satisfied by the variance–covariance matrix of one of the two observation vectors. Two other conditions (i.e., compound symmetry for one of the two observation vectors; absence of correlation between the components of one observation vector, combined with a particular pattern of joint heteroscedasticity in the two observation vectors) are also considered and discussed. When both observation vectors possess the same variance–covariance matrix up to a positive multiplicative constant, the circularity condition is shown to be necessary and sufficient. “Observation vectors” may designate partial realizations of temporal or spatial stochastic processes as well as profile vectors of repeated measures. From the proof, it follows that an effective sample size appropriately defined can measure the discrepancy from the more general sufficient condition for valid unmodified t testing in correlation analysis with autocorrelated and heteroscedastic sample data. The proof is complemented by a simulation study. Finally, the differences between the role of the circularity condition in the correlation analysis and its role in the repeated measures ANOVA (i.e., where it was first introduced) are scrutinized, and the link between the circular variance–covariance structure and the centering of observations with respect to the sample mean is emphasized.     

On Comparing Cumulative Incidence Functions Using an Empirical Likelihood Ratio Type Test

We consider the competing risks problem with two risks and develop empirical likelihood ratio type tests for testing the null hypothesis that the cumulative incidence functions corresponding to these two risks are equal against the alternatives: (a) they are not equal and (b) they are linearly ordered. The asymptotic null distributions of the proposed test statistics are shown to have simple distribution-free representations in terms of a standard Brownian motion process. The results of a simulation study indicate that the proposed test for testing for the presence of the linear order is more powerful than a test designed for the same situation in Aly et al. (1994 Aly, E.A.A., Kochar, S.C., McKeague, I.W. (1994). Some tests for comparing cumulative incidence functions and cause-specific hazard rates. J. Am. Stat. Assoc. 89:994–999.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). To illustrate the theoretical results, we discuss an example involving survival times of mice exposed to radiation.     

Convergence of Sample Eigenvectors of Spiked Population Model

In this paper, for the general non Gaussian spiked population model, where a few fixed eigenvalues of the population covariance matrix are separated from others, we investigate the convergence properties of the eigenvectors of sample covariance matrices corresponding to the spiked population eigenvalues and angle between the population eigenvectors and sample eigenvectors as both the sample size and population size are large.