Statistical consulting

The paper and the special issue focus on the activity of statistical consulting and its varieties. This includes academic consulting, consulting to and in industry as well as statistics in public media.     

Multi-phase sampling under cost constraints

Multi-phase sampling (M-PhS) scheme is useful when the interest is in the estimation of the population mean of an expensive variable strictly connected with other cheaper (auxiliary) variables. The MSE is an accuracy measure of an estimator. Usually it decreases as the sample size increases. In practice the sample size cannot become arbitrarily large for possible cost constraints. From a practical point of view it would be useful to know the sample sizes which guarantee the best accuracy of the estimates for fixed costs. These “optimum” sample sizes can be, in some cases, computable but not admissible. In other cases, they can be neither admissible nor computable. The main goal of this paper is to propose a solution for both these situations. It will be clear that in both situations the solution is to consider a M-PhS scheme with one or more phases less.     

H-likelihood: problems and solutions

In recent issues of this journal it has been asserted in two papers that the use of h-likelihood is wrong, in the sense of giving unsatisfactory estimates of some parameters for binary data (Kuk and Cheng, 1999; Waddington and Thompson, 2004) or theoretically unsound (Kuk and Cheng, 1999). We wish to refute both these assertions.     

Nonconglomerative coherent conditional probabilities in statistical inference

Noteworthy connections among conglomerability, countable additivity and coherence are discussed in detail, reaching the conclusion that nonconglomerable conditional probabilities must not be doomed and play a significant role in statistical inference. Extended and updated version of a contributed paper presented at the International Conference on “Information Processing and Management of Uncertainty in knowledge-based systems”, IPMU 2004, Perugia, Italy.     

The multivariate alpha-power model

The main object of this paper is to propose a multivariate extension to the alpha-power model which is an alternative to the multivariate skew-normal model (Arellano-Valle and Azzalini, 2008). It also extends the power-normal model discussed in Gupta and Gupta (2008) by making it more flexible. Inference is dealt with by using the likelihood approach and a pseudo-likelihood approach based on conditional distributions which, although slightly less efficient, is simpler to implement. An application to a real data set is used to demonstrate the usefulness of the extension.     

Skew-normal alpha-power model

In this paper, authors study properties and inference for the newly introduced skew-normal alpha-power model, generalizing both, the power-normal and skew-normal models. Inference is approached via maximum likelihood. Fisher information matrix is derived and shown to be nonsingular at the whole parametric space. Special emphasis is placed on the special case of the power–skew-normal model. Studies with real data illustrate the fact that the model can be very useful in applications, being able to overfit less general models entertained in the literature.     

Asymptotic relative efficiency of wald tests in measurement error models

In this paper, asymptotic relative efficiency (ARE) of Wald tests for the Tweedie class of models with log-linear mean, is considered when the aux¬iliary variable is measured with error. Wald test statistics based on the naive maximum likelihood estimator and on a consistent estimator which is obtained by using Nakarnura's (1990) corrected score function approach are defined. As shown analytically, the Wald statistics based on the naive and corrected score function estimators are asymptotically equivalents in terms of ARE. On the other hand, the asymptotic relative efficiency of the naive and corrected Wald statistic with respect to the Wald statistic based on the true covariate equals to the square of the correlation between the unobserved and the observed co-variate. A small scale numerical Monte Carlo study and an example illustrate the small sample size situation.     

Likelihood-Based Inference for Multivariate Skew-Normal Regression Models

In this article, we present EM algorithms for performing maximum likelihood estimation for three multivariate skew-normal regression models of considerable practical interest. We also consider the restricted estimation of the parameters of certain important special cases of two models. The methodology developed is applied in the analysis of longitudinal data on dental plaque and cholesterol levels.     

Equivariant prediction in finite populations

A theory of equivariant prediction is developed for predicting the population total in finite populations. Minimum risk equivariant predictors (MREP) are derived under the location, scale and locationscale superpopulation models. Under the general linear model, it is shown that the best(linear) unbiased predictor (B(L)UP) is an MREP.     

Finite population prediction under a linear function superpopulation model:a bayesian perspective

Exact and approximate Bayesian inference is developed for the prediction problem in finite populations under a linear functional superpopulation model. The models considered are the usual regression models involving two variables, X and Y, where the independent variable X is measured with error. The approach is based on the conditional distribution of Y given X and our predictor is the posterior mean of the quantity of interest (population total and population variance) given the observed data. Empirical investigations about optimal purposive samples and possible model misspecifications based on comparisons with the corresponding models where X is measured without error are also reported.