Modeling and Anomalous Cluster Detection for Point Processes Using Process Convolutions

We present a model using process convolutions, which describes spatial and temporal variations of the intensity of events that occur at random geographical locations. An inhomogeneous Poisson process is used to model the intensity over a spatial region with multiplicative spatial and temporal covariate effects. Temporal variation in the structure of the intensity is obtained by employing a time-varying process for the convolution. Use of a compactly supported kernel in the convolution improves the computational efficiency. Additionally, anomalous cluster detection in the event rates is developed based on exceedance probabilities. The methods are demonstrated on data of major crimes in Cincinnati during 2006. Supplementary materials for this article are available online.     

Fixed Rank Filtering for Spatio-Temporal Data

Datasets from remote-sensing platforms and sensor networks are often spatial, temporal, and very large. Processing massive amounts of data to provide current estimates of the (hidden) state from current and past data is challenging, even for the Kalman filter. A large number of spatial locations observed through time can quickly lead to an overwhelmingly high-dimensional statistical model. Dimension reduction without sacrificing complexity is our goal in this article. We demonstrate how a Spatio-Temporal Random Effects (STRE) component of a statistical model reduces the problem to one of fixed dimension with a very fast statistical solution, a methodology we call Fixed Rank Filtering (FRF). This is compared in a simulation experiment to successive, spatial-only predictions based on an analogous Spatial Random Effects (SRE) model, and the value of incorporating temporal dependence is quantified. A remote-sensing dataset of aerosol optical depth (AOD), from the Multi-angle Imaging SpectroRadiometer (MISR) instrument on the Terra satellite, is used to compare spatio-temporal FRF with spatial-only prediction. FRF achieves rapid production of optimally filtered AOD predictions, along with their prediction standard errors. In our case, over 100,000 spatio-temporal data were processed: Parameter estimation took 64.4 seconds and optimal predictions and their standard errors took 77.3 seconds to compute. Supplemental materials giving complete details on the design and analysis of a simulation experiment, the simulation code, and the MISR data used are available on-line.     

Classification of Mixtures of Spatial Point Processes via Partial Bayes Factors

Motivated by the problem of minefield detection, we investigate the problem of classifying mixtures of spatial point processes. In particular we are interested in testing the hypothesis that a given dataset was generated by a Poisson process versus a mixture of a Poisson process and a hard-core Strauss process. We propose testing this hypothesis by comparing the evidence for each model by using partial Bayes factors. We use the term partial Bayes factor to describe a Bayes factor, a ratio of integrated likelihoods, based on only part of the available information, namely that information contained in a small number of functionals of the data. We applied our method to both real and simulated data, and considering the difficulty of classifying these point patterns by eye, our approach overall produced good results.     

Zero Expectile Processes and Bayesian Spatial Regression

We introduce new classes of stationary spatial processes with asymmetric, sub-Gaussian marginal distributions using the idea of expectiles. We derive theoretical properties of the proposed processes. Moreover, we use the proposed spatial processes to formulate a spatial regression model for point-referenced data where the spatially correlated errors have skewed marginal distribution. We introduce a Bayesian computational procedure for model fitting and inference for this class of spatial regression models. We compare the performance of the proposed method with the traditional Gaussian process-based spatial regression through simulation studies and by applying it to a dataset on air pollution in California.     

Analysis of Data from a Series of Events by a Geometric Process Model

Geometric process was first introduced by Lam.A stochastic process {X_i,i=1,2,...} iscalled a geometric process (GP) if,for some a>0,{a~(i-1)X_i,i=1,2,...} forms a renewal process.In thispaper,the GP is used to analyze the data from a series of events.A nonparametric method is introduced forthe estimation of the three parameters in the GP.The limiting distributions of the three estimators are studied.Through the analysis of some real data sets,the GP model is compared with other three homogeneous andnonhomogeneous Poisson models.It seems that on average the GP model is the best model among these fourmodels in analyzing the data from a series of events.     

The product representation of a locally dependent random graph

Directed graphs with random black and white colourings of edges such that the colours of edges from different vertices are mutually independent are called locally dependent random graphs. Two random graphs are equivalent if they cannot be distinguished from percolation processes on them if only the vertices are seen. A necessary and sufficient condition is given for when a locally dependent random graph is equivalent to a product random graph; that is one in which the edges can be grouped in such a way that within each group the colours of the edges are equivalent and between groups they are independent. As an application the random graph corresponding to a spatial general epidemic model is considered.     

