Dear all,

I use aov() with weighted means, then I use TukeyHSD() and it shows the following message:

Error in result[, i] <- prj[, select, drop = FALSE] %*% rep.int(1, df[i]) :
number of items to replace is not a multiple of replacement length

How do I fix this error?

Thanks a lot