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
>