A new generalized two-sided class of distributions with an emphasis on two-sided generalized normal distribution

The ordinary-G class of distributions is defined to have the cumulative distribution function (cdf) as the value of the cdf of the ordinary distribution F whose range is the unit interval at G, that is, F(G), and it generalizes the ordinary distribution. In this work, we consider the standard two-sided power distribution to define other classes like the beta-G and the Kumaraswamy-G classes. We extend the idea of two-sidedness to other ordinary distributions like normal. After studying the basic properties of the new class in general setting, we consider the two-sided generalized normal distribution with maximum likelihood estimation procedure.     

Disaster response on September 11, 2001 through the lens of statistical network analysis

The rescue and relief operations triggered by the September 11, 2001 attacks on the World Trade Center in New York City demanded collaboration among hundreds of organisations. To shed light on the response to the September 11, 2001 attacks and help to plan and prepare the response to future disasters, we study the inter-organisational network that emerged in response to the attacks. Studying the inter-organisational network can help to shed light on (1) whether some organisations dominated the inter-organisational network and facilitated communication and coordination of the disaster response; (2) whether the dominating organisations were supposed to coordinate disaster response or emerged as coordinators in the wake of the disaster; and (3) the degree of network redundancy and sensitivity of the inter-organisational network to disturbances following the initial disaster. We introduce a Bayesian framework which can answer the substantive questions of interest while being as simple and parsimonious as possible. The framework allows organisations to have varying propensities to collaborate, while taking covariates into account, and allows to assess whether the inter-organisational network had network redundancy—in the form of transitivity—by using a test which may be regarded as a Bayesian score test. We discuss implications in terms of disaster management.     

A Bayesian approach to non-parametric monotone function estimation

Summary.  The paper proposes two Bayesian approaches to non-parametric monotone function estimation. The first approach uses a hierarchical Bayes framework and a characterization of smooth monotone functions given by Ramsay that allows unconstrained estimation. The second approach uses a Bayesian regression spline model of Smith and Kohn with a mixture distribution of constrained normal distributions as the prior for the regression coefficients to ensure the monotonicity of the resulting function estimate. The small sample properties of the two function estimators across a range of functions are provided via simulation and compared with existing methods. Asymptotic results are also given that show that Bayesian methods provide consistent function estimators for a large class of smooth functions. An example is provided involving economic demand functions that illustrates the application of the constrained regression spline estimator in the context of a multiple-regression model where two functions are constrained to be monotone.     

In Search of the Optimum Harvey Wallbanger Recipe via Mixture Experiment Techniques

A Harvey Wallbanger is an alcoholic drink composed of three ingredients: orange juice, vodka, and Galliano. One of the authors (Sahrmann), a connoisseur of the libation, was perplexed upon discovering that different bartender handbooks listed different recipes for the drink. Determined to shed some light on this irksome problem, he issued invitations for an evening workshop session (party) to former and current attendees of an experimental-design training course taught within his company. A mixture experiment was designed and data collected using the class participants. This article describes the experimental design and the analysis of the data. The results of the experiment indicate that reducing the proportion of vodka produces better-flavored drinks. The results of the experiment were not entirely conclusive, however, so recommendations are made for future experimentation.     

Model-Based Classification via Mixtures of Multivariate t-Factor Analyzers

A model-based classification technique is developed, based on mixtures of multivariate t-factor analyzers. Specifically, two related mixture models are developed and their classification efficacy studied. An AECM algorithm is used for parameter estimation, and convergence of these algorithms is determined using Aitken's acceleration. Two different techniques are proposed for model selection: the BIC and the ICL. Our classification technique is applied to data on red wine samples from Italy and to fatty acid measurements on Italian olive oils. These results are discussed and compared to more established classification techniques; under this comparison, our mixture models give excellent classification performance.     

Examples of Uncorrelated Dependent Random Variables Using a Bivariate Mixture

A particular mixture of bivariate distributions is used to present examples of dependent uncorrelated random variables and independent random variables. A necessary and sufficient condition for the independence for such a bivariate distribution is given.     

Global Robustness with Respect to the Loss Function and the Prior

We propose a class [I,S] of loss functions for modeling the imprecise preferences of the decision maker in Bayesian Decision Theory. This class is built upon two extreme loss functions I and S which reflect the limited information about the loss function. We give an approximation of the set of Bayes actions for every loss function in [I,S] and every prior in a mixture class; if the decision space is a subset of , we obtain the exact set.     

Estimation of scale parameters in mixture distributions

Simultaneous estimation of scale parameters is considered in mixture distributions under squared-error loss. A general class of estimators is obtained which dominates the componentwise best multiple estimators and the moment estimators. As special cases, improved estimators are obtained for the multivariate t-distribution and the p-variate Lomax distribution.     

Model-based clustering for social networks

Summary.  Network models are widely used to represent relations between interacting units or actors. Network data often exhibit transitivity, meaning that two actors that have ties to a third actor are more likely to be tied than actors that do not, homophily by attributes of the actors or dyads, and clustering. Interest often focuses on finding clusters of actors or ties, and the number of groups in the data is typically unknown. We propose a new model, the latent position cluster model , under which the probability of a tie between two actors depends on the distance between them in an unobserved Euclidean 'social space', and the actors' locations in the latent social space arise from a mixture of distributions, each corresponding to a cluster. We propose two estimation methods: a two-stage maximum likelihood method and a fully Bayesian method that uses Markov chain Monte Carlo sampling. The former is quicker and simpler, but the latter performs better. We also propose a Bayesian way of determining the number of clusters that are present by using approximate conditional Bayes factors. Our model represents transitivity, homophily by attributes and clustering simultaneously and does not require the number of clusters to be known. The model makes it easy to simulate realistic networks with clustering, which are potentially useful as inputs to models of more complex systems of which the network is part, such as epidemic models of infectious disease. We apply the model to two networks of social relations. A free software package in the R statistical language, latentnet, is available to analyse data by using the model.     

Bayesian inference and model comparison for asymmetric smooth transition heteroskedastic models

Inference, quantile forecasting and model comparison for an asymmetric double smooth transition heteroskedastic model is investigated. A Bayesian framework in employed and an adaptive Markov chain Monte Carlo scheme is designed. A mixture prior is proposed that alleviates the usual identifiability problem as the speed of transition parameter tends to zero, and an informative prior for this parameter is suggested, that allows for reliable inference and a proper posterior, despite the non-integrability of the likelihood function. A formal Bayesian posterior model comparison procedure is employed to compare the proposed model with its two limiting cases: the double threshold GARCH and symmetric ARX GARCH models. The proposed methods are illustrated using both simulated and international stock market return series. Some illustrations of the advantages of an adaptive sampling scheme for these models are also provided. Finally, Bayesian forecasting methods are employed in a Value-at-Risk study of the international return series. The results generally favour the proposed smooth transition model and highlight explosive and smooth nonlinear behaviour in financial markets.