Noninformative priors for linear combinations of normal means with unequal variances

For normal populations with unequal variances, we develop matching priors and reference priors for a linear combination of the means. Here, we find three second-order matching priors: a highest posterior density (HPD) matching prior, a cumulative distribution function (CDF) matching prior, and a likelihood ratio (LR) matching prior. Furthermore, we show that the reference priors are all first-order matching priors, but that they do not satisfy the second-order matching criterion that establishes the symmetry and the unimodality of the posterior under the developed priors. The results of a simulation indicate that the second-order matching prior outperforms the reference priors in terms of matching the target coverage probabilities, in a frequentist sense. Finally, we compare the Bayesian credible intervals based on the developed priors with the confidence intervals derived from real data.     

Choice of the ridge factor from the correlation matrix determinant

Ridge regression is the alternative method to ordinary least squares, which is mostly applied when a multiple linear regression model presents a worrying degree of collinearity. A relevant topic in ridge regression is the selection of the ridge parameter, and different proposals have been presented in the scientific literature. Since the ridge estimator is biased, its estimation is normally based on the calculation of the mean square error (MSE) without considering (to the best of our knowledge) whether the proposed value for the ridge parameter really mitigates the collinearity. With this goal and different simulations, this paper proposes to estimate the ridge parameter from the determinant of the matrix of correlation of the data, which verifies that the variance inflation factor (VIF) is lower than the traditionally established threshold. The possible relation between the VIF and the determinant of the matrix of correlation is also analysed. Finally, the contribution is illustrated with three real examples.     

A statistical procedure for representing state fragility and transition paths

State fragility is a concept that entered the political discourse in the last decades producing remarkable implications for aid allocation and international policies. The operationalization of this concept has generated a number of composite indices to produce rankings of fragile states. However, the temporal dimension of the driving forces leading to fragility has been rather neglected. This article discusses a statistical procedure that helps to represent the global fragility of a country and the path that a country has followed or will follow in the future when possibly entering into (or escaping from) a fragility condition. Specifically, multiple factor analysis is applied to depict vulnerable and weak countries, and to identify the fundamental forces that determine their overall fragility. Moreover, the trajectories of countries along the years are estimated using partial factor scores. Finally, the path of each country is predicted by means of parsimonious regression models, based on a reduced set of explanatory variables, and according to scenarios elaborated from available international outlooks.     

Mixtures of Gaussian copula factor analyzers for clustering high dimensional data

Mixtures of factor analyzers is a useful model-based clustering method which can avoid the curse of dimensionality in high-dimensional clustering. However, this approach is sensitive to both diverse non-normalities of marginal variables and outliers, which are commonly observed in multivariate experiments. We propose mixtures of Gaussian copula factor analyzers (MGCFA) for clustering high-dimensional clustering. This model has two advantages; (1) it allows different marginal distributions to facilitate fitting flexibility of the mixture model, (2) it can avoid the curse of dimensionality by embedding the factor-analytic structure in the component-correlation matrices of the mixture distribution.An EM algorithm is developed for the fitting of MGCFA. The proposed method is free of the curse of dimensionality and allows any parametric marginal distribution which fits best to the data. It is applied to both synthetic data and a microarray gene expression data for clustering and shows its better performance over several existing methods.     

A Note on Estimation of a Distribution Function in a Nonparametric Set-up Using Stein’s Shrinkage Estimation Technique

Based on Stein’s famous shrinkage estimation of a multivariate normal distribution, we propose a new type of estimators of the distribution function of a random variable in a nonparametric setup. The proposed estimators are then compared with the empirical distribution function, which is the best equivariant estimator under a well-known loss function. Our extensive simulation study shows that our proposed estimators can perform better for moderate to large sample sizes.     

Book Reviews

Books reviewed:
David Griffiths, W. Douglas Stirling, and K. Laurence Weldon, Understanding Data: Principles and Practice of Statistics
Ingwer Borg and Patrick Groenen, Modern Multidimensional Scaling: Theory and Applications
Jeffrey H. Dorfman, Bayesian Economics Through Numerical Methods: A Guide to Econometrics and Decision-making with Prior Information
Marek Musiela and Marek Rutkowski, Martingale Methods in Financial Modelling: Theory and Applications
Aad W. van der Vaart and Jon A. Wellner, Weak Convergence and Empirical Processes     

Inference for general parametric functions in box‐Cox‐type transformation models

The authors propose a simple but general method of inference for a parametric function of the Box‐Cox‐type transformation model. Their approach is built upon the classical normal theory but takes parameter estimation into account. It quickly leads to test statistics and confidence intervals for a linear combination of scaled or unsealed regression coefficients, as well as for the survivor function and marginal effects on the median or other quantité functions of an original response. The authors show through simulations that the finite‐sample performance of their method is often superior to the delta method, and that their approach is robust to mild departures from normality of error distributions. They illustrate their approach with a numerical example.     

Model-based measurement of latent risk in time series with applications

Summary.  Risk is at the centre of many policy decisions in companies, governments and other institutions. The risk of road fatalities concerns local governments in planning countermeasures, the risk and severity of counterparty default concerns bank risk managers daily and the risk of infection has actuarial and epidemiological consequences. However, risk cannot be observed directly and it usually varies over time. We introduce a general multivariate time series model for the analysis of risk based on latent processes for the exposure to an event, the risk of that event occurring and the severity of the event. Linear state space methods can be used for the statistical treatment of the model. The new framework is illustrated for time series of insurance claims, credit card purchases and road safety. It is shown that the general methodology can be effectively used in the assessment of risk.     

Sex hormones and choice under risk

We study the correlation of choice under risk in Holt–Laury lotteries for gains and losses with gender, the use of hormonal contraceptives, menstrual cycle information, salivary testosterone, estradiol, progesterone, and cortisol as well as the digit ratio (2D:4D; length of the index finger to the ring finger of the right hand) in more than 200 subjects (45% females). In males, salivary testosterone is negatively correlated with risk aversion for gains only. In females, salivary cortisol is positively correlated with risk aversion for gains only. No other significant correlations between risk preferences and salivary hormones are observed. No significant correlations between risk preferences and the menstrual cycle are observed in naturally cycling females. No significant correlations between risk preferences and the digit ratio are observed in either gender and/or race.