A Bayesian model for longitudinal circular data

-

English
20 Pages
Read an excerpt
Gain access to the library to view online
Learn more

Description


The analysis of short longitudinal series of circular data may be problematic and to some extent has not been completely developed. In this paper we present a Bayesian analysis of a model for such data. The model is based on a radial projection onto the circle of a particular bivariate normal distribution. Inferences about the parameters of the model are based on samples from the corresponding joint posterior density which are obtained using a Metropolis-within-Gibbs scheme after the introduction of suitable latent variables. The procedure is illustrated both using a simulated data set and a realdata set previously analyzed in the literature.

Subjects

Informations

Published by
Published 01 September 2011
Reads 35
Language English
Report a problem
 Working Paper 11-27 Statistics and Econometrics Series 020 September 2011   Departamento de Estadística  Universidad Carlos III de Madrid Calle Madrid, 126 28903 Getafe (Spain) Fax (34) 91 624-98-49 A BAYESIAN MODEL FOR LONGITUDINAL CIRCULAR DATA  Gabriel Núñez-Antonio* and Eduardo Gutiérrez-Peña**   Abstract  The analysis of short longitudinal series of circular data may be problematic and to some extent has not been completely developed. In this paper we present a Bayesian analysis of a model for such data. The model is based on a radial projection onto the circle of a particular bivariate normal distribution. Inferences about the parameters of the model are based on samples from the corresponding joint posterior density which are obtained using a Metropolis-within-Gibbs scheme after the introduction of suitable latent variables. The procedure is illustrated both using a simulated data set and a real-data set previously analyzed in the literature.   Keywords: Circular data; longitudinal data; Gibbs sampler; latent variables; mixed-effects linear models; projected normal distribution.   * Department of Statistics, Universidad Carlos III de Madrid, Spain. E-mail address: gabriel.nunez@uc3m.es. **Department of Probability and Statistics, IIMAS-UNAM, Mexico. E-mail address: eduardo@sigma.iimas.unam.mx.    
1IntroductionSeveralapproacheshavebeenproposedforanalyzinglongitudinaldata.Forareviewthereaderisreferred,forexample,toDiggleetal.(2002),Fitzmauriceetal.(2004),HedekerandGibbons(2006),andGelmanandHill(2007).Thesebooksalldiscusslongitudinalmodelsfor‘scalar’(i.e.linear)responsesasopposedtocircular.Incontrast,methodologicalproposalstodescriberelation-shipswithinrepeatedmeasurementsofcirculardataareratherlimited.Thismaybeduetothedifficultiesinworkingwithprobabilitydistributionscommonlyassociatedwithdirectionaldataandtotheintrinsicdependencyinherenttolongitudinalstructures.Circulardataareaparticularcaseofdirectionaldata.Specifically,circulardatarepresentdi-rectionsintwodimensions.ForasurveythereaderisreferredtoFisheretal.(1987),Fisher(1993),MardiaandJupp(2000),andJammalamadakaandSenGupta(2001).SeealsoArnoldandSenGupta(2006)foranoverviewoftheapplicationsofcirculardataanalysisinecologicalandenvironmentalsciences.Fromatheoreticalpointofview,therearethreebasicapproachestodirectionalstatistics,whichmaybecalledtheembedding,wrappingandintrinsicapproaches;see,MardiaandJupp(2000).Consequently,thereareseveralwaysofgeneratingprobabilitydistributionsforcirculardata.Onerelativalystraightforwardwayistoradiallyprojectontheunitcircleprobabilitydistributionsoriginallydefinedontheplane.Inthegeneralcase,letYbeaq-dimensionalrandomvectorsuchthatPr(Y=0)=0.ThenU=||Y||1Yisarandompointontheq-dimensionalunitsphere.Itsmeandirectionistheunitvectorη=E(U),whereρ=||E(U)||,0ρ1;hereE()representstheusualexpectationforrandomvectors,and||∙||representstheusualEuclideannorm.Theparameterρiscalledthemeanresultantlengthandrepresentsameasureofconcentrationfordirectionaldistributions(see,forexample,MardiaandJupp,2000,andPresnelletal.,1998).AnimportantinstanceisthatinwhichYhasaq-variateNormaldistribution,Nq(∙|µ,Λ),withmeanvectorµ=E(Y)andprecisionmatrixΛ=Var(Y)1.InthiscaseUissaidtohaveaq-dimensionalprojectednormaldistribution,heredenotedbyPN(∙|µ,Λ).Inthecircularcase,2