Applied Statistics : Practical 11
|
|
- Jasper Holland
- 6 years ago
- Views:
Transcription
1 Applied Statistics : Practical 11 This practical will introduce basic tools for geostatistics in R. You may need first to install and load a few packages. The packages sp and lattice contain useful function and structures for the management of spatially distributed data. The package gstat provides tools for the analysis of geostatistical data. 1. Exploratory data analysis of geostatistical data In this item we consider a classical dataset in geostatistics which is avaliable in the sp package. The data set consists of 155 samples of top soil heavy metal concentrations (ppm), along with a number of soil and landscape variables. The samples were collected in a flood plain of the river Meuse, near the village Stein (The Netherlands). > library(sp) > data(meuse) > head(meuse) ## x y cadmium copper lead zinc elev dist om ffreq soil ## ## ## ## ## ## ## lime landuse dist.m ## 1 1 Ah 50 ## 2 1 Ah 30 ## 3 1 Ah 150 ## 4 0 Ga 270 ## 5 0 Ah 380 ## 6 0 Ga 470 You can see that the dataset reports the geographical coordinates (x and y) as well as the measurements associated to each data point. We need first to instruct R about which column correspond to the coordinates: > coordinates(meuse) <- c('x','y') This changes the nature of the dataset, which is now treated as a SpatialPointsDataFrame, i.e. a data frame which associated spatial location. This is one of the data structure that are available in R for spatial data and it is provided by the package sp. If we now try to plot the data, R gives as the locations of the data point: > plot(meuse) Other types of plots are avalable to explore the other variables in the data set, for example the quantity of zinc in the soil: > spplot(meuse,'zinc',do.log=true) > bubble(meuse,'zinc',do.log=true) Can you interpret the output of this plot? (Check the help if needed) In these plots, each data point is associated to an interval of zinc quantity and the coding is based either on color (spplot) 1
2 or the dimension of the dot (bubble). A logarithmic scale is used to determine the intervals in which to partition the range of quantity of zinc. While the SpatialPointsDataFrame structure is usually preferred for geostatistical data, other types of structure are available in the sp package. For example, the river borders can be described (and plotted) using SpatialPolygons > data(meuse.riv) > meuse.lst <- list(polygons(list(polygon(meuse.riv)), "meuse.riv")) > meuse.sr <- SpatialPolygons(meuse.lst) > plot(meuse.sr, col = "grey") > plot(meuse, add = TRUE) Looking at the geography, can you suggest any interpretation for the variation in the quantity of zinc? The quantity of zinc appears to be larger the closer the data points are to the river. To explore the spatial dependence, we can first plot the semivariogram cloud, i.e. the empirical semivarogram for all the distances observed in the dataset (we are transforming the quantity of zinc on the log scale first): > library(gstat) > cld <- variogram(log(zinc) ~ 1, data=meuse, cloud = TRUE) > plot(cld, main = 'Semivariogram cloud') 2
3 Semivariogram cloud 3 semivariance 2 1 or better the binned semivariogram distance > svgm <- variogram(log(zinc) ~ 1, width=100, data=meuse) > plot(svgm, main = 'Binned Semivariogram',pch=19) 3
4 Binned Semivariogram 0.6 semivariance distance The parameter width controls the bandwidth, try to change the value and see what happens. It is also possible to include covariates in the formula (in place of 1), to account for a drift term in the model. Does the plot of the binned semivariogram suggests the presence of spatial dependence? If yes, does the process appear to be (second order) stationary? The binned semivariogram suggests the presence of spatial dependence beacuse its values increase with the spatial distance. The process appears to be stationary beacuse the binned semivariogram does not diverge to plus infinity for large distances. We can then fit a model semivariogram to the binned semivariogram, via weighted least squares, using the vgm and fit.variogram function. vgm provides the expressions for a galley of model semivariograms, you can type vgm() for a complete list. Let us now fit a spherical semivariogram: > sph.model<-fit.variogram(svgm, vgm(psill=0.6, "Sph", range=800, nugget=0.2)) > plot(svgm,sph.model) > sph.model We needed to specify a starting point for the model parameters in the optimization algorithm: the partial sill (0.6), the range (800) and the nugget (0.2). Reasonable choices for these starting values can be obtained looking at the binned variogram plot. In particular, the algorithm may fail if you select completely unreasonable choices for the (effective) range. Try now to fit an exponential semivariogram model, what are the estimated nugget, sill and effective range? Which one best fits the binned semivariogram? 4
5 > exp.model<-fit.variogram(svgm, vgm(0.6, "Exp", 800, 0.2)) > plot(svgm,exp.model) 0.6 semivariance distance > exp.model τ 2 = , σ 2 = = and R = = The spherical semivariogram is a best fit for the data, because the fitted exponential model has a range larger than the one observed in the data. 2. Simulated examples In this item we consider a few simulated examples to explore the empirical and fitted variogram in known situations. To simulate a spatially dependent random field, we can use the gstat package. > library(gstat) > # define first a grid for simulations > xy <- expand.grid(1:100, 1:100) > names(xy)<-c("x","y") > > # spherical example > sph.field<-gstat(formula=z~ 1,locations=~x+y, dummy=t, beta = 1, + model=vgm(psill=0.025, range=10, model='sph'), nmax=20) > sph.sim<-predict(sph.field,newdata=xy,nsim=1) ## [using unconditional Gaussian simulation] 5
6 > # to plot the simulated field > library(sp) > gridded(sph.sim) = ~x+y > spplot(sph.sim[1]) Try changing the sill and the range to appreciate how they influence the realization of the random field. Estimate the binned semivariogram and fit a spherical semivariogram model. > sph.vario <- variogram( sim1~ 1, width=5, data=sph.sim) > sph.model<-fit.variogram(sph.vario, vgm(0.02, "Sph", 10, 0.005)) > plot(sph.vario,sph.model) semivariance distance Simulate then a new random field with a pure nugget (i.e. spatial white noise), by chosing the Nug model. What happens if you try to fit a spherical semivariogram in this case? > # nugget simulation > noise.field<-gstat(formula=z~ 1,locations=~x+y, dummy=t, beta = 1, model=vgm(psill=0,ran > noise.sim<-predict(noise.field,newdata=xy,nsim=1) ## [using unconditional Gaussian simulation] > # to plot the simulated field > gridded(noise.sim) = ~x+y 6
7 > spplot(obj=noise.sim[1]) > svgm <- variogram(sim1 ~ 1, width=5,data=noise.sim) > sph.model<-fit.variogram(svgm, vgm(0.02, "Sph", 10, 0.005)) ## Warning in fit.variogram(svgm, vgm(0.02, "Sph", 10, 0.005)): singular model in variogram fit > plot(svgm,sph.model) R gives a warning because the model is degenerate. However, the fit may be not too bad, because if the sill and the nugget almost coincide whatever the range the semivariogram is essentially constant. This estimates is also very dependent on the width of the bins in the binned semivariogram. 3. Kriging prediction Let us consider again the Meuse dataset, now with the aim of reconstructing a smooth surface of the (logarithm of) lead concentration. > library(gstat) > library(sp) > data(meuse) > coordinates(meuse) <- c('x','y') > # river meuse > data(meuse.riv) > meuse.lst <- list(polygons(list(polygon(meuse.riv)), "meuse.riv")) > meuse.sr <- SpatialPolygons(meuse.lst) > > # grid for prediction > data(meuse.grid) > coordinates(meuse.grid) <- c('x','y') > meuse.grid <- as(meuse.grid, 'SpatialPixelsDataFrame') Let us start by assuming the process as stationary and thus estimating the semivariogram first: > v <- variogram(log(lead) ~ 1, meuse) > plot(v,pch=19) > #spherical > v.fit <- fit.variogram(v, vgm(1, "Sph", 800, 1)) > plot(v, v.fit, pch = 19) > #exponential > v.fit <- fit.variogram(v, vgm(1, "Exp", 800, 1)) > plot(v, v.fit, pch = 19) > # gaussian > v.fit <- fit.variogram(v, vgm(1, "Gau", 800, 1)) > plot(v, v.fit, pch = 19) Which model provides the best fit? semivariogram. The spherical model is the one that fits best the binned > v.fit <- fit.variogram(v, vgm(1, "Sph", 800, 1)) 7
8 There are two alternative ways to obtain the kriging prediction, by using the krige function or by fitting a model with gstat and then use the predict option for that model. Both are doing the same mathematical operations (actually, krige is just a wrapper for gstat and predict ). Let us consider the simple kriging prediction first (we to pretend to know the true mean, even if we estimate it from the data). > # > beta.hat <- mean(log(meuse$lead)) # assumed known > #simple kriging prediction > lz.sk <- krige(log(lead)~1, meuse, meuse.grid, v.fit, beta = beta.hat) ## [using simple kriging] > plot(lz.sk) > # or > lz.mod<-gstat(formula=log(lead)~1, data=meuse, model=v.fit, beta = beta.hat) > lz.sk2<-predict(lz.mod,meuse.grid) ## [using simple kriging] > plot(lz.sk2) > > # the spplot gives you prediction > # variance as well (although the colorscale is not great) > spplot(lz.sk, col.regions=bpy.colors(n = 100, cutoff.tails = 0.2), + main = 'Simple Kriging') The ordinary kriging prediction can be easily obtained by not providing a mean value: > ## ordinary kriging > lz.ok <- krige(log(lead)~1, meuse, meuse.grid, v.fit) ## [using ordinary kriging] Looking at the concentration of lead and the position of the observation with respect to the river, > data(meuse.riv) > meuse.lst <- list(polygons(list(polygon(meuse.riv)), "meuse.riv")) > meuse.sr <- SpatialPolygons(meuse.lst) > plot(meuse.sr, col = "grey") > plot(meuse, add = TRUE) > spplot(meuse,'lead',do.log=true) we can see that the lead concentration appears to change with the distance from the river (which is one of the parameter in the dataset, dist). Let us consider first the scatterplot of log(lead) and dist. > plot(meuse$dist,log(meuse$lead)) 8
9 log(meuse$lead) meuse$dist This suggests to fit a universal kriging model where the mean is a function of the distance from the river (possibly a square root, looking at the scatterplot). > lead.gstat <- gstat(id = 'lead', formula = log(lead) ~ sqrt(dist), + data = meuse, model=v.fit) > lead.gstat ## data: ## lead : formula = log(lead)`~`sqrt(dist) ; data dim = 155 x 12 ## variograms: ## model psill range ## lead[1] Nug ## lead[2] Sph > lead.uk<-predict(lead.gstat,newdata = meuse.grid) ## [using universal kriging] > plot(lead.uk,main="universal kriging") The gstat function automatically (and iteratively) estimates the drift, the residuals and the residual variogram, then the predict function compute the universal kriging predictor for the new locations. It is also possible to get an evaluation of the drift using the option BLUE=TRUE: 9
10 > lead.trend<-predict(lead.gstat,newdata = meuse.grid,blue=true) ## [generalized least squares trend estimation] > plot(lead.trend,main="drift") You can also get the predicted value of the field or the estimated trend in a specific new observation (but note that you need to provide also the correspondent value for the predictors in the non-stationary mean model, in this case the distance from the river) > new_obs<-data.frame(x=179660,y=331860,dist= ) > coordinates(new_obs)<-c('x','y') > predict(lead.gstat,newdata=new_obs) ## [using universal kriging] ## coordinates lead.pred lead.var ## 1 (179660, ) > predict(lead.gstat,newdata=new_obs,blue=true) ## [generalized least squares trend estimation] ## coordinates lead.pred lead.var ## 1 (179660, )
The meuse data set: a brief tutorial for the gstat R package
The meuse data set: a brief tutorial for the gstat R package Edzer Pebesma March 11, 2017 1 Introduction The meuse data set provided by package sp is a data set comprising of four heavy metals measured
More informationPackage automap. May 30, 2011
Package automap May 30, 2011 Version 1.0-10 Date 2010/30/05 Title Automatic interpolation package Author Paul Hiemstra Maintainer Paul Hiemstra Description
More informationUniversity of California, Los Angeles Department of Statistics
University of California, Los Angeles Department of Statistics Statistics C173/C273 Instructor: Nicolas Christou Ordinary kriging using geor and gstat In this document we will discuss kriging using the
More informationAnalysing Spatial Data in R: Representing Spatial Data
Analysing Spatial Data in R: Representing Spatial Data Roger Bivand Department of Economics Norwegian School of Economics and Business Administration Bergen, Norway 31 August 2007 Object framework To begin
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 Tobler s Law All places are related, but nearby places are related more than distant places Corollary:
More informationSpatial Interpolation - Geostatistics 4/3/2018
Spatial Interpolation - Geostatistics 4/3/201 (Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Distance between pairs of points Lag Mean Tobler s Law All places are related, but nearby places
More informationTowards meaningful spatial statistics
Towards meaningful spatial statistics ue IfGI-Logo 1.6 Logovarianten Edzer Pebesma, Christoph Stasch, Simon Scheider, Werner Kuhn Einsatz in internationalen bzw. achigen Präsentationen. iche: Briefbogen,
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 11/3/2016 GEO327G/386G, UT Austin 1 Tobler s Law All places are related, but nearby places are related
More informationKriging Peter Claussen 9/5/2017
Kriging Peter Claussen 9/5/2017 Libraries automap : Automatic interpolation package library(automap) ## Loading required package: sp ## Warning: package 'sp' was built under R version 3.3.2 library(gstat)
More informationUniversity of California, Los Angeles Department of Statistics
Statistics C173/C273 University of California, Los Angeles Department of Statistics Cross Validation Instructor: Nicolas Christou Cross validation is a technique that allows us to compare values with true
More information1 Layering, and Alternatives to Layering
Lattice With Layers layer() and as.layer() (from latticeextra) 1 J H Maindonald http://www.maths.anu.edu.au/~johnm Centre for Mathematics and Its Applications Australian National University. Contents 1
More informationSoftware Tutorial Session Universal Kriging
Software Tutorial Session Universal Kriging The example session with PG2000 which is described in this and Part 1 is intended as an example run to familiarise the user with the package. This documented
More informationSome methods for the quantification of prediction uncertainties for digital soil mapping: Universal kriging prediction variance.
Some methods for the quantification of prediction uncertainties for digital soil mapping: Universal kriging prediction variance. Soil Security Laboratory 2018 1 Universal kriging prediction variance In
More informationGeo372 Vertiefung GIScience. Spatial Interpolation
Geo372 Vertiefung GIScience Spatial Interpolation Herbstsemester Ross Purves Last week We looked at data quality and integration We saw how uncertainty could be introduced by the concepts we chose, by
More informationGridding and Contouring in Geochemistry for ArcGIS
Gridding and Contouring in Geochemistry for ArcGIS The Geochemsitry for ArcGIS extension includes three gridding options: Minimum Curvature Gridding, Kriging and a new Inverse Distance Weighting algorithm.
More informationPackage rtop. July 2, 2014
Package rtop July 2, 2014 Type Package Title Interpolation of data with variable spatial support Version 0.3-45 Date 2014-02-06 Author Maintainer Depends sp Imports gstat Suggests
More informationAnalysing Spatial Data in R: Vizualising Spatial Data
Analysing Spatial Data in R: Vizualising Spatial Data Roger Bivand Department of Economics Norwegian School of Economics and Business Administration Bergen, Norway 31 August 2007 Vizualising Spatial Data
More informationLab 2 - Simulation Study
Lab 2 - Simulation Study Matthew T. Pratola 1, Jervyn Ang 1 and Andrew MacDougal 2 February 5, 2009 1 Department of Statistics and Actuarial Science, Simon Fraser University, Bunaby, BC, Canada 2 Department
More informationASSIGNMENT 3 Cobalt data T:\sys502\arcview\projects\Cobalt_2 T:\sys502\arcview\ projects\cobalt_2\cobalt_2.mxd T:\sys502\arcview\projects\Cobalt_2
ESE 502 Tony E. Smith ASSIGNMENT 3 (1) In this study you will use some of the Cobalt data from Vancouver Island Geochemistry data set in B&G (pp.150, 202) to carry out essentially the same type of analysis
More informationTutorial Session -- Semi-variograms
Tutorial Session -- Semi-variograms The example session with PG2000 which is described below is intended as an example run to familiarise the user with the package. This documented example illustrates
More informationGEO 580 Lab 4 Geostatistical Analysis
GEO 580 Lab 4 Geostatistical Analysis In this lab, you will move beyond basic spatial analysis (querying, buffering, clipping, joining, etc.) to learn the application of some powerful statistical techniques
More informationPackage geofd. R topics documented: October 5, Version 1.0
Version 1.0 Package geofd October 5, 2015 Author Ramon Giraldo , Pedro Delicado , Jorge Mateu Maintainer Pedro Delicado
More informationPractical 11: Interpolating Point Data in R
Practical 11: Interpolating Point Data in R This practical provides an introduction to some techniques which are useful for interpolating point data across space in R. Interpolation describes a means of
More informationspsann Optimization of Sample Configurations Using Spatial Simulated Annealing
spsann Optimization of Sample Configurations Using Spatial Simulated Annealing Alessandro Samuel-Rosa June 23, 2017 Contents 1 Introduction 1 2 Getting Started 2 2.1 Package structure...........................
More informationThe sgeostat Package
The sgeostat Package April 11, 2007 Title An Object-oriented Framework for Geostatistical Modeling in S+ Version 1.0-21 Author S original by James J. Majure Iowa State University R
More informationUniversity of California, Los Angeles Department of Statistics. Spatial prediction
University of California, Los Angeles Department of Statistics Statistics C173/C273 Instructor: Nicolas Christou Spatial prediction Spatial prediction: a. One of our goals in geostatistics is to predict
More informationASSIGNMENT 4. (a) First you will construct a visual representation of the data a follows:
ESE 502 Tony E. Smith ASSIGNMENT 4 This study is designed as an extension of Assignment 3. Hence there is no need to repeat everything you said (or perhaps didn't say!) in that assignment. One strategy
More informationThe gstat Package. March 12, 2008
The gstat Package March 12, 2008 Version 0.9-44 Date 2008/03/11 Title geostatistical modelling, prediction and simulation Author Edzer J. Pebesma and others Maintainer Edzer J. Pebesma
More informationGeostatistics 2D GMS 7.0 TUTORIALS. 1 Introduction. 1.1 Contents
GMS 7.0 TUTORIALS 1 Introduction Two-dimensional geostatistics (interpolation) can be performed in GMS using the 2D Scatter Point module. The module is used to interpolate from sets of 2D scatter points
More informationThe gstat Package. September 21, 2007
The gstat Package September 21, 2007 Version 0.9-40 Date 2007/09/03 Title geostatistical modelling, prediction and simulation Author Edzer J. Pebesma and others Maintainer Edzer J.
More informationPackage gstat. July 2, 2014
Package gstat July 2, 2014 Version 1.0-19 Date 2014/04/01 Title spatial and spatio-temporal geostatistical modelling, prediction and simulation variogram modelling; simple, ordinary and universal point
More informationWatershed Sciences 4930 & 6920 GEOGRAPHIC INFORMATION SYSTEMS
HOUSEKEEPING Watershed Sciences 4930 & 6920 GEOGRAPHIC INFORMATION SYSTEMS Quizzes Lab 8? WEEK EIGHT Lecture INTERPOLATION & SPATIAL ESTIMATION Joe Wheaton READING FOR TODAY WHAT CAN WE COLLECT AT POINTS?
More informationGIS Exercise - Spring, 2011
GIS Exercise - Spring, 2011 Maria Antonia Brovelli Laura Carcano, Marco Minghini ArcGIS exercise 3 - Global trend removal Introduction: Besides exploring statistical characteristics and distributional
More informationNon Stationary Variograms Based on Continuously Varying Weights
Non Stationary Variograms Based on Continuously Varying Weights David F. Machuca-Mory and Clayton V. Deutsch Centre for Computational Geostatistics Department of Civil & Environmental Engineering University
More informationUniversity of California, Los Angeles Department of Statistics. Assign NA values outside the area defined by the observed data points
University of California, Los Angeles Department of Statistics Statistics C173/C273 Instructor: Nicolas Christou Assign NA values outside the area defined by the observed data points In this document we
More informationFlexible Lag Definition for Experimental Variogram Calculation
Flexible Lag Definition for Experimental Variogram Calculation Yupeng Li and Miguel Cuba The inference of the experimental variogram in geostatistics commonly relies on the method of moments approach.
More informationNon-linearity and spatial correlation in landslide susceptibility mapping
Non-linearity and spatial correlation in landslide susceptibility mapping C. Ballabio, J. Blahut, S. Sterlacchini University of Milano-Bicocca GIT 2009 15/09/2009 1 Summary Landslide susceptibility modeling
More informationUsing Kriging Methods to Estimate Damage Distribution
Lehigh University Lehigh Preserve Theses and Dissertations 2011 Using Kriging Methods to Estimate Damage Distribution Bing Wang Lehigh University Follow this and additional works at: http://preserve.lehigh.edu/etd
More informationSurface Smoothing Using Kriging
1 AutoCAD Civil 3D has the ability to add data points to a surface based on mathematical criteria. This gives you the ability to strengthen the surface definition in areas where data may be sparse or where
More informationFurther Simulation Results on Resampling Confidence Intervals for Empirical Variograms
University of Wollongong Research Online Centre for Statistical & Survey Methodology Working Paper Series Faculty of Engineering and Information Sciences 2010 Further Simulation Results on Resampling Confidence
More informationWhat can we represent as a Surface?
Geography 38/42:376 GIS II Topic 7: Surface Representation and Analysis (Chang: Chapters 13 & 15) DeMers: Chapter 10 What can we represent as a Surface? Surfaces can be used to represent: Continuously
More informationExercise 3-1: Soil property mapping
Exercise 3-1: Soil property mapping Mapping objectives: Create a soil analysis point layer from tabular data in QGIS Create a continuous surface soil property map using Kriging Interpolation in SAGA GIS
More informationPackage plotgooglemaps
Package plotgooglemaps February 15, 2013 Type Package Title Plot SP data as HTML map mashup over Google Maps Version 1.3 Date 2012-07-06 Author Milan Kilibarda Maintainer Milan Kilibarda
More informationDijkstra's Algorithm
Shortest Path Algorithm Dijkstra's Algorithm To find the shortest path from the origin node to the destination node No matrix calculation Floyd s Algorithm To find all the shortest paths from the nodes
More informationMethods to define confidence intervals for kriged values: Application on Precision Viticulture data.
Methods to define confidence intervals for kriged values: Application on Precision Viticulture data. J-N. Paoli 1, B. Tisseyre 1, O. Strauss 2, J-M. Roger 3 1 Agricultural Engineering University of Montpellier,
More informationAn Automated Method for Hot-to-Cold Geometry Mapping
Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2015-05-01 An Automated Method for Hot-to-Cold Geometry Mapping Brandon Levi Doolin Brigham Young University - Provo Follow this
More informationCombining punctual and ordinal contour data for accurate floodplain topography mapping
Proceedings of Spatial Accuracy 2016 [ISBN: 978-2-95-45-5] Combining punctual and ordinal contour data for accurate floodplain topography mapping Carole Delenne *1,3, Jean-Stéphane Bailly 2, Mathieu Dartevelle
More informationThe mapmisc package. Patrick Brown. August 28, 2017
The mapmisc package Patrick Brown August 28, 2017 This document will be incomplete if rgdal is unavailable or there is on internet connection when this document is compiled. The full document is at http://diseasemapping.
More informationgstat user s manual ln_zinc Nug(0) Sph( ) distance
gstat user s manual 0.7 0.6 543 500 564 589 452 477 533 574 457 415 semivariance 0.5 0.4 0.3 419 457 547 0.2 299 57 0.1 ln_zinc 0.0536187 Nug(0) + 0.581735 Sph(892.729) 0 0 200 400 600 800 1000 1200 1400
More informationPackage gstat. March 12, 2017
Version 1.1-5 Package gstat March 12, 2017 Title Spatial and Spatio-Temporal Geostatistical Modelling, Prediction and Simulation Variogram modelling; simple, ordinary and universal point or block (co)kriging;
More informationPackage UncerIn2. November 24, 2015
Version 2.0 Date 2015-11-10 Type Package Package UncerIn2 November 24, 2015 Title Implements Models of Uncertainty into the Interpolation Functions Author Tomas Burian Maintainer Tomas
More informationPackage Grid2Polygons
Package Grid2Polygons February 15, 2013 Version 0.1-2 Date 2013-01-28 Title Convert Spatial Grids to Polygons Author Jason C. Fisher Maintainer Jason C. Fisher Depends R (>= 2.15.0),
More informationMINI-PAPER A Gentle Introduction to the Analysis of Sequential Data
MINI-PAPER by Rong Pan, Ph.D., Assistant Professor of Industrial Engineering, Arizona State University We, applied statisticians and manufacturing engineers, often need to deal with sequential data, which
More informationA Comparison of Spatial Prediction Techniques Using Both Hard and Soft Data
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Dissertations and Theses in Statistics Statistics, Department of 5-2011 A Comparison of Spatial Prediction Techniques Using
More informationSurface Creation & Analysis with 3D Analyst
Esri International User Conference July 23 27 San Diego Convention Center Surface Creation & Analysis with 3D Analyst Khalid Duri Surface Basics Defining the surface Representation of any continuous measurement
More informationPackage intamap. R topics documented: May 2, Version Date Title Procedures for Automated Interpolation
Version 1.4-9 Date 2018-05-02 Title Procedures for Automated Interpolation Package intamap May 2, 2018 Author Edzer Pebesma , Jon Olav Skoien with
More informationUniversity of California, Los Angeles Department of Statistics
Statistics C173/C273 University of California, Los Angeles Department of Statistics Instructor: Nicolas Christou Computing the variogram using the geor package in R Spatial statistics computations can
More informationVariogram Inversion and Uncertainty Using Dynamic Data. Simultaneouos Inversion with Variogram Updating
Variogram Inversion and Uncertainty Using Dynamic Data Z. A. Reza (zreza@ualberta.ca) and C. V. Deutsch (cdeutsch@civil.ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta
More informationAirborne discrete return LiDAR data was collected on September 3-4, 2007 by
SUPPLEMENTAL MATERIAL 2 LiDAR Specifications Airborne discrete return LiDAR data was collected on September 3-4, 2007 by Watershed Sciences, Inc. (Corvallis, Oregon USA). LiDAR was collected approximately
More informationLocal spatial-predictor selection
University of Wollongong Research Online Centre for Statistical & Survey Methodology Working Paper Series Faculty of Engineering and Information Sciences 2013 Local spatial-predictor selection Jonathan
More informationOasis montaj Advanced Mapping
Oasis montaj Advanced Mapping As more information becomes available in digital format, Earth science specialists are recognizing that it is essential to work with a variety of data from different sources.
More information8: Tree-based regression
8: Tree-based regression John H Maindonald June 18, 2018 Ideas and issues illustrated by the graphs in this vignette The fitting of a tree proceeds by making a succession of splits on the x-variable or
More informationR Package CircSpatial for the Imaging - Kriging - Simulation. of Circular-Spatial Data
R Package CircSpatial for the Imaging - Kriging - Simulation 2 4-20 -10 0 y 2 4-20 -10 0 10 of Circular-Spatial Data Bill Morphet PhD Advisor Juergen Symanzik April, 2008 1 Circular Random Variable (CRV)
More informationModeling Multiple Rock Types with Distance Functions: Methodology and Software
Modeling Multiple Rock Types with Distance Functions: Methodology and Software Daniel A. Silva and Clayton V. Deutsch The sub division of the deposit into estimation domains that are internally consistent
More informationApplied Geostatistics with JMP Mauromoustakos A. and K. Thompson, U of Arkansas, Fayetteville, AR
Paper 213-27 Applied Geostatistics with JMP Mauromoustakos A. and K. Thompson, U of Arkansas, Fayetteville, AR ABSTRACT Presentation of basic visual geostatistical techniques using JMP V4 software will
More informationYou will begin by exploring the locations of the long term care facilities in Massachusetts using descriptive statistics.
Getting Started 1. Create a folder on the desktop and call it your last name. 2. Copy and paste the data you will need to your folder from the folder specified by the instructor. Exercise 1: Explore the
More informationFathom Dynamic Data TM Version 2 Specifications
Data Sources Fathom Dynamic Data TM Version 2 Specifications Use data from one of the many sample documents that come with Fathom. Enter your own data by typing into a case table. Paste data from other
More informationGeostatistical Reservoir Characterization of McMurray Formation by 2-D Modeling
Geostatistical Reservoir Characterization of McMurray Formation by 2-D Modeling Weishan Ren, Oy Leuangthong and Clayton V. Deutsch Department of Civil & Environmental Engineering, University of Alberta
More informationLesson 5 overview. Concepts. Interpolators. Assessing accuracy Exercise 5
Interpolation Tools Lesson 5 overview Concepts Sampling methods Creating continuous surfaces Interpolation Density surfaces in GIS Interpolators IDW, Spline,Trend, Kriging,Natural neighbors TopoToRaster
More information2D Geostatistical Modeling and Volume Estimation of an Important Part of Western Onland Oil Field, India.
and Volume Estimation of an Important Part of Western Onland Oil Field, India Summary Satyajit Mondal*, Liendon Ziete, and B.S.Bisht ( GEOPIC, ONGC) M.A.Z.Mallik (E&D, Directorate, ONGC) Email: mondal_satyajit@ongc.co.in
More informationBMEGUI Tutorial 1 Spatial kriging
BMEGUI Tutorial 1 Spatial kriging 1. Objective The primary objective of this exercise is to get used to the basic operations of BMEGUI using a purely spatial dataset. The analysis will consist in an exploratory
More informationFMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu
FMA901F: Machine Learning Lecture 3: Linear Models for Regression Cristian Sminchisescu Machine Learning: Frequentist vs. Bayesian In the frequentist setting, we seek a fixed parameter (vector), with value(s)
More informationApplied Statistics : Practical 9
Applied Statistics : Practical 9 This practical explores nonparametric regression and shows how to fit a simple additive model. The first item introduces the necessary R commands for nonparametric regression
More informationCPSC 695. Methods for interpolation and analysis of continuing surfaces in GIS Dr. M. Gavrilova
CPSC 695 Methods for interpolation and analysis of continuing surfaces in GIS Dr. M. Gavrilova Overview Data sampling for continuous surfaces Interpolation methods Global interpolation Local interpolation
More informationLab 12: Sampling and Interpolation
Lab 12: Sampling and Interpolation What You ll Learn: -Systematic and random sampling -Majority filtering -Stratified sampling -A few basic interpolation methods Data for the exercise are found in the
More informationPart I, Chapters 4 & 5. Data Tables and Data Analysis Statistics and Figures
Part I, Chapters 4 & 5 Data Tables and Data Analysis Statistics and Figures Descriptive Statistics 1 Are data points clumped? (order variable / exp. variable) Concentrated around one value? Concentrated
More informationA Geostatistical and Flow Simulation Study on a Real Training Image
A Geostatistical and Flow Simulation Study on a Real Training Image Weishan Ren (wren@ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta Abstract A 12 cm by 18 cm slab
More informationUniversity of California, Los Angeles Department of Statistics
University of California, Los Angeles Department of Statistics Statistics 12 Instructor: Nicolas Christou Data analysis with R - Some simple commands When you are in R, the command line begins with > To
More informationAdditional File 5: Step-by-step guide to using Rwui to create webapps that display R script results on Google dynamic maps
Additional File 5: Step-by-step guide to using Rwui to create webapps that display R script results on Google dynamic maps Richard Newton, Andrew Deonarine and Lorenz Wernisch July 7, 2012 1 Simple example
More informationMap overlay and spatial aggregation in sp
Map overlay and spatial aggregation in sp Edzer Pebesma Nov 2010 Contents 1 Introduction 1 2 Geometry overlays 2 3 Extraction, aggregation of attributes 5 3.1 Extracting attribute values.....................
More informationSpatial and multi-scale data assimilation in EO-LDAS. Technical Note for EO-LDAS project/nceo. P. Lewis, UCL NERC NCEO
Spatial and multi-scale data assimilation in EO-LDAS Technical Note for EO-LDAS project/nceo P. Lewis, UCL NERC NCEO Abstract Email: p.lewis@ucl.ac.uk 2 May 2012 In this technical note, spatial data assimilation
More informationMarkov Bayes Simulation for Structural Uncertainty Estimation
P - 200 Markov Bayes Simulation for Structural Uncertainty Estimation Samik Sil*, Sanjay Srinivasan and Mrinal K Sen. University of Texas at Austin, samiksil@gmail.com Summary Reservoir models are built
More informationClustering Lecture 5: Mixture Model
Clustering Lecture 5: Mixture Model Jing Gao SUNY Buffalo 1 Outline Basics Motivation, definition, evaluation Methods Partitional Hierarchical Density-based Mixture model Spectral methods Advanced topics
More informationUniversity of California, Los Angeles Department of Statistics
University of California, Los Angeles Department of Statistics Statistics 13 Instructor: Nicolas Christou Data analysis with R - Some simple commands When you are in R, the command line begins with > To
More informationFurther Maths Notes. Common Mistakes. Read the bold words in the exam! Always check data entry. Write equations in terms of variables
Further Maths Notes Common Mistakes Read the bold words in the exam! Always check data entry Remember to interpret data with the multipliers specified (e.g. in thousands) Write equations in terms of variables
More informationIn SigmaPlot, a restricted but useful version of the global fit problem can be solved using the Global Fit Wizard.
Global Fitting Problem description Global fitting (related to Multi-response onlinear Regression) is a process for fitting one or more fit models to one or more data sets simultaneously. It is a generalization
More informationDirect Sequential Co-simulation with Joint Probability Distributions
Math Geosci (2010) 42: 269 292 DOI 10.1007/s11004-010-9265-x Direct Sequential Co-simulation with Joint Probability Distributions Ana Horta Amílcar Soares Received: 13 May 2009 / Accepted: 3 January 2010
More informationResampling Methods for Dependent Data
S.N. Lahiri Resampling Methods for Dependent Data With 25 Illustrations Springer Contents 1 Scope of Resampling Methods for Dependent Data 1 1.1 The Bootstrap Principle 1 1.2 Examples 7 1.3 Concluding
More informationVariography Setting up the Parameters When a data file has been loaded, and the Variography tab selected, the following screen will be displayed.
Variography - Introduction The variogram (or semi-variogram) is a graph relating the variance of the difference in value of a variable at pairs of sample points to the separation distance between those
More informationsurface but these local maxima may not be optimal to the objective function. In this paper, we propose a combination of heuristic methods: first, addi
MetaHeuristics for a Non-Linear Spatial Sampling Problem Eric M. Delmelle Department of Geography and Earth Sciences University of North Carolina at Charlotte eric.delmelle@uncc.edu 1 Introduction In spatial
More informationBeyond The Vector Data Model - Part Two
Beyond The Vector Data Model - Part Two Introduction Spatial Analyst Extension (Spatial Analysis) What is your question? Selecting a method of analysis Map Algebra Who is the audience? What is Spatial
More informationLearner Expectations UNIT 1: GRAPICAL AND NUMERIC REPRESENTATIONS OF DATA. Sept. Fathom Lab: Distributions and Best Methods of Display
CURRICULUM MAP TEMPLATE Priority Standards = Approximately 70% Supporting Standards = Approximately 20% Additional Standards = Approximately 10% HONORS PROBABILITY AND STATISTICS Essential Questions &
More informationIntroduction Basics Simple Statistics Graphics. Using R for Data Analysis and Graphics. 4. Graphics
Using R for Data Analysis and Graphics 4. Graphics Overview 4.1 Overview Several R graphics functions have been presented so far: > plot(d.sport[,"kugel"], d.sport[,"speer"], + xlab="ball push", ylab="javelin",
More informationModeling Uncertainty in the Earth Sciences Jef Caers Stanford University
Modeling spatial continuity Modeling Uncertainty in the Earth Sciences Jef Caers Stanford University Motivation uncertain uncertain certain or uncertain uncertain Spatial Input parameters Spatial Stochastic
More informationInstitute for Statics und Dynamics of Structures Fuzzy Time Series
Institute for Statics und Dynamics of Structures Fuzzy Time Series Bernd Möller 1 Description of fuzzy time series 2 3 4 5 Conclusions Folie 2 von slide422 1 Description of fuzzy time series 2 3 4 5 Conclusions
More informationPackage SSN. July 2, 2014
Package SSN July 2, 2014 Type Package Title Spatial Modeling on Stream Networks Version 1.1.3 Date 2014-5-27 Depends R (>= 3.0.2), RSQLite, sp Imports MASS, igraph (>= 0.6), maptools, lattice, methods
More informationDISTRIBUTION STATEMENT A Approved for public release: distribution unlimited.
AVIA Test Selection through Spatial Variance Bounding Method for Autonomy Under Test By Miles Thompson Senior Research Engineer Aerospace, Transportation, and Advanced Systems Lab DISTRIBUTION STATEMENT
More informationCreating Surfaces. Steve Kopp Steve Lynch
Steve Kopp Steve Lynch Overview Learn the types of surfaces and the data structures used to store them Emphasis on surface interpolation Learn the interpolation workflow Understand how interpolators work
More informationUsing Excel for Graphical Analysis of Data
Using Excel for Graphical Analysis of Data Introduction In several upcoming labs, a primary goal will be to determine the mathematical relationship between two variable physical parameters. Graphs are
More informationHierarchical Modelling for Large Spatial Datasets
Hierarchical Modelling for Large Spatial Datasets Sudipto Banerjee 1 and Andrew O. Finley 2 1 Department of Forestry & Department of Geography, Michigan State University, Lansing Michigan, U.S.A. 2 Biostatistics,
More information