Introduction. Semivariogram Cloud
|
|
- Griffin Gaines
- 5 years ago
- Views:
Transcription
1 Introduction Data: set of n attribute measurements {z(s i ), i = 1,, n}, available at n sample locations {s i, i = 1,, n} Objectives: Slide 1 quantify spatial auto-correlation, or attribute dissimilarity typically expressed as: 1 2 [z(s i) z(s j )] 2 as a function of separation distance between sample pairs s i and s j introduce the sample semivariogram, its characteristics, and provide some examples NOTE: Spatial auto-correlation is a second-order characteristic of spatial variation, and hence the sample semivariogram should be computed from data whose spatial variation is not explained by first-order effects justify the need of going beyond the sample semivariogram to a semivariogram model introduce parametric functions of distance that can be used as formal theoretical semivariogram models discuss issues of fitting semivariogram models to sample semivariogram values Semivariogram Cloud Definition: A scatter-plot of attribute squared semidifferences between all possible pairs of samples measured at different locations, versus their separation distance Computational procedure: Slide 2 1 construct Euclidean distance matrix D = [d ij, i = 1,, n, j = 1,, n] between all n 2 pairs of data locations, where d ij is defined as: d ij = h ij = s i s j 2 construct squared semidifference matrix E = [e ij, i = 1,, n, j = 1,, n] between all n 2 pairs of attribute values, where e ij is defined as: e ij = 1 2 [z(s i) z(s j )] 2 3 plot each distance value d ij against the corresponding squared semidifference e ij ; in other words, plot e = vec(e) versus d = vec(d) The plot of all pairs {d ij, e ij } is termed a semivariogram cloud
2 Semivariogram Cloud Example Bay Area rain gauge precipitation Semivariogram cloud semidifferences squared Slide distance A measure of dissimilarity between attribute values measured at different locations, ie, a spatial measure of attribute dissimilarity Expected graph pattern: As the distance d ij between sample pairs increases, the corresponding squared semidifference e ij should also increase Difficult to interpret, so we consider groups of sample pairs separated by similar distances ie, average squared semidifferences within distance classes (x-axis bins in the right graph above) Semivariogram Cloud Versus Plot 9 Semivariogram cloud 12 Sample semivariogram 8 7 squared semidifference semivariance values Slide distance Going from the first to the second: distance class or midpoint define a set of L distance classes; the l-th class has limits: (d l t l, d l + t l ], where d l is the class midpoint and t l is half the class width (or distance tolerance) for a given distance class (d l t l, d l + t l ], the semivariogram value ˆγ(d l ) is the average of n(d l ) << n 2 squared attribute semidifferences computed from sample pairs whose inter-distances d ij satisfy: d l t l < d ij d l + t l in other words, the semivariogram plot can be regarded as a summary of the semivariogram cloud, according to some distance-based grouping of samples
3 Computing Sample Semivariograms 1 compute distance matrix D = [d ij, i = 1,, n, j = 1,, n] and squared semidifference matrix E = [e ij, i = 1,, n, j = 1,, n] between n 2 data pairs D = d 12 d 13 d 14 d 15 d 12 d 23 d 24 d 25 d 13 d 23 d 34 d 35 d 14 d 24 d 34 d 45 d 15 d 25 d 35 d 45 E = e 12 e 13 e 14 e 15 e 12 e 23 e 24 e 25 e 13 e 23 e 34 e 35 e 14 e 24 e 34 e 45 e 15 e 25 e 35 e 45 Slide 5 2 for a given distance class (d l t l, d l + t l ], find entries of E that correspond to entries of D falling in that distance class, eg: d 12 d 13 d 14 d 15 d 12 d 23 d 24 d 25 d 13 d 23 d 34 d 35 d 14 d 24 d 34 d 45 d 15 d 25 d 35 d 45 e 12 e 13 e 14 e 15 e 12 e 23 e 24 e 25 e 13 e 23 e 34 d 35 e 14 e 24 e 34 e 45 e 15 e 25 e 35 e 45 3 sample semivariogram ˆγ(d l ) for that class is the average of the n(d l ) squared semidifferences, e-values, whose corresponding distances, d-values, fall in class (d l t l, d l + t l ]; ie, the mean of all e-values in boxes in the matrix on the right above Examples of Semivariogram Computation Locations separated by distance class (5 15] Histogram of squared semi differences (5 15] 7 y coordinates frequency # of data 68 mean 262 std dev 475 minimum 1 q(25) 1 median 39 q(75) 22 maximum Slide x coordinates Locations separated by distance class (15 25] squared semi-difference eij Histogram of squared semi differences (15 25] 7 y coordinates frequency # of data 16 mean 386 std dev 733 minimum 1 q(25) 3 median 129 q(75) 49 maximum x coordinates squared semi-difference eij ˆγ((5 15]) = 262, ˆγ((15 25]) = 386 = averages of values displayed in histograms Map views linking sample pairs that contribute to such histograms are extremely informative
4 Sample Semivariogram Plots Consider a set of L distance classes with midpoints {d l, l = 1,, L} and tolerances {t l, l = 1,, L} The plot of semivariance values {ˆγ(d l ), l = 1,, L} versus the average sample inter-distance for each class is called a sample semivariogram ˆγ(d l ) = 1 n(d l ) e c = 1 n(d l ) 2n(d l ) c=1 n(d l ) d ij (d l t l,d l +t l ] [z(s i ) z(s j )] 2 Slide 7 Bay Area rain gauge precipitation Sample semivariogram semivariance γ ( h) lag distance h numbers above bullets denote # of sample pairs contributing to ˆγ(d l ) at each lag distance could also graph variances of e-values within the distance classes; ˆγ() =, always Semivariogram Characteristics 2 Variogram, γ(h) Increasing Variability Slide Lag Distance (h) sill: limit semivariogram value (plateau) is approximately equal to sample variance (for representative sample) range: distance at which semivariogram reaches (or starts oscillating around) sill = distance of influence of any datum on another nugget effect: discontinuity at origin (ˆγ() > ); sum of measurement error and micro-structures (variability at scales smaller than sampling interval) watch out for sparse data, outliers and positional or attribute errors transformation of Euclidean distance into statistical distance bearing imprint of specific phenomenon
5 Sample Semivariogram Shape & Interpretation (1) Quadratic shape near origin: Image Semivariogram: quadratic shape at origin North γ 8 Slide East Distance Interpretation: highly continuous (extremely smooth) spatial attribute variability spatial attribute is differentiable typical variables: elevation, temperature, Sample Semivariogram Shape & Interpretation (2) Linear shape near origin: Image Semivariogram: linear shape at origin North γ 8 Slide East Distance Interpretation: continuous variability (not extremely smooth) of spatial attribute attribute is not differentiable typical variables: ore grades,
6 Sample Semivariogram Shape & Interpretation (3) Discontinuous near origin: Image Semivariogram: discontinuous at origin Slide 11 North -1 γ East Distance Interpretation: highly irregular (quasi-random) spatial variability at small scales typical variables: precipitation, Sample Semivariogram Shape & Interpretation (4) Oscillating (around sill): 2 DN (image intensity) values Semivariogram of DN values North 2 15 γ 8 4 Slide East 2 Distance Interpretation: periodic variability of spatial attribute yields sinusoidal semivariogram semivariogram shape possibly due to limited sampling need to provide physical evidence for periodicity frequently encountered in time series
7 The Need for Semivariogram Models Problems: (i) sill, range, and relative nugget, cannot be determined directly from the sample semivariogram plot, (ii) a continuum of semivariogram values γ(d) for any distance vector d is required in interpolation, but sample semivariogram values {ˆγ(d l ), l = 1,, L} are typically calculated only for few (L) distances {d l, l = 1,, L} Slide 13 Semivariogram model definition: parametric function γ(d; θ) fitted to sample semivariogram values {ˆγ(d l ), l = 1,, L}; θ denotes parameter vector with, eg, range, and sill (for a given semivariogram function) Sample semivariogram of precipitation Sample and model semivariogram of precipitation γ 8 γ sample variogram model variogram Distance (degrees) Distance (degrees) semivariogram modeling is more than a curve fitting exercise; Warning: cannot use any curve as semivariogram model!!! Valid Semivariogram Models: Pure Nugget Effect 12 Nugget effect variogram model semivariance γ(d) Slide lag distance d, if d = γ(d; θ) = σ, if d > θ = [σ], where σ denotes attribute variance indicates complete absence of spatial correlation could occur due to measurement error and microstructure, ie, features occurring at scales smaller than sampling interval
8 Valid Semivariogram Models: Spherical 15 Spherical variogram model (sill=, range=3) semivariance γ(d) 5 Slide lag distance d 3 σ d 2 r 1 d 3, if d < r 2 r γ(d; θ) = σ, if d r θ = [σ r], where r is the model range linear behavior at origin clearly defined range parameter r Valid Semivariogram Models: Exponential 15 Exponential variogram model (sill=, range=3) semivariance γ(d) 5 Slide lag distance d γ(d; θ) = σ 1 exp 3d r θ = [σ r] linear behavior at origin; rises faster than spherical; reaches sill asymptotically effective range parameter r; distance at which 95% of sill reached
9 Valid Semivariogram Models: Gaussian 15 Gaussian variogram model (sill=, range=3) semivariance γ(d) 5 Slide lag distance d γ(d; θ) = σ 1 exp 3d2 r 2 θ = [σ r] quadratic behavior at origin; implies smooth spatial variability of attribute values; reaches sill asymptotically effective range parameter r; distance at which 95% of sill reached Valid Semivariogram Models: Nugget + Exponential 1 Standardized nugget+exponential variogram model 8 γ Slide 18 Distance, if d = γ(d; θ) = a + ([σ a][1 exp( 3d )]), if d r θ = [σ a r] discontinuous at origin; reaches sill asymptotically practical range parameter r; distance at which 95% of sill reached a/σ = relative nugget contribution = proportion (to total sill) of purely random spatial variability more complex models can be built by adding or multiplying valid models
10 Fitting Semivariogram Models to Sample Data Or fitting valid semivariogram functions (curves) to sample semivariogram values Manual fitting: select number of semivariograms, their type (functional form), sill, and range model behavior at origin (nugget effect, shape of semivariogram at distances smaller than first lag) using prior knowledge about phenomenon Slide 19 Automatic fitting: least squares fit (ordinary, generalized, weighted): choose semivariogram model parameters (typically iteratively) so as to minimize discrepancy between model and sample semivariogram values over all lags; other methods also available treat with caution, especially with sparse data and outliers Cross-validation: given a proposed parameter set, ie, a semivariogram model, perform cross-validation using geostatistical interpolation, and record resulting error statistics repeat with different model parameters, and select as optimal model the one whose parameters yield best cross-validation error statistics Summary Spatial auto-correlation can be quantified by looking at attribute dissimilarity as a function of separation distance The semivariogram cloud is too cloudy for detecting meaningful patterns Slide 2 The semivariogram plot is constructed by averaging squared semidifferences within distance bins to smooth out the variability in the semivariogram cloud NOTE: Watch out for trends (first-order effects) in the data; a sample semivariogram quantifies second-order effects and might be contaminated by variations due to trends/drifts A quantitative way to encapsulate a sample semivariogram is through a parametric semivariogram model Fitting procedures exist for estimating the parameters of semivariogram models, ie, for fitting model semivariograms to sample semivariograms The final semivariogram model can be used for simulation (pattern generation) and geostatistical interpolation NOTE: A semivariogram model is a spatial process model, whose parameters are inferred from the sample data through the sample semivariogram
11 Introduction Data: set of n attribute measurements {z(s i ), i = 1,, n}, available at n sample locations {s i, i = 1,, n} Objectives: (i) predict or interpolate unknown attribute value z(s p ) at location s p from the n sample data, and (ii) assess reliability of predicted value Slide 1 Geostatistical spatial interpolation: predicted attribute value = weighted linear combination of sample data values + attribute mean, if known (non-linear methods also exist) a semivariogram model is used to determine the weights, which account for: spatial auto-correlation between sample data and unknown value spatial auto-correlation between sample data themselves (data redundancy) in addition, and contrary to most interpolation algorithms, geostatistics offers a measure of reliability (prediction error variance) regarding the attribute prediction SK prediction: ẑ(s p ) = m + Simple Kriging (SK) n w p (s i )[z(s i ) m] = wp T r i=1 w p = [w p (s i ), i = 1,, n] T : (n 1) vector of SK-weights assigned to n sample data for prediction at location s p ; superscript T denotes transposition r = [z(s i ) m, i = 1,, n] T : (n 1) vector of residual data from known mean m Slide 2 z(s 1 ) m ẑ(s p ) = m + w p (s 1 ) w p (s i ) w p (s n ) z(s i ) m wp T z(s n ) m r use semivariogram model to determine weights at each prediction location; typically, it is the covariogram model that is used due to computational reasons Lecture Notes Introduction to Geostatistical Spatial Interpolation total # of slides = 22
12 Semivariogram / Covariogram / Correlogram Model 15 Semivariogram model 15 Covariogram model 15 Correlogram model variance γ( ) = σ() σ() variance ρ() 1 unit correlation semivariance γ(d) 5 covariance σ(d) 5 correlation ρ(d) 5 Slide 3 range lag distance d range lag distance d range lag distance d Conversion between models, with σ() = γ( ) being the sill of the semivariogram model: Semivariogram covariogram: σ(d) = σ() γ(d) Covariogram correlogram: Semivariogram correlogram: Covariogram semivariogram: ρ(d) = σ(d) σ() ρ(d) = 1 γ(d) σ() γ(d) = σ() σ(d) Requisites for Geostatistical Interpolation I Slide 4 Data-to-data and data-to-unknown distances: d 1j d 1n D = d i1 d in d n1 d nj and d p = d 1p d ip d np Comments: as any other interpolation method, one accounts for the proximity of the n sample locations to the prediction location s p Note: Vector d p changes from one prediction location s p to another, hence the subscript p unlike other interpolation methods, one also accounts for the proximity between sample locations themselves (sample configuration or data layout) Note: Matrix D of sample-to-sample distances is the same for all prediction locations Lecture Notes Introduction to Geostatistical Spatial Interpolation total # of slides = 22
13 Requisites for Geostatistical Interpolation II From distance matrices to model covariance matrices: Take any distance value d ij and d ip, ie, any entry in D and d p, and transform it, via the covariogram model, to a covariance value σ(d ij ) and σ(d ip ) Slide 5 Data-to-data and data-to-unknown model covariances: σ() σ(d 1j ) σ(d 1n ) Σ = σ(d i1 ) σ() σ(d in ) and σ p = σ(d n1 ) σ(d nj ) σ() σ(d 1p ) σ(d ip ) σ(d np ) data-to-data covariance matrix Σ: (n n) matrix with model covariance values σ(d ij ) between any two sample locations separated by distance d ij data-to-unknown covariance vector σ p : (n 1) vector with model covariance values σ(d ip ) between the n sample locations and the prediction location s p Note: Vector σ p changes from one prediction location s p to another, hence the subscript p Requisites for Geostatistical Interpolation III Slide 6 Data-to-data and data-to-unknown model covariances: Comments: Σ = σ() σ(d 1j ) σ(d 1n ) σ(d i1 ) σ() σ(d in ) σ(d n1 ) σ(d nj ) σ() and σ p = σ(d 1p ) σ(d ip ) σ(d np ) data-to-data covariance matrix Σ: encapsulates the redundancy between the sample data; for positive spatial auto-correlation, the more clustered is the sample layout, the more redundant are the sample data (less information content); a clustered sample layout typically translates into larger entries in Σ data-to-unknown covariance vector σ p : encapsulates the statistical proximity (correlation) between the sample data and the unknown attribute value z(s p ) at the prediction location s p ; that correlation is a function of distance between sample and prediction locations, not of the actual (unknown) value z(s p ); The larger the entries of vector σ p, the stronger the predictive power of sample data Lecture Notes Introduction to Geostatistical Spatial Interpolation total # of slides = 22
14 Simple Kriging (SK) System & Weights σ() σ(d 1n ) σ(d n1 ) σ() w p (s 1 ) w p (s n ) = σ(d 1p ) σ(d np ) Σw p = σ p Comments: Slide 7 the SK system is a (disguised) version of the normal equations for the case of regression with no intercept term: X T Xb = X T y, where X is the design matrix and y is the vector of data on the dependent variable; in regression, the data-to-data covariance is estimated as X T X/n, and the data-to-unknown covariance as X T y/n the weights vector w p is obtained by solving the SK system, as w p = Σ 1 σ p, anew at each prediction location s p since the entries of σ p change entries of w p do not depend on data values or on sill, σ(), of covariogram model: σ() ρ() ρ(d 1n ) ρ(d n1 ) 1 w p (s 1 ) w p (s n ) = σ() ρ(d 1p ) ρ(d np ) w p (s 1 ) w p (s n ) Interpreting the Simple Kriging Weights 1 ρ(d 1n ) = 1 σ() ρ(d n1 ) 1 1 ρ(d 1p ) σ() w p = Σ 1 σ p ρ(d np ) Slide 8 if sample interdistances d ij are larger than correlogram range, then ρ(d ij ) =, and Σ = σ()i, the (n n) identity matrix; this entails that w p (s i ) = ρ(d ip ), ie, weights are equal to correlogram values but in general, Σ = σ()i, ie, sample interdistances are within correlation range, in which case Σ 1 modulates σ p : influence of samples in clusters is downplayed the closer the sample data to the prediction location, and the more spread out the data over the study region, the better the SK prediction is expected to be for sample data far away (beyond correlation range) from the prediction location s p, ρ(d ip ) = and w p (s i ) = : all weighs are equal to for prediction at a sample location s p s i, data-to-unknown covariance vector σ p = σ i is same as i-th column of Σ; this yields w p (s i ) = 1 if s i = s p, otherwise: only sample co-located with prediction location receives non-zero (= 1) weight Lecture Notes Introduction to Geostatistical Spatial Interpolation total # of slides = 22
15 Simple Kriging Prediction and Error Variance Once the SK weights are computed as w p = Σ 1 σ p, they are substituted in the following equations to compute the SK prediction ẑ(s p ) and associated error variance ˆσ(s p ) Slide 9 SK prediction does not depend on sill σ() of covariogram model: z(s 1 ) m ẑ(s p ) = m+wp T r = m+[w p (s 1 ) w p (s n )] = m+ n w p (s i )[z(s i ) m] i=1 z(s n ) m SK prediction error variance does depend on covariogram model sill σ(): σ(d 1p ) ˆσ(s p ) = σ() wp T σ p = σ() [w p (s 1 ) w p (s n )] = σ() n w p (s i )σ(d ip ) n=1 σ(d np ) which can also be written as: ˆσ(s p ) = σ() 1 n i=1 w p(s i )ρ(d ip ) Interpreting the SK Prediction and Error Variance ẑ(s p ) = m + Comments: n w p (s i )[z(s i ) m] i=1 n ˆσ(s p ) = σ() w p (s i )σ(d ip ) i=1 Slide for sample data far away (beyond correlation range) from the prediction location s p, w p (s i ) =, i: all weighs are equal to In this case, the SK prediction equals the known mean m and the SK error variance equals the known covariogram sill: ẑ(s p ) = m and ˆσ(s p ) = σ(); away from the sample data, SK yields back the (assumed known) attribute overall mean and variance for prediction at a sample location s p s i, w p (s i ) = 1 if s i = s p, otherwise: the SK prediction identifies the known sample datum and the SK error variance is zero: ẑ(s i ) = z(s i ) and ˆσ(s i ) = ; SK is an exact interpolation algorithm for all other prediction locations, the SK predictions depend on the sample data configuration and their values, while the SK error variances depend only on the sample data configuration; both SK predictions and error variances depend on the covariogram model σ(d) adopted Lecture Notes Introduction to Geostatistical Spatial Interpolation total # of slides = 22
16 Determining the SK Weights: Step Local data configuration (1) (2) (5) (3) (4) 646 (6) 783 (7) Slide (n x n) matrix of data-to-data inter-distances: D = i, j-th element of D: d ij = s i s j (n x 1) vector of prediction-to-data-location distances: T d p = i-th element of d p : d ip = s i s p Determining the SK Weights: Step Local data configuration 1 Correlogram model (1) (2) (5) (3) 8 6 ρ (d) ρ(d) = exp( 3d ) (4) 646 (6) 783 (7) 2 Slide d 361 exp( 3 361/) exp( 3 447/) exp( 3 671/) exp( 3 86/) = exp( 3 894/) exp( 3 949/) exp( /) 2 d p σ p =sill exp( 3d p /range) These would be the weights if one ignored auto-correlation between sample data Lecture Notes Introduction to Geostatistical Spatial Interpolation total # of slides = 22
17 Determining the SK Weights: Step Local data configuration 1 Correlogram model (1) (2) (5) (3) 8 6 ρ (d) ρ(d) = exp( 3d ) (4) 646 (6) 783 (7) 2 Slide SK system: Σ d w p (s 1 ) 34 w p (s 2 ) 26 w p (s 3 ) 13 w p (s 4 ) = 9 w p (s 5 ) 7 w p (s 6 ) 6 w p (s 7 ) 2 w p σ p i, j-th element of matrix Σ: σ ij = 1 exp( 3 d ij /) Determining the SK Weights: Step 4 Slide 14 w p (s 1 ) w p (s 2 ) w p (s 3 ) w p (s 4 ) w p (s 5 ) w p (s 6 ) w p (s 7 ) wp = Σ 1 σp 145 SK weights (1) (2) -1 2 (5) (3) prediction = variance = (4) 28 7 (6) (7) original weights vector (w p = σ p ) modified by Σ 1 to account for sample redundancy; eg, w p (s 1 ) = 27 instead of ρ(d 1p ) = 34 Lecture Notes Introduction to Geostatistical Spatial Interpolation total # of slides = 22
11/8/2018. Spatial Interpolation & Geostatistics. Kriging Step 1
(Z i Z j ) 2 / 2 (Z i Zj) 2 / 2 Semivariance y 11/8/2018 Spatial Interpolation & Geostatistics Kriging Step 1 Describe spatial variation with Semivariogram Lag Distance between pairs of points Lag Mean
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 y Kriging Step 1 Describe spatial variation with Semivariogram (Z i Z j ) 2 / 2 Point cloud Map 3
More informationExploring the World of Ordinary Kriging. Dennis J. J. Walvoort. Wageningen University & Research Center Wageningen, The Netherlands
Exploring the World of Ordinary Kriging Wageningen University & Research Center Wageningen, The Netherlands July 2004 (version 0.2) What is? What is it about? Potential Users a computer program for exploring
More informationInvestigation of Monthly Pan Evaporation in Turkey with Geostatistical Technique
Investigation of Monthly Pan Evaporation in Turkey with Geostatistical Technique Hatice Çitakoğlu 1, Murat Çobaner 1, Tefaruk Haktanir 1, 1 Department of Civil Engineering, Erciyes University, Kayseri,
More informationTypes of Spatial Data
Spatial Data Types of Spatial Data Point pattern Point referenced geostatistical Block referenced Raster / lattice / grid Vector / polygon Point Pattern Data Interested in the location of points, not their
More informationBasics of Point-Referenced Data Models
Basics of Point-Referenced Data Models Basic tool is a spatial process, {Y (s), s D}, where D R r Chapter 2: Basics of Point-Referenced Data Models p. 1/45 Basics of Point-Referenced Data Models Basic
More informationPRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH
PRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH SURESH TRIPATHI Geostatistical Society of India Assumptions and Geostatistical Variogram
More informationLecture 5 Geostatistics
Lecture 5 Geostatistics Lecture Outline Spatial Estimation Spatial Interpolation Spatial Prediction Sampling Spatial Interpolation Methods Spatial Prediction Methods Interpolating Raster Surfaces with
More informationChapter 4 - Fundamentals of spatial processes Lecture notes
TK4150 - Intro 1 Chapter 4 - Fundamentals of spatial processes Lecture notes Odd Kolbjørnsen and Geir Storvik January 30, 2017 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites
More informationIntroduction to Geostatistics
Introduction to Geostatistics Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, 2017 1 Department of Biostatistics, Bloomberg School of Public Health, Johns Hopkins University, Baltimore,
More informationKriging Luc Anselin, All Rights Reserved
Kriging Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline Principles Kriging Models Spatial Interpolation
More informationPoint-Referenced Data Models
Point-Referenced Data Models Jamie Monogan University of Georgia Spring 2013 Jamie Monogan (UGA) Point-Referenced Data Models Spring 2013 1 / 19 Objectives By the end of these meetings, participants should
More informationA robust statistically based approach to estimating the probability of contamination occurring between sampling locations
A robust statistically based approach to estimating the probability of contamination occurring between sampling locations Peter Beck Principal Environmental Scientist Image placeholder Image placeholder
More informationInfluence of parameter estimation uncertainty in Kriging: Part 2 Test and case study applications
Hydrology and Earth System Influence Sciences, of 5(), parameter 5 3 estimation (1) uncertainty EGS in Kriging: Part Test and case study applications Influence of parameter estimation uncertainty in Kriging:
More informationSpatial statistics, addition to Part I. Parameter estimation and kriging for Gaussian random fields
Spatial statistics, addition to Part I. Parameter estimation and kriging for Gaussian random fields 1 Introduction Jo Eidsvik Department of Mathematical Sciences, NTNU, Norway. (joeid@math.ntnu.no) February
More informationIndex. Geostatistics for Environmental Scientists, 2nd Edition R. Webster and M. A. Oliver 2007 John Wiley & Sons, Ltd. ISBN:
Index Akaike information criterion (AIC) 105, 290 analysis of variance 35, 44, 127 132 angular transformation 22 anisotropy 59, 99 affine or geometric 59, 100 101 anisotropy ratio 101 exploring and displaying
More informationUniversity of California, Los Angeles Department of Statistics. Effect of variogram parameters on kriging weights
University of California, Los Angeles Department of Statistics Statistics C173/C273 Instructor: Nicolas Christou Effect of variogram parameters on kriging weights We will explore in this document how the
More informationGridding of precipitation and air temperature observations in Belgium. Michel Journée Royal Meteorological Institute of Belgium (RMI)
Gridding of precipitation and air temperature observations in Belgium Michel Journée Royal Meteorological Institute of Belgium (RMI) Gridding of meteorological data A variety of hydrologic, ecological,
More informationData Break 8: Kriging the Meuse RiverBIOS 737 Spring 2004 p.1/27
Data Break 8: Kriging the Meuse River BIOS 737 Spring 2004 Data Break 8: Kriging the Meuse RiverBIOS 737 Spring 2004 p.1/27 Meuse River: Reminder library(gstat) Data included in gstat library. data(meuse)
More informationThe ProbForecastGOP Package
The ProbForecastGOP Package April 24, 2006 Title Probabilistic Weather Field Forecast using the GOP method Version 1.3 Author Yulia Gel, Adrian E. Raftery, Tilmann Gneiting, Veronica J. Berrocal Description
More informationGEOSTATISTICS. Dr. Spyros Fountas
GEOSTATISTICS Dr. Spyros Fountas Northing (m) 140550 140450 140350 Trent field Disturbed area Andover 140250 Panholes 436950 437050 437150 437250 437350 Easting (m) Trent Field Westover Farm (Blackmore,
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry
More informationLecture 9: Introduction to Kriging
Lecture 9: Introduction to Kriging Math 586 Beginning remarks Kriging is a commonly used method of interpolation (prediction) for spatial data. The data are a set of observations of some variable(s) of
More informationI don t have much to say here: data are often sampled this way but we more typically model them in continuous space, or on a graph
Spatial analysis Huge topic! Key references Diggle (point patterns); Cressie (everything); Diggle and Ribeiro (geostatistics); Dormann et al (GLMMs for species presence/abundance); Haining; (Pinheiro and
More informationGeostatistics: Kriging
Geostatistics: Kriging 8.10.2015 Konetekniikka 1, Otakaari 4, 150 10-12 Rangsima Sunila, D.Sc. Background What is Geostatitics Concepts Variogram: experimental, theoretical Anisotropy, Isotropy Lag, Sill,
More informationPracticum : Spatial Regression
: Alexandra M. Schmidt Instituto de Matemática UFRJ - www.dme.ufrj.br/ alex 2014 Búzios, RJ, www.dme.ufrj.br Exploratory (Spatial) Data Analysis 1. Non-spatial summaries Numerical summaries: Mean, median,
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models 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 informationFluvial Variography: Characterizing Spatial Dependence on Stream Networks. Dale Zimmerman University of Iowa (joint work with Jay Ver Hoef, NOAA)
Fluvial Variography: Characterizing Spatial Dependence on Stream Networks Dale Zimmerman University of Iowa (joint work with Jay Ver Hoef, NOAA) March 5, 2015 Stream network data Flow Legend o 4.40-5.80
More informationA kernel indicator variogram and its application to groundwater pollution
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session IPS101) p.1514 A kernel indicator variogram and its application to groundwater pollution data Menezes, Raquel University
More informationØ Set of mutually exclusive categories. Ø Classify or categorize subject. Ø No meaningful order to categorization.
Statistical Tools in Evaluation HPS 41 Dr. Joe G. Schmalfeldt Types of Scores Continuous Scores scores with a potentially infinite number of values. Discrete Scores scores limited to a specific number
More informationCopyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
Chapter 15. SPATIAL INTERPOLATION 15.1 Elements of Spatial Interpolation 15.1.1 Control Points 15.1.2 Type of Spatial Interpolation 15.2 Global Methods 15.2.1 Trend Surface Models Box 15.1 A Worked Example
More informationOn dealing with spatially correlated residuals in remote sensing and GIS
On dealing with spatially correlated residuals in remote sensing and GIS Nicholas A. S. Hamm 1, Peter M. Atkinson and Edward J. Milton 3 School of Geography University of Southampton Southampton SO17 3AT
More informationSpatial Analysis II. Spatial data analysis Spatial analysis and inference
Spatial Analysis II Spatial data analysis Spatial analysis and inference Roadmap Spatial Analysis I Outline: What is spatial analysis? Spatial Joins Step 1: Analysis of attributes Step 2: Preparing for
More informationENGRG Introduction to GIS
ENGRG 59910 Introduction to GIS Michael Piasecki October 13, 2017 Lecture 06: Spatial Analysis Outline Today Concepts What is spatial interpolation Why is necessary Sample of interpolation (size and pattern)
More informationSpatial Data Mining. Regression and Classification Techniques
Spatial Data Mining Regression and Classification Techniques 1 Spatial Regression and Classisfication Discrete class labels (left) vs. continues quantities (right) measured at locations (2D for geographic
More informationSpatial Statistics with Image Analysis. Lecture L02. Computer exercise 0 Daily Temperature. Lecture 2. Johan Lindström.
C Stochastic fields Covariance Spatial Statistics with Image Analysis Lecture 2 Johan Lindström November 4, 26 Lecture L2 Johan Lindström - johanl@maths.lth.se FMSN2/MASM2 L /2 C Stochastic fields Covariance
More informationChapter 1. Summer School GEOSTAT 2014, Spatio-Temporal Geostatistics,
Chapter 1 Summer School GEOSTAT 2014, Geostatistics, 2014-06-19 sum- http://ifgi.de/graeler Institute for Geoinformatics University of Muenster 1.1 Spatial Data From a purely statistical perspective, spatial
More informationSpace-time data. Simple space-time analyses. PM10 in space. PM10 in time
Space-time data Observations taken over space and over time Z(s, t): indexed by space, s, and time, t Here, consider geostatistical/time data Z(s, t) exists for all locations and all times May consider
More informationModels for spatial data (cont d) Types of spatial data. Types of spatial data (cont d) Hierarchical models for spatial data
Hierarchical models for spatial data Based on the book by Banerjee, Carlin and Gelfand Hierarchical Modeling and Analysis for Spatial Data, 2004. We focus on Chapters 1, 2 and 5. Geo-referenced data arise
More information11. Kriging. ACE 492 SA - Spatial Analysis Fall 2003
11. Kriging ACE 492 SA - Spatial Analysis Fall 2003 c 2003 by Luc Anselin, All Rights Reserved 1 Objectives The goal of this lab is to further familiarize yourself with ESRI s Geostatistical Analyst, extending
More informationCOMPARISON OF DIGITAL ELEVATION MODELLING METHODS FOR URBAN ENVIRONMENT
COMPARISON OF DIGITAL ELEVATION MODELLING METHODS FOR URBAN ENVIRONMENT Cahyono Susetyo Department of Urban and Regional Planning, Institut Teknologi Sepuluh Nopember, Indonesia Gedung PWK, Kampus ITS,
More informationPAPER 206 APPLIED STATISTICS
MATHEMATICAL TRIPOS Part III Thursday, 1 June, 2017 9:00 am to 12:00 pm PAPER 206 APPLIED STATISTICS Attempt no more than FOUR questions. There are SIX questions in total. The questions carry equal weight.
More informationAn Introduction to Spatial Statistics. Chunfeng Huang Department of Statistics, Indiana University
An Introduction to Spatial Statistics Chunfeng Huang Department of Statistics, Indiana University Microwave Sounding Unit (MSU) Anomalies (Monthly): 1979-2006. Iron Ore (Cressie, 1986) Raw percent data
More informationStatistícal Methods for Spatial Data Analysis
Texts in Statistícal Science Statistícal Methods for Spatial Data Analysis V- Oliver Schabenberger Carol A. Gotway PCT CHAPMAN & K Contents Preface xv 1 Introduction 1 1.1 The Need for Spatial Analysis
More informationGIST 4302/5302: Spatial Analysis and Modeling
GIST 4302/5302: Spatial Analysis and Modeling Basics of Statistics Guofeng Cao www.myweb.ttu.edu/gucao Department of Geosciences Texas Tech University guofeng.cao@ttu.edu Spring 2015 Outline of This Week
More informationHandbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp
Handbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp Marcela Alfaro Córdoba August 25, 2016 NCSU Department of Statistics Continuous Parameter
More informationLab #3 Background Material Quantifying Point and Gradient Patterns
Lab #3 Background Material Quantifying Point and Gradient Patterns Dispersion metrics Dispersion indices that measure the degree of non-randomness Plot-based metrics Distance-based metrics First-order
More informationTricks to Creating a Resource Block Model. St John s, Newfoundland and Labrador November 4, 2015
Tricks to Creating a Resource Block Model St John s, Newfoundland and Labrador November 4, 2015 Agenda 2 Domain Selection Top Cut (Grade Capping) Compositing Specific Gravity Variograms Block Size Search
More informationE(x i ) = µ i. 2 d. + sin 1 d θ 2. for d < θ 2 0 for d θ 2
1 Gaussian Processes Definition 1.1 A Gaussian process { i } over sites i is defined by its mean function and its covariance function E( i ) = µ i c ij = Cov( i, j ) plus joint normality of the finite
More informationSpatial Cross-correlation Models for Vector Intensity Measures (PGA, Ia, PGV and Sa s) Considering Regional Site Conditions
Spatial Cross-correlation Models for Vector Intensity Measures (PGA, Ia, PGV and Sa s) Considering Regional Site Conditions Gang Wang and Wenqi Du Department of Civil and Environmental Engineering Hong
More informationSoil Moisture Modeling using Geostatistical Techniques at the O Neal Ecological Reserve, Idaho
Final Report: Forecasting Rangeland Condition with GIS in Southeastern Idaho Soil Moisture Modeling using Geostatistical Techniques at the O Neal Ecological Reserve, Idaho Jacob T. Tibbitts, Idaho State
More informationIntroduction. Spatial Processes & Spatial Patterns
Introduction Spatial data: set of geo-referenced attribute measurements: each measurement is associated with a location (point) or an entity (area/region/object) in geographical (or other) space; the domain
More informationAn Introduction to Pattern Statistics
An Introduction to Pattern Statistics Nearest Neighbors The CSR hypothesis Clark/Evans and modification Cuzick and Edwards and controls All events k function Weighted k function Comparative k functions
More informationInterpolation {x,y} Data with Suavity. Peter K. Ott Forest Analysis & Inventory Branch BC Ministry of FLNRO Victoria, BC
Interpolation {x,y} Data with Suavity Peter K. Ott Forest Analysis & Inventory Branch BC Ministry of FLNRO Victoria, BC Peter.Ott@gov.bc.ca 1 The Goal Given a set of points: x i, y i, i = 1,2,, n find
More informationGeostatistics in Hydrology: Kriging interpolation
Chapter Geostatistics in Hydrology: Kriging interpolation Hydrologic properties, such as rainfall, aquifer characteristics (porosity, hydraulic conductivity, transmissivity, storage coefficient, etc.),
More informationWhat s for today. All about Variogram Nugget effect. Mikyoung Jun (Texas A&M) stat647 lecture 4 September 6, / 17
What s for today All about Variogram Nugget effect Mikyoung Jun (Texas A&M) stat647 lecture 4 September 6, 2012 1 / 17 What is the variogram? Let us consider a stationary (or isotropic) random field Z
More informationAdaptive Sampling of Clouds with a Fleet of UAVs: Improving Gaussian Process Regression by Including Prior Knowledge
Master s Thesis Presentation Adaptive Sampling of Clouds with a Fleet of UAVs: Improving Gaussian Process Regression by Including Prior Knowledge Diego Selle (RIS @ LAAS-CNRS, RT-TUM) Master s Thesis Presentation
More informationPackage ProbForecastGOP
Type Package Package ProbForecastGOP February 19, 2015 Title Probabilistic weather forecast using the GOP method Version 1.3.2 Date 2010-05-31 Author Veronica J. Berrocal , Yulia
More information7 Geostatistics. Figure 7.1 Focus of geostatistics
7 Geostatistics 7.1 Introduction Geostatistics is the part of statistics that is concerned with geo-referenced data, i.e. data that are linked to spatial coordinates. To describe the spatial variation
More informationUmeå University Sara Sjöstedt-de Luna Time series analysis and spatial statistics
Umeå University 01-05-5 Sara Sjöstedt-de Luna Time series analysis and spatial statistics Laboration in ArcGIS Geostatistical Analyst These exercises are aiming at helping you understand ArcGIS Geostatistical
More informationØ Set of mutually exclusive categories. Ø Classify or categorize subject. Ø No meaningful order to categorization.
Statistical Tools in Evaluation HPS 41 Fall 213 Dr. Joe G. Schmalfeldt Types of Scores Continuous Scores scores with a potentially infinite number of values. Discrete Scores scores limited to a specific
More informationLognormal Measurement Error in Air Pollution Health Effect Studies
Lognormal Measurement Error in Air Pollution Health Effect Studies Richard L. Smith Department of Statistics and Operations Research University of North Carolina, Chapel Hill rls@email.unc.edu Presentation
More informationAn Introduction to Spatial Autocorrelation and Kriging
An Introduction to Spatial Autocorrelation and Kriging Matt Robinson and Sebastian Dietrich RenR 690 Spring 2016 Tobler and Spatial Relationships Tobler s 1 st Law of Geography: Everything is related to
More informationMultivariate Geostatistics
Hans Wackernagel Multivariate Geostatistics An Introduction with Applications Third, completely revised edition with 117 Figures and 7 Tables Springer Contents 1 Introduction A From Statistics to Geostatistics
More informationUNIVERSITY OF MASSACHUSETTS. Department of Mathematics and Statistics. Basic Exam - Applied Statistics. Tuesday, January 17, 2017
UNIVERSITY OF MASSACHUSETTS Department of Mathematics and Statistics Basic Exam - Applied Statistics Tuesday, January 17, 2017 Work all problems 60 points are needed to pass at the Masters Level and 75
More informationStatistics for analyzing and modeling precipitation isotope ratios in IsoMAP
Statistics for analyzing and modeling precipitation isotope ratios in IsoMAP The IsoMAP uses the multiple linear regression and geostatistical methods to analyze isotope data Suppose the response variable
More informationA MultiGaussian Approach to Assess Block Grade Uncertainty
A MultiGaussian Approach to Assess Block Grade Uncertainty Julián M. Ortiz 1, Oy Leuangthong 2, and Clayton V. Deutsch 2 1 Department of Mining Engineering, University of Chile 2 Department of Civil &
More informationRegression revisited (again)
by I. Clark* Synopsis One of the seminal pioneering papers in reserve evaluation was published by Danie Krige in 1951. In that paper he introduced the concept of regression techniques in providing better
More informationCorrecting Variogram Reproduction of P-Field Simulation
Correcting Variogram Reproduction of P-Field Simulation Julián M. Ortiz (jmo1@ualberta.ca) Department of Civil & Environmental Engineering University of Alberta Abstract Probability field simulation is
More informationBuilding Blocks for Direct Sequential Simulation on Unstructured Grids
Building Blocks for Direct Sequential Simulation on Unstructured Grids Abstract M. J. Pyrcz (mpyrcz@ualberta.ca) and C. V. Deutsch (cdeutsch@ualberta.ca) University of Alberta, Edmonton, Alberta, CANADA
More informationReport on Kriging in Interpolation
Tabor Reedy ENVS421 3/12/15 Report on Kriging in Interpolation In this project I explored use of the geostatistical analyst extension and toolbar in the process of creating an interpolated surface through
More informationImproving Spatial Data Interoperability
Improving Spatial Data Interoperability A Framework for Geostatistical Support-To To-Support Interpolation Michael F. Goodchild, Phaedon C. Kyriakidis, Philipp Schneider, Matt Rice, Qingfeng Guan, Jordan
More informationCBMS Lecture 1. Alan E. Gelfand Duke University
CBMS Lecture 1 Alan E. Gelfand Duke University Introduction to spatial data and models Researchers in diverse areas such as climatology, ecology, environmental exposure, public health, and real estate
More informationCREATION OF DEM BY KRIGING METHOD AND EVALUATION OF THE RESULTS
CREATION OF DEM BY KRIGING METHOD AND EVALUATION OF THE RESULTS JANA SVOBODOVÁ, PAVEL TUČEK* Jana Svobodová, Pavel Tuček: Creation of DEM by kriging method and evaluation of the results. Geomorphologia
More informationDescriptive Univariate Statistics and Bivariate Correlation
ESC 100 Exploring Engineering Descriptive Univariate Statistics and Bivariate Correlation Instructor: Sudhir Khetan, Ph.D. Wednesday/Friday, October 17/19, 2012 The Central Dogma of Statistics used to
More informationSimilarity and Dissimilarity
1//015 Similarity and Dissimilarity COMP 465 Data Mining Similarity of Data Data Preprocessing Slides Adapted From : Jiawei Han, Micheline Kamber & Jian Pei Data Mining: Concepts and Techniques, 3 rd ed.
More informationComparison of rainfall distribution method
Team 6 Comparison of rainfall distribution method In this section different methods of rainfall distribution are compared. METEO-France is the French meteorological agency, a public administrative institution
More informationNon-Ergodic Probabilistic Seismic Hazard Analyses
Non-Ergodic Probabilistic Seismic Hazard Analyses M.A. Walling Lettis Consultants International, INC N.A. Abrahamson University of California, Berkeley SUMMARY A method is developed that relaxes the ergodic
More informationOverview of Statistical Analysis of Spatial Data
Overview of Statistical Analysis of Spatial Data Geog 2C Introduction to Spatial Data Analysis Phaedon C. Kyriakidis www.geog.ucsb.edu/ phaedon Department of Geography University of California Santa Barbara
More informationIntensity Analysis of Spatial Point Patterns Geog 210C Introduction to Spatial Data Analysis
Intensity Analysis of Spatial Point Patterns Geog 210C Introduction to Spatial Data Analysis Chris Funk Lecture 5 Topic Overview 1) Introduction/Unvariate Statistics 2) Bootstrapping/Monte Carlo Simulation/Kernel
More informationAdvances in Locally Varying Anisotropy With MDS
Paper 102, CCG Annual Report 11, 2009 ( 2009) Advances in Locally Varying Anisotropy With MDS J.B. Boisvert and C. V. Deutsch Often, geology displays non-linear features such as veins, channels or folds/faults
More informationA Parametric Spatial Bootstrap
A Parametric Spatial Bootstrap Liansheng Tang a ; William R. Schucany b ; Wayne A. Woodward b ; and Richard F. Gunst b July 17, 2006 a Department of Biostatistics, University of Washington, P.O. Box 357232,
More informationA Short Note on the Proportional Effect and Direct Sequential Simulation
A Short Note on the Proportional Effect and Direct Sequential Simulation Abstract B. Oz (boz@ualberta.ca) and C. V. Deutsch (cdeutsch@ualberta.ca) University of Alberta, Edmonton, Alberta, CANADA Direct
More informationRegression Revisited (again) Isobel Clark Geostokos Limited, Scotland. Abstract
Regression Revisited (again) Isobel Clark Geostokos Limited, Scotland Abstract One of the seminal pioneering papers in reserve evaluation was published by Danie Krige in 1951. In this paper he introduced
More informationENVIRONMENTAL MONITORING Vol. II - Geostatistical Analysis of Monitoring Data - Mark Dowdall, John O Dea GEOSTATISTICAL ANALYSIS OF MONITORING DATA
GEOSTATISTICAL ANALYSIS OF MONITORING DATA Mark Dowdall Norwegian Radiation Protection Authority, Environmental Protection Unit, Polar Environmental Centre, Tromso, Norway John O Dea Institute of Technology,
More informationAreal Unit Data Regular or Irregular Grids or Lattices Large Point-referenced Datasets
Areal Unit Data Regular or Irregular Grids or Lattices Large Point-referenced Datasets Is there spatial pattern? Chapter 3: Basics of Areal Data Models p. 1/18 Areal Unit Data Regular or Irregular Grids
More informationRegression with correlation for the Sales Data
Regression with correlation for the Sales Data Scatter with Loess Curve Time Series Plot Sales 30 35 40 45 Sales 30 35 40 45 0 10 20 30 40 50 Week 0 10 20 30 40 50 Week Sales Data What is our goal with
More informationGeog 210C Spring 2011 Lab 6. Geostatistics in ArcMap
Geog 210C Spring 2011 Lab 6. Geostatistics in ArcMap Overview In this lab you will think critically about the functionality of spatial interpolation, improve your kriging skills, and learn how to use several
More informationGlobal Spatial Autocorrelation Clustering
Global Spatial Autocorrelation Clustering Luc Anselin http://spatial.uchicago.edu join count statistics Moran s I Moran scatter plot non-parametric spatial autocorrelation Join Count Statistics Recap -
More informationSpatial Analysis 2. Spatial Autocorrelation
Spatial Analysis 2 Spatial Autocorrelation Spatial Autocorrelation a relationship between nearby spatial units of the same variable If, for every pair of subareas i and j in the study region, the drawings
More informationCorrelation model for spatially distributed ground-motion intensities
EARTHQUAKE ENGINEERING AND STRUCTURAL DYNAMICS Earthquake Engng Struct. Dyn. 2009; 38:1687 1708 Published online 28 April 2009 in Wiley InterScience (www.interscience.wiley.com)..922 Correlation model
More informationSpatio-temporal precipitation modeling based on time-varying regressions
Spatio-temporal precipitation modeling based on time-varying regressions Oleg Makhnin Department of Mathematics New Mexico Tech Socorro, NM 87801 January 19, 2007 1 Abstract: A time-varying regression
More informationIs there still room for new developments in geostatistics?
Is there still room for new developments in geostatistics? Jean-Paul Chilès MINES ParisTech, Fontainebleau, France, IAMG 34th IGC, Brisbane, 8 August 2012 Matheron: books and monographs 1962-1963: Treatise
More informationWhat s for today. Introduction to Space-time models. c Mikyoung Jun (Texas A&M) Stat647 Lecture 14 October 16, / 19
What s for today Introduction to Space-time models c Mikyoung Jun (Texas A&M) Stat647 Lecture 14 October 16, 2012 1 / 19 Space-time Data So far we looked at the data that vary over space Now we add another
More informationAccelerated Advanced Algebra. Chapter 1 Patterns and Recursion Homework List and Objectives
Chapter 1 Patterns and Recursion Use recursive formulas for generating arithmetic, geometric, and shifted geometric sequences and be able to identify each type from their equations and graphs Write and
More informationLinear Regression. In this problem sheet, we consider the problem of linear regression with p predictors and one intercept,
Linear Regression In this problem sheet, we consider the problem of linear regression with p predictors and one intercept, y = Xβ + ɛ, where y t = (y 1,..., y n ) is the column vector of target values,
More informationToward an automatic real-time mapping system for radiation hazards
Toward an automatic real-time mapping system for radiation hazards Paul H. Hiemstra 1, Edzer J. Pebesma 2, Chris J.W. Twenhöfel 3, Gerard B.M. Heuvelink 4 1 Faculty of Geosciences / University of Utrecht
More informationSpatial analysis is the quantitative study of phenomena that are located in space.
c HYON-JUNG KIM, 2016 1 Introduction Spatial analysis is the quantitative study of phenomena that are located in space. Spatial data analysis usually refers to an analysis of the observations in which
More informationSPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA
SPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA D. Pokrajac Center for Information Science and Technology Temple University Philadelphia, Pennsylvania A. Lazarevic Computer
More informationBAYESIAN MODEL FOR SPATIAL DEPENDANCE AND PREDICTION OF TUBERCULOSIS
BAYESIAN MODEL FOR SPATIAL DEPENDANCE AND PREDICTION OF TUBERCULOSIS Srinivasan R and Venkatesan P Dept. of Statistics, National Institute for Research Tuberculosis, (Indian Council of Medical Research),
More information