[R] TukeyHSD fails on my data

[Gav Wood]

Howdo folks,

So I have my data (attached). There are two columns I'm interested in; 
"algname" and "dur". I'd like to know how dur changes with algname. 
algname is nominal and there are 7 possibilities. There are two more 
nominal independents, "task" and "id", so my model is:

dur ~ algname+task+id

 From the research I've done, a TukeyHSD seems to be what I need, so I do:

TukeyHSD(aov(dur ~ algname+task+id, cstats), ordered=T)

But I get this back:

Error in rep.int(n, length(means)) : unimplemented type 'NULL' in 'rep'
In addition: Warning messages:
1: non-factors ignored: task in: replications(paste("~", xx), data = mf)
2: non-factors ignored: id in: replications(paste("~", xx), data = mf)

Can anyone shed any light on the situation?

Gav