# Re: [R] CircStats and Anova

From: Spencer Graves <spencer.graves_at_pdf.com>
Date: Mon 16 May 2005 - 05:09:13 EST

It's not clear to me what you are asking. Have you considered "aov.circular" in package "circular"? The following example comes from the help page for that function:

``` >      x <- c(rvonmises(50, 0, 1), rvonmises(100, pi/3, 10))
>      group <- c(rep(0, 50), rep(1, 100))
>
>      aov.circular(x, group)

```

Call:
aov.circular(x = x, group = group)

Circular Analysis of Variance: High Concentration F-Test

df SS MS F p

```Between   1  5.901 5.9009 34.15 3.135e-08
Within  148 29.431 0.1989    NA        NA
Total   149 35.332 0.2371    NA        NA

>      aov.circular(x, group, method="LRT")

```

Call:
aov.circular(x = x, group = group, method = "LRT")

Circular Analysis of Variance: Likelihood Ratio Test

df: 1
ChiSq: 29.31
p.value: 6.163e-08

Does this provide what you want? If no, I suggest you follow the procedure in the posting guide,
"http://www.R-project.org/posting-guide.html". Some people answer their own questions following that process. If they still have questions, they increase their chances of getting useful replies in part because their subsequent posts provides more of the information someone else needs to understand your concern and reply in a way that you will understand. For me, one of the most useful things is a simple example that someone else can easily copy from an email and test in R, craft a reply and send it in a few seconds.

```	  hope this helps.
spencer graves

```

> Hi,
> If I have two sets of directional data (in radians) and want to compare them with a multifactorial anova.
>
> Is it even legitimate to compare circular data with an anova? The books I've picked up from the library don't really say, but it looks unlikely.
>
> If it is allowable, is my having stored the data as circular (X = as.circular(A)) something the aov() function will take into account, assuming the anova would be somewhat different for directional information?
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help