Outline. Point Pattern Analysis Part I. Revisit IRP/CSR

Similar documents
CHAPTER VI Statistical Analysis of Experimental Data

Summary of the lecture in Biostatistics

KLT Tracker. Alignment. 1. Detect Harris corners in the first frame. 2. For each Harris corner compute motion between consecutive frames

Lecture Notes Types of economic variables

Simple Linear Regression

Chapter 4 (Part 1): Non-Parametric Classification (Sections ) Pattern Classification 4.3) Announcements

Introduction to local (nonparametric) density estimation. methods

Lecture 3. Sampling, sampling distributions, and parameter estimation

LECTURE - 4 SIMPLE RANDOM SAMPLING DR. SHALABH DEPARTMENT OF MATHEMATICS AND STATISTICS INDIAN INSTITUTE OF TECHNOLOGY KANPUR

Chapter 14 Logistic Regression Models

To use adaptive cluster sampling we must first make some definitions of the sampling universe:

Chapter 3 Sampling For Proportions and Percentages

Chapter 5 Properties of a Random Sample

STATISTICAL PROPERTIES OF LEAST SQUARES ESTIMATORS. x, where. = y - ˆ " 1

Statistics Descriptive and Inferential Statistics. Instructor: Daisuke Nagakura

Lecture 3 Probability review (cont d)

6. Nonparametric techniques

A Combination of Adaptive and Line Intercept Sampling Applicable in Agricultural and Environmental Studies

MEASURES OF DISPERSION

ENGI 4421 Joint Probability Distributions Page Joint Probability Distributions [Navidi sections 2.5 and 2.6; Devore sections

f f... f 1 n n (ii) Median : It is the value of the middle-most observation(s).

Sampling Theory MODULE V LECTURE - 14 RATIO AND PRODUCT METHODS OF ESTIMATION

hp calculators HP 30S Statistics Averages and Standard Deviations Average and Standard Deviation Practice Finding Averages and Standard Deviations

1. The weight of six Golden Retrievers is 66, 61, 70, 67, 92 and 66 pounds. The weight of six Labrador Retrievers is 54, 60, 72, 78, 84 and 67.

Chapter Business Statistics: A First Course Fifth Edition. Learning Objectives. Correlation vs. Regression. In this chapter, you learn:

2.28 The Wall Street Journal is probably referring to the average number of cubes used per glass measured for some population that they have chosen.

Idea is to sample from a different distribution that picks points in important regions of the sample space. Want ( ) ( ) ( ) E f X = f x g x dx

= 1. UCLA STAT 13 Introduction to Statistical Methods for the Life and Health Sciences. Parameters and Statistics. Measures of Centrality

Functions of Random Variables

Simulation Output Analysis

Centroids & Moments of Inertia of Beam Sections

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

Arithmetic Mean Suppose there is only a finite number N of items in the system of interest. Then the population arithmetic mean is

Chapter -2 Simple Random Sampling

Multiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades

Chapter 13, Part A Analysis of Variance and Experimental Design. Introduction to Analysis of Variance. Introduction to Analysis of Variance

Objectives of Multiple Regression

Econometric Methods. Review of Estimation

n -dimensional vectors follow naturally from the one

Chapter -2 Simple Random Sampling

X ε ) = 0, or equivalently, lim

Analysis of System Performance IN2072 Chapter 5 Analysis of Non Markov Systems

EP2200 Queueing theory and teletraffic systems. Queueing networks. Viktoria Fodor KTH EES/LCN KTH EES/LCN

Block-Based Compact Thermal Modeling of Semiconductor Integrated Circuits

L5 Polynomial / Spline Curves

GOALS The Samples Why Sample the Population? What is a Probability Sample? Four Most Commonly Used Probability Sampling Methods


Special Instructions / Useful Data

Chapter 13 Student Lecture Notes 13-1

ENGI 3423 Simple Linear Regression Page 12-01

ECON 482 / WH Hong The Simple Regression Model 1. Definition of the Simple Regression Model

Bayesian belief networks

is the score of the 1 st student, x

The TDT. (Transmission Disequilibrium Test) (Qualitative and quantitative traits) D M D 1 M 1 D 2 M 2 M 2D1 M 1

Lecture 2 - What are component and system reliability and how it can be improved?

Estimation of Stress- Strength Reliability model using finite mixture of exponential distributions

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

Evaluation of uncertainty in measurements

Unsupervised Learning and Other Neural Networks

Continuous Distributions

Lecture 1 Review of Fundamental Statistical Concepts

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Parameter, Statistic and Random Samples

The Mathematical Appendix

Midterm Exam 1, section 1 (Solution) Thursday, February hour, 15 minutes

Chapter 11 Systematic Sampling

Mean is only appropriate for interval or ratio scales, not ordinal or nominal.

Module 7: Probability and Statistics

Simple Linear Regression

Centroids Method of Composite Areas

Bayes (Naïve or not) Classifiers: Generative Approach

