[R] Post Hoc Groupings

[Jarrett Byrnes]

Indeed, that does it.  Odd.  I guess as Slip was a number, it needed to 
be categorized.  Interesting...

Not if only I can figure out how to get the TukeyHSD grouping table....

On Oct 26, 2005, at 7:35 PM, Mark Lyman wrote:

> Looking at the errors your code produces, it looks like you need to 
> make Dock and Slip factors.
>
> dock_2004_data$Dockf<-factor(dock_2004_data$Dock)
>
> dock_2004_data$Slipf<-factor(dock_2004_data$Slip)
>
> rich.aov <- aov(X.open ~ Dockf*Slipf, data=dock_2004_data)
>
> TukeyHSD(rich.aov, c("Dockf", "Slipf"))
>
>
> Jarrett Byrnes wrote:
>
>> Indeed, the following works as well
>> On Oct 26, 2005, at 5:23 PM, P Ehlers wrote:
>>
>>
>>> fm1 <- aov(breaks ~ wool*tension, data = warpbreaks)
>>> TukeyHSD(fm1, c("wool","tension", "wool:tension"))
>>>
>>
>> However, when working with my own dataset,  I get the following 
>> errors.  I have some inkling this may be due to a slightly unbalanced 
>> sample size, but am uncertain of this.
>>
>> > rich.aov <- aov(X.open ~ Dock*Slip, data=dock_2004_data)
>> > TukeyHSD(rich.aov, c("Dock", "Slip"))
>>
>> Error in rep.int(n, length(means)) : unimplemented type 'NULL' in 
>> 'rep'
>> In addition: Warning messages:
>> 1: non-factors ignored: Slip in: replications(paste("~", xx), data = 
>> mf)
>> 2: non-factors ignored: Dock, Slip in: replications(paste("~", xx), 
>> data = mf)
>>
>> I am pleased to know that these errors are not quite the
>>
>> ______________________________________________
>> 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
>>
>>
>>
>