Bayesian analysis of circular data using wrapped distributions

被引:29
作者
Ravindran P. [1 ]
Ghosh S.K. [2 ]
机构
[1] Hoffmann-La Roche Inc, Nutley
[2] Department of Statistics, NC State University, Raleigh
关键词
Bayesian inference; Markov chain Monte Carlo (MCMC); Wrapped Cauchy; Wrapped Double Exponential; Wrapped Normal;
D O I
10.1080/15598608.2011.10483731
中图分类号
学科分类号
摘要
Circular data arise in a number of different areas such as geological, meteorological, biological and industrial sciences. Standard statistical techniques can not be used to model circular data due to the circular geometry of the sample space. One of the common methods to analyze circular data is known as the wrapping approach. This approach is based on a simple fact that a probability distribution on a circle can be obtained by wrapping a probability distribution defined on the real line. A large class of probability distributions that are flexible to account for different features of circular data can be obtained by the aforementioned approach. However, the likelihood-based inference for wrapped distributions can be very complicated and computationally intensive. The EM algorithm to compute the MLE is feasible, but is computationally unsatisfactory. A data augmentation method using slice sampling is proposed to overcome such computational difficulties. The proposed method turns out to be flexible and computationally efficient to fit a wide class of wrapped distributions. In addition, a new model selection criteria for circular data is developed. Results from an extensive simulation study are presented to validate the performance of the proposed estimation method and the model selection criteria. Application to a real data set is also presented and parameter estimates are compared to those that are available in the literature. © 2011 Copyright Taylor and Francis Group, LLC.
引用
收藏
页码:547 / 561
页数:14
相关论文
共 23 条
[1]  
Best D.J., Fisher N.I., Efficient simulation of the von Mises distribution, Appl. Statist., 28, pp. 152-157, (1978)
[2]  
Coles S., Inference for circular distributions and processes, Statistics and Computing, 8, 2, pp. 105-113, (1998)
[3]  
Damien P., Wakefield J., Walker S., Gibbs sampling for Bayesian non-conjugate hierarchical models by using auxiliary variables, J. R. Statist. Soc. B, 61, pp. 331-344, (1999)
[4]  
Damien P., Walker S., A full Bayesian analysis of circular data using von Mises distribution, The Canadian Journal of Statistics, 27, pp. 291-298, (1999)
[5]  
Fisher N.I., Statistical Analysis of Circular Data, (1993)
[6]  
Fisher N.I., Lee A.J., Time series analysis of circular data, J. R. Statist. Soc. B, 56, pp. 327-339, (1994)
[7]  
Gelfand A.E., Ghosh S.K., Model choice: A minimum posterior predictive loss function, Biometrika, 85, pp. 1-11, (1998)
[8]  
Higdon D.M., Auxiliary variable methods for Markov Chain Monte Carlo with applications, J. Am. Statist. Assoc., 93, pp. 585-595, (1998)
[9]  
Jammalamadaka S.R., Kozubowski T.J., New families of wrapped distributions for modeling skew circular data, Communications in Statistics-Theory and Methods, 33, pp. 2059-2074, (2004)
[10]  
Jammalamadaka S.R., Sengupta A., Topics in Circular Statistics, (2001)