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.     

Measuring the complexity of generalized linear hierarchical models

Measuring a statistical model's complexity is important for model criticism and comparison. However, it is unclear how to do this for hierarchical models due to uncertainty about how to count the random effects. The authors develop a complexity measure for generalized linear hierarchical models based on linear model theory. They demonstrate the new measure for binomial and Poisson observables modeled using various hierarchical structures, including a longitudinal model and an areal‐data model having both spatial clustering and pure heterogeneity random effects. They compare their new measure to a Bayesian index of model complexity, the effective number p_D of parameters (Spiegelhalter, Best, Carlin & van der Linde 2002); the comparisons are made in the binomial and Poisson cases via simulation and two real data examples. The two measures are usually close, but differ markedly in some instances where p_D is arguably inappropriate. Finally, the authors show how the new measure can be used to approach the difficult task of specifying prior distributions for variance components, and in the process cast further doubt on the commonly‐used vague inverse gamma prior.     

Using Related Samples in Assessing Conformance to Safety Goals: A Nuclear Reactor Safety Application

Quantifying safety goals is a key to the regulation of activities which are beneficial on the whole but entail some risks in being performed. Determining compliance with safety goals involves dealing with uncertainties. A recent article by Bier^(I) describes some of the difficulties encountered using measures with uncertainty to determine compliance with safety goals for nuclear reactors. This paper uses a hierarchical Bayes approach to address two practical modeling problems in determining safety goal compliance under uncertainty: (1) allowing some modeling assumptions to be relaxed, and (2) allowing data from previous related samples to be included in the analysis. The two issues effect each other to the extent that relaxing some assumptions allows the use of a broader range of data. The usefulness of these changes and their impact on assessing safety compliance for nuclear reactors is shown.     

Describing the Dynamics of Attention to TV Commercials: A Hierarchical Bayes Analysis of the Time to Zap an Ad

This paper provides insights into the dynamics of attention to TV commercials via an analysis of the length of time that commercials are viewed before being 'zapped'. The model, which incorporates a flexible baseline hazard rate and captures unobserved heterogeneity across both consumers and commercials using a hierarchical Bayes approach, is estimated on two datasets in which commercial viewing is captured by a passive online device that continually monitors a household's TV viewing. Consistent with previous findings in psychology about the nature of attentional engagement in TV viewing, baseline hazard rates are found to be non-monotonic. In addition, the data show considerable ad-to-ad and household-to-household heterogeneity in zapping behavior. While one of the datasets contains some information on characteristics of the ads, these data do not reveal any firm links between the ad heterogeneity and the ad characteristics. A number of methodological and computational issues arise in the hierarchical Bayes analysis.     

减少商品混凝土现浇楼板早期收缩裂缝的构造措施

主要分析近年来一些工程中,现浇楼板出现早期收缩裂缝的比例有较大幅度增加的主要原因,从设计角度提出减少混凝土早期收缩裂缝的构造技术措施.     

Linear Discriminant Analysis of Multivariate Spatial–Temporal Regressions

Abstract.  We consider classification of the realization of a multivariate spatial–temporal Gaussian random field into one of two populations with different regression mean models and factorized covariance matrices. Unknown means and common feature vector covariance matrix are estimated from training samples with observations correlated in space and time, assuming spatial–temporal correlations to be known. We present the first-order asymptotic expansion of the expected error rate associated with a linear plug-in discriminant function. Our results are applied to ecological data collected from the Lithuanian Economic Zone in the Baltic Sea.     

Empirical Bayes estimation of parameters in Markov transition probability matrix with computational methods

Empirical Bayes estimator for the transition probability matrix is worked out in the cases where we have belief regarding the parameters, For example, where the states seem to be equal or not. In both cases, priors are in accordance with our beliefs. Using EM algorithm, computational methods for different hyperparameters of the empirical Bayes are described. Also, robustness of empirical Bayes procedure is investigated.     

Wavelet Shrinkage with Double Weibull Prior

In this article, we propose a denoising methodology in the wavelet domain based on a Bayesian hierarchical model using Double Weibull prior. We propose two estimators, one based on posterior mean (Double Weibull Wavelet Shrinker, DWWS) and the other based on larger posterior mode (DWWS-LPM), and show how to calculate them efficiently. Traditionally, mixture priors have been used for modeling sparse wavelet coefficients. The interesting feature of this article is the use of non-mixture prior. We show that the methodology provides good denoising performance, comparable even to state-of-the-art methods that use mixture priors and empirical Bayes setting of hyperparameters, which is demonstrated by extensive simulations on standardly used test functions. An application to real-word dataset is also considered.     

Probability matching priors for two-sided tolerance intervals in balanced one-way and two-way nested random effects models

We consider two-sided Bayesian tolerance intervals, with approximate frequentist validity, for a future observation in balanced one-way and two-way nested random effects models. Probability matching conditions, specific to this problem, are derived in either case via a technique that involves inversion of approximate posterior characteristic functions. In addition to yielding probability matching priors for the present problem, these conditions are useful in evaluating certain other priors that have received attention in the literature.     

Estimation in an empirical bayes model for longitudinal and cross‐sectionally clustered binary data

Some studies generate data that can be grouped into clusters in more than one way. Consider for instance a smoking prevention study in which responses on smoking status are collected over several years in a cohort of students from a number of different schools. This yields longitudinal data, also cross‐sectionaliy clustered in schools. The authors present a model for analyzing binary data of this type, combining generalized estimating equations and estimation of random effects to address the longitudinal and cross‐sectional dependence, respectively. The estimation procedure for this model is discussed, as are the results of a simulation study used to investigate the properties of its estimates. An illustration using data from a smoking prevention trial is given.