The method of approximate particular solutions to simulate an anomalous mobile-immobile transport process

In the current work, a generalized mathematical model based on the Coimbra time fractional derivative of variable order, which describes an anomalous mobile-immobile transport process in complex systems is investigated numerically. A robust numerical technique based on the meshfree strong form method combined with an efficient time-stepping scheme is performed to compute the approximate solution of the problem with high accuracy. For this purpose, firstly, an effective implicit time discretization approach is used for discretizing the variable-order time fractional problem in the time direction. Then a global meshless technique based on the method of approximate particular solutions is performed to fully discretize the model in the spatial domain. The validity and performance of the procedure to numerically simulate the proposed generalized solute transport model on regular and irregular domains are demonstrated through some numerical examples.     

A dynamical model for the overland flow on an infiltrating catchment

An approximate lumped parameter model for the surface runoff phenomenon on a catchment is presented. It consists of two submodels which are spatial discretizations of the basic partial differential equations of overland flow and infiltration. The aim of the model is to describe the dynamical input-output relationship between the rainfall rate and the surface runoff from the catchment. The boundary conditions of the infiltration process are modelled in an approximate way so that the number of state variables in the model can be reduced. With the aid of a single model parameter it is possible to describe simple catchment shapes like linearly converging and diverging surfaces. The model structure is flexible so that it can also be applied to more complex catchment configurations. The simulation results show that the phenomenon under consideration can be properly described by the model structure presented.     

Stochastic modelling of tropical cyclone tracks

A stochastic model for the tracks of tropical cyclones that allows for the computerised generation of a large number of synthetic cyclone tracks is introduced. This will provide a larger dataset than previously available for the assessment of risks in areas affected by tropical cyclones. To improve homogeneity, the historical tracks are first split into six classes. The points of cyclone genesis are modelled as a spatial Poisson point process, the intensity of which is estimated using a generalised version of a kernel estimator. For these points, initial values of direction, translation speed, and wind speed are drawn from histograms of the historical values of these variables observed in the neighbourhood of the respective points, thereby generating a first 6-h segment of a track. The subsequent segments are then generated by drawing changes in theses variables from histograms of the historical data available near the cyclone’s current location. A termination probability for the track is determined after each segment as a function of wind speed and location. In the present paper, the model is applied to historical cyclone data from the western North Pacific, but it is general enough to be transferred to other ocean basins with only minor adjustments. A version for the North Atlantic is currently under preparation.     

Return and Value at Risk using the Dirichlet Process

There exists a wide variety of models for return, and the chosen model determines the tool required to calculate the value at risk (VaR). This paper introduces an alternative methodology to model‐based simulation by using a Monte Carlo simulation of the Dirichlet process. The model is constructed in a Bayesian framework, using properties initially described by Ferguson. A notable advantage of this model is that, on average, the random draws are sampled from a mixed distribution that consists of a prior guess by an expert and the empirical process based on a random sample of historical asset returns. The method is relatively automatic and similar to machine learning tools, e.g. the estimate is updated as new data arrive.     

Hierarchical Bayesian autoregressive models for large space–time data with applications to ozone concentration modelling

Increasingly large volumes of space–time data are collected everywhere by mobile computing applications, and in many of these cases, temporal data are obtained by registering events, for example, telecommunication or Web traffic data. Having both the spatial and temporal dimensions adds substantial complexity to data analysis and inference tasks. The computational complexity increases rapidly for fitting Bayesian hierarchical models, as such a task involves repeated inversion of large matrices. The primary focus of this paper is on developing space–time autoregressive models under the hierarchical Bayesian setup. To handle large data sets, a recently developed Gaussian predictive process approximation method is extended to include autoregressive terms of latent space–time processes. Specifically, a space–time autoregressive process, supported on a set of a smaller number of knot locations, is spatially interpolated to approximate the original space–time process. The resulting model is specified within a hierarchical Bayesian framework, and Markov chain Monte Carlo techniques are used to make inference. The proposed model is applied for analysing the daily maximum 8‐h average ground level ozone concentration data from 1997 to 2006 from a large study region in the Eastern United States. The developed methods allow accurate spatial prediction of a temporally aggregated ozone summary, known as the primary ozone standard, along with its uncertainty, at any unmonitored location during the study period. Trends in spatial patterns of many features of the posterior predictive distribution of the primary standard, such as the probability of noncompliance with respect to the standard, are obtained and illustrated. Copyright © 2012 John Wiley & Sons, Ltd.     

