Estimating functions for inhomogeneous spatial point processes with incomplete covariate data
|
|
- Megan Cain
- 5 years ago
- Views:
Transcription
1 Estimating functions for inhomogeneous spatial point processes with incomplete covariate data Rasmus aagepetersen Department of Mathematics Aalborg University Denmark August 15, / 23
2 Data (Barro Colorado Island Forest Dynamics Plot) Observation window: S = [0, 1000] [0, 500]m 2 Beilschmiedia Ocotea Elevation Gradient norm (steepness) Question: tree intensities related to elevation and gradient? 2 / 23
3 Log-linear models for intensity function z(u) = ( z 1 (u),...,z p (u) ) vector of covariates for each location u in observation window. E.g. z(u) = (1, z elev (u), z grad (u)) for rain forest example. Log-linear model for intensity function: Interpretation: λ(u; β) = exp(z(u)β T ) λ(u; β)dn u = P(point occurs in neighbourhood N u around u) 3 / 23
4 Estimation of β X: point pattern observed in Poisson process log likelihood function: l(β) = z(u)β T u X Poisson score estimating function: u(β) = z(u) u X λ(u; β)du z(u)λ(u; β)du Also applicable for non-poisson point processes with intensity function λ( ; β) (Schoenberg, 2005, aagepetersen, 2007) 4 / 23
5 Missing covariate data Elevation covariate interpolated from elevation observations on grid. However, evaluating u(β) = z(u) z(u)λ(u; β)du u X requires z(u) observed for any u! 5 / 23
6 Approximations of log likelihood I Suppose z(u) observed at finite set of locations Q. Rathbun (1996) approximate z(u)λ(u; β)du z(u)λ(u; β)du where z(u)λ(u; β) unbiased prediction of z(u)λ(u; β), u given z(u), u Q. 6 / 23
7 Approximations of log likelihood I Suppose z(u) observed at finite set of locations Q. Rathbun (1996) approximate z(u)λ(u; β)du z(u)λ(u; β)du where z(u)λ(u; β) unbiased prediction of z(u)λ(u; β), u given z(u), u Q. Riemann approximation: z(u)λ(u; β)du w(u)z(u)λ(u; β) u Q where w(u) quadrature weight for u Q. 7 / 23
8 Approximations of log likelihood II: spatstat Approximation of score function used in R package spatstat (Baddeley and Turner) u(β) u spat (β) = u X z(u) u Q w(u)z(u)λ(u; β) but now Q = X D includes observed points in addition to dummy points D. Note: X events of point pattern (supplied by nature ) D dummy points controlled by scientist. Two types of weights: grid or dirichlet 8 / 23
9 Quadrature schemes in spatstat Grid Dirichlet * * * * * * * * * * * * * w(u) = C v #(X C v )1, u C v where = v D C v w(u) area of Dirichlet cell for u in Dirichlet tesselation generated by Q. 9 / 23
10 spatstat: relation to generalized linear models and iterative weighted least squares Estimating function u spat (β) = u X z(u) u X D w(u)z(u)λ(u; β) formally equivalent to score function of weighted Poisson regression (log link): u spat (β) = u X D w(u)z(u) [ y u λ(u; β) ] with weights w(u) and observations y u = 1[u X]/w(u). Hence implementation straightforward using e.g. glm() in R (Poisson family, log link). 10 / 23
11 Distribution of parameter estimates from approximate score functions? Problem: hard to obtain distribution of parameter estimates from approximate score functions 11 / 23
12 Monte Carlo approximation of integral Consider M random uniform dummy points D on (simple random dummy points). Assume wlog = 1. Rathbun et al. (2006): Monte Carlo approx. of integral: u rath (β) = z(u) 1 z(u)λ(u; β) M u X u D = u(β) [ 1 z(u)λ(u; β) M u D z(u)λ(u; β)du ] 12 / 23
13 Monte Carlo approximation of integral Consider M random uniform dummy points D on (simple random dummy points). Assume wlog = 1. Rathbun et al. (2006): Monte Carlo approx. of integral: u rath (β) = z(u) 1 z(u)λ(u; β) M u X u D = u(β) [ 1 z(u)λ(u; β) M u D z(u)λ(u; β)du ] CLT for Monte Carlo approximation: M 1/2[ 1 z(u)λ(u; β) z(u)λ(u; β)du] d N(0, G) M u D where G asymptotic covariance matrix. 13 / 23
14 Distribution of parameter estimate (Poisson process case) Poisson score asymptotically normal where i(β) Fisher information. u(β) N ( 0, i(β) ) Suppose ˆβ solution of u rath (β) = 0 where u rath (β) = u(β) [ 1 z(u)λ(u; β) M u D z(u)λ(u; β)du ] 14 / 23
15 Distribution of parameter estimate (Poisson process case) Poisson score asymptotically normal where i(β) Fisher information. u(β) N ( 0, i(β) ) Suppose ˆβ solution of u rath (β) = 0 where u rath (β) = u(β) [ 1 z(u)λ(u; β) M u D z(u)λ(u; β)du ] Then ˆβ N ( β, i(β) 1 i(β) 1 G i(β) 1) NB: increasing intensity asymptotics: N = E#X! (fixed observation window ) 15 / 23
16 Stratified dummy points Alternative: one uniformly sampled dummy point in each of M cells (stratified dummy points) Suppose f continuously differentiable. Then CLT [ 1 ] d M f (u) f (u)du N(0, G) M u D (faster rate of convergence). Asymptotic covariance matrix G depends on partial derivatives f i u j of f. 16 / 23
17 Monte Carlo versions of spatstat (aagepetersen, 2007) D point process of dummy points of intensity M. Monte Carlo version of dirichlet u dir (β) = u X z(u) u X D (either simple or stratified dummy points) λ(u; β) z(u) λ(u; β) M 17 / 23
18 Monte Carlo versions of spatstat (aagepetersen, 2007) D point process of dummy points of intensity M. Monte Carlo version of dirichlet u dir (β) = u X z(u) u X D (either simple or stratified dummy points) λ(u; β) z(u) λ(u; β) M Monte Carlo version of grid: (stratified dummy points) u grid (β) = u X z(u) 1 M u X D 1 z(u)λ(u; β) #(X C u ) 1 Advantage: implementation using glm() just as for usual spatstat. 18 / 23
19 Distribution of parameter estimates (Poisson process case) Grid: u grid (β) and u rath (β) with stratified dummy points same asymptotic distribution - hence same asymptotic distribution for parameter estimates. Dirichlet: more subtle behaviour stratified dummy points: u dir (β) less efficient than u grid (β)/u rath (β) simple dummy points: asymptotic covariance matrix tends to that of MLE when q = M/N increases (N = E#X). 19 / 23
20 Numerical example: Poisson process Case of Poisson process with covariate vector (1, z elev ) β elev = 0.1. Ratios of asymptotic standard errors for estimating function estimate ˆβ elev relative to asymptotic standard error for MLE. simple stratified q u rath u grid u dir u dir q = M/N proportion of dummy points. u dir better than Rathbuns Monte Carlo approximation u rath but not useful in case of stratified dummy points. 20 / 23
21 Clustered point processes Consider cluster process X = c Y X c where Y stationary Poisson point process of intensity κ > 0. Given Y, clusters X c are independent Poisson processes with intensity functions λ c (u) = α exp(z 2:p (u)β T 2:p)h(u c; ω) Intensity function of X is then of log-linear form exp(z(u)β T ) where β 1 = log(κα) and z 1 (u) = 1. Estimate ˆβ using u(β), u dir (β) or u grid (β) still asymptotically normal but covariance matrix now depends on ψ = (κ, ω). 21 / 23
22 Minimum contrast estimation of ψ Closed form expression K(t; ψ) and estimate ˆK(t) = u,v X 1[ v u t] λ(u; ˆβ)λ(v; ˆβ) v u available for K-function (extension to inhomogeneous case). 22 / 23
23 Minimum contrast estimation of ψ Closed form expression K(t; ψ) and estimate ˆK(t) = u,v X 1[ v u t] λ(u; ˆβ)λ(v; ˆβ) v u available for K-function (extension to inhomogeneous case). Second step for estimation of ψ: ˆψ = arg min ψ r 0 [ˆK(t) c K(t; ψ) c ] 2 dt NB: only requires z(u) for u X. 23 / 23
24 Inference regarding clustering So far focused on intensity parameter β. Sometimes clustering of main interest but still need to estimate β to adjust for aggregation due to covariates. Joint asymptotic normality of (ˆβ, ˆψ) established in & Guan (2007) for wide class of mixing point processes including cluster processes log Gaussian Cox processes (ψ covariance parameter of Gaussian field). 24 / 23
25 References Baddeley, A., Møller, J., and aagepetersen, R. (2000). Non- and semi-parametric estimation of interaction in inhomogeneous point processes, Statistica Neerlandica, 54, aagepetersen, R. (2007). An estimating function approach to inference for inhomogeneous Neyman-Scott processes, Biometrics, 63, aagepetersen, R. (2007) Estimating functions for inhomogeneous spatial point processes with incomplete covariate data, submitted. aagepetersen, R. and Guan, Y. (2007). Two-step estimation for inhomogeneous spatial point processes, in preparation. 25 / 23
Estimating 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 informationSpatial analysis of tropical rain forest plot data
Spatial analysis of tropical rain forest plot data Rasmus Waagepetersen Department of Mathematical Sciences Aalborg University December 11, 2010 1/45 Tropical rain forest ecology Fundamental questions:
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 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 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 November 8, 2010 1/34
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationSpatial point processes: Theory and practice illustrated with R
Spatial point processes: Theory and practice illustrated with R Department of Mathematical Sciences Aalborg University Lecture IV, February 24, 2011 1/28 Contents of lecture IV Papangelou conditional intensity.
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 informationP n. This is called the law of large numbers but it comes in two forms: Strong and Weak.
Large Sample Theory Large Sample Theory is a name given to the search for approximations to the behaviour of statistical procedures which are derived by computing limits as the sample size, n, tends to
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 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 informationLikelihood-Based Methods
Likelihood-Based Methods Handbook of Spatial Statistics, Chapter 4 Susheela Singh September 22, 2016 OVERVIEW INTRODUCTION MAXIMUM LIKELIHOOD ESTIMATION (ML) RESTRICTED MAXIMUM LIKELIHOOD ESTIMATION (REML)
More information1 Local Asymptotic Normality of Ranks and Covariates in Transformation Models
Draft: February 17, 1998 1 Local Asymptotic Normality of Ranks and Covariates in Transformation Models P.J. Bickel 1 and Y. Ritov 2 1.1 Introduction Le Cam and Yang (1988) addressed broadly the following
More informationChapter 1: A Brief Review of Maximum Likelihood, GMM, and Numerical Tools. Joan Llull. Microeconometrics IDEA PhD Program
Chapter 1: A Brief Review of Maximum Likelihood, GMM, and Numerical Tools Joan Llull Microeconometrics IDEA PhD Program Maximum Likelihood Chapter 1. A Brief Review of Maximum Likelihood, GMM, and Numerical
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 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 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 informationStatistics 222, Spatial Statistics. Outline for the day: 1. Problems and code from last lecture. 2. Likelihood. 3. MLE. 4. Simulation.
Statistics 222, Spatial Statistics. Outline for the day: 1. Problems and code from last lecture. 2. Likelihood. 3. MLE. 4. Simulation. 1. Questions and code from last time. The difference between ETAS
More informationSpatial modelling using INLA
Spatial modelling using INLA August 30, 2012 Bayesian Research and Analysis Group Queensland University of Technology Log Gaussian Cox process Hierarchical Poisson process λ(s) = exp(z(s)) for Z(s) a Gaussian
More informationFor iid Y i the stronger conclusion holds; for our heuristics ignore differences between these notions.
Large Sample Theory Study approximate behaviour of ˆθ by studying the function U. Notice U is sum of independent random variables. Theorem: If Y 1, Y 2,... are iid with mean µ then Yi n µ Called law of
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 informationPh.D. Qualifying Exam Friday Saturday, January 3 4, 2014
Ph.D. Qualifying Exam Friday Saturday, January 3 4, 2014 Put your solution to each problem on a separate sheet of paper. Problem 1. (5166) Assume that two random samples {x i } and {y i } are independently
More informationSpring 2017 Econ 574 Roger Koenker. Lecture 14 GEE-GMM
University of Illinois Department of Economics Spring 2017 Econ 574 Roger Koenker Lecture 14 GEE-GMM Throughout the course we have emphasized methods of estimation and inference based on the principle
More informationGibbs point processes : modelling and inference
Gibbs point processes : modelling and inference J.-F. Coeurjolly (Grenoble University) et J.-M Billiot, D. Dereudre, R. Drouilhet, F. Lavancier 03/09/2010 J.-F. Coeurjolly () Gibbs models 03/09/2010 1
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 informationLinear Methods for Prediction
Chapter 5 Linear Methods for Prediction 5.1 Introduction We now revisit the classification problem and focus on linear methods. Since our prediction Ĝ(x) will always take values in the discrete set G we
More informationGeneralized Linear Models. Kurt Hornik
Generalized Linear Models Kurt Hornik Motivation Assuming normality, the linear model y = Xβ + e has y = β + ε, ε N(0, σ 2 ) such that y N(μ, σ 2 ), E(y ) = μ = β. Various generalizations, including general
More informationAn adapted intensity estimator for linear networks with an application to modelling anti-social behaviour in an urban environment
An adapted intensity estimator for linear networks with an application to modelling anti-social behaviour in an urban environment M. M. Moradi 1,2,, F. J. Rodríguez-Cortés 2 and J. Mateu 2 1 Institute
More informationSTAT 6350 Analysis of Lifetime Data. Failure-time Regression Analysis
STAT 6350 Analysis of Lifetime Data Failure-time Regression Analysis Explanatory Variables for Failure Times Usually explanatory variables explain/predict why some units fail quickly and some units survive
More informationPOLI 8501 Introduction to Maximum Likelihood Estimation
POLI 8501 Introduction to Maximum Likelihood Estimation Maximum Likelihood Intuition Consider a model that looks like this: Y i N(µ, σ 2 ) So: E(Y ) = µ V ar(y ) = σ 2 Suppose you have some data on Y,
More informationGibbs models estimation of galaxy point processes with the ABC Shadow algorithm
Gibbs models estimation of galaxy point processes with the ABC Shadow algorithm Lluı s Hurtado Gil Universidad CEU San Pablo (lluis.hurtadogil@ceu.es) Radu Stoica Universite de Lorraine, IECL (radu-stefan.stoica@univlorraine.fr)
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 informationNormalising constants and maximum likelihood inference
Normalising constants and maximum likelihood inference Jakob G. Rasmussen Department of Mathematics Aalborg University Denmark March 9, 2011 1/14 Today Normalising constants Approximation of normalising
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 informationVariations. ECE 6540, Lecture 10 Maximum Likelihood Estimation
Variations ECE 6540, Lecture 10 Last Time BLUE (Best Linear Unbiased Estimator) Formulation Advantages Disadvantages 2 The BLUE A simplification Assume the estimator is a linear system For a single parameter
More informationMachine learning: Hypothesis testing. Anders Hildeman
Location of trees 0 Observed trees 50 100 150 200 250 300 350 400 450 500 0 100 200 300 400 500 600 700 800 900 1000 Figur: Observed points pattern of the tree specie Beilschmiedia pendula. Location of
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 informationif n is large, Z i are weakly dependent 0-1-variables, p i = P(Z i = 1) small, and Then n approx i=1 i=1 n i=1
Count models A classical, theoretical argument for the Poisson distribution is the approximation Binom(n, p) Pois(λ) for large n and small p and λ = np. This can be extended considerably to n approx Z
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 informationEstimation of Dynamic Regression Models
University of Pavia 2007 Estimation of Dynamic Regression Models Eduardo Rossi University of Pavia Factorization of the density DGP: D t (x t χ t 1, d t ; Ψ) x t represent all the variables in the economy.
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 informationJean-Michel Billiot, Jean-François Coeurjolly and Rémy Drouilhet
Colloque International de Statistique Appliquée pour le Développement en Afrique International Conference on Applied Statistics for Development in Africa Sada 07 nn, 1?? 007) MAXIMUM PSEUDO-LIKELIHOOD
More informationFinal Exam. 1. (6 points) True/False. Please read the statements carefully, as no partial credit will be given.
1. (6 points) True/False. Please read the statements carefully, as no partial credit will be given. (a) If X and Y are independent, Corr(X, Y ) = 0. (b) (c) (d) (e) A consistent estimator must be asymptotically
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 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 informationGraduate Econometrics I: Maximum Likelihood II
Graduate Econometrics I: Maximum Likelihood II Yves Dominicy Université libre de Bruxelles Solvay Brussels School of Economics and Management ECARES Yves Dominicy Graduate Econometrics I: Maximum Likelihood
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 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 informationExploring persistence in the forward rate series
From the SelectedWorks of Derek Bond February, 2009 Exploring persistence in the forward rate series Derek Bond, University of Ulster Available at: https://works.bepress.com/derek_bond/28/ Exploring persistence
More informationSpatial Models: A Quick Overview Astrostatistics Summer School, 2017
Spatial Models: A Quick Overview Astrostatistics Summer School, 2017 Murali Haran Department of Statistics Penn State University April 10, 2017 1 What this tutorial will cover I will explain why spatial
More informationPoisson regression 1/15
Poisson regression 1/15 2/15 Counts data Examples of counts data: Number of hospitalizations over a period of time Number of passengers in a bus station Blood cells number in a blood sample Number of typos
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 informationThe Poisson transform for unnormalised statistical models. Nicolas Chopin (ENSAE) joint work with Simon Barthelmé (CNRS, Gipsa-LAB)
The Poisson transform for unnormalised statistical models Nicolas Chopin (ENSAE) joint work with Simon Barthelmé (CNRS, Gipsa-LAB) Part I Unnormalised statistical models Unnormalised statistical models
More informationStatistical study of spatial dependences in long memory random fields on a lattice, point processes and random geometry.
Statistical study of spatial dependences in long memory random fields on a lattice, point processes and random geometry. Frédéric Lavancier, Laboratoire de Mathématiques Jean Leray, Nantes 9 décembre 2011.
More informationA NOTE ON ROBUST ESTIMATION IN LOGISTIC REGRESSION MODEL
Discussiones Mathematicae Probability and Statistics 36 206 43 5 doi:0.75/dmps.80 A NOTE ON ROBUST ESTIMATION IN LOGISTIC REGRESSION MODEL Tadeusz Bednarski Wroclaw University e-mail: t.bednarski@prawo.uni.wroc.pl
More informationLecture 9 Point Processes
Lecture 9 Point Processes Dennis Sun Stats 253 July 21, 2014 Outline of Lecture 1 Last Words about the Frequency Domain 2 Point Processes in Time and Space 3 Inhomogeneous Poisson Processes 4 Second-Order
More informationVarious types of likelihood
Various types of likelihood 1. likelihood, marginal likelihood, conditional likelihood, profile likelihood, adjusted profile likelihood 2. semi-parametric likelihood, partial likelihood 3. empirical likelihood,
More informationLecture 3 September 1
STAT 383C: Statistical Modeling I Fall 2016 Lecture 3 September 1 Lecturer: Purnamrita Sarkar Scribe: Giorgio Paulon, Carlos Zanini Disclaimer: These scribe notes have been slightly proofread and may have
More informationNow consider the case where E(Y) = µ = Xβ and V (Y) = σ 2 G, where G is diagonal, but unknown.
Weighting We have seen that if E(Y) = Xβ and V (Y) = σ 2 G, where G is known, the model can be rewritten as a linear model. This is known as generalized least squares or, if G is diagonal, with trace(g)
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 informationNon-maximum likelihood estimation and statistical inference for linear and nonlinear mixed models
Optimum Design for Mixed Effects Non-Linear and generalized Linear Models Cambridge, August 9-12, 2011 Non-maximum likelihood estimation and statistical inference for linear and nonlinear mixed models
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2009 Prof. Gesine Reinert Our standard situation is that we have data x = x 1, x 2,..., x n, which we view as realisations of random
More informationInference For High Dimensional M-estimates. Fixed Design Results
: Fixed Design Results Lihua Lei Advisors: Peter J. Bickel, Michael I. Jordan joint work with Peter J. Bickel and Noureddine El Karoui Dec. 8, 2016 1/57 Table of Contents 1 Background 2 Main Results and
More informationLinear Methods for Prediction
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike License. Your use of this material constitutes acceptance of that license and the conditions of use of materials on this
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 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 informationFunctional Principal Component Analysis of Spatio-Temporal Point Processes with Applications in Disease Surveillance
Functional Principal Component Analysis of Spatio-emporal Point Processes with Applications in Disease Surveillance Yehua Li Department of Statistics and Statistical Laboratory, Iowa State University,
More informationSTAT331. Cox s Proportional Hazards Model
STAT331 Cox s Proportional Hazards Model In this unit we introduce Cox s proportional hazards (Cox s PH) model, give a heuristic development of the partial likelihood function, and discuss adaptations
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 informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationSparse Nonparametric Density Estimation in High Dimensions Using the Rodeo
Outline in High Dimensions Using the Rodeo Han Liu 1,2 John Lafferty 2,3 Larry Wasserman 1,2 1 Statistics Department, 2 Machine Learning Department, 3 Computer Science Department, Carnegie Mellon University
More informationPOISSON POINT PROCESS MODELS SOLVE THE PSEUDO-ABSENCE PROBLEM FOR PRESENCE-ONLY DATA IN ECOLOGY 1
The Annals of Applied Statistics 2010, Vol. 4, No. 3, 1383 1402 DOI: 10.1214/10-AOAS331 Institute of Mathematical Statistics, 2010 POISSON POINT PROCESS MODELS SOLVE THE PSEUDO-ABSENCE PROBLEM FOR PRESENCE-ONLY
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 informationGARCH Models Estimation and Inference
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 1 Likelihood function The procedure most often used in estimating θ 0 in
More informationUsing Estimating Equations for Spatially Correlated A
Using Estimating Equations for Spatially Correlated Areal Data December 8, 2009 Introduction GEEs Spatial Estimating Equations Implementation Simulation Conclusion Typical Problem Assess the relationship
More informationEstimation Theory Fredrik Rusek. Chapters 6-7
Estimation Theory Fredrik Rusek Chapters 6-7 All estimation problems Summary All estimation problems Summary Efficient estimator exists All estimation problems Summary MVU estimator exists Efficient estimator
More informationParametric Inference Maximum Likelihood Inference Exponential Families Expectation Maximization (EM) Bayesian Inference Statistical Decison Theory
Statistical Inference Parametric Inference Maximum Likelihood Inference Exponential Families Expectation Maximization (EM) Bayesian Inference Statistical Decison Theory IP, José Bioucas Dias, IST, 2007
More informationEM Algorithm II. September 11, 2018
EM Algorithm II September 11, 2018 Review EM 1/27 (Y obs, Y mis ) f (y obs, y mis θ), we observe Y obs but not Y mis Complete-data log likelihood: l C (θ Y obs, Y mis ) = log { f (Y obs, Y mis θ) Observed-data
More informationQuasi-likelihood Scan Statistics for Detection of
for Quasi-likelihood for Division of Biostatistics and Bioinformatics, National Health Research Institutes & Department of Mathematics, National Chung Cheng University 17 December 2011 1 / 25 Outline for
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 informationOther Survival Models. (1) Non-PH models. We briefly discussed the non-proportional hazards (non-ph) model
Other Survival Models (1) Non-PH models We briefly discussed the non-proportional hazards (non-ph) model λ(t Z) = λ 0 (t) exp{β(t) Z}, where β(t) can be estimated by: piecewise constants (recall how);
More informationA general mixed model approach for spatio-temporal regression data
A general mixed model approach for spatio-temporal regression data Thomas Kneib, Ludwig Fahrmeir & Stefan Lang Department of Statistics, Ludwig-Maximilians-University Munich 1. Spatio-temporal regression
More informationST745: Survival Analysis: Cox-PH!
ST745: Survival Analysis: Cox-PH! Eric B. Laber Department of Statistics, North Carolina State University April 20, 2015 Rien n est plus dangereux qu une idee, quand on n a qu une idee. (Nothing is more
More informationGraduate Econometrics I: Maximum Likelihood I
Graduate Econometrics I: Maximum Likelihood I Yves Dominicy Université libre de Bruxelles Solvay Brussels School of Economics and Management ECARES Yves Dominicy Graduate Econometrics I: Maximum Likelihood
More informationUNIVERSITY OF CALIFORNIA, SAN DIEGO
UNIVERSITY OF CALIFORNIA, SAN DIEGO Estimation of the primary hazard ratio in the presence of a secondary covariate with non-proportional hazards An undergraduate honors thesis submitted to the Department
More information