Superfície de Resposta
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
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
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
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
Mauro, ?I() Att Em 10 de janeiro de 2012 21:54, Mauro Sznelwar <sznelwar@uol.com.br> escreveu:
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
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
as is ou seja o argumento é interpretado como é por exemplo x1 * x2 é interpretado como uma interacao em uma formula (se forem continuas, 3 parametros) I(x1*x2) é um unico termo, com um parametro, dado pelo produto das variaveis mais em na sessao 23 Fórmulas e especificação de modelos em http://leg.ufpr.br/~paulojus/embrapa/Rembrapa/ On Tue, 10 Jan 2012, Mauro Sznelwar wrote:
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
esta é uma funcao de 2 variaveis nao precisa pacote especial nenhum basta estipular o grid de valores para x e z, calcular y e plotar com image() persp() contour() ou do lattice wireframe() levelplot() countourplot() os helps destas funcoes contem exemplos assim como demo("graphics") On Tue, 10 Jan 2012, Evelyn Alecrim wrote:
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
participantes (4)
-
Eder David Borges da Silva -
Evelyn Alecrim -
Mauro Sznelwar -
Paulo Justiniano