截面数据插值的组合模型及其实证

采用Moran指标分析截面数据空间自相关性.针对具有空间自相关截面数据缺失插值问题,分别建立一阶空间自回归插值模型与克立格方法插值模型,在此基础上建立截面数据的组合插值模型,并用信息熵法确定组合插值模型的加权系数.用2003年福建部分市县城镇化水平的截面数据建立插值模型实证研究结果表明:组合插值模型的效果优于单项插值模型的效果.     

Influence of measurement data metrics on the identification of Interval Fields for the representation of spatial variability in finite element models

Due to the exponential growth in computing power, numerical modelling techniques method have gained an increasing amount of interest for engineering and design applications. Nowadays, the deterministic finite element (FE) method, an efficient tool to accurately solve the Partial Differential Equations (PDE) that govern most real-world problems, has become an indispensable tool for an engineer in various design stages. A more recent trend herein is to use the ever increasing computing power incorporate uncertainty and variability, which is omnipresent is all real-live applications, into these FE models. Several advanced techniques for incorporating either variability between nominally identical parts or spatial variability within one part into the FE models, have been introduced in this context. For the representation of spatial variability on the parameters of an FE model in a possibilistic context, the theory of Interval Fields (IF) was proven to show promising results. Following this approach, variability in the input FE model is introduced as the superposition of base vectors, depicting the spatial ‘patterns’, which are scaled by interval factors, which represent the actual variability. Application of this concept, however, requires identification of the governing parameters of these interval fields, i. e. the base vectors and interval scalars. Recent work of the authors therefore was focussed on finding a solution to the inverse problem, where the spatial uncertainty on the output side of the model is known from measurement data, but the spatial variability on the input parameters is unknown. Based on an a priori knowledge on the constituting base vectors of the interval field, the simulated output of the IFFEM computation is compared to measured data, and the input parameters are iteratively adjusted in order to minimize the discrepancy between the variability in simulation and measurement data. This discrepancy is defined based on geometric properties of the convex sets of both measurement and simulation data. However, the robustness of this methodology with respect to the size of the measurement data set that is used for the identification, as yet remains unclear. This paper therefore is focussed on the investigation of this robustness, by performing the identification on different measurement sets, depicting the same variability in the dynamic response of a simple FE model, which contain a decreasing amount of measurement replica. It was found that accurate identification remains feasible, even under a limited amount of measurement replica, which is highly relevant in the context of a non-probabilistic representation of variability in the FE model parameters. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Adaptive Bayesian Nonstationary Modeling for Large Spatial Datasets Using Covariance Approximations

Gaussian process models have been widely used in spatial statistics but face tremendous modeling and computational challenges for very large nonstationary spatial datasets. To address these challenges, we develop a Bayesian modeling approach using a nonstationary covariance function constructed based on adaptively selected partitions. The partitioned nonstationary class allows one to knit together local covariance parameters into a valid global nonstationary covariance for prediction, where the local covariance parameters are allowed to be estimated within each partition to reduce computational cost. To further facilitate the computations in local covariance estimation and global prediction, we use the full-scale covariance approximation (FSA) approach for the Bayesian inference of our model. One of our contributions is to model the partitions stochastically by embedding a modified treed partitioning process into the hierarchical models that leads to automated partitioning and substantial computational benefits. We illustrate the utility of our method with simulation studies and the global Total Ozone Matrix Spectrometer (TOMS) data. Supplementary materials for this article are available online.     

Modeling and Exploring Multivariate Spatial Variation: A Test Procedure for Isotropy of Multivariate Spatial Data

In this paper an exploratory technique based on the diagonalization of cross-variogram matrices is described. Through the definition of a model for the analysis and simulation of multivariate spatial data, a test procedure for the assumption of isotropy of multivariate spatial data is proposed. Applications to simulated and real data are reported.     

The distance-decay function of geographical gravity model: Power law or exponential law?

