Designing and integrating composite networks for monitoring multivariate gaussian pollution fields

Networks of ambient monitoring stations are used to monitor environmental pollution fields such as those for acid rain and air pollution. Such stations provide regular measurements of pollutant concentrations. The networks are established for a variety of purposes at various times so often several stations measuring different subsets of pollutant concentrations can be found in compact geographical regions. The problem of statistically combining these disparate information sources into a single 'network' then arises. Capitalizing on the efficiencies so achieved can then lead to the secondary problem of extending this network. The subject of this paper is a set of 31 air pollution monitoring stations in southern Ontario. Each of these regularly measures a particular subset of ionic sulphate, sulphite, nitrite and ozone. However, this subset varies from station to station. For example only two stations measure all four. Some measure just one. We describe a Bayesian framework for integrating the measurements of these stations to yield a spatial predictive distribution for unmonitored sites and unmeasured concentrations at existing stations. Furthermore we show how this network can be extended by using an entropy maximization criterion. The methods assume that the multivariate response field being measured has a joint Gaussian distribution conditional on its mean and covariance function. A conjugate prior is used for these parameters, some of its hyperparameters being fitted empirically.     

一种基于加权F-范数的半正定矩阵的逼近方法

正定性是许多金融预测模型的重要假设前提，然而从实际样本中得到的相关系数矩阵并不能保证其正定性。为此在介绍如何根据样本设定相关系数矩阵以及范数逼近原理的基础上，如何根据该原理找到与之最接近的相关系数矩阵，即最接近的单位对角半正定对称矩阵。通过实证，验证了其方法的有效性。     

A permutation-based Bayesian approach for inverse covariance estimation

Abstract

Covariance estimation and selection for multivariate datasets in a high-dimensional regime is a fundamental problem in modern statistics. Gaussian graphical models are a popular class of models used for this purpose. Current Bayesian methods for inverse covariance matrix estimation under Gaussian graphical models require the underlying graph and hence the ordering of variables to be known. However, in practice, such information on the true underlying model is often unavailable. We therefore propose a novel permutation-based Bayesian approach to tackle the unknown variable ordering issue. In particular, we utilize multiple maximum a posteriori estimates under the DAG-Wishart prior for each permutation, and subsequently construct the final estimate of the inverse covariance matrix. The proposed estimator has smaller variability and yields order-invariant property. We establish posterior convergence rates under mild assumptions and illustrate that our method outperforms existing approaches in estimating the inverse covariance matrices via simulation studies.     

相对论性理想流体熵密度的计算

求解自引力非辐射体系的熵结构，这个问题很有意义．本文一般地导出相对论性流体热力学第二定律的协变表达式．利用理想流体的熵流守恒和局域能量守恒等基本守恒律推导出熵密度的计算公式．     

Spatio-temporal change-point modeling

There is by now a substantial literature on spatio-temporal modeling. However, to date, there exists essentially no literature which addresses the issue of process change from a certain time. In fact, if we look at change points for purely time series data, the customary form is to propose a model involving a mean or level shift. We see little attempting to capture a change in association structure. Part of the concern is how to specify flexible ways to bridge the association across the time point and still ensure that a proper joint distribution has been defined for all of the data. Introducing a spatial component evidently adds further complication. We want to allow for a change-point reflecting change in both temporal and spatial association. In this paper we propose a constructive, flexible model formulation through additive specifications. We also demonstrate how computational concerns benefit from the availability of temporal order. Finally, we illustrate with several simulated datasets to examine the capability of the model to detect different types of structural changes.     

Some approaches to statistical analysis of a weather modification experiment

Data from a weather modification experiment are examined and a number of statistical analyses reported. The validity of earlier inferences is studied as are the utilities of various statistical methods. The experiment is described. The original analysis of North American Weather Consultants, who conducted the experiment, is reviewed. Data summarization is reported. A major approach to analysis is through the use of cloud-physics covari-ates in regression analyses. Finally, a multivariate analysis is discussed. It appears that the covariates may have been affected by treatment (cloud seeding) and that their use is invalid, not only reducing error variances but removing treatment effect. Some recommendations for improved design of similar future experiments are given in a concluding section, including preliminary trial use of blocking by storms.     

Approximate Bayesian Inference in Spatial Generalized Linear Mixed Models

Abstract.  In this paper we propose fast approximate methods for computing posterior marginals in spatial generalized linear mixed models. We consider the common geostatistical case with a high dimensional latent spatial variable and observations at known registration sites. The methods of inference are deterministic, using no simulation-based inference. The first proposed approximation is fast to compute and is 'practically sufficient', meaning that results do not show any bias or dispersion effects that might affect decision making. Our second approximation, an improvement of the first version, is 'practically exact', meaning that one would have to run MCMC simulations for very much longer than is typically done to detect any indication of error in the approximate results. For small-count data the approximations are slightly worse, but still very accurate. Our methods are limited to likelihood functions that give unimodal full conditionals for the latent variable. The methods help to expand the future scope of non-Gaussian geostatistical models as illustrated by applications of model choice, outlier detection and sampling design. The approximations take seconds or minutes of CPU time, in sharp contrast to overnight MCMC runs for solving such problems.     

The maximum likelihood estimators in the growth curve model with serial covariance structure

The maximum likelihood estimators of unknown parameters in the growth curve model with serial covariance structure under some conditions are derived in the paper.