Australian & New Zealand Journal of Statistics最新文献

Various non-parametric goodness of fit tests have already been investigated in the literature. However, those tests are rarely used in the case of missing observations. We here study the goodness of fit test for missing data based on L_p distances along with Kolmogorov–Smirnov and Cramer–von-Mises distances when missingness occurs at random. The asymptotic distributions of the proposed test statistics have been derived under contiguous alternatives that enable us to investigate the asymptotic local power of the tests. We also study the performance of the tests for finite samples using simulation, and the tests perform well for those cases. The usefulness of the tests is illustrated on three real data sets.

In this paper, we analysed classification rules under Bayesian decision theory. The setup we considered here is fairly general, which can represent all possible parametric models. The Bayes classification rule we investigated minimises the Bayes risk under general loss functions. Among the existing literatures, the 0-1 loss function appears most frequently, under which the Bayes classification rule is determined by the posterior predictive densities. Theoretically, we extended the Bernstein–von Mises theorem to the multiple-sample case. On this basis, the oracle property of Bayes classification rule has been discussed in detail, which refers to the convergence of the Bayes classification rule to the one built from the true distributions, as the sample size tends to infinity. Simulations show that the Bayes classification rules do have some advantages over the traditional classifiers, especially when the number of features approaches the sample size.

We introduce new estimators of the inhomogeneous K-function and the pair correlation function of a spatial point process as well as the cross K-function and the cross pair correlation function of a bivariate spatial point process under the assumption of second-order intensity-reweighted stationarity. These estimators rely on a ‘global’ normalisation factor which depends on an aggregation of the intensity function, while the existing estimators depend ‘locally’ on the intensity function at the individual observed points. The advantages of our new global estimators over the existing local estimators are demonstrated by theoretical considerations and a simulation study.

In any other circumstance, it might make sense to define the extent of the terrain (Data Science) first, and then locate and describe the landmarks (Principles). But this data revolution we are experiencing defies a cadastral survey. Areas are continually being annexed into Data Science. For example, biometrics was traditionally statistics for agriculture in all its forms but now, in Data Science, it means the study of characteristics that can be used to identify an individual. Examples of non-intrusive measurements include height, weight, fingerprints, retina scan, voice, photograph/video (facial landmarks and facial expressions) and gait. A multivariate analysis of such data would be a complex project for a statistician, but a software engineer might appear to have no trouble with it at all. In any applied-statistics project, the statistician worries about uncertainty and quantifies it by modelling data as realisations generated from a probability space. Another approach to uncertainty quantification is to find similar data sets, and then use the variability of results between these data sets to capture the uncertainty. Both approaches allow ‘error bars’ to be put on estimates obtained from the original data set, although the interpretations are different. A third approach, that concentrates on giving a single answer and gives up on uncertainty quantification, could be considered as Data Engineering, although it has staked a claim in the Data Science terrain. This article presents a few (actually nine) statistical principles for data scientists that have helped me, and continue to help me, when I work on complex interdisciplinary projects.

The theoretical foundation for a number of model selection criteria is established in the context of inhomogeneous point processes and under various asymptotic settings: infill, increasing domain and combinations of these. For inhomogeneous Poisson processes we consider Akaike's information criterion and the Bayesian information criterion, and in particular we identify the point process analogue of ‘sample size’ needed for the Bayesian information criterion. Considering general inhomogeneous point processes we derive new composite likelihood and composite Bayesian information criteria for selecting a regression model for the intensity function. The proposed model selection criteria are evaluated using simulations of Poisson processes and cluster point processes.

We suggest several constructions suitable to define the depth of set-valued observations with respect to a sample of convex sets or with respect to the distribution of a random closed convex set. With the concept of a depth, it is possible to determine if a given convex set should be regarded an outlier with respect to a sample of convex closed sets. Some of our constructions are motivated by the known concepts of half-space depth and band depth for function-valued data. A novel construction derives the depth from a family of non-linear expectations of random sets. Furthermore, we address the role of positions of sets for evaluation of their depth. Two case studies concern interval regression for Greek wine data and detection of outliers in a sample of particles.

In survey sampling studies, statistical inference can be constructed either using design based randomisation or super population model. Design-based inference using judgement post-stratified (JPS) sampling is available in the literature. This paper develops statistical inference based on super population model in a finite population setting using JPS sampling design. For a JPS sample, first a simple random sample (SRS) is constructed without replacement. The sample units in this SRS are then stratified based on judgement ranking in a small comparison set to induce a data structure in the sample. The paper shows that the mean of a JPS sample is model unbiased and has smaller mean square prediction error (MSPE) than the MSPE of a simple random sample mean. Using an unbiased estimator of the MSPE, the paper also constructs prediction confidence interval for the population mean. A small-scale empirical study shows that the JPS sample predictor performs better than an SRS predictor when the quality of ranking information in JPS sampling is not poor. The paper also shows that the coverage probabilities of prediction intervals are very close to the nominal coverage probability. Proposed inferential procedure is applied to a real data set obtained from an agricultural research farm.

For modelling the location of pyramidal cells in the human cerebral cortex, we suggest a hierarchical point process in that exhibits anisotropy in the form of cylinders extending along the z-axis. The model consists first of a generalised shot noise Cox process for the xy-coordinates, providing cylindrical clusters, and next of a Markov random field model for the z-coordinates conditioned on the xy-coordinates, providing either repulsion, aggregation or both within specified areas of interaction. Several cases of these hierarchical point processes are fitted to two pyramidal cell data sets, and of these a final model allowing for both repulsion and attraction between the points seem adequate. We discuss how the final model relates to the so-called minicolumn hypothesis in neuroscience.

It was a major breakthrough when design-based stereological methods for vertical sections were developed by Adrian Baddeley and coworkers in the 1980s. Most importantly, it was shown how to estimate in a design-based fashion surface area from observations in random vertical sections with uniform position and uniform rotation around the vertical axis. The great practical importance of these developments is due to the fact that some biostructures can only be recognised on vertical sections. Later, local design-based estimation of mean particle volume from vertical sections was developed. In the present paper, we review these important advances in stereology. Quite recently, vertical sections have gained renewed interest, since it has been shown that mean particle shape can be estimated from such sections. These new developments are also reviewed in the present paper.