Qual o significado deste I em I(x^2)+I(z^2) Evelyn, Sempre poste um CMR, ajuda um monte... uma opção pode ser: x <- rnorm(50) z <- rnorm(50) y <- x^2+z^2 m <- lm(y~I(x^2)+I(z^2)) summary(m) Ngrid<-100 xx<-seq(-2,2,l=Ngrid) zz<-seq(-2,2,l=Ngrid) dat <- expand.grid(xx=xx,zz=zz) dat$ypre <- predict(m,newdata=list(x=dat$xx,z=dat$zz)) require(lattice) contourplot(ypre~xx*zz,dat) wireframe(ypre~xx*zz,dat) Veja o help para detalhes de cores, rotação, etc.. Att Em 10 de janeiro de 2012 13:37, Evelyn Alecrim <evelyn.alecrim@gmail.com> escreveu:
Olá Pessoal,
Estou precisando produzir um gráfico de superfície para o seguinte modelo : lm(y=a+bx+cz+dx*z) Alguém poderia me dizer que pacote poderia usar e qual o código? Obrigada...
-- Evelyn da Fonseca Alecrim Bragion Mestranda em Ecologia Aplicada - UFLA Lavras - MG/Brasil