Inference in cylindrical models having latent Markovian classes - with an application to ocean current data

Abstract

Spatial direction vector data can be represented cylindrically by linear magnitudes and circular angles. We analyze such data by using a hierarchical Markov random field model with latent discrete classes and conditionally independent cylindrical data given the classes. The structure of a Potts model segments the spatial domain, and each class defines a cylindrical density that represents a specific structure. We consider two types of cylindrical distributions; the Weibull sine-skewed von Mises distribution, which is skewed in the circular part, and the generalized Pareto-type wrapped Cauchy distribution, which is heavy-tailed in the linear part. In this setting, we develop a statistically efficient block composite likelihood method for parameter estimation. The method is shown to provide much faster convergence than an expectation–maximization approach. However, the convergence is less stable for the block composite likelihood method, and we suggest a hybrid estimation approach for practical use. We apply the approach to study ocean surface currents in the Norwegian Sea. The models are able to describe the currents in terms of interpretable local regimes of two–three classes. Scoring rules are used to measure predictive performance of the two cylindrical densities. Results indicate that there is clearly skew angular components, and possibly also some heavy tails in magnitude.