If I have a matrix x: >x slug surgery swat prey predator 1 1 22 2 91 2 2 40 8 115 3 3 48 3 110
slug = individual is tested in each swat, prey and predator odour treatments surgery = different surgical treatment on slug values in swat, prey, predator columns = average headings of slug 2 levels of treatments: surgery and odour (swat, prey, predator) how do I test (appropriate code) a MANOVA repeated measures with an interaction between surgery and odour (swat,prey,predator)? do I need to re-arrange the matrix? the solution might be: y <- cbind(swat,prey,predator) fit <- manova(y ~ surgery * swat,prey,predator+ Error(slug), data = x) ??? Greg ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.