In this thesis we study the problem of estimation of parametric covariance and variogram functions for spatial and spatio- temporal random processes. It has the following principal parts. Variogram Estimation: We consider the "weighted" least squares criterion of fitting a parametric variogram function to second order stationary geo-statistical processes. Two new weight functions are investigated as alternative to the commonly used weight function proposed by Cressie (1985). We discuss asymptotic convergence properties of the sample variogram estimator and estimators of unknown parameters of parametric variogram functions, under a "mixed increasing domain" sampling design as proposed by Lahiriet al. While empirical results of Mean Square Errors, for parameter estimation, obtained using both the proposed functions are found to be comparatively better, we also theoretically establish that under general conditions one of the proposed weight functions give estimates with better asymptotic effciency. Spatio-Temporal Covariance Estimation: Over the past decade, there have been some important advances in methods for constructing valid spatiotemporal covariance functions; but not much attention has been given - so far - on methods of parameter estimation. In this thesis we propose a new frequency domain approach to estimating parameters of spatio-temporal covariance functions. We derive asymptotic strong consistency properties of the estimators using the concept of stochastic equicontinuity. The theory is illustrated with a simulation. Non-Linearity of Geostatistical Data: Linear prediction theory for spatial data is well established and substantial literature is available on the subject. Relatively less is known about non-linearity. In our final and ongoing, research problem we propose a non-linear predictor for geostatistical data. We demonstrate that the predictor is a function of higher order moments. This leads us to construct spatial bispectra for parametric third order moments.

We analyze by simulation the properties of two time domain and two frequency domain estimators for low order autoregressive fractionally integrated moving average Gaussian models, ARFIMA (p,d,q). The estimators considered are the exact maximum likelihood for demeaned data, EML, the associated modified profile likelihood, MPL, and the Whittle estimator with, WLT, and without tapered data, WL. Length of the series is 100. The estimators are compared in terms of pile-up effect, mean square error, bias, and empirical confidence level. The tapered version of the Whittle likelihood turns out to be a reliable estimator for ARMA and ARFIMA models. Its small losses in performance in case of ``well-behaved" models are compensated sufficiently in more ``difficult" models. The modified profile likelihood is an alternative to the WLT but is computationally more demanding. It is either equivalent to the EML or more favorable than the EML. For fractionally integrated models, particularly, it dominates clearly the EML. The WL has serious deficiencies for large ranges of parameters, and so cannot be recommended in general. The EML, on the other hand, should only be used with care for fractionally integrated models due to its potential large negative bias of the fractional integration parameter. In general, one should proceed with caution for ARMA(1,1) models with almost canceling roots, and, in particular, in case of the EML and the MPL for inference in the vicinity of a moving average root of +1. (author's abstract) / Series: Preprint Series / Department of Applied Statistics and Data Processing

No description available.