<div dir="ltr"><div><div>Prezados Edson e Felipe,<br><br></div>obrigado pela pronta ajuda.<br><br></div>Att.,<br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><b><br></b><div style="text-align:left"><i><b>Emerson</b></i><br></div></div></div></div>
<br><div class="gmail_quote">Em 23 de agosto de 2016 13:57, Felipe <span dir="ltr"><<a href="mailto:felipe.e.barletta@gmail.com" target="_blank">felipe.e.barletta@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Segue um exemplo simples com dados fictícios:<br>
##### Exemplo com x1:Tempo e x2:temperatura e y:reposta<br>
x1 <- c(-1,-1,1,1,0,0,0,0,0,sqrt(2),<wbr>-sqrt(2),0,0)<br>
x2 <- c(-1,1,-1,1,0,0,0,0,0,0,0,sqrt<wbr>(2),-sqrt(2))<br>
y2 <- c(76.5,77,78,79.5,79.9,80.3,80<wbr>,79.7,79.8,78.4,75.6,78.5,77)<br>
mod <- lm(y2 ~ x1*x2+I(x1^2)+I(x2^2))<br>
<br>
summary(mod)<br>
<br>
alr3::pureErrorAnova(mod)<br>
<br>
<br>
### Para construir a superficie de resposta, inicialmente, deve ser gerada<br>
### uma sequencia valores dos niveis dos fatores<br>
x <- seq(-sqrt(2),sqrt(2),l=30)<br>
y <- seq(-sqrt(2),sqrt(2),l=30)<br>
<br>
z <- function(x=x,y=y){<br>
    y <- mod$coef[1]+mod$coef[2]*x+mod$<wbr>coef[3]*y+mod$coef[4]*x^2+<br>
         mod$coef[5]*y^2+mod$coef[6]*<wbr>x*y<br>
    }<br>
mod<br>
<br>
## Comando outer prepara os dados para serem utilizados pela função persp<br>
z <- outer(x,y,z)<br>
z<br>
persp(x,y,z,theta=-35,phi=5,ex<wbr>pand=0.5,xlab='Tempo',ticktype<wbr>='detailed',<br>
      ylab='Temperatura',zlab='Respo<wbr>sta',scale=T,col='gray80')<br>
<br>
<br>
image(x,y,z,col=heat.colors(15<wbr>5))<br>
contour(x,y,z,add=T)<br>
<br>
# Incluindo legenda<br>
fields::image.plot(x,y,z,col=r<wbr>ev(heat.colors(155)))<br>
<br>
<br>
<br>
######### Ponto estacionario<br>
b <- matrix(c(mod$coef[2],mod$coef[<wbr>3]));b<br>
B <- matrix(c(mod$coef[4],mod$coef[<wbr>6]/2,mod$coef[6]/2,mod$coef[5]<wbr>),ncol=2);B<br>
<br>
x0 <- -.5*solve(B)%*%b;x0<br>
<br>
tempo <- x0[1]*5+85<br>
temperatura <- x0[2]*5+175<br>
<br>
image(x,y,z,col=heat.colors(11<wbr>5))<br>
contour(x,y,z,add=T)<br>
points(x0[1],x0[2],col='darkgr<wbr>een',pch=19)<br>
<br>
# Ponto de maximo<br>
y0 <- mod$coef[1]+1/2%*%t(x0)%*%b<br>
y0<span class="HOEnZb"><font color="#888888"><br>
<br>
-- <br>
Atenciosamente<br>
Felipe E. Barletta Mendes<br>
Estatístico(UFPR) - Conre3 9766-A<br>
Mestrando em Bioestatística(UEM)<br>
<a href="tel:%2B55%20%2841%29-92077191" value="+554192077191" target="_blank">+55 (41)-92077191</a><br>
<a href="tel:%2B55%20%2841%29-33287216" value="+554133287216" target="_blank">+55 (41)-33287216</a><br>
<br>
</font></span></blockquote></div><br></div>