Naïve Bayes MIT Course Notes Cynthia Rudin

Nonparametric Techniques

Module 7. Lecture 7: Statistical parameter estimation

12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model

STRATIFIED SAMPLING IN AGRICULTURAL SURVEYS

Lecture 2: The Simple Regression Model

Chapter 4: Elements of Statistics

(b) By independence, the probability that the string 1011 is received correctly is

Chapter Two. An Introduction to Regression ( )

Feature Selection: Part 2. 1 Greedy Algorithms (continued from the last lecture)

A NEW LOG-NORMAL DISTRIBUTION

9 U-STATISTICS. Eh =(m!) 1 Eh(X (1),..., X (m ) ) i.i.d

Chapter 8: Statistical Analysis of Simulated Data

STA302/1001-Fall 2008 Midterm Test October 21, 2008

1 Onto functions and bijections Applications to Counting

Discrete Mathematics and Probability Theory Fall 2016 Seshia and Walrand DIS 10b

9.1 Introduction to the probit and logit models

Lecture 07: Poles and Zeros

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

THE ROYAL STATISTICAL SOCIETY HIGHER CERTIFICATE

Centers of Gravity - Centroids

Statistics. Correlational. Dr. Ayman Eldeib. Simple Linear Regression and Correlation. SBE 304: Linear Regression & Correlation 1/3/2018

STA 108 Applied Linear Models: Regression Analysis Spring Solution for Homework #1

ENGI 4421 Propagation of Error Page 8-01

Analysis of Variance with Weibull Data

(Monte Carlo) Resampling Technique in Validity Testing and Reliability Testing

Convergence of the Desroziers scheme and its relation to the lag innovation diagnostic

Chapter Statistics Background of Regression Analysis

Statistical modelling and latent variables (2)

Transcription:

Pot Patter Aalyss Part I Outle Revst IRP/CSR, frst- ad secod order effects What s pot patter aalyss (PPA)? Desty-based pot patter measures Dstace-based pot patter measures Revst IRP/CSR Equal probablty: ay pot has equal probablty of beg ay posto or, equvaletly, each small subarea of the map has a equal chace of recevg a pot Idepedece: the postog of ay pot s depedet of the postog of ay other pot 9 8 7 6 5 4 3 2 2 4 6 8 k x P( k,, x) = k x x k λ e P( k) = k! λ k ad λ = x

Frst- & Secod Order Effects IRP/CSR s ot realstc geography Frst-Order effect: varatos the desty of a process across space due to varatos evromet propertes: No equal probablty Secod-Order effect: teracto betwee locatos: NO depedece What Is Pot Patter Aalyss (PPA) Pot patters, where the oly data are the locatos of a set of pot objects Represet the smplest possble spatal data Examples Hot-spot detecto (crme, dsease, geologcal actvtes) Dsease & evrometal relatos Freeway accdets Vegetato, archaeologcal studes What Is PPA (Cot.) I a pot process the basc propertes of the process are set by a sgle parameter, the probablty that ay small area wll receve a pot

What Is PPA (Cot.) Requremets for a set of evets to costtute a pot patter Mapable patter o the plae Objectvely determed study area A eumerato or cesus of all ettes of terest, ot a sample A oe-to-oe correspodece betwee objects the study area ad evets the patter Evet locatos must be proper (ot be the cetrods of polygos) Descrbg a Pot Patter Pot desty (frst-order or secod-order?) Pot separato (frst-order or secodorder?) Whe frst-order effects are marked, absolute locato s a mportat determat of observatos, ad a pot patter has clear varatos across space the umber of evets per ut area are observed Whe secod-order effects are strog, there s teracto betwee locatos, depedg o the dstace betwee them, ad thus relatve locato s mportat Descrbg a Pot Patter I practce t s close to mpossble to dstgush from varato the evromet or teracto Low desty larger separato Hgh desty smaller separato Frst-order or secod order?

Some Notatos s (x, y ) A set of locatos S wth evets S = { s, s2,..., s,..., s} Each evet s has two coordates: s = ( x, y ) The study rego A has a area a Descrbg a Pot Patter Mea Ceter: the mea X coordate ad the mea Y coordate x y = = s = ( µ = x, µ y ), Stadard Dstace: a measure of how dspersed the evets are aroud ther mea ceter d = = 2 2 [( x µ ) + ( y µ ) ] x y Descrbg a Pot Patter A summary crcle ca the be plotted for the pot patter, cetered at the mea ceter wth radus d (SD) If the stadard dstace s computed separately for each axs, a summary ellpse ca be obtaed (2 ds) Summary crcle Summary ellpse

Measurg Pot Patter Two basc famles Desty based measures: st order effect Dstace based measures: 2 d order effect Desty-Based Measures Frst-order effect Three methods Crude desty: Overall Quadrat cout: couts of evets each quadrat Desty estmato: cotuous Sestve to the defto of the study area: MAUP & Edge Effect Crude Desty Overall testy of evets the study area λ = o.( S A = ) a a Descrptve statstcs for summary purpose too smple

