Re: [R-br] Superfície de Resposta

[Emerson Cotta Bodevan via R-br]

Prezados Edson e Felipe,

obrigado pela pronta ajuda.

Att.,


*Emerson*

Em 23 de agosto de 2016 13:57, Felipe 
escreveu:

> Segue um exemplo simples com dados fictícios:
> # Exemplo com x1:Tempo e x2:temperatura e y:reposta
> x1 <- c(-1,-1,1,1,0,0,0,0,0,sqrt(2),-sqrt(2),0,0)
> x2 <- c(-1,1,-1,1,0,0,0,0,0,0,0,sqrt(2),-sqrt(2))
> y2 <- c(76.5,77,78,79.5,79.9,80.3,80,79.7,79.8,78.4,75.6,78.5,77)
> mod <- lm(y2 ~ x1*x2+I(x1^2)+I(x2^2))
>
> summary(mod)
>
> alr3::pureErrorAnova(mod)
>
>
> ### Para construir a superficie de resposta, inicialmente, deve ser gerada
> ### uma sequencia valores dos niveis dos fatores
> x <- seq(-sqrt(2),sqrt(2),l=30)
> y <- seq(-sqrt(2),sqrt(2),l=30)
>
> z <- function(x=x,y=y){
> y <- mod$coef[1]+mod$coef[2]*x+mod$coef[3]*y+mod$coef[4]*x^2+
>  mod$coef[5]*y^2+mod$coef[6]*x*y
> }
> mod
>
> ## Comando outer prepara os dados para serem utilizados pela função persp
> z <- outer(x,y,z)
> z
> persp(x,y,z,theta=-35,phi=5,expand=0.5,xlab='Tempo',ticktype='detailed',
>   ylab='Temperatura',zlab='Resposta',scale=T,col='gray80')
>
>
> image(x,y,z,col=heat.colors(155))
> contour(x,y,z,add=T)
>
> # Incluindo legenda
> fields::image.plot(x,y,z,col=rev(heat.colors(155)))
>
>
>
> # Ponto estacionario
> b <- matrix(c(mod$coef[2],mod$coef[3]));b
> B <- matrix(c(mod$coef[4],mod$coef[6]/2,mod$coef[6]/2,mod$coef[5]
> ),ncol=2);B
>
> x0 <- -.5*solve(B)%*%b;x0
>
> tempo <- x0[1]*5+85
> temperatura <- x0[2]*5+175
>
> image(x,y,z,col=heat.colors(115))
> contour(x,y,z,add=T)
> points(x0[1],x0[2],col='darkgreen',pch=19)
>
> # Ponto de maximo
> y0 <- mod$coef[1]+1/2%*%t(x0)%*%b
> y0
>
> --
> Atenciosamente
> Felipe E. Barletta Mendes
> Estatístico(UFPR) - Conre3 9766-A
> Mestrando em Bioestatística(UEM)
> +55 (41)-92077191
> +55 (41)-33287216
>
>
___
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.

Re: [R-br] Superfície de Resposta

[Felipe via R-br]

Segue um exemplo simples com dados fictícios:
# Exemplo com x1:Tempo e x2:temperatura e y:reposta
x1 <- c(-1,-1,1,1,0,0,0,0,0,sqrt(2),-sqrt(2),0,0)
x2 <- c(-1,1,-1,1,0,0,0,0,0,0,0,sqrt(2),-sqrt(2))
y2 <- c(76.5,77,78,79.5,79.9,80.3,80,79.7,79.8,78.4,75.6,78.5,77)
mod <- lm(y2 ~ x1*x2+I(x1^2)+I(x2^2))

summary(mod)

alr3::pureErrorAnova(mod)


### Para construir a superficie de resposta, inicialmente, deve ser gerada
### uma sequencia valores dos niveis dos fatores
x <- seq(-sqrt(2),sqrt(2),l=30)
y <- seq(-sqrt(2),sqrt(2),l=30)

z <- function(x=x,y=y){
y <- mod$coef[1]+mod$coef[2]*x+mod$coef[3]*y+mod$coef[4]*x^2+
 mod$coef[5]*y^2+mod$coef[6]*x*y
}
mod

## Comando outer prepara os dados para serem utilizados pela função persp
z <- outer(x,y,z)
z
persp(x,y,z,theta=-35,phi=5,expand=0.5,xlab='Tempo',ticktype='detailed',
  ylab='Temperatura',zlab='Resposta',scale=T,col='gray80')


image(x,y,z,col=heat.colors(155))
contour(x,y,z,add=T)

# Incluindo legenda
fields::image.plot(x,y,z,col=rev(heat.colors(155)))



# Ponto estacionario
b <- matrix(c(mod$coef[2],mod$coef[3]));b
B <- matrix(c(mod$coef[4],mod$coef[6]/2,mod$coef[6]/2,mod$coef[5]),ncol=2);B

x0 <- -.5*solve(B)%*%b;x0

tempo <- x0[1]*5+85
temperatura <- x0[2]*5+175

image(x,y,z,col=heat.colors(115))
contour(x,y,z,add=T)
points(x0[1],x0[2],col='darkgreen',pch=19)

# Ponto de maximo
y0 <- mod$coef[1]+1/2%*%t(x0)%*%b
y0

--
Atenciosamente
Felipe E. Barletta Mendes
Estatístico(UFPR) - Conre3 9766-A
Mestrando em Bioestatística(UEM)
+55 (41)-92077191
+55 (41)-33287216

___
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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] Superfície de Resposta

[Edson Lira via R-br]

Veja se ajuda
https://www.r-project.org/conferences/useR-2011/TalkSlides/Invited/Gromping-Design_of_Experiments.pdfhttps://cran.r-project.org/doc/contrib/Vikneswaran-ED_companion.pdfhttps://cran.r-project.org/web/views/ExperimentalDesign.html
 [  ]'s
Edson Lira 
Estatístico
Manaus-Amazonas 

Em Terça-feira, 23 de Agosto de 2016 8:39, Emerson Cotta Bodevan via R-br 
 escreveu:
 

 Caros colegas, bom dia.

Estou estudando sobre superfície de resposta. Alguém poderia me indicar pacotes 
sobre esse assunto no R?

Também, sendo possível, pacotes para planejamento e análises de experimentos.

Agradeço qualquer sugestão e indicação.

Att.,

Emerson

___
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� c�igo m�imo 
reproduz�el.

  ___
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.

[R-br] Superfície de Resposta

[Emerson Cotta Bodevan via R-br]

Caros colegas, bom dia.

Estou estudando sobre superfície de resposta. Alguém poderia me indicar
pacotes sobre esse assunto no R?

Também, sendo possível, pacotes para planejamento e análises de
experimentos.

Agradeço qualquer sugestão e indicação.

Att.,

*Emerson*
___
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.

Re: [R-br] Superfície de Resposta

[Paulo Justiniano]

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.alec...@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.

[R-br] Superfície de Resposta

[Evelyn Alecrim]

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.

Re: [R-br] Superfície de Resposta

[Eder David Borges da Silva]

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.alec...@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.
___
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.

Re: [R-br] Superfície de Resposta

[Mauro Sznelwar]

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.alec...@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.

Re: [R-br] Superfície de Resposta

[Eder David Borges da Silva]

Mauro,
?I()
Att

Em 10 de janeiro de 2012 21:54, Mauro Sznelwar sznel...@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.alec...@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.
___
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.