CircStats and Anova
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
Mueller, Adrienne wrote:
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? Thanks in advance, Adrienne [[alternative HTML version deleted]]
______________________________________________ R-help at stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html