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