Specify model with polynomial interaction terms up to degree n

Hello,

Em 08-07-2012 03:00, Peter Ehlers escreveu:
On 2012-07-07 14:56, YTP wrote:
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
Are you sure that you submitted the poly() exactly as
you show above? I can get your error with raw=FALSE or
NULL or NA or 0 but not with raw=TRUE (or any nonzero
number).

This is with R version 2.15.1 Patched (2012-06-27 r59661).

Peter Ehlers

It works with me too.
This has some days already, but I remember the error message and that I 
had tested it before posting the code above.

sessionInfo()
R version 2.15.0 (2012-03-30)
Platform: i386-pc-mingw32/i386 (32-bit)

locale:
[1] LC_COLLATE=Portuguese_Portugal.1252 
LC_CTYPE=Portuguese_Portugal.1252
[3] LC_MONETARY=Portuguese_Portugal.1252 LC_NUMERIC=C 

[5] LC_TIME=Portuguese_Portugal.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

Rui Barradas
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 at 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.

Specify model with polynomial interaction terms up to degree n

Thread (15 messages)