[R-sig-ME] Modeling polar coordinates

David Duffy davidD at qimr.edu.au
Wed Aug 25 01:31:49 CEST 2010


On Tue, 24 Aug 2010, Titus von der Malsburg wrote:

> Your [Luca's] suggestion about using sin and cos combinations sounds
> interesting, but your description is a bit too terse for me.
>
>  mer.cos <- lmer(cos(x) ~ 1 + factor + (1|subject), my.data.frame)
>  mer.sin <- lmer(sin(x) ~ 1 + factor + (1|subject), my.data.frame)
>
> This is of course modulo some transformation that makes the residuals
> being normally distributed.  (Which would that be?)

I can't see it being that simple. AFAICT (Song, Correlated Data Analysis), 
the Fisher and Lee model uses a tan(z/2) link function for the mean, *and* 
1/2 sec^2(z/2) for the dispersion.  You might be able to directly maximize 
the wrapped normal likelihood (from the circular package) for your model.

Just 2c, David Duffy.

-- 
| David Duffy (MBBS PhD)                                         ,-_|\
| email: davidD at qimr.edu.au  ph: INT+61+7+3362-0217 fax: -0101  /     *
| Epidemiology Unit, Queensland Institute of Medical Research   \_,-._/
| 300 Herston Rd, Brisbane, Queensland 4029, Australia  GPG 4D0B994A v




More information about the R-sig-mixed-models mailing list