<p>Previous papers have elaborated formal gradient analysis for spatial processes, focusing on the distribution theory for directional derivatives associated with a response variable assumed to follow a Gaussian process model. In the current work, these ideas are extended to additionally accommodate one or more continuous covariate(s) whose directional derivatives are of interest and to relate the behavior of the directional derivatives of the response surface to those of the covariate surface(s). It is of interest to assess whether, in some sense, the gradients of the response follow those of the explanatory variable(s), thereby gaining insight into the local relationships between the variables. The joint Gaussian structure of the spatial random effects and associated directional derivatives allows for explicit distribution theory and, hence, kriging across the spatial region using multivariate normal theory. The gradient analysis is illustrated for bivariate and multivariate spatial models, non-Gaussian responses such as presence-absence and point patterns, and outlined for several additional spatial modeling frameworks that commonly arise in the literature. Working within a hierarchical modeling framework, posterior samples enable all gradient analyses to occur as post model fitting procedures.</p> / Dissertation

No description available.

This thesis concerns the development, estimation and investigation of a general anisotropic spatial correlation function, within model-based geostatistics, expressed as a Gaussian linear mixed model, and estimated using residual maximum likelihood (REML). The Matern correlation function is attractive because of its parameter which controls smoothness of the spatial process, and which can be estimated from the data. This function is combined with geometric anisotropy, with an extension permitting different distance metrics, forming a flexible spatial covariance model which incorporates as special cases many infinite- range spatial covariance functions in common use. Derivatives of the residual log-likelihood with respect to the four correlation-model parameters are derived, and the REML algorithm coded in Splus for testing and refinement as a precursor to its implementation into the software ASReml, with additional generality of linear mixed models. Suggestions are given regarding initial values for the estimation. A residual likelihood ratio test for anisotropy is also developed and investigated. Application to three soil-based examples reveals that anisotropy does occur in practice, and that this technique is able to fit covariance models previously unavailable or inaccessible. Simulations of isotropic and anisotropic data with and without a nugget effect reveal the following principal points. Inclusion of some closely-spaced locations greatly improves estimation, particularly of the Matern smoothness parameter, and of the nugget variance when present. The presence of geometric anisotropy does not adversely affect parameter estimation. Presence of a nugget effect introduces greater uncertainty into the parameter estimates, most dramatically for the smoothness parameter, and also increases the chance of non-convergence and decreases the power of the test for anisotropy. Estimation is more difficult with very “unsmooth" processes (Matern smoothness parameter 0.1 or 0.25) | non- convergence is more likely and estimates are less precise and/or more biased. However it is still often possible to fit the full model including both anisotropy and nugget effect using REML with as few as 100 observations. Additional simulations involving model misspecification reveal that ignoring anisotropy when it is present can substantially increase the mean squared error of prediction, but overfitting by attempting to model anisotropy when it is absent is less damaging. Further, plug-in estimates of prediction error variance are reasonable estimates of the actual mean squared error of prediction, regardless of the model fitted, weakening the argument requiring Bayesian approaches to properly allow for uncertainty in the parameter estimates when estimating prediction error variance. The most valuable outcome of this research is the implementation of an anisotropic Matern correlation function in ASReml, including the full generality of Gaussian linear mixed models which permits additional fixed and random effects, making publicly available the facility to fit, via REML estimation, a much wider range of variance models than has previously been readily accessible. This greatly increases the probability and ease with which a well-fitting covariance model can be found for a spatial data set, thus contributing to improved geostatistical spatial analysis. / http://proxy.library.adelaide.edu.au/login?url= http://library.adelaide.edu.au/cgi-bin/Pwebrecon.cgi?BBID=1297562 / Thesis (Ph.D.) -- University of Adelaide, School of Agriculture, Food and Wine, 2007