Hi Rui, Thanks for responding. I did not write "raw=raw", and I'm not sure why R would return such a misleading error message. Indeed, the same error message comes up when I run the 2nd part of your code:
> m <- matrix(1:6, ncol=2) > p6 <- poly(m, degree=6, raw=TRUE) Error in poly(dots[[1L]], degree, raw = raw) : 'degree' must be less than number of unique points However I can think of a workaround for this, and I get the main idea from your example, which does exactly what I was hoping for, thanks very much! -- View this message in context: http://r.789695.n4.nabble.com/Specify-model-with-polynomial-interaction-terms-up-to-degree-n-tp4635130p4635732.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.
