×

Selection of the neighborhood structure for space-time Markov random field models. (English) Zbl 1145.62376

Summary: A space-time, univariate dataset is assumed to have been sampled from a 3-dimensional Markov random field where the data dependence structure is modeled through pairwise interaction parameters. The likelihood function depends upon (1) an undirected, 3-dimensional graph, where edges connect observation points, and (2) the parameter dimension that captures possible space-time anisotropy of data interaction. Automatic model selection to discriminate both the graph and the model dimension is suggested on the basis of a penalized Pseudo-likelihood function. In most cases, the procedure can be implemented using standard statistical packages capable of GLM estimation. Weak consistency of the criterion is shown to hold under mild and easily verifiable sufficient conditions. Its performance in small samples is studied providing simulation results

MSC:

62M40 Random fields; image analysis
62M05 Markov processes: estimation; hidden Markov models
Full Text: DOI

References:

[1] Akaike H (1969) Fitting Autoregressive models for prediction. Annals of the Institute of Mathematical Statistics21, 243–247 · Zbl 0202.17301 · doi:10.1007/BF02532251
[2] Barnett S (1990) Matrices: methods and applications. New York: Oxford University Press · Zbl 0706.15001
[3] Besag JN (1974) Spatial interactions and the statistical analysis of lattice data. Journal of the Royal Statistical Society B36, 192–225
[4] Besag JN, York J, Mollié A (1991) Bayesian image restoration, with two applications in spatial statistics. Annals of the Institute of Statistical Mathematics43, 1–59 · Zbl 0760.62029 · doi:10.1007/BF00116466
[5] Burnham KP, Anderson DR (1998) Model selection and inference. New York: Springer Verlag
[6] Chiristakos G (1992) Random field models in earth sciences. New York: Academic Press
[7] Cressie N (1993) Statistics for spatial data (revised edition). New York: Wiley
[8] Dobrushin RL (1968) The description of a random field by mean of conditional probabilities and conditions of its regularity. Theory of Probability and its Applications13, 197–224 · doi:10.1137/1113026
[9] Geman D, Geman S (1984) Stochastic relaxation, Gibbs distributions and the Bayesian restoration of images. IEEE PAMI6, 721–741 · Zbl 0573.62030
[10] Geyer CJ, Thompson EA (1992) Constrained Monte Carlo maximum likelihood for dependent data. Journal of the Royal Statistical Society B54, 657–699
[11] Griffith DA, Lagona F (1988) On the quality of likelihood-based estimators in spatial autoregressive models when the data dependence structure is misspecified. Journal of Statistical Planning and Inference69, 153–174 · Zbl 0953.62102 · doi:10.1016/S0378-3758(97)00156-0
[12] Gumpertz ML, Graham JM, Ristaino JB (1997) Autologistic model of spatial patterns of phytophtora epidemic in bell pepper: effect of soil variables on disease presence. Journal of Agricultural, Biological and Environmental Statistics2(2), 131–156 · doi:10.2307/1400400
[13] Guyon X (1995) Random fields on a network. New York: Springer-Verlag · Zbl 0839.60003
[14] Guyon X, Yao JF (1999) On the underfitting and overfitting sets of models chosen by order selection criteria. Journal of Multivariate Analysis70, 221–249 · Zbl 1070.62516 · doi:10.1006/jmva.1999.1828
[15] Hoeffding W (1963) Probability inequality for sums of bounded random variables. Journal of the American Statistical Association58, 12–30 · Zbl 0127.10602 · doi:10.2307/2282952
[16] Huffer FW, Wu H (1998) Markov chain Monte Carlo for autologistic regression models with application to the distribution of plant species. Biometrics54, 509–524 · Zbl 1058.62677 · doi:10.2307/3109759
[17] Hughes JP, Guttorp P, Charles SP (1999) A non-homogeneous hidden Markov model for precipitation occurrence. Applied Statistics48, 15–30 · Zbl 0920.62141
[18] Ji C, Seymour L (1998) A consistent model selection procedure for Markov random fields based on penalized likelihood. Annals of Applied Probability6, 423–443 · Zbl 0856.62082
[19] Jona Lasinio G (2001) Modeling and exploring spatial variation. Journal of Multivariate Analysis77, 295–317 · Zbl 0997.62043 · doi:10.1006/jmva.2000.1938
[20] Kashyap R, Chellappa R (1983) Estimation and choice of neighbours in spatial interaction models of images. IEEE Transactions on Information Theory29, 60–72 · Zbl 0507.62082 · doi:10.1109/TIT.1983.1056610
[21] Knorr-Held L, Besag J (1998) Modelling risk in time and space. Statistics in Medicine17, 2045–2060 · doi:10.1002/(SICI)1097-0258(19980930)17:18<2045::AID-SIM943>3.0.CO;2-P
[22] Lagona F (2002) Adjacency selection in Markov random fields for high resolution, hyperspectral data. Journal of Geographical Systems-Special Issue on High Spatial Resolution Hyperspectral Imagery4(1), 53–68
[23] McCullagh PA, Nelder JA (1989) Generalized linear models 2nd edn. London: Chapman and Hall
[24] Onsanger L (1944) Cristal statistics I: a two-dimensional model with order-disorder transition. Phisical Review65, 117–149 · Zbl 0060.46001 · doi:10.1103/PhysRev.65.117
[25] Petrov VV (1995) Limit theorems in probability theory. Oxford: Clarendon Press
[26] Pickard DK (1987) Inference for discrete Markov fields: the simplest nontrivial case. Journal of the American Statistical Association82, 90–96 · Zbl 0621.62091 · doi:10.2307/2289128
[27] Schwartz G (1978) Estimating the dimension of a model. Annals of Statistics6, 461–464 · Zbl 0379.62005 · doi:10.1214/aos/1176344136
[28] Taxt T, Lundervold A, Angelsen B (1990) Noise Reduction and segmentation in time-varying ultrasound images. In: The Tenth International Conference on Pattern Recognition Atlanta City NJ, pp 591–596
[29] Tonellato SF (2001) A multivariate time series model for the analysis and prediction of carbon monoxide atmospheric concentrations. Applied Statistics50, 187–200
[30] Waller LA, Carlin BP, Xia H, Gelfand A (1997) Hierarchical spatio-temporal mapping of disease rates. Journal of the American Statistical Association92, 607–617 · Zbl 0889.62094 · doi:10.2307/2965708
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.