Quadrat Couts Cout the umbers of evets a set of equal-area quadrats ad record these couts as a frequecy dstrbuto Equally-spaced patters wll have most quadrats wth smlar couts Clustered patter wll have a few hgh cout quadrates ad may empty oes Two approaches: Exhaustve cesus Radom samplg Quadrat Couts Approach I: Exhaustve cesus of quadrats Completely fll the study rego wth o overlaps The choce of org, quadrat oretato, ad quadrat sze affects the observed frequecy dstrbuto MAUP Quadrat Couts (Cot.) Approach II: Radom samplg approach Possble to crease the sample sze smply by addg more quadrats (for sparse patters) May descrbe a pot patter wthout havg complete data o the etre patter

Quadrat Couts (Cot.) Other shapes of quadrats Desty Estmato Assume: the patter has a desty aywhere the study area ot just locatos where there s a evet Estmated by coutg the umber of evets a rego, or kerel, cetered at the locato where the estmate s to be made Varato of Quadrat Couts Two approaches Smple desty estmato Kerel desty estmato Smple Desty Estmato Smply cout the umber of evets of a patter a rego (usually, crcle) the calculate the desty by dvdg t wth the area of the rego Badwdth r Too large or too small? what s the proper sze? Desty does ot smoothly vary across space ) o.[ S C( p, r)] λ p = 2 πr C(p,r) s a crcle of radus e cetered at the locato of terest p

Kerel Desty Estmato Itroduce a dstace effect wth a kerel fucto Weght earby evets more heavly tha dstat oes estmatg the local desty e.g. a pot at half of badwdth r oly couts for ½ pot Kerel Desty Estmato (Cot.) Badwdth r mpact the smoothess of the resultg desty feld depeds Narrow kerels produce bumpy surfaces Wde kerels produce smooth surfaces Applcatos of Desty-Based Methods Vsualze a pot patter to detect hot spots Check whether or ot that process s frst-order statoary from the local testy varatos Lk pot objects to other geographc data (e.g. dsease ad polluto)

Summary Pot patter s the oly patter st order vs. 2 d order SV Descrptve summary of pot patter Mea ceter ad stadard dstace Desty measures Quadrat Couts Desty Estmato Dstace-Based Measures Secod-order effect Look at patters wth the dstaces amog pot evets Two approaches: Nearest Neghbor Dstace (NND) Dstace fuctos G fucto F fucto K fucto Nearest Neghbor Dstace NND: the dstace from a evet to the earest evet Usually Eucldea dstace Pythagorea Theorem: d( s, sj) = ( x x j ) + ( y y j ) d ( s ) = m( d( s, s )) m j Mea NND: Summarzes all the earest-eghbor dstaces by a sgle mea value: d s d = = m( ) m Way too smple lose a lot useful formato

Expected Behavor of Mea NND Clustered: All NND are short small mea NND Equally-spaced: Mmum dstaces are loger larger mea NND G Fucto G(d) exames the cumulatve frequecy dstrbuto of the NNDs btw the evet pots What fracto of all the NNDs the patter are less tha a specfed dstace (d) Smplest dstace fucto Illustrato

G fucto (radom) Expected Behavor of G(d).75 Clustered: Rapd.5 crease, starts at shorter dstace Equally-spaced: Also rapd crease but starts at a.25 much loger dstace 2 3 4 G fucto (clustered) G fucto (eve-spacg) Dstace.75.75.5.5.25.25 2 3 4 Dstace 2 3 4 Dstace F Fucto F(d) s also the cumulatve frequecy dstrbuto of the certa NNDs BUT, NNDs are calculated btw a group of radom pots the study area ad evet pots Icreased sample sze smoother curve G(d) vs. F(d) Both use the dea of NNDs Oly NNDs are used F(d) s smoother tha G(d) Combe both more revealg. G ad F fuctos 9. 8. 7. 6. 5. 4. 3. 2.... 2. 4. 6. 8...75.5.25 2 3 4 Dstace F fucto G fucto

G ad F fuctos (eve spacg).75.5 F fucto G fucto G ad F fuctos (radom).25.75 2 3 4 Dstace.5.75 G ad F fuctos (clustered) F fucto G fucto.25.5 2 3 4 Dstace.25 F fucto G fucto 2 3 4 Dstace Rpley s K Fucto K(d) makes use of all ter-evet dstace Makes t a complex fucto to calculate But, t cludes all possble dstace formato about the patter stead of NNDs Calculatg the K(d) For each dstace d:. Determe how may evets le wth that dstace (d) of each evet 2. Calculate a average umber of evets at the dstace d 3. Dvde by overall study area evet desty

Iterpretato of K(d) Revew (PPT Statstcs) Desty-based : Quadrat Couts Kerel desty estmato Dstace-based: Mea NND G & F fuctos K fucto