Spatial analysis of tropical rain forest plot data
|
|
- Sheila Candice Jackson
- 5 years ago
- Views:
Transcription
1 Spatial analysis of tropical rain forest plot data Rasmus Waagepetersen Department of Mathematical Sciences Aalborg University December 11, /45
2 Tropical rain forest ecology Fundamental questions: which factors influence the spatial distribution of rain forest trees and what is the reason for the high biodiversity of rain forests? Key factors: environment: topography, soil composition,... seed dispersal limitation: by wind, birds or mammals... 2/45
3 Tropical rain forest ecology Fundamental questions: which factors influence the spatial distribution of rain forest trees and what is the reason for the high biodiversity of rain forests? Key factors: environment: topography, soil composition,... seed dispersal limitation: by wind, birds or mammals... Outline: data examples introduction to spatial point processes applications to tropical rain forest data 3/45
4 Example: Capparis Frondosa and environment observation window W = 1000 m 500 m seed dispersal clustering environment inhomogeneity Elevation Potassium content in soil. Quantify dependence on environmental variables and seed dispersal using statistics for spatial point processes. 4/45
5 Example: modes of seed dispersal and clustering Three species with different modes of seed dispersal: Acalypha Diversifolia explosive capsules Loncocharpus Heptaphyllus wind Capparis Frondosa bird/mammal Is degree of clustering related to mode of seed dispersal? 5/45
6 Spatial point process Spatial point process: random collection of points (finite number of points in bounded sets) 6/45
7 Intensity function and product density X: spatial point process. A and B small subregions. B A 7/45
8 Intensity function and product density X: spatial point process. A and B small subregions. Intensity function of point process X ρ(u) A P(X has a point in A), u A A B 8/45
9 Intensity function and product density X: spatial point process. A and B small subregions. Intensity function of point process X ρ(u) A P(X has a point in A), u A A B Second order product density ρ (2) (u,v) A B P(X has a point in each of A and B) u A, v B 9/45
10 Pair correlation and K-function Pair correlation function g(u,v) = ρ(2) (u,v) ρ(u)ρ(v) NB: independent points ρ (2) (u,v) = ρ(u)ρ(v) g(u,v) = 1 10/45
11 Pair correlation and K-function Pair correlation function g(u,v) = ρ(2) (u,v) ρ(u)ρ(v) NB: independent points ρ (2) (u,v) = ρ(u)ρ(v) g(u,v) = 1 K-function K(t) = h t g(h)dh (provided g(u,v) = g(u v) i.e. X second-order reweighted stationary, Baddeley, Møller, Waagepetersen, 2000) 11/45
12 Pair correlation and K-function Pair correlation function g(u,v) = ρ(2) (u,v) ρ(u)ρ(v) NB: independent points ρ (2) (u,v) = ρ(u)ρ(v) g(u,v) = 1 K-function K(t) = h t g(h)dh (provided g(u,v) = g(u v) i.e. X second-order reweighted stationary, Baddeley, Møller, Waagepetersen, 2000) Examples of pair correlation and K-functions: g(t) K(t) t t 12/45
13 Mean and covariances of counts A and B subsets of the plane. N(A) and N(B) random numbers/counts of points in A and B. E[N(A)] = µ(a) = A ρ(u)du A B Cov[N(A),N(B)] = A B ρ(u)du+ ρ(u)ρ(v)[g(u, v) 1]dudv A B NB: can compute means and covariances for any sets A and B! (in contrast to quadrat count methods) 13/45
14 The Poisson process X is a Poisson process with intensity function ρ( ) if for any bounded region B: 1. N(B) is Poisson distributed with mean µ(b) = B ρ(u)du 2. Given N(B) = n, the n points are independent and identically distributed with density proportional to intensity function ρ( ). Homogeneous: ρ = 150/0.7 Inhomogeneous: ρ(x,y) e 10.6y 14/45
15 Back to rain forest: parametric models for intensity and pair correlation Study influence of covariates Z(u) = (Z 1 (u),...,z p (u)) using log-linear model for intensity function: logρ(u;β) = βz(u) T ρ(u;β) = exp(βz(u) T ) where βz(u) T = β 1 Z 1 (u)+β 2 Z 2 (u)+...+β p Z p (u) 15/45
16 Capparis Frondosa and Poisson process? Fit model with covariates elevation and Potassium. Fitted intensity function ρ(u; ˆβ) = exp(ˆβ 0 +ˆβ 1 Elev(u)+ˆβ 2 K(u)) : Estimated K-function and K(t) = πt 2 -function for Poisson process: K(t) Not Poisson process - aggregation due to unobserved factors (e.g. seed dispersal) t 16/45
17 Cox processes X is a Cox process driven by the random intensity function Λ if, conditional on Λ = λ, X is a Poisson process with intensity function λ. Example: log Gaussian Cox process (Møller, Syversveen, W, 1998) logλ(u) = βz(u) T +Y(u) where {Y(u)} Gaussian random field. βz(u) T Y(u) βz(u) T +Y(u) /45
18 Intensity and pair correlation function for log Gaussian Cox processes Log linear intensity logρ(u;β) = µ+z(u)β T 18/45
19 Intensity and pair correlation function for log Gaussian Cox processes Log linear intensity Pair correlation function: logρ(u;β) = µ+z(u)β T g(u v;ψ) = exp[c(u v;σ 2,α)], ψ = (σ 2,α) where σ 2 variance of Gaussian field and c( ;α) covariance function. 19/45
20 Intensity and pair correlation function for log Gaussian Cox processes Log linear intensity Pair correlation function: logρ(u;β) = µ+z(u)β T g(u v;ψ) = exp[c(u v;σ 2,α)], ψ = (σ 2,α) where σ 2 variance of Gaussian field and c( ;α) covariance function. Examples: σ 2 exp( u v /α) and σ 2 exp( u v 2 /α) exp. covariance Gau. covariance t t 20/45
21 Cluster process: Inhomogeneous Thomas process (W, 2007) Parents stationary Poisson point process intensity κ Offspring distributed around mothers according to Gaussian density with standard deviation ω Inhomogeneity: offspring survive according to probability p(u) exp(z(u)β T ) depending on covariates (independent thinning) /45
22 Cluster process as Cox process Inhomogeneous Thomas is a Cox process with intensity function Λ(u) = Λ 0 (u)exp[z(u)β T ] where Λ 0 (u) = v C k(u v) and k Gaussian density with standard deviation ω. More generally, k can be any bivariate probability density function - but not all choices lead to explicit pair correlation function. 22/45
23 Intensity and pair correlation function for cluster-cox processes Log linear intensity logρ(u;β) = µ+z(u)β T 23/45
24 Intensity and pair correlation function for cluster-cox processes Log linear intensity logρ(u;β) = µ+z(u)β T Pair correlation function for inhomogeneous Thomas: g(u v;ψ) = 1+exp( u v 2 /(4ω) 2 )/(4ω 2 κπ) = 1+σ 2 exp( u v 2 /α), ψ = (κ,ω) or ψ = (σ 2,α) 24/45
25 Intensity and pair correlation function for cluster-cox processes Log linear intensity logρ(u;β) = µ+z(u)β T Pair correlation function for inhomogeneous Thomas: g(u v;ψ) = 1+exp( u v 2 /(4ω) 2 )/(4ω 2 κπ) = 1+σ 2 exp( u v 2 /α), ψ = (κ,ω) or ψ = (σ 2,α) Pair correlation function for more general Bessel cluster processes: g(u v;ψ) = 1+σ 2( u v /α)ν K ν ( u v /α) 2 ν 1 Γ(ν) Special case ν = 1/2: g(u v;ψ) = 1+σ 2 exp( u v /α) 25/45
26 Parameter estimation Possibilities: 1. Maximum likelihood estimation (Monte Carlo computation of likelihood function) 2. Simple estimating functions based on intensity function and pair correlation function - inspired by methods for count variables: least squares, composite likelihood, quasi-likelihood,... 26/45
27 Example: composite likelihood I (Schoenberg, 2005; W, 2007) Consider indicators X i = 1[N i > 0] for presence of points in cells C i. P(X i = 1) = ρ β (u i ) C i. 27/45
28 Example: composite likelihood I (Schoenberg, 2005; W, 2007) Consider indicators X i = 1[N i > 0] for presence of points in cells C i. P(X i = 1) = ρ β (u i ) C i. Composite Bernouilli likelihood n (P(X i = 1)) X i (1 P(X i = 1)) 1 X i i=1 has limit ( C i 0) L(β) = [ u X W n ρ β (u i ) X i (1 ρ β (u i ) C i ) 1 X i i=1 ρ(u; β)] exp( ρ(u;β)du) W 28/45
29 Example: composite likelihood I (Schoenberg, 2005; W, 2007) Consider indicators X i = 1[N i > 0] for presence of points in cells C i. P(X i = 1) = ρ β (u i ) C i. Composite Bernouilli likelihood n (P(X i = 1)) X i (1 P(X i = 1)) 1 X i i=1 has limit ( C i 0) L(β) = [ u X W Estimate ˆβ maximizes L(β). n ρ β (u i ) X i (1 ρ β (u i ) C i ) 1 X i i=1 ρ(u; β)] exp( ρ(u;β)du) W NB: L(β) formally equivalent to likelihood function of a Poisson process with intensity function ρ β ( ). 29/45
30 Example: minimum contrast estimation for ψ Computationally easy approach if X second-order reweighted stationary so that K-function well-defined. Estimate of K-function (Baddeley, Møller and W, 2000): ˆK β (t) = 1[0 < u v t] e u,v ρ(u;β)ρ(v;β) u,v X W Unbiased if β true regression parameter. Minimum contrast estimation: minimize squared distance between theoretical K and ˆK: ˆψ = argmin ψ r 0 (ˆKˆβ(t) K(t;ψ) ) 2 dt K(t) K estim. Inhom. Thomas t 30/45
31 Two-step estimation Obtain estimates (ˆβ, ˆψ) in two steps 1. obtain ˆβ using composite likelihood 2. obtain ˆψ using minimum contrast 31/45
32 Clustering and mode of seed dispersal Fit Thomas cluster process with log linear model for intensity function. Acalypha and Capparis: positive dependence on elevation and potassium (significantly positive coefficients ˆβ = (0.02, 0.005) and ˆβ = (0.03,0.004)). Loncocharpus: negative dependence on nitrogen and phosphorous (ˆβ = ( 0.03, 0.16)). 32/45
33 Recall ω = width of clusters. Estimates of ω for explosive, wind and bird/mammal: Estimates of K-functions for bird/mammal dispersed species K(t) K inhom. Thomas K hom. K Poisson t Triangles: model without covariates. 33/45
34 Decomposition of variance for rain forest tree point patterns (Shen, Jalilian, W, in progress) Question: how much of the spatial variation for rain forest trees is due to environment? 34/45
35 Decomposition of variance for rain forest tree point patterns (Shen, Jalilian, W, in progress) Question: how much of the spatial variation for rain forest trees is due to environment? Variance of a count N(B) (number of points in region B) for a stationary Cox process (constant intensity ρ): 35/45
36 Decomposition of variance for rain forest tree point patterns (Shen, Jalilian, W, in progress) Question: how much of the spatial variation for rain forest trees is due to environment? Variance of a count N(B) (number of points in region B) for a stationary Cox process (constant intensity ρ): VarN(B) = ρdu Variance=Poisson variance B 36/45
37 Decomposition of variance for rain forest tree point patterns (Shen, Jalilian, W, in progress) Question: how much of the spatial variation for rain forest trees is due to environment? Variance of a count N(B) (number of points in region B) for a stationary Cox process (constant intensity ρ): VarN(B) = ρdu + ρ 2 [g(u,v) 1]dudv B B Variance=Poisson variance+extra variance due to random intensity B 37/45
38 Decomposition of variance for log linear random intensity: 38/45
39 Decomposition of variance for log linear random intensity: VarlogΛ(u) = VarβZ(u) T Variance=Environment 39/45
40 Decomposition of variance for log linear random intensity: VarlogΛ(u) = VarβZ(u) T +VarY(u) = σ 2 Z +σ2 Y Variance=Environment+Seed dispersal 40/45
41 Decomposition of variance for log linear random intensity: VarlogΛ(u) = VarβZ(u) T +VarY(u) = σ 2 Z +σ2 Y Variance=Environment+Seed dispersal Note Z(u) = βz(u) T regarded as stationary random process. 41/45
42 Decomposition of variance for log linear random intensity: VarlogΛ(u) = VarβZ(u) T +VarY(u) = σz 2 +σ2 Y Variance=Environment+Seed dispersal Note Z(u) = βz(u) T regarded as stationary random process. Estimate β and σy 2 using two-step approach. Simple empirical estimate of σz 2: ˆσ Z 2 = 1 ( Z(u) Z) 2 n G u G Compute σ 2 Z σ 2 Z +σ2 Y and σ 2 Y σ 2 Z +σ2 Y 42/45
43 Decomposition of variance for log linear random intensity: VarlogΛ(u) = VarβZ(u) T +VarY(u) = σz 2 +σ2 Y Variance=Environment+Seed dispersal Note Z(u) = βz(u) T regarded as stationary random process. Estimate β and σy 2 using two-step approach. Simple empirical estimate of σz 2: ˆσ Z 2 = 1 ( Z(u) Z) 2 n G u G Compute σ 2 Z σ 2 Z +σ2 Y and σ 2 Y σ 2 Z +σ2 Y Can also define closely related R 2 summarizing how much of variation in Λ is due to Z. 43/45
44 Additive model for random intensity function (Jalilian, Guan, W, in progress) Alternative to log linear model: Λ(u) = βz(u) T +Y(u) Cox process superposition of point processes with (random) intensity functions Z(u) = βz(u) T and Y(u) Straightforward variance decomposition for Λ: VarΛ(u) = Var Z(u)+VarY(u) = σ 2 Z +σ2 Y R 2 = σ 2 Z σ 2 Z +σ2 Y 44/45
45 Results Consider pair correlation functions of the form g(u v;σ 2,α) = 1+σ 2 exp( u v δ /α) δ = 1 or δ = 2 Species Λ δ R 2 Goodness of fit ( AIC ) Acalypha log linear log linear additive additive Lonchocarpus log linear log linear additive additive Capparis log linear log linear additive additive Best fit with log linear model and δ = 1. Largest R 2 for bird/mammal dispersion. Smallest for explosive capsules. 45/45
46 Fitted pair correlation functions Plots show g(u v) 1: c(r) Exp Gau LG Exp LG Gau c(r) Exp Gau LG Exp LG Gau c(r) Exp Gau LG Exp LG Gau r r r 46/45
47 Final remarks Log-linear model gives better fit than additive Better fit with exponential covariance function than with Gaussian Gaussian covariance function (Thomas process) tails decay too fast value of R 2 related to mode of seed dispersal? 47/45
48 Thank you for your attention! 48/45
Decomposition of variance for spatial Cox processes
Decomposition of variance for spatial Cox processes Rasmus Waagepetersen Department of Mathematical Sciences Aalborg University Joint work with Abdollah Jalilian and Yongtao Guan November 8, 2010 1/34
More informationDecomposition of variance for spatial Cox processes
Decomposition of variance for spatial Cox processes Rasmus Waagepetersen Department of Mathematical Sciences Aalborg University Joint work with Abdollah Jalilian and Yongtao Guan December 13, 2010 1/25
More informationEstimating functions for inhomogeneous spatial point processes with incomplete covariate data
Estimating functions for inhomogeneous spatial point processes with incomplete covariate data Rasmus aagepetersen Department of Mathematics Aalborg University Denmark August 15, 2007 1 / 23 Data (Barro
More informationESTIMATING FUNCTIONS FOR INHOMOGENEOUS COX PROCESSES
ESTIMATING FUNCTIONS FOR INHOMOGENEOUS COX PROCESSES Rasmus Waagepetersen Department of Mathematics, Aalborg University, Fredrik Bajersvej 7G, DK-9220 Aalborg, Denmark (rw@math.aau.dk) Abstract. Estimation
More informationDecomposition of variance for spatial Cox processes Jalilian, Abdollah; Guan, Yongtao; Waagepetersen, Rasmus Plenge
Aalborg Universitet Decomposition of variance for spatial Cox processes Jalilian, Abdollah; Guan, Yongtao; Waagepetersen, Rasmus Plenge Publication date: 2011 Document Version Early version, also known
More informationQuasi-likelihood for Spatial Point Processes
Quasi-likelihood for Spatial Point Processes Yongtao Guan Miami, USA Abdollah Jalilian Kermanshah, Iran Rasmus aagepetersen Aalborg, Denmark Summary. Fitting regression models for intensity functions of
More informationLecture 2: Poisson point processes: properties and statistical inference
Lecture 2: Poisson point processes: properties and statistical inference Jean-François Coeurjolly http://www-ljk.imag.fr/membres/jean-francois.coeurjolly/ 1 / 20 Definition, properties and simulation Statistical
More informationTwo step estimation for Neyman-Scott point process with inhomogeneous cluster centers. May 2012
Two step estimation for Neyman-Scott point process with inhomogeneous cluster centers Tomáš Mrkvička, Milan Muška, Jan Kubečka May 2012 Motivation Study of the influence of covariates on the occurrence
More informationChapter 2. Poisson point processes
Chapter 2. Poisson point processes Jean-François Coeurjolly http://www-ljk.imag.fr/membres/jean-francois.coeurjolly/ Laboratoire Jean Kuntzmann (LJK), Grenoble University Setting for this chapter To ease
More informationEstimating functions for inhomogeneous spatial point processes with incomplete covariate data
Estimating functions for inhomogeneous spatial point processes with incomplete covariate data Rasmus aagepetersen Department of Mathematical Sciences, Aalborg University Fredrik Bajersvej 7G, DK-9220 Aalborg
More informationOn Model Fitting Procedures for Inhomogeneous Neyman-Scott Processes
On Model Fitting Procedures for Inhomogeneous Neyman-Scott Processes Yongtao Guan July 31, 2006 ABSTRACT In this paper we study computationally efficient procedures to estimate the second-order parameters
More informationAALBORG UNIVERSITY. An estimating function approach to inference for inhomogeneous Neyman-Scott processes. Rasmus Plenge Waagepetersen
AALBORG UNIVERITY An estimating function approach to inference for inhomogeneous Neyman-cott processes by Rasmus Plenge Waagepetersen R-2005-30 eptember 2005 Department of Mathematical ciences Aalborg
More informationSecond-Order Analysis of Spatial Point Processes
Title Second-Order Analysis of Spatial Point Process Tonglin Zhang Outline Outline Spatial Point Processes Intensity Functions Mean and Variance Pair Correlation Functions Stationarity K-functions Some
More informationPoint Processes. Tonglin Zhang. Spatial Statistics for Point and Lattice Data (Part II)
Title: patial tatistics for Point Processes and Lattice Data (Part II) Point Processes Tonglin Zhang Outline Outline imulated Examples Interesting Problems Analysis under tationarity Analysis under Nonstationarity
More informationRESEARCH REPORT. Optimal Estimation of the Intensity Function of a Spatial Point Process.
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING RESEARCH REPORT www.csgb.dk 2011 Yongtao Guan, Abdollah Jalilian and Rasmus aagepetersen Optimal Estimation of the Intensity Function of a Spatial
More informationA Thinned Block Bootstrap Variance Estimation. Procedure for Inhomogeneous Spatial Point Patterns
A Thinned Block Bootstrap Variance Estimation Procedure for Inhomogeneous Spatial Point Patterns May 22, 2007 Abstract When modeling inhomogeneous spatial point patterns, it is of interest to fit a parametric
More informationAn introduction to spatial point processes
An introduction to spatial point processes Jean-François Coeurjolly 1 Examples of spatial data 2 Intensities and Poisson p.p. 3 Summary statistics 4 Models for point processes A very very brief introduction...
More informationProperties of spatial Cox process models
Properties of spatial Cox process models Jesper Møller Department of Mathematical Sciences, Aalborg University Abstract: Probabilistic properties of Cox processes of relevance for statistical modelling
More informationRESEARCH REPORT. A note on gaps in proofs of central limit theorems. Christophe A.N. Biscio, Arnaud Poinas and Rasmus Waagepetersen
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING 2017 www.csgb.dk RESEARCH REPORT Christophe A.N. Biscio, Arnaud Poinas and Rasmus Waagepetersen A note on gaps in proofs of central limit theorems
More informationSecond-order pseudo-stationary random fields and point processes on graphs and their edges
Second-order pseudo-stationary random fields and point processes on graphs and their edges Jesper Møller (in collaboration with Ethan Anderes and Jakob G. Rasmussen) Aalborg University Jesper Møller (Aalborg
More informationSummary statistics for inhomogeneous spatio-temporal marked point patterns
Summary statistics for inhomogeneous spatio-temporal marked point patterns Marie-Colette van Lieshout CWI Amsterdam The Netherlands Joint work with Ottmar Cronie Summary statistics for inhomogeneous spatio-temporal
More informationVariance Estimation for Statistics Computed from. Inhomogeneous Spatial Point Processes
Variance Estimation for Statistics Computed from Inhomogeneous Spatial Point Processes Yongtao Guan April 14, 2007 Abstract This paper introduces a new approach to estimate the variance of statistics that
More informationPoint process models for earthquakes with applications to Groningen and Kashmir data
Point process models for earthquakes with applications to Groningen and Kashmir data Marie-Colette van Lieshout colette@cwi.nl CWI & Twente The Netherlands Point process models for earthquakes with applications
More informationModern statistics for spatial point processes
Modern statistics for spatial point processes June 21, 2006 Jesper Møller and Rasmus P Waagepetersen Department of Mathematical Sciences, Aalborg University Abstract: We summarize and discuss the current
More informationPrediction. is a weighted least squares estimate since it minimizes. Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark
Prediction Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark March 22, 2017 WLS and BLUE (prelude to BLUP) Suppose that Y has mean β and known covariance matrix V (but Y need not
More informationNegative Association, Ordering and Convergence of Resampling Methods
Negative Association, Ordering and Convergence of Resampling Methods Nicolas Chopin ENSAE, Paristech (Joint work with Mathieu Gerber and Nick Whiteley, University of Bristol) Resampling schemes: Informal
More informationRESEARCH REPORT. Log Gaussian Cox processes on the sphere. Francisco Cuevas-Pacheco and Jesper Møller. No.
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING 2018 www.csgb.dk RESEARCH REPORT Francisco Cuevas-Pacheco and Jesper Møller Log Gaussian Cox processes on the sphere No. 04, March 2018 Log Gaussian
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 informationStructured Spatio-Temporal Shot-Noise Cox Point Process Models, with a View to Modelling Forest Fires
Scandinavian Journal of Statistics, Vol. 37: 2 25, 2010 doi: 10.1111/j.1467-9469.2009.00670.x Published by Blackwell Publishing Ltd. Structured Spatio-Temporal Shot-Noise Cox Point Process Models, with
More informationLinear models. Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark. October 5, 2016
Linear models Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark October 5, 2016 1 / 16 Outline for today linear models least squares estimation orthogonal projections estimation
More informationRESEARCH REPORT. Multivariate Product-Shot-noise Cox Point Process Models.
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING 2015 www.csgb.dk RESEARCH REPORT Abdollah Jalilian, Yongtao Guan, Jorge Mateu and Rasmus Waagepetersen Multivariate Product-Shot-noise Cox Point Process
More informationDoubly Inhomogeneous Cluster Point Processes
Doubly Inhomogeneous Cluster Point Processes Tomáš Mrkvička, Samuel Soubeyrand May 2016 Abscissa Motivation - Group dispersal model It is dispersal model, where particles are released in groups by a single
More informationAalborg Universitet. Modern statistics for spatial point processes Møller, Jesper; Waagepetersen, Rasmus Plenge. Publication date: 2006
Aalborg Universitet Modern statistics for spatial point processes Møller, Jesper; Waagepetersen, Rasmus Plenge Publication date: 2006 Document Version Accepted author manuscript, peer reviewed version
More informationWLS and BLUE (prelude to BLUP) Prediction
WLS and BLUE (prelude to BLUP) Prediction Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark April 21, 2018 Suppose that Y has mean X β and known covariance matrix V (but Y need
More informationGeometric anisotropic spatial point pattern analysis and Cox processes
Downloaded from vbn.aau.dk on: January 29, 2019 Aalborg Universitet Geometric anisotropic spatial point pattern analysis and Cox processes Møller, Jesper; Toftaker, Håkon Publication date: 2012 Document
More informationREGULARIZED ESTIMATING EQUATIONS FOR MODEL SELECTION OF CLUSTERED SPATIAL POINT PROCESSES
Statistica Sinica 25 (2015), 173-188 doi:http://dx.doi.org/10.5705/ss.2013.208w REGULARIZED ESTIMATING EQUATIONS FOR MODEL SELECTION OF CLUSTERED SPATIAL POINT PROCESSES Andrew L. Thurman 1, Rao Fu 2,
More informationReparametrization of COM-Poisson Regression Models with Applications in the Analysis of Experimental Count Data
Reparametrization of COM-Poisson Regression Models with Applications in the Analysis of Experimental Count Data Eduardo Elias Ribeiro Junior 1 2 Walmes Marques Zeviani 1 Wagner Hugo Bonat 1 Clarice Garcia
More informationModel Selection and Geometry
Model Selection and Geometry Pascal Massart Université Paris-Sud, Orsay Leipzig, February Purpose of the talk! Concentration of measure plays a fundamental role in the theory of model selection! Model
More informationHandbook of Spatial Statistics. Jesper Møller
Handbook of Spatial Statistics Jesper Møller May 7, 2008 ii Contents 1 1 2 3 3 5 4 Parametric methods 7 4.1 Introduction.............................. 7 4.2 Setting and notation.........................
More informationCox s proportional hazards model and Cox s partial likelihood
Cox s proportional hazards model and Cox s partial likelihood Rasmus Waagepetersen October 12, 2018 1 / 27 Non-parametric vs. parametric Suppose we want to estimate unknown function, e.g. survival function.
More informationPh.D. Qualifying Exam Friday Saturday, January 6 7, 2017
Ph.D. Qualifying Exam Friday Saturday, January 6 7, 2017 Put your solution to each problem on a separate sheet of paper. Problem 1. (5106) Let X 1, X 2,, X n be a sequence of i.i.d. observations from a
More informationHierarchical Modeling for Univariate Spatial Data
Hierarchical Modeling for Univariate Spatial Data Geography 890, Hierarchical Bayesian Models for Environmental Spatial Data Analysis February 15, 2011 1 Spatial Domain 2 Geography 890 Spatial Domain This
More informationCox s proportional hazards/regression model - model assessment
Cox s proportional hazards/regression model - model assessment Rasmus Waagepetersen September 27, 2017 Topics: Plots based on estimated cumulative hazards Cox-Snell residuals: overall check of fit Martingale
More informationInference Methods for the Conditional Logistic Regression Model with Longitudinal Data Arising from Animal Habitat Selection Studies
Inference Methods for the Conditional Logistic Regression Model with Longitudinal Data Arising from Animal Habitat Selection Studies Thierry Duchesne 1 (Thierry.Duchesne@mat.ulaval.ca) with Radu Craiu,
More informationStochastic Quasi-likelihood for Case-Control Point Pattern Data
Stochastic Quasi-likelihood for Case-Control Point Pattern Data Ganggang Xu, Rasmus Waagepetersen and Yongtao Guan January 8, 2018 Abstract We propose a novel stochastic quasi-likelihood estimation procedure
More informationWrapped Gaussian processes: a short review and some new results
Wrapped Gaussian processes: a short review and some new results Giovanna Jona Lasinio 1, Gianluca Mastrantonio 2 and Alan Gelfand 3 1-Università Sapienza di Roma 2- Università RomaTRE 3- Duke University
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 informationREGRESSION WITH SPATIALLY MISALIGNED DATA. Lisa Madsen Oregon State University David Ruppert Cornell University
REGRESSION ITH SPATIALL MISALIGNED DATA Lisa Madsen Oregon State University David Ruppert Cornell University SPATIALL MISALIGNED DATA 10 X X X X X X X X 5 X X X X X 0 X 0 5 10 OUTLINE 1. Introduction 2.
More informationSpatial point processes
Mathematical sciences Chalmers University of Technology and University of Gothenburg Gothenburg, Sweden June 25, 2014 Definition A point process N is a stochastic mechanism or rule to produce point patterns
More informationRESEARCH REPORT. Inhomogeneous spatial point processes with hidden second-order stationarity. Ute Hahn and Eva B.
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING www.csgb.dk RESEARCH REPORT 2013 Ute Hahn and Eva B. Vedel Jensen Inhomogeneous spatial point processes with hidden second-order stationarity No.
More informationIntegrated Likelihood Estimation in Semiparametric Regression Models. Thomas A. Severini Department of Statistics Northwestern University
Integrated Likelihood Estimation in Semiparametric Regression Models Thomas A. Severini Department of Statistics Northwestern University Joint work with Heping He, University of York Introduction Let Y
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationTail negative dependence and its applications for aggregate loss modeling
Tail negative dependence and its applications for aggregate loss modeling Lei Hua Division of Statistics Oct 20, 2014, ISU L. Hua (NIU) 1/35 1 Motivation 2 Tail order Elliptical copula Extreme value copula
More informationLecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE. Rick Katz
1 Lecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home
More informationGeneralized Linear Models
Generalized Linear Models Advanced Methods for Data Analysis (36-402/36-608 Spring 2014 1 Generalized linear models 1.1 Introduction: two regressions So far we ve seen two canonical settings for regression.
More informationAdvanced Statistics I : Gaussian Linear Model (and beyond)
Advanced Statistics I : Gaussian Linear Model (and beyond) Aurélien Garivier CNRS / Telecom ParisTech Centrale Outline One and Two-Sample Statistics Linear Gaussian Model Model Reduction and model Selection
More informationBayesian Hierarchical Models
Bayesian Hierarchical Models Gavin Shaddick, Millie Green, Matthew Thomas University of Bath 6 th - 9 th December 2016 1/ 34 APPLICATIONS OF BAYESIAN HIERARCHICAL MODELS 2/ 34 OUTLINE Spatial epidemiology
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 informationESTIMATING THE MEAN LEVEL OF FINE PARTICULATE MATTER: AN APPLICATION OF SPATIAL STATISTICS
ESTIMATING THE MEAN LEVEL OF FINE PARTICULATE MATTER: AN APPLICATION OF SPATIAL STATISTICS Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, N.C.,
More informationGaussian with mean ( µ ) and standard deviation ( σ)
Slide from Pieter Abbeel Gaussian with mean ( µ ) and standard deviation ( σ) 10/6/16 CSE-571: Robotics X ~ N( µ, σ ) Y ~ N( aµ + b, a σ ) Y = ax + b + + + + 1 1 1 1 1 1 1 1 1 1, ~ ) ( ) ( ), ( ~ ), (
More informationModel Selection for Geostatistical Models
Model Selection for Geostatistical Models Richard A. Davis Colorado State University http://www.stat.colostate.edu/~rdavis/lectures Joint work with: Jennifer A. Hoeting, Colorado State University Andrew
More informationMA 575 Linear Models: Cedric E. Ginestet, Boston University Non-parametric Inference, Polynomial Regression Week 9, Lecture 2
MA 575 Linear Models: Cedric E. Ginestet, Boston University Non-parametric Inference, Polynomial Regression Week 9, Lecture 2 1 Bootstrapped Bias and CIs Given a multiple regression model with mean and
More informationA General Unified Niche-Assembly/Dispersal-Assembly Theory of Forest Species Biodiversity
A General Unified Niche-Assembly/Dispersal-Assembly Theory of Forest Species Biodiversity Keith Rennolls CMS, University of Greenwich, Park Row, London SE10 9LS k.rennolls@gre.ac.uk Abstract: A generalised
More informationAalborg Universitet. Parametric methods for spatial point processes Møller, Jesper. Publication date: 2008
Aalborg Universitet Parametric methods for spatial point processes Møller, Jesper Publication date: 2008 Document Version Publisher's PDF, also known as Version of record Link to publication from Aalborg
More informationOutline. Statistical inference for linear mixed models. One-way ANOVA in matrix-vector form
Outline Statistical inference for linear mixed models Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark general form of linear mixed models examples of analyses using linear mixed
More informationarxiv: v2 [math.st] 5 May 2018
Log Gaussian Cox processes on the sphere Francisco Cuevas-Pacheco a, Jesper Møller a, a Department of Mathematical Sciences, Aalborg University, Denmark arxiv:1803.03051v2 [math.st] 5 May 2018 Abstract
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 informationA Spatio-Temporal Point Process Model for Firemen Demand in Twente
University of Twente A Spatio-Temporal Point Process Model for Firemen Demand in Twente Bachelor Thesis Author: Mike Wendels Supervisor: prof. dr. M.N.M. van Lieshout Stochastic Operations Research Applied
More informationWeighted Least Squares
Weighted Least Squares The standard linear model assumes that Var(ε i ) = σ 2 for i = 1,..., n. As we have seen, however, there are instances where Var(Y X = x i ) = Var(ε i ) = σ2 w i. Here w 1,..., w
More informationEstimation and Model Selection in Mixed Effects Models Part I. Adeline Samson 1
Estimation and Model Selection in Mixed Effects Models Part I Adeline Samson 1 1 University Paris Descartes Summer school 2009 - Lipari, Italy These slides are based on Marc Lavielle s slides Outline 1
More informationInvestigating seed dispersal and natural. moment methods
Investigating seed dispersal and natural enemy attack with wavelet ltvariances and moment methods Helene C. Muller Landau and Matteo Detto Smithsonian Tropical Research Institute Heterogeneous spatial
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 informationHierarchical Modelling for Univariate Spatial Data
Hierarchical Modelling for Univariate Spatial 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
More informationOutline of GLMs. Definitions
Outline of GLMs Definitions This is a short outline of GLM details, adapted from the book Nonparametric Regression and Generalized Linear Models, by Green and Silverman. The responses Y i have density
More informationCovariance function estimation in Gaussian process regression
Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationMS&E 226: Small Data. Lecture 11: Maximum likelihood (v2) Ramesh Johari
MS&E 226: Small Data Lecture 11: Maximum likelihood (v2) Ramesh Johari ramesh.johari@stanford.edu 1 / 18 The likelihood function 2 / 18 Estimating the parameter This lecture develops the methodology behind
More informationParametric Modelling of Over-dispersed Count Data. Part III / MMath (Applied Statistics) 1
Parametric Modelling of Over-dispersed Count Data Part III / MMath (Applied Statistics) 1 Introduction Poisson regression is the de facto approach for handling count data What happens then when Poisson
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationRESEARCH REPORT. A logistic regression estimating function for spatial Gibbs point processes.
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING www.csgb.dk RESEARCH REPORT 2013 Adrian Baddeley, Jean-François Coeurjolly, Ege Rubak and Rasmus Waagepetersen A logistic regression estimating function
More informationLikelihood Inference in the Presence of Nuisance Parameters
PHYSTAT2003, SLAC, September 8-11, 2003 1 Likelihood Inference in the Presence of Nuance Parameters N. Reid, D.A.S. Fraser Department of Stattics, University of Toronto, Toronto Canada M5S 3G3 We describe
More informationStatistical Analysis of Spatio-temporal Point Process Data. Peter J Diggle
Statistical Analysis of Spatio-temporal Point Process Data Peter J Diggle Department of Medicine, Lancaster University and Department of Biostatistics, Johns Hopkins University School of Public Health
More informationThinning-stable point processes as a model for bursty spatial data
Thinning-stable point processes as a model for bursty spatial data Chalmers University of Technology, Gothenburg, Sweden Paris, Jan 14th 2015 Communications Science. XXth Century Fixed line telephony Scientific
More informationIntroduction to General and Generalized Linear Models
Introduction to General and Generalized Linear Models Mixed effects models - II Henrik Madsen, Jan Kloppenborg Møller, Anders Nielsen April 16, 2012 H. Madsen, JK. Møller, A. Nielsen () Chapman & Hall
More informationIE 303 Discrete-Event Simulation
IE 303 Discrete-Event Simulation 1 L E C T U R E 5 : P R O B A B I L I T Y R E V I E W Review of the Last Lecture Random Variables Probability Density (Mass) Functions Cumulative Density Function Discrete
More informationInformation geometry for bivariate distribution control
Information geometry for bivariate distribution control C.T.J.Dodson + Hong Wang Mathematics + Control Systems Centre, University of Manchester Institute of Science and Technology Optimal control of stochastic
More informationOutline for today. Computation of the likelihood function for GLMMs. Likelihood for generalized linear mixed model
Outline for today Computation of the likelihood function for GLMMs asmus Waagepetersen Department of Mathematics Aalborg University Denmark Computation of likelihood function. Gaussian quadrature 2. Monte
More informationSpatial cluster point processes related to Poisson-Voronoi tessellations Møller, Jesper; Rasmussen, Jakob Gulddahl
Aalborg Universitet Spatial cluster point processes related to Poisson-Voronoi tessellations Møller, Jesper; Rasmussen, Jakob Gulddahl Publication date: 2013 Document Version Early version, also known
More informationCSCI-567: Machine Learning (Spring 2019)
CSCI-567: Machine Learning (Spring 2019) Prof. Victor Adamchik U of Southern California Mar. 19, 2019 March 19, 2019 1 / 43 Administration March 19, 2019 2 / 43 Administration TA3 is due this week March
More informationFundamentals of Data Assimila1on
014 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO July 14-16, 014 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
More informationChapter 4: Factor Analysis
Chapter 4: Factor Analysis In many studies, we may not be able to measure directly the variables of interest. We can merely collect data on other variables which may be related to the variables of interest.
More informationIntroduction to Bayes and non-bayes spatial statistics
Introduction to Bayes and non-bayes spatial statistics Gabriel Huerta Department of Mathematics and Statistics University of New Mexico http://www.stat.unm.edu/ ghuerta/georcode.txt General Concepts Spatial
More informationLecture 10 Spatio-Temporal Point Processes
Lecture 10 Spatio-Temporal Point Processes Dennis Sun Stats 253 July 23, 2014 Outline of Lecture 1 Review of Last Lecture 2 Spatio-temporal Point Processes 3 The Spatio-temporal Poisson Process 4 Modeling
More information( t) Cox regression part 2. Outline: Recapitulation. Estimation of cumulative hazards and survival probabilites. Ørnulf Borgan
Outline: Cox regression part 2 Ørnulf Borgan Department of Mathematics University of Oslo Recapitulation Estimation of cumulative hazards and survival probabilites Assumptions for Cox regression and check
More informationGIST 4302/5302: Spatial Analysis and Modeling Point Pattern Analysis
GIST 4302/5302: Spatial Analysis and Modeling Point Pattern Analysis Guofeng Cao www.spatial.ttu.edu Department of Geosciences Texas Tech University guofeng.cao@ttu.edu Fall 2018 Spatial Point Patterns
More informationLong-Run Covariability
Long-Run Covariability Ulrich K. Müller and Mark W. Watson Princeton University October 2016 Motivation Study the long-run covariability/relationship between economic variables great ratios, long-run Phillips
More informationVariable Selection and Model Choice in Survival Models with Time-Varying Effects
Variable Selection and Model Choice in Survival Models with Time-Varying Effects Boosting Survival Models Benjamin Hofner 1 Department of Medical Informatics, Biometry and Epidemiology (IMBE) Friedrich-Alexander-Universität
More informationIntroduction to Statistical Analysis
Introduction to Statistical Analysis Changyu Shen Richard A. and Susan F. Smith Center for Outcomes Research in Cardiology Beth Israel Deaconess Medical Center Harvard Medical School Objectives Descriptive
More informationComparing Non-informative Priors for Estimation and Prediction in Spatial Models
Environmentrics 00, 1 12 DOI: 10.1002/env.XXXX Comparing Non-informative Priors for Estimation and Prediction in Spatial Models Regina Wu a and Cari G. Kaufman a Summary: Fitting a Bayesian model to spatial
More information