Areal Unit Data Regular or Irregular Grids or Lattices Large Point-referenced Datasets
|
|
- Ronald Fields
- 6 years ago
- Views:
Transcription
1 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
2 Areal Unit Data Regular or Irregular Grids or Lattices Large Point-referenced Datasets Is there spatial pattern? Do we want to smooth the data? Chapter 3: Basics of Areal Data Models p. 1/18
3 Areal Unit Data Regular or Irregular Grids or Lattices Large Point-referenced Datasets Is there spatial pattern? Do we want to smooth the data? Inference for new areal units? Chapter 3: Basics of Areal Data Models p. 1/18
4 Areal Unit Data Regular or Irregular Grids or Lattices Large Point-referenced Datasets Is there spatial pattern? Do we want to smooth the data? Inference for new areal units? Descriptive/algorithmic vs. Model-based Chapter 3: Basics of Areal Data Models p. 1/18
5 Areal unit data Actual Transformed SIDS Rates < >3.5 Chapter 3: Basics of Areal Data Models p. 2/18
6 Proximity matrices W, entries w ij (with w ii = 0). Choices for w ij : w ij = 1 if i,j share a common boundary (possibly a common vertex) w ij is an inverse distance between units w ij = 1 if distance between units is K w ij = 1 for m nearest neighbors. Chapter 3: Basics of Areal Data Models p. 3/18
7 Proximity matrices W, entries w ij (with w ii = 0). Choices for w ij : w ij = 1 if i,j share a common boundary (possibly a common vertex) w ij is an inverse distance between units w ij = 1 if distance between units is K w ij = 1 for m nearest neighbors. W is typically symmetric, but need not be Chapter 3: Basics of Areal Data Models p. 3/18
8 Proximity matrices W, entries w ij (with w ii = 0). Choices for w ij : w ij = 1 if i,j share a common boundary (possibly a common vertex) w ij is an inverse distance between units w ij = 1 if distance between units is K w ij = 1 for m nearest neighbors. W is typically symmetric, but need not be W : standardize row i by w i+ = j w ij (so matrix is now row stochastic, but probably no longer symmetric). Chapter 3: Basics of Areal Data Models p. 3/18
9 Proximity matrices W, entries w ij (with w ii = 0). Choices for w ij : w ij = 1 if i,j share a common boundary (possibly a common vertex) w ij is an inverse distance between units w ij = 1 if distance between units is K w ij = 1 for m nearest neighbors. W is typically symmetric, but need not be W : standardize row i by w i+ = j w ij (so matrix is now row stochastic, but probably no longer symmetric). W elements often called weights ; interpretation Chapter 3: Basics of Areal Data Models p. 3/18
10 Proximity matrices W, entries w ij (with w ii = 0). Choices for w ij : w ij = 1 if i,j share a common boundary (possibly a common vertex) w ij is an inverse distance between units w ij = 1 if distance between units is K w ij = 1 for m nearest neighbors. W is typically symmetric, but need not be W : standardize row i by w i+ = j w ij (so matrix is now row stochastic, but probably no longer symmetric). W elements often called weights ; interpretation Could also define first-order neighbors W (1), second-order neighbors W (2), etc. Chapter 3: Basics of Areal Data Models p. 3/18
11 Measures of spatial association Moran s I: essentially an areal covariogram" I = n i j w ij(y i Ȳ )(Y j Ȳ ) ( i j w ij) i (Y i Ȳ )2 Chapter 3: Basics of Areal Data Models p. 4/18
12 Measures of spatial association Moran s I: essentially an areal covariogram" I = n i j w ij(y i Ȳ )(Y j Ȳ ) ( i j w ij) i (Y i Ȳ )2 Geary s C: essentially an areal variogram" C = (n 1) i j w ij(y i Y j ) 2 ( i j w ij) i (Y i Ȳ )2 Chapter 3: Basics of Areal Data Models p. 4/18
13 Measures of spatial association Moran s I: essentially an areal covariogram" I = n i j w ij(y i Ȳ )(Y j Ȳ ) ( i j w ij) i (Y i Ȳ )2 Geary s C: essentially an areal variogram" C = (n 1) i j w ij(y i Y j ) 2 ( i j w ij) i (Y i Ȳ )2 Both are asymptotically normal if Y i are i.i.d.; Moran has mean 1/(n 1) 0, Geary has mean 1 Chapter 3: Basics of Areal Data Models p. 4/18
14 Measures of spatial association Moran s I: essentially an areal covariogram" I = n i j w ij(y i Ȳ )(Y j Ȳ ) ( i j w ij) i (Y i Ȳ )2 Geary s C: essentially an areal variogram" C = (n 1) i j w ij(y i Y j ) 2 ( i j w ij) i (Y i Ȳ )2 Both are asymptotically normal if Y i are i.i.d.; Moran has mean 1/(n 1) 0, Geary has mean 1 Significance testing by comparing to a collection of say 1000 random permutations of the Y i Chapter 3: Basics of Areal Data Models p. 4/18
15 Measures of spatial association (cont d) < >563 Figure 1: Choropleth map of 1999 average verbal SAT scores, lower 48 U.S. states. Chapter 3: Basics of Areal Data Models p. 5/18
16 Measures of spatial association (cont d) For these data, the spatial.cor function in S+SpatialStats gives a Moran s I of , with associated standard error estimate very strong evidence against H 0 : no spatial correlation Chapter 3: Basics of Areal Data Models p. 6/18
17 Measures of spatial association (cont d) For these data, the spatial.cor function in S+SpatialStats gives a Moran s I of , with associated standard error estimate very strong evidence against H 0 : no spatial correlation spatial.cor also gives a Geary s C of , with associated standard error estimate again, very strong evidence against H 0 (departure from 1) Chapter 3: Basics of Areal Data Models p. 6/18
18 Measures of spatial association (cont d) For these data, the spatial.cor function in S+SpatialStats gives a Moran s I of , with associated standard error estimate very strong evidence against H 0 : no spatial correlation spatial.cor also gives a Geary s C of , with associated standard error estimate again, very strong evidence against H 0 (departure from 1) Warning: These data have not been adjusted for covariates, such as the proportion of students who take the exam (Midwestern colleges have historically relied on the ACT, not the SAT; only the best and brightest students in these states would bother taking the SAT) Chapter 3: Basics of Areal Data Models p. 6/18
19 Measures of spatial association (cont d) For these data, the spatial.cor function in S+SpatialStats gives a Moran s I of , with associated standard error estimate very strong evidence against H 0 : no spatial correlation spatial.cor also gives a Geary s C of , with associated standard error estimate again, very strong evidence against H 0 (departure from 1) Warning: These data have not been adjusted for covariates, such as the proportion of students who take the exam (Midwestern colleges have historically relied on the ACT, not the SAT; only the best and brightest students in these states would bother taking the SAT) the map, I, and C all motivate the search for spatial covariates! Chapter 3: Basics of Areal Data Models p. 6/18
20 Correlogram (via Moran s I) rho(d) distance Replace w ij with w (1) ij taken from W (1) I (1) Chapter 3: Basics of Areal Data Models p. 7/18
21 Correlogram (via Moran s I) rho(d) distance Replace w ij with w (1) ij taken from W (1) I (1) Replace w ij with w (2) ij taken from W (2) I (2), etc. Chapter 3: Basics of Areal Data Models p. 7/18
22 Correlogram (via Moran s I) rho(d) distance Replace w ij with w (1) ij taken from W (1) I (1) Replace w ij with w (2) ij taken from W (2) I (2), etc. Plot I (r) vs. r; an initial decline in r followed by variation around 0 spatial pattern! Chapter 3: Basics of Areal Data Models p. 7/18
23 Correlogram (via Moran s I) rho(d) distance Replace w ij with w (1) ij taken from W (1) I (1) Replace w ij with w (2) ij taken from W (2) I (2), etc. Plot I (r) vs. r; an initial decline in r followed by variation around 0 spatial pattern! spatial analogue of the temporal lag autocorrelation plot Chapter 3: Basics of Areal Data Models p. 7/18
24 Rasterized binary data map NORTH SOUTH land use classification non-forest forest Chapter 3: Basics of Areal Data Models p. 8/18
25 Binary data correlogram logodds N NW distance NE E A version for a binary map, using two-way tables and log odds ratios at the pixel level Chapter 3: Basics of Areal Data Models p. 9/18
26 Binary data correlogram logodds N NW distance NE E A version for a binary map, using two-way tables and log odds ratios at the pixel level Note strongest pattern is to the north (N), but in no direction are the values 0 even at 40 km Chapter 3: Basics of Areal Data Models p. 9/18
27 Spatial smoothers To smooth Y i, replace with Ŷi = i w ijy j w i+ Chapter 3: Basics of Areal Data Models p. 10/18
28 Spatial smoothers To smooth Y i, replace with Ŷi = i w ijy j w i+ More generally, we could include the value actually observed for unit i, and revise our smoother to (1 α)y i + αŷi For 0 < α < 1, this is a linear (convex) combination in shrinkage" form Chapter 3: Basics of Areal Data Models p. 10/18
29 Spatial smoothers To smooth Y i, replace with Ŷi = i w ijy j w i+ More generally, we could include the value actually observed for unit i, and revise our smoother to (1 α)y i + αŷi For 0 < α < 1, this is a linear (convex) combination in shrinkage" form Finally, we could try model-based smoothing, i.e., based on E(Y i Data), i.e., the mean of the predictive distribution. Smoothers then emerge as byproducts of the hierarchical spatial models we use to explain the Y i s Chapter 3: Basics of Areal Data Models p. 10/18
30 Markov random fields Consider Y = (Y 1,Y 2,...,Y n ) and the set {p(y i y j, j i)} Chapter 3: Basics of Areal Data Models p. 11/18
31 Markov random fields Consider Y = (Y 1,Y 2,...,Y n ) and the set {p(y i y j, j i)} We know p(y 1,y 2,...y n ) determines {p(y i y j,j i)} (the set of full conditional distributions) Chapter 3: Basics of Areal Data Models p. 11/18
32 Markov random fields Consider Y = (Y 1,Y 2,...,Y n ) and the set {p(y i y j, j i)} We know p(y 1,y 2,...y n ) determines {p(y i y j,j i)} (the set of full conditional distributions) Does {p(y i y j,j i)} determine p(y 1,y 2,...y n )??? Chapter 3: Basics of Areal Data Models p. 11/18
33 Markov random fields Consider Y = (Y 1,Y 2,...,Y n ) and the set {p(y i y j, j i)} We know p(y 1,y 2,...y n ) determines {p(y i y j,j i)} (the set of full conditional distributions) Does {p(y i y j,j i)} determine p(y 1,y 2,...y n )??? If the full conditionals are compatible, then by Brook s Lemma, p(y 1,...,y n ) = p(y 1 y 2,...,y n ) p(y 10 y 2,...,y n ) p(y 2 y 10,y 3,...,y n ) p(y 20 y 10,y 3,...,y n )... p(y n y 10,...,y n 1,0 ) p(y n0 y 10,...,y n 1,0 ) p(y 10,...,y n0 ) If left side is proper, the fact that it integrates to 1 determines the normalizing constant! Chapter 3: Basics of Areal Data Models p. 11/18
34 Local modeling Suppose we specify the full conditionals such that p(y i y j,j i) = p(y i y j i ), where i is the set of neighbors of cell (region) i. When does {p(y i y j i )} determine p(y 1,y 2,...y n )? Chapter 3: Basics of Areal Data Models p. 12/18
35 Local modeling Suppose we specify the full conditionals such that p(y i y j,j i) = p(y i y j i ), where i is the set of neighbors of cell (region) i. When does {p(y i y j i )} determine p(y 1,y 2,...y n )? Def n: a clique is a set of cells such that each element is a neighbor of every other element Chapter 3: Basics of Areal Data Models p. 12/18
36 Local modeling Suppose we specify the full conditionals such that p(y i y j,j i) = p(y i y j i ), where i is the set of neighbors of cell (region) i. When does {p(y i y j i )} determine p(y 1,y 2,...y n )? Def n: a clique is a set of cells such that each element is a neighbor of every other element Def n: a potential function of order k is a function of k arguments that is exchangeable in these arguments Chapter 3: Basics of Areal Data Models p. 12/18
37 Local modeling Suppose we specify the full conditionals such that p(y i y j,j i) = p(y i y j i ), where i is the set of neighbors of cell (region) i. When does {p(y i y j i )} determine p(y 1,y 2,...y n )? Def n: a clique is a set of cells such that each element is a neighbor of every other element Def n: a potential function of order k is a function of k arguments that is exchangeable in these arguments Def n: p(y 1,...,y n ) is a Gibbs distribution if, as a function of the y i, it is a product of potentials on cliques Chapter 3: Basics of Areal Data Models p. 12/18
38 Local modeling Suppose we specify the full conditionals such that p(y i y j,j i) = p(y i y j i ), where i is the set of neighbors of cell (region) i. When does {p(y i y j i )} determine p(y 1,y 2,...y n )? Def n: a clique is a set of cells such that each element is a neighbor of every other element Def n: a potential function of order k is a function of k arguments that is exchangeable in these arguments Def n: p(y 1,...,y n ) is a Gibbs distribution if, as a function of the y i, it is a product of potentials on cliques Example: For binary data, k = 2, we might take Q(y i,y j ) = I(y i = y j ) = y i y j + (1 y i )(1 y j ) Chapter 3: Basics of Areal Data Models p. 12/18
39 Local modeling Cliques of size 1 independence Chapter 3: Basics of Areal Data Models p. 13/18
40 Local modeling Cliques of size 1 independence Cliques of size 2 pairwise difference form p(y 1,y 2,...y n ) exp 1 2τ 2 (y i y j ) 2 I(i j) and therefore p(y i y j,j i) = N( j i y j /m i,τ 2 /m i ), where m i is the number of neighbors of i i,j Chapter 3: Basics of Areal Data Models p. 13/18
41 Local modeling Cliques of size 1 independence Cliques of size 2 pairwise difference form p(y 1,y 2,...y n ) exp 1 2τ 2 (y i y j ) 2 I(i j) and therefore p(y i y j,j i) = N( j i y j /m i,τ 2 /m i ), where m i is the number of neighbors of i Hammersley-Clifford Theorem: If we have a Markov Random Field (i.e., {p(y i y j i )} uniquely determine p(y 1,y 2,...y n )), then the latter is a Gibbs distribution i,j Chapter 3: Basics of Areal Data Models p. 13/18
42 Local modeling Cliques of size 1 independence Cliques of size 2 pairwise difference form p(y 1,y 2,...y n ) exp 1 2τ 2 (y i y j ) 2 I(i j) and therefore p(y i y j,j i) = N( j i y j /m i,τ 2 /m i ), where m i is the number of neighbors of i Hammersley-Clifford Theorem: If we have a Markov Random Field (i.e., {p(y i y j i )} uniquely determine p(y 1,y 2,...y n )), then the latter is a Gibbs distribution Geman and Geman result : If we have a joint Gibbs distribution, then we have a Markov Random Field i,j Chapter 3: Basics of Areal Data Models p. 13/18
43 Conditional autoregressive (CAR) model Gaussian (autonormal) case p(y i y j,j i) = N j b ij y j,τ 2 i Using Brook s Lemma we can obtain { p(y 1,y 2,...y n ) exp 1 } 2 y D 1 (I B)y where B = {b ij } and D is diagonal with D ii = τ 2 i. Chapter 3: Basics of Areal Data Models p. 14/18
44 Conditional autoregressive (CAR) model Gaussian (autonormal) case p(y i y j,j i) = N j b ij y j,τ 2 i Using Brook s Lemma we can obtain { p(y 1,y 2,...y n ) exp 1 } 2 y D 1 (I B)y where B = {b ij } and D is diagonal with D ii = τ 2 i. suggests a multivariate normal distribution with µ Y = 0 and Σ Y = (I B) 1 D Chapter 3: Basics of Areal Data Models p. 14/18
45 Conditional autoregressive (CAR) model Gaussian (autonormal) case p(y i y j,j i) = N j b ij y j,τ 2 i Using Brook s Lemma we can obtain { p(y 1,y 2,...y n ) exp 1 } 2 y D 1 (I B)y where B = {b ij } and D is diagonal with D ii = τ 2 i. suggests a multivariate normal distribution with µ Y = 0 and Σ Y = (I B) 1 D D 1 (I B) symmetric requires b ij τ 2 i = b ji τ 2 j for all i,j Chapter 3: Basics of Areal Data Models p. 14/18
46 CAR Model (cont d) Returning to W, let b ij = w ij /w i+ and τ 2 i = τ 2 /w i+, so p(y 1,y 2,...y n ) exp { 1 } 2τ 2y (D w W)y where D w is diagonal with (D w ) ii = w i+ and thus p(y 1,y 2,...y n ) exp 1 2τ 2 w ij (y i y j ) 2 i j Intrinsic autoregressive (IAR) model! Chapter 3: Basics of Areal Data Models p. 15/18
47 CAR Model (cont d) Returning to W, let b ij = w ij /w i+ and τ 2 i = τ 2 /w i+, so p(y 1,y 2,...y n ) exp { 1 } 2τ 2y (D w W)y where D w is diagonal with (D w ) ii = w i+ and thus p(y 1,y 2,...y n ) exp 1 2τ 2 w ij (y i y j ) 2 i j Intrinsic autoregressive (IAR) model! Joint distribution is improper since (D w W)1 = 0, so requires a constraint say, i y i = 0 Chapter 3: Basics of Areal Data Models p. 15/18
48 CAR Model (cont d) Returning to W, let b ij = w ij /w i+ and τ 2 i = τ 2 /w i+, so p(y 1,y 2,...y n ) exp { 1 } 2τ 2y (D w W)y where D w is diagonal with (D w ) ii = w i+ and thus p(y 1,y 2,...y n ) exp 1 2τ 2 w ij (y i y j ) 2 i j Intrinsic autoregressive (IAR) model! Joint distribution is improper since (D w W)1 = 0, so requires a constraint say, i y i = 0 Not a data model a random effects model! Chapter 3: Basics of Areal Data Models p. 15/18
49 CAR Model Issues With τ 2 unknown, appropriate exponent for τ 2? Hodges et al.: n I, where I is the # of islands Chapter 3: Basics of Areal Data Models p. 16/18
50 CAR Model Issues With τ 2 unknown, appropriate exponent for τ 2? Hodges et al.: n I, where I is the # of islands Proper version: replace D w W by D w ρw, and choose ρ so that Σ y = (D w ρw) 1 exists! This in turn implies Y i Y j i N(ρ j w ijy j,τ 2 /m i ) Chapter 3: Basics of Areal Data Models p. 16/18
51 CAR Model Issues With τ 2 unknown, appropriate exponent for τ 2? Hodges et al.: n I, where I is the # of islands Proper version: replace D w W by D w ρw, and choose ρ so that Σ y = (D w ρw) 1 exists! This in turn implies Y i Y j i N(ρ j w ijy j,τ 2 /m i ) To ρ or not to ρ? (that is the question) Chapter 3: Basics of Areal Data Models p. 16/18
52 CAR Model Issues With τ 2 unknown, appropriate exponent for τ 2? Hodges et al.: n I, where I is the # of islands Proper version: replace D w W by D w ρw, and choose ρ so that Σ y = (D w ρw) 1 exists! This in turn implies Y i Y j i N(ρ j w ijy j,τ 2 /m i ) To ρ or not to ρ? (that is the question) Advantages: Chapter 3: Basics of Areal Data Models p. 16/18
53 CAR Model Issues With τ 2 unknown, appropriate exponent for τ 2? Hodges et al.: n I, where I is the # of islands Proper version: replace D w W by D w ρw, and choose ρ so that Σ y = (D w ρw) 1 exists! This in turn implies Y i Y j i N(ρ j w ijy j,τ 2 /m i ) To ρ or not to ρ? (that is the question) Advantages: makes distribution proper Chapter 3: Basics of Areal Data Models p. 16/18
54 CAR Model Issues With τ 2 unknown, appropriate exponent for τ 2? Hodges et al.: n I, where I is the # of islands Proper version: replace D w W by D w ρw, and choose ρ so that Σ y = (D w ρw) 1 exists! This in turn implies Y i Y j i N(ρ j w ijy j,τ 2 /m i ) To ρ or not to ρ? (that is the question) Advantages: makes distribution proper adds parametric flexibility Chapter 3: Basics of Areal Data Models p. 16/18
55 CAR Model Issues With τ 2 unknown, appropriate exponent for τ 2? Hodges et al.: n I, where I is the # of islands Proper version: replace D w W by D w ρw, and choose ρ so that Σ y = (D w ρw) 1 exists! This in turn implies Y i Y j i N(ρ j w ijy j,τ 2 /m i ) To ρ or not to ρ? (that is the question) Advantages: makes distribution proper adds parametric flexibility ρ = 0 interpretable as independence Chapter 3: Basics of Areal Data Models p. 16/18
56 CAR Models with ρ parameter Disadvantages: why should we expect Y i to be a proportion of average of neighbors a sensible spatial interpretation? Chapter 3: Basics of Areal Data Models p. 17/18
57 CAR Models with ρ parameter Disadvantages: why should we expect Y i to be a proportion of average of neighbors a sensible spatial interpretation? calibration of ρ as a correlation, e.g., ρ = 0.80 yields 0.1 Moran s I 0.15, ρ = 0.90 yields 0.2 Moran s I 0.25, ρ = 0.99 yields Moran s I 0.5 Chapter 3: Basics of Areal Data Models p. 17/18
58 CAR Models with ρ parameter Disadvantages: why should we expect Y i to be a proportion of average of neighbors a sensible spatial interpretation? calibration of ρ as a correlation, e.g., ρ = 0.80 yields 0.1 Moran s I 0.15, ρ = 0.90 yields 0.2 Moran s I 0.25, ρ = 0.99 yields Moran s I 0.5 So, used with random effects, scope of spatial pattern may be limited Chapter 3: Basics of Areal Data Models p. 17/18
59 Comments The CAR specifies Σ 1 Y, not Σ Y (as in the point-level modeling case), so does not directly model association Chapter 3: Basics of Areal Data Models p. 18/18
60 Comments The CAR specifies Σ 1 Y, not Σ Y (as in the point-level modeling case), so does not directly model association (Σ 1 Y ) ii = 1/τ 2 i ; (Σ 1 Y ) ij = 0 cond l independence Chapter 3: Basics of Areal Data Models p. 18/18
61 Comments The CAR specifies Σ 1 Y, not Σ Y (as in the point-level modeling case), so does not directly model association (Σ 1 Y ) ii = 1/τ 2 i ; (Σ 1 Y ) ij = 0 cond l independence Prediction at new sites is ad hoc, in that if p(y 0 y 1,y 2,...y n ) = N( j w 0j y j /w 0+,τ 2 /w 0+ ) then p(y 0,y 1,...y n ) well-defined but not CAR Chapter 3: Basics of Areal Data Models p. 18/18
62 Comments The CAR specifies Σ 1 Y, not Σ Y (as in the point-level modeling case), so does not directly model association (Σ 1 Y ) ii = 1/τ 2 i ; (Σ 1 Y ) ij = 0 cond l independence Prediction at new sites is ad hoc, in that if p(y 0 y 1,y 2,...y n ) = N( j w 0j y j /w 0+,τ 2 /w 0+ ) then p(y 0,y 1,...y n ) well-defined but not CAR Non-Gaussian case: For binary data, the autologistic: p(y i y j,j i) exp φ w ij I(y i = y j ) j Chapter 3: Basics of Areal Data Models p. 18/18
Areal data models. Spatial smoothers. Brook s Lemma and Gibbs distribution. CAR models Gaussian case Non-Gaussian case
Areal data models Spatial smoothers Brook s Lemma and Gibbs distribution CAR models Gaussian case Non-Gaussian case SAR models Gaussian case Non-Gaussian case CAR vs. SAR STAR models Inference for areal
More informationLecture 18. Models for areal data. Colin Rundel 03/22/2017
Lecture 18 Models for areal data Colin Rundel 03/22/2017 1 areal / lattice data 2 Example - NC SIDS SID79 3 EDA - Moran s I If we have observations at n spatial locations s 1,... s n ) I = n i=1 n n j=1
More informationTechnical Vignette 5: Understanding intrinsic Gaussian Markov random field spatial models, including intrinsic conditional autoregressive models
Technical Vignette 5: Understanding intrinsic Gaussian Markov random field spatial models, including intrinsic conditional autoregressive models Christopher Paciorek, Department of Statistics, University
More informationCSC 412 (Lecture 4): Undirected Graphical Models
CSC 412 (Lecture 4): Undirected Graphical Models Raquel Urtasun University of Toronto Feb 2, 2016 R Urtasun (UofT) CSC 412 Feb 2, 2016 1 / 37 Today Undirected Graphical Models: Semantics of the graph:
More informationAsymptotic standard errors of MLE
Asymptotic standard errors of MLE Suppose, in the previous example of Carbon and Nitrogen in soil data, that we get the parameter estimates For maximum likelihood estimation, we can use Hessian matrix
More informationSummary STK 4150/9150
STK4150 - Intro 1 Summary STK 4150/9150 Odd Kolbjørnsen May 22 2017 Scope You are expected to know and be able to use basic concepts introduced in the book. You knowledge is expected to be larger than
More informationSpatial Autocorrelation (2) Spatial Weights
Spatial Autocorrelation (2) Spatial Weights Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline
More informationMarkov Random Fields
Markov Random Fields Umamahesh Srinivas ipal Group Meeting February 25, 2011 Outline 1 Basic graph-theoretic concepts 2 Markov chain 3 Markov random field (MRF) 4 Gauss-Markov random field (GMRF), and
More informationWeb Appendices: Hierarchical and Joint Site-Edge Methods for Medicare Hospice Service Region Boundary Analysis
Web Appendices: Hierarchical and Joint Site-Edge Methods for Medicare Hospice Service Region Boundary Analysis Haijun Ma, Bradley P. Carlin and Sudipto Banerjee December 8, 2008 Web Appendix A: Selecting
More informationBasics of Geographic Analysis in R
Basics of Geographic Analysis in R Spatial Autocorrelation and Spatial Weights Yuri M. Zhukov GOV 2525: Political Geography February 25, 2013 Outline 1. Introduction 2. Spatial Data and Basic Visualization
More informationMarkov Random Fields
Markov Random Fields 1. Markov property The Markov property of a stochastic sequence {X n } n 0 implies that for all n 1, X n is independent of (X k : k / {n 1, n, n + 1}), given (X n 1, X n+1 ). Another
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 informationLecture 7 Autoregressive Processes in Space
Lecture 7 Autoregressive Processes in Space Dennis Sun Stanford University Stats 253 July 8, 2015 1 Last Time 2 Autoregressive Processes in Space 3 Estimating Parameters 4 Testing for Spatial Autocorrelation
More informationPart 6: Multivariate Normal and Linear Models
Part 6: Multivariate Normal and Linear Models 1 Multiple measurements Up until now all of our statistical models have been univariate models models for a single measurement on each member of a sample of
More informationAreal data. Infant mortality, Auckland NZ districts. Number of plant species in 20cm x 20 cm patches of alpine tundra. Wheat yield
Areal data Reminder about types of data Geostatistical data: Z(s) exists everyhere, varies continuously Can accommodate sudden changes by a model for the mean E.g., soil ph, two soil types with different
More informationHierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets
Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets Abhirup Datta 1 Sudipto Banerjee 1 Andrew O. Finley 2 Alan E. Gelfand 3 1 University of Minnesota, Minneapolis,
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 informationIntroduction to Graphical Models
Introduction to Graphical Models STA 345: Multivariate Analysis Department of Statistical Science Duke University, Durham, NC, USA Robert L. Wolpert 1 Conditional Dependence Two real-valued or vector-valued
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 informationBayesian SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling
Bayesian SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling Jon Wakefield Departments of Statistics and Biostatistics University of Washington 1 / 37 Lecture Content Motivation
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 informationClassification. Chapter Introduction. 6.2 The Bayes classifier
Chapter 6 Classification 6.1 Introduction Often encountered in applications is the situation where the response variable Y takes values in a finite set of labels. For example, the response Y could encode
More informationMultivariate spatial modeling
Multivariate spatial modeling Point-referenced spatial data often come as multivariate measurements at each location Chapter 7: Multivariate Spatial Modeling p. 1/21 Multivariate spatial modeling Point-referenced
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 informationMATH 829: Introduction to Data Mining and Analysis Graphical Models I
MATH 829: Introduction to Data Mining and Analysis Graphical Models I Dominique Guillot Departments of Mathematical Sciences University of Delaware May 2, 2016 1/12 Independence and conditional independence:
More informationMarkov random fields. The Markov property
Markov random fields The Markov property Discrete time: (X k X k!1,x k!2,... = (X k X k!1 A time symmetric version: (X k! X!k = (X k X k!1,x k+1 A more general version: Let A be a set of indices >k, B
More informationStatistics 352: Spatial statistics. Jonathan Taylor. Department of Statistics. Models for discrete data. Stanford University.
352: 352: Models for discrete data April 28, 2009 1 / 33 Models for discrete data 352: Outline Dependent discrete data. Image data (binary). Ising models. Simulation: Gibbs sampling. Denoising. 2 / 33
More informationChapter 3 - Temporal processes
STK4150 - Intro 1 Chapter 3 - Temporal processes Odd Kolbjørnsen and Geir Storvik January 23 2017 STK4150 - Intro 2 Temporal processes Data collected over time Past, present, future, change Temporal aspect
More informationOutline. Overview of Issues. Spatial Regression. Luc Anselin
Spatial Regression Luc Anselin University of Illinois, Urbana-Champaign http://www.spacestat.com Outline Overview of Issues Spatial Regression Specifications Space-Time Models Spatial Latent Variable Models
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 informationMAP Examples. Sargur Srihari
MAP Examples Sargur srihari@cedar.buffalo.edu 1 Potts Model CRF for OCR Topics Image segmentation based on energy minimization 2 Examples of MAP Many interesting examples of MAP inference are instances
More informationReview: Directed Models (Bayes Nets)
X Review: Directed Models (Bayes Nets) Lecture 3: Undirected Graphical Models Sam Roweis January 2, 24 Semantics: x y z if z d-separates x and y d-separation: z d-separates x from y if along every undirected
More informationGauge Plots. Gauge Plots JAPANESE BEETLE DATA MAXIMUM LIKELIHOOD FOR SPATIALLY CORRELATED DISCRETE DATA JAPANESE BEETLE DATA
JAPANESE BEETLE DATA 6 MAXIMUM LIKELIHOOD FOR SPATIALLY CORRELATED DISCRETE DATA Gauge Plots TuscaroraLisa Central Madsen Fairways, 996 January 9, 7 Grubs Adult Activity Grub Counts 6 8 Organic Matter
More informationSimultaneous Multi-frame MAP Super-Resolution Video Enhancement using Spatio-temporal Priors
Simultaneous Multi-frame MAP Super-Resolution Video Enhancement using Spatio-temporal Priors Sean Borman and Robert L. Stevenson Department of Electrical Engineering, University of Notre Dame Notre Dame,
More informationApproaches for Multiple Disease Mapping: MCAR and SANOVA
Approaches for Multiple Disease Mapping: MCAR and SANOVA Dipankar Bandyopadhyay Division of Biostatistics, University of Minnesota SPH April 22, 2015 1 Adapted from Sudipto Banerjee s notes SANOVA vs MCAR
More informationProbabilistic Graphical Models
2016 Robert Nowak Probabilistic Graphical Models 1 Introduction We have focused mainly on linear models for signals, in particular the subspace model x = Uθ, where U is a n k matrix and θ R k is a vector
More informationSpatial and Environmental Statistics
Spatial and Environmental Statistics Dale Zimmerman Department of Statistics and Actuarial Science University of Iowa January 17, 2019 Dale Zimmerman (UIOWA) Spatial and Environmental Statistics January
More informationBayesian Image Segmentation Using MRF s Combined with Hierarchical Prior Models
Bayesian Image Segmentation Using MRF s Combined with Hierarchical Prior Models Kohta Aoki 1 and Hiroshi Nagahashi 2 1 Interdisciplinary Graduate School of Science and Engineering, Tokyo Institute of Technology
More informationSTAT 518 Intro Student Presentation
STAT 518 Intro Student Presentation Wen Wei Loh April 11, 2013 Title of paper Radford M. Neal [1999] Bayesian Statistics, 6: 475-501, 1999 What the paper is about Regression and Classification Flexible
More informationPart 8: GLMs and Hierarchical LMs and GLMs
Part 8: GLMs and Hierarchical LMs and GLMs 1 Example: Song sparrow reproductive success Arcese et al., (1992) provide data on a sample from a population of 52 female song sparrows studied over the course
More informationEXAM IN STATISTICAL MACHINE LEARNING STATISTISK MASKININLÄRNING
EXAM IN STATISTICAL MACHINE LEARNING STATISTISK MASKININLÄRNING DATE AND TIME: August 30, 2018, 14.00 19.00 RESPONSIBLE TEACHER: Niklas Wahlström NUMBER OF PROBLEMS: 5 AIDING MATERIAL: Calculator, mathematical
More informationIntroduction. Semivariogram Cloud
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
More informationEconometría 2: Análisis de series de Tiempo
Econometría 2: Análisis de series de Tiempo Karoll GOMEZ kgomezp@unal.edu.co http://karollgomez.wordpress.com Segundo semestre 2016 II. Basic definitions A time series is a set of observations X t, each
More informationBayesian spatial hierarchical modeling for temperature extremes
Bayesian spatial hierarchical modeling for temperature extremes Indriati Bisono Dr. Andrew Robinson Dr. Aloke Phatak Mathematics and Statistics Department The University of Melbourne Maths, Informatics
More informationNearest Neighbor Gaussian Processes for Large Spatial Data
Nearest Neighbor Gaussian Processes for Large Spatial Data Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, 2017 1 Department of Biostatistics, Bloomberg School of Public Health, Johns
More informationMCMC Sampling for Bayesian Inference using L1-type Priors
MÜNSTER MCMC Sampling for Bayesian Inference using L1-type Priors (what I do whenever the ill-posedness of EEG/MEG is just not frustrating enough!) AG Imaging Seminar Felix Lucka 26.06.2012 , MÜNSTER Sampling
More informationBayesian inference for general Gaussian graphical models with application to multivariate lattice data
Bayesian inference for general Gaussian graphical models with application to multivariate lattice data Adrian Dobra, Alex Lenkoski and Abel odriguez arxiv:1005.4094v1 [stat.me] 21 May 2010 Abstract We
More informationA spatio-temporal model for extreme precipitation simulated by a climate model
A spatio-temporal model for extreme precipitation simulated by a climate model Jonathan Jalbert Postdoctoral fellow at McGill University, Montréal Anne-Catherine Favre, Claude Bélisle and Jean-François
More informationUndirected Graphical Models: Markov Random Fields
Undirected Graphical Models: Markov Random Fields 40-956 Advanced Topics in AI: Probabilistic Graphical Models Sharif University of Technology Soleymani Spring 2015 Markov Random Field Structure: undirected
More informationSpatial Statistics with Image Analysis. Outline. A Statistical Approach. Johan Lindström 1. Lund October 6, 2016
Spatial Statistics Spatial Examples More Spatial Statistics with Image Analysis Johan Lindström 1 1 Mathematical Statistics Centre for Mathematical Sciences Lund University Lund October 6, 2016 Johan Lindström
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 informationSpatial autocorrelation: robustness of measures and tests
Spatial autocorrelation: robustness of measures and tests Marie Ernst and Gentiane Haesbroeck University of Liege London, December 14, 2015 Spatial Data Spatial data : geographical positions non spatial
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science Algorithms For Inference Fall 2014
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.438 Algorithms For Inference Fall 2014 Recitation 3 1 Gaussian Graphical Models: Schur s Complement Consider
More informationLayered Synthesis of Latent Gaussian Trees
Layered Synthesis of Latent Gaussian Trees Ali Moharrer, Shuangqing Wei, George T. Amariucai, and Jing Deng arxiv:1608.04484v2 [cs.it] 7 May 2017 Abstract A new synthesis scheme is proposed to generate
More informationStat260: Bayesian Modeling and Inference Lecture Date: February 10th, Jeffreys priors. exp 1 ) p 2
Stat260: Bayesian Modeling and Inference Lecture Date: February 10th, 2010 Jeffreys priors Lecturer: Michael I. Jordan Scribe: Timothy Hunter 1 Priors for the multivariate Gaussian Consider a multivariate
More informationChris Bishop s PRML Ch. 8: Graphical Models
Chris Bishop s PRML Ch. 8: Graphical Models January 24, 2008 Introduction Visualize the structure of a probabilistic model Design and motivate new models Insights into the model s properties, in particular
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationFrailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Mela. P.
Frailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Melanie M. Wall, Bradley P. Carlin November 24, 2014 Outlines of the talk
More informationSpatial Smoothing in Stan: Conditional Auto-Regressive Models
Spatial Smoothing in Stan: Conditional Auto-Regressive Models Charles DiMaggio, PhD, NYU School of Medicine Stephen J. Mooney, PhD, University of Washington Mitzi Morris, Columbia University Dan Simpson,
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 informationBayesian Areal Wombling for Geographic Boundary Analysis
Bayesian Areal Wombling for Geographic Boundary Analysis Haolan Lu, Haijun Ma, and Bradley P. Carlin haolanl@biostat.umn.edu, haijunma@biostat.umn.edu, and brad@biostat.umn.edu Division of Biostatistics
More informationGaussian Processes 1. Schedule
1 Schedule 17 Jan: Gaussian processes (Jo Eidsvik) 24 Jan: Hands-on project on Gaussian processes (Team effort, work in groups) 31 Jan: Latent Gaussian models and INLA (Jo Eidsvik) 7 Feb: Hands-on project
More informationThe Nature of Geographic Data
4 The Nature of Geographic Data OVERVIEW Elaborates on the spatial is special theme Focuses on how phenomena vary across space and the general nature of geographic variation Describes the main principles
More informationDale L. Zimmerman Department of Statistics and Actuarial Science, University of Iowa, USA
SPATIAL STATISTICS Dale L. Zimmerman Department of Statistics and Actuarial Science, University of Iowa, USA Keywords: Geostatistics, Isotropy, Kriging, Lattice Data, Spatial point patterns, Stationarity
More informationSpatial inference. Spatial inference. Accounting for spatial correlation. Multivariate normal distributions
Spatial inference I will start with a simple model, using species diversity data Strong spatial dependence, Î = 0.79 what is the mean diversity? How precise is our estimate? Sampling discussion: The 64
More informationA fast sampler for data simulation from spatial, and other, Markov random fields
A fast sampler for data simulation from spatial, and other, Markov random fields Andee Kaplan Iowa State University ajkaplan@iastate.edu June 22, 2017 Slides available at http://bit.ly/kaplan-phd Joint
More informationA = {(x, u) : 0 u f(x)},
Draw x uniformly from the region {x : f(x) u }. Markov Chain Monte Carlo Lecture 5 Slice sampler: Suppose that one is interested in sampling from a density f(x), x X. Recall that sampling x f(x) is equivalent
More informationTemporal vs. Spatial Data
Temporal vs. Spatial Data Temporal 1 dimensional Units: day, week, month Lag: t, t-1, t-2 Durbin-Watson Spatial 2-3 dimensional Units: county, mile, region Lag: near neighbor, networks (?) Moran s I Differencing
More informationECE 592 Topics in Data Science
ECE 592 Topics in Data Science Final Fall 2017 December 11, 2017 Please remember to justify your answers carefully, and to staple your test sheet and answers together before submitting. Name: Student ID:
More informationPattern Theory and Its Applications
Pattern Theory and Its Applications Rudolf Kulhavý Institute of Information Theory and Automation Academy of Sciences of the Czech Republic, Prague Ulf Grenander A Swedish mathematician, since 1966 with
More informationMachine Learning, Fall 2009: Midterm
10-601 Machine Learning, Fall 009: Midterm Monday, November nd hours 1. Personal info: Name: Andrew account: E-mail address:. You are permitted two pages of notes and a calculator. Please turn off all
More informationCMPE 58K Bayesian Statistics and Machine Learning Lecture 5
CMPE 58K Bayesian Statistics and Machine Learning Lecture 5 Multivariate distributions: Gaussian, Bernoulli, Probability tables Department of Computer Engineering, Boğaziçi University, Istanbul, Turkey
More informationImage segmentation combining Markov Random Fields and Dirichlet Processes
Image segmentation combining Markov Random Fields and Dirichlet Processes Jessica SODJO IMS, Groupe Signal Image, Talence Encadrants : A. Giremus, J.-F. Giovannelli, F. Caron, N. Dobigeon Jessica SODJO
More informationPermutation-invariant regularization of large covariance matrices. Liza Levina
Liza Levina Permutation-invariant covariance regularization 1/42 Permutation-invariant regularization of large covariance matrices Liza Levina Department of Statistics University of Michigan Joint work
More informationProbabilistic Graphical Models
Probabilistic Graphical Models David Sontag New York University Lecture 4, February 16, 2012 David Sontag (NYU) Graphical Models Lecture 4, February 16, 2012 1 / 27 Undirected graphical models Reminder
More informationA Bayesian Perspective on Residential Demand Response Using Smart Meter Data
A Bayesian Perspective on Residential Demand Response Using Smart Meter Data Datong-Paul Zhou, Maximilian Balandat, and Claire Tomlin University of California, Berkeley [datong.zhou, balandat, tomlin]@eecs.berkeley.edu
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 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 informationHierarchical Modeling for Spatio-temporal Data
Hierarchical Modeling for Spatio-temporal Data 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
More informationAn EM algorithm for Gaussian Markov Random Fields
An EM algorithm for Gaussian Markov Random Fields Will Penny, Wellcome Department of Imaging Neuroscience, University College, London WC1N 3BG. wpenny@fil.ion.ucl.ac.uk October 28, 2002 Abstract Lavine
More informationHidden Markov Models for precipitation
Hidden Markov Models for precipitation Pierre Ailliot Université de Brest Joint work with Peter Thomson Statistics Research Associates (NZ) Page 1 Context Part of the project Climate-related risks for
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 informationChapter 4 - Fundamentals of spatial processes Lecture notes
Chapter 4 - Fundamentals of spatial processes Lecture notes Geir Storvik January 21, 2013 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites Mostly positive correlation Negative
More informationSpatial Misalignment
Spatial Misalignment Jamie Monogan University of Georgia Spring 2013 Jamie Monogan (UGA) Spatial Misalignment Spring 2013 1 / 28 Objectives By the end of today s meeting, participants should be able to:
More informationProbabilistic Graphical Models Lecture Notes Fall 2009
Probabilistic Graphical Models Lecture Notes Fall 2009 October 28, 2009 Byoung-Tak Zhang School of omputer Science and Engineering & ognitive Science, Brain Science, and Bioinformatics Seoul National University
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 informationModelling geoadditive survival data
Modelling geoadditive survival data Thomas Kneib & Ludwig Fahrmeir Department of Statistics, Ludwig-Maximilians-University Munich 1. Leukemia survival data 2. Structured hazard regression 3. Mixed model
More informationMultinomial Data. f(y θ) θ y i. where θ i is the probability that a given trial results in category i, i = 1,..., k. The parameter space is
Multinomial Data The multinomial distribution is a generalization of the binomial for the situation in which each trial results in one and only one of several categories, as opposed to just two, as in
More informationGIST 4302/5302: Spatial Analysis and Modeling
GIST 4302/5302: Spatial Analysis and Modeling Review Guofeng Cao www.gis.ttu.edu/starlab Department of Geosciences Texas Tech University guofeng.cao@ttu.edu Spring 2016 Course Outlines Spatial Point Pattern
More informationPart 7: Hierarchical Modeling
Part 7: Hierarchical Modeling!1 Nested data It is common for data to be nested: i.e., observations on subjects are organized by a hierarchy Such data are often called hierarchical or multilevel For example,
More informationResearch Division Federal Reserve Bank of St. Louis Working Paper Series
Research Division Federal Reserve Bank of St Louis Working Paper Series Kalman Filtering with Truncated Normal State Variables for Bayesian Estimation of Macroeconomic Models Michael Dueker Working Paper
More informationGraphical Models and Independence Models
Graphical Models and Independence Models Yunshu Liu ASPITRG Research Group 2014-03-04 References: [1]. Steffen Lauritzen, Graphical Models, Oxford University Press, 1996 [2]. Christopher M. Bishop, Pattern
More informationIntroduction to Spatial Statistics and Modeling for Regional Analysis
Introduction to Spatial Statistics and Modeling for Regional Analysis Dr. Xinyue Ye, Assistant Professor Center for Regional Development (Department of Commerce EDA University Center) & School of Earth,
More informationOn prediction and density estimation Peter McCullagh University of Chicago December 2004
On prediction and density estimation Peter McCullagh University of Chicago December 2004 Summary Having observed the initial segment of a random sequence, subsequent values may be predicted by calculating
More informationGaussian predictive process models for large spatial data sets.
Gaussian predictive process models for large spatial data sets. Sudipto Banerjee, Alan E. Gelfand, Andrew O. Finley, and Huiyan Sang Presenters: Halley Brantley and Chris Krut September 28, 2015 Overview
More informationSwitching Regime Estimation
Switching Regime Estimation Series de Tiempo BIrkbeck March 2013 Martin Sola (FE) Markov Switching models 01/13 1 / 52 The economy (the time series) often behaves very different in periods such as booms
More informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
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 informationNonresponse weighting adjustment using estimated response probability
Nonresponse weighting adjustment using estimated response probability Jae-kwang Kim Yonsei University, Seoul, Korea December 26, 2006 Introduction Nonresponse Unit nonresponse Item nonresponse Basic strategy
More informationMarkov Chain Monte Carlo
Markov Chain Monte Carlo Recall: To compute the expectation E ( h(y ) ) we use the approximation E(h(Y )) 1 n n h(y ) t=1 with Y (1),..., Y (n) h(y). Thus our aim is to sample Y (1),..., Y (n) from f(y).
More information