The distance-decay function of the geographical gravity model is originally an inverse power law, which suggests a scaling process in spatial interaction. However, the distance exponent of the model cannot be reasonably explained with the ideas from Euclidean geometry. This results in a dimension dilemma in geographical analysis. Consequently, a negative exponential function was used to replace the inverse power function to serve for a distance-decay function. But a new puzzle arose that the exponential-based gravity model goes against the first law of geography. This paper is devoted for solving these kinds of problems by mathematical reasoning and empirical analysis. New findings are as follows. First, the distance exponent of the gravity model is demonstrated to be a fractal dimension using the geometric measure relation. Second, the similarities and differences between the gravity models and spatial interaction models are revealed using allometric relations. Third, a four-parameter gravity model possesses a symmetrical expression, and we need dual gravity models to describe spatial flows. The observational data of China's cities and regions (29 elements indicative of 841 data points) in 2010 are employed to verify the theoretical inferences. A conclusion can be reached that the geographical gravity model based on power-law decay is more suitable for analyzing large, complex, and scale-free regional and urban systems. This study lends further support to the suggestion that the underlying rationale of fractal structure is entropy maximization. Moreover, it suggests that many dimensional dilemmas of spatial modeling can be solved using the concepts from fractal geometry.     

Spatial Regression Models,Response Surfaces,and Process Optimization

Abstract

Spatial regression models are developed as a complementary alternative to second-order polynomial response surfaces in the context of process optimization. These models provide estimates of design variable effects and smooth, data-faithful approximations to the unknown response function over the design space. The predicted response surfaces are driven by the covariance structures of the models. Several structures, isotropic and anisotropic, are considered and connections with thin plate splines are reviewed. Estimation of covariance parameters is achieved via maximum likelihood and residual maximum likelihood. A feature of the spatial regression approach is the visually appealing graphical summaries that are produced. These allow rapid and intuitive identification of process windows on the design space for which the response achieves target performance. Relevant design issues are briefly discussed and spatial designs, such as the packing designs available in Gosset, are suggested as a suitable design complement. The spatial regression models also perform well with no global design, for example with data obtained from series of designs on the same space of design variables. The approach is illustrated with an example involving the optimization of components in a DNA amplification assay. A Monte Carlo comparison of the spatial models with both thin plate splines and second-order polynomial response surfaces for a scenario motivated by the example is also given. This shows superior performance of the spatial models to the second-order polynomials with respect to both prediction over the complete design space and for cross-validation prediction error in the region of the optimum. An anisotropic spatial regression model performs best for a high noise case and both this model and the thin plate spline for a low noise case. Spatial regression is recommended for construction of response surfaces in all process optimization applications.     

High-Dimensional Multivariate Time Series With Additional Structure

High-dimensional multivariate time series are challenging due to the dependent and high-dimensional nature of the data, but in many applications there is additional structure that can be exploited to reduce computing time along with statistical error. We consider high-dimensional vector autoregressive processes with spatial structure, a simple and common form of additional structure. We propose novel high-dimensional methods that take advantage of such structure without making model assumptions about how distance affects dependence. We provide nonasymptotic bounds on the statistical error of parameter estimators in high-dimensional settings and show that the proposed approach reduces the statistical error. An application to air pollution in the USA demonstrates that the estimation approach reduces both computing time and prediction error and gives rise to results that are meaningful from a scientific point of view, in contrast to high-dimensional methods that ignore spatial structure. In practice, these high-dimensional methods can be used to decompose high-dimensional multivariate time series into lower-dimensional multivariate time series that can be studied by other methods in more depth. Supplementary materials for this article are available online.     

空间变系数模型的统计诊断

空间变系数模型作为一类有效的空间数据分析方法已经得到了广泛的应用.本文主要研究该模型的统计诊断与影响分析方法。首先我们基于数据删除模型定义了Cook统计量,其次我们基于均值漂移模型讨论了异常点的检验问题。     

Modelling of spatial variability of soil undrained shear strength by conditional random fields for slope reliability analysis

Conditional random field model can make best use of limited site investigation data to properly characterize the spatial variation of soil properties. This paper aims to propose a simplified approach for generating conditional random fields of soil undrained shear strength. A numerical method is adopted to validate the effectiveness of the proposed approach. With the proposed approach, the analytical posterior statistics of spatially varying undrained shear strength conditioned on the known values at measurement locations can be obtained. The conditional random field model of undrained shear strength is constructed using the field vane shear test data at a site of the west side highway in New York and the probability of slope failure is estimated by subset simulation. A clay slope under undrained conditions is investigated as an example to illustrate the proposed approach. The effects of borehole location and borehole layout scheme on the slope reliability are addressed. The results indicate that the proposed approach not only can well incorporate the limited site investigation data into modelling of the actual spatial variation of soil parameters by conditional random fields, but also can capture the depth-dependent nature of soil properties. The realizations of conditional random fields generated by the proposed approach can be well constrained to the site investigation data.