Sequential, Bayesian Geostatistics: A Principled Method for Large Data Sets

The principled statistical application of Gaussian random field models used in geostatistics has historically been limited to data sets of a small size. This limitation is imposed by the requirement to store and invert the covariance matrix of all the samples to obtain a predictive distribution at unsampled locations, or to use likelihood-based covariance estimation. Various ad hoc approaches to solve this problem have been adopted, such as selecting a neighborhood region and/or a small number of observations to use in the kriging process, but these have no sound theoretical basis and it is unclear what information is being lost. In this article, we present a Bayesian method for estimating the posterior mean and covariance structures of a Gaussian random field using a sequential estimation algorithm. By imposing sparsity in a well-defined framework, the algorithm retains a subset of " basis vectors " that best represent the " true " posterior Gaussian random field model in the relative entropy sense. This allows a principled treatment of Gaussian random field models on very large data sets. The method is particularly appropriate when the Gaussian random field model is regarded as a latent variable model, which may be nonlinearly related to the observations. We show the application of the sequential, sparse Bayesian estimation in Gaussian random field models and discuss its merits and drawbacks.     

Geoarchaeological studies on Roman time harbour sediments in Cologne — comparison of different OSL dating techniques

Due to the construction of a new North-South subway in Cologne, Roman time harbour sediments were exposed and were sampled for luminescence dating. A very good independent age control was given by the precise knowledge of the chronology of Roman activity and by radiocarbon ages of charcoal samples. Hence, different methodological approaches within luminescence dating were applied for Holocene heterogeneously bleached fluvial samples and were compared to the known ages. For one sample, optically stimulated luminescence (OSL) dating was applied to coarse-grained quartz using a single aliquot regenerative-dose (SAR) protocol. After D_e-measurements, different statistical approaches were tested (i.e. arithmetic mean, median, minimum age model, finite mixture model, leading edge method and the Fuchs and Lang approach). It is demonstrated that the Fuchs and Lang approach along with the leading edge method yielded the best matching OSL ages with respect to the known ages.     

Theory of Comparative Advantage: Do Transportation Costs Matter?

In this paper we present a formal analysis that incorporates returns to transportation into a Ricardian framework to predict trade patterns. The important point gained from this analysis is that increasing returns to transportation, coupled with appropriate distances between trading partners, can be shown to reverse Ricardian predictions even when there are no international differences in tastes, technology, or factor endowments. Additional gains from trade may emerge from reductions in aggregate delivery costs owing to scale economies.     

Vedisch kapana

Ohne Zusammenfassung     

ARTIFICIAL NEURAL NETWORKS: A NEW APPROACH TO MODELING INTERREGIONAL TELECOMMUNICATION FLOWS*

ABSTRACT. During the last thirty years there has been much research effort in regional science devoted to modeling interactions over geographic space. Theoretical approaches for studying these phenomena have been modified considerably. This paper suggests a new modeling approach, based upon a general nested sigmoid neural network model. Its feasibility is illustrated in the context of modeling interregional telecommunication traffic in Austria, and its performance is evaluated in comparison with the classical regression approach of the gravity type. The application of this neural network approach may be viewed as a three-stage process. The first stage refers to the identification of an appropriate network from the family of two-layered feedforward networks with 3 input nodes, one layer of (sigmoidal) intermediate nodes and one (sigmoidal) output node (logistic activation function). There is no general procedure to address this problem. We solved this issue experimentally. The input-output dimensions have been chosen in order to make the comparison with the gravity model as close as possible. The second stage involves the estimation of the network parameters of the selected neural network model. This is performed via the adaptive setting of the network parameters (training, estimation) by means of the application of a least mean squared error goal and the error back propagating technique, a recursive learning procedure using a gradient search to minimize the error goal. Particular emphasis is laid on the sensitivity of the network performance to the choice of the initial network parameters, as well as on the problem of overfitting. The final stage of applying the neural network approach refers to the testing of the interregional teletraffic flows predicted. Prediction quality is analyzed by means of two performance measures, average relative variance and the coefficient of determination, as well as by the use of residual analysis. The analysis shows that the neural network model approach outperforms the classical regression approach to modeling telecommunication traffic in Austria.