Estimation of conditional auto-regressive models

<p>Conditional auto-regressive (CAR) models are frequently used with spatial data. However, the likelihood of such a model is expensive to compute even for a moderately sized data set of around 1000 sites. For models involving latent variables, the likelihood is not usually available in closed...

Full description

Bibliographic Details
Main Author: Sha, Z
Other Authors: Ripley, B
Format: Thesis
Published: 2016
Subjects:
_version_ 1826277725890936832
author Sha, Z
author2 Ripley, B
author_facet Ripley, B
Sha, Z
author_sort Sha, Z
collection OXFORD
description <p>Conditional auto-regressive (CAR) models are frequently used with spatial data. However, the likelihood of such a model is expensive to compute even for a moderately sized data set of around 1000 sites. For models involving latent variables, the likelihood is not usually available in closed form.</p> <p>In this thesis we use a Monte Carlo approximation to the likelihood (extending the approach of Geyer and Thompson (1992)), and develop two strategies for maximising this. One strategy is to limit the step size by defining an experimental region using a Monte Carlo approximation to the variance of the estimates. The other is to use response surface methodology. The iterative procedures are fully automatic, with user-specified options to control the simulation and convergence criteria. Both strategies are implemented in our R package mclcar.</p> <p>We demonstrate aspects of the algorithms on simulated data on a torus, and achieve similar results to others in a short computational time on two datasets from the literature. We then use the methods on a challenging problem concerning forest restoration with data from around 7000 trees arranged in transects within study plots. We modelled the growth rate of the trees by a linear mixed effects model with CAR spatial error and CAR random e ects for study plots in an acceptable computational time.</p> <p>Our proposed methods can be used for similar models to provide a clearly defined framework for maximising Monte Carlo approximations to likelihoods and reconstructing likelihood surfaces near the maximum.</p>
first_indexed 2024-03-06T23:33:14Z
format Thesis
id oxford-uuid:6cc56943-2b4d-4931-895a-f3ab67e48e3a
institution University of Oxford
last_indexed 2024-03-06T23:33:14Z
publishDate 2016
record_format dspace
spelling oxford-uuid:6cc56943-2b4d-4931-895a-f3ab67e48e3a2022-03-26T19:13:17ZEstimation of conditional auto-regressive modelsThesishttp://purl.org/coar/resource_type/c_db06uuid:6cc56943-2b4d-4931-895a-f3ab67e48e3aStatisticsSpatial analysis (Statistics)Monte Carlo methodORA Deposit2016Sha, ZRipley, B<p>Conditional auto-regressive (CAR) models are frequently used with spatial data. However, the likelihood of such a model is expensive to compute even for a moderately sized data set of around 1000 sites. For models involving latent variables, the likelihood is not usually available in closed form.</p> <p>In this thesis we use a Monte Carlo approximation to the likelihood (extending the approach of Geyer and Thompson (1992)), and develop two strategies for maximising this. One strategy is to limit the step size by defining an experimental region using a Monte Carlo approximation to the variance of the estimates. The other is to use response surface methodology. The iterative procedures are fully automatic, with user-specified options to control the simulation and convergence criteria. Both strategies are implemented in our R package mclcar.</p> <p>We demonstrate aspects of the algorithms on simulated data on a torus, and achieve similar results to others in a short computational time on two datasets from the literature. We then use the methods on a challenging problem concerning forest restoration with data from around 7000 trees arranged in transects within study plots. We modelled the growth rate of the trees by a linear mixed effects model with CAR spatial error and CAR random e ects for study plots in an acceptable computational time.</p> <p>Our proposed methods can be used for similar models to provide a clearly defined framework for maximising Monte Carlo approximations to likelihoods and reconstructing likelihood surfaces near the maximum.</p>
spellingShingle Statistics
Spatial analysis (Statistics)
Monte Carlo method
Sha, Z
Estimation of conditional auto-regressive models
title Estimation of conditional auto-regressive models
title_full Estimation of conditional auto-regressive models
title_fullStr Estimation of conditional auto-regressive models
title_full_unstemmed Estimation of conditional auto-regressive models
title_short Estimation of conditional auto-regressive models
title_sort estimation of conditional auto regressive models
topic Statistics
Spatial analysis (Statistics)
Monte Carlo method
work_keys_str_mv AT shaz estimationofconditionalautoregressivemodels