<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;">Caros amigos usuários do RStudio</span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;"><br>
</span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;">Preciso gerar uma superfície de resposta pelo RStudio. Tenho as superfícies geradas pelo Design Expert (bem mais simples de mexer, mas é um software pago e tenho somente a versão de teste)
e comparo as equações geradas com as obtidas pelo R. </span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;">Não tenho experiência com o RStudio, por isso tenho insegurança se o que fiz está ou não correto. </span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;">Poderiam me dar a opinião de vocês? Caso tenham também algum argumento que melhore/otimize meus ajustes, ficaria imensamente grata.</span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;"><br>
</span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;">Meu planejamento é um 2^2 com triplicata do ponto central. Anotei as dúvidas nos scripts.</span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;"><br>
</span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;">Dados:</span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;"><span>x1 <- c(-1,+1,-1,+1,0,0,0)<br>
</span>
<div>x2 <- c(-1,-1,+1,+1,0,0,0)<br>
</div>
<div>y <- c(6.07,6.67,6.25,6.7,6.43,6.47,6.44)</div>
<div><br>
</div>
<div>dados1<-data.frame(x1,x2,y)<br>
</div>
<span></span></span></font></div>
<div style=""><font color="#201f1e"><span style="font-size: 14.6667px;"><span>dados1</span><br>
</span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="color: rgb(32, 31, 30); font-family: "Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 14.6667px; background-color: rgb(255, 255, 255); display: inline !important"><br>
</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
modelo<-lm(y~x1*x2)<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="color: rgb(32, 31, 30); font-family: "Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 14.6667px; background-color: rgb(255, 255, 255); display: inline !important">
<div>summary(modelo)<br>
</div>
<div><br>
</div>
<div><b>O código abaixo fornece o mesmo resultado, não entendo o porquê. Deve ser facultativo usar qualquer um dos scripts.</b></div>
<div><span style="color: rgb(237, 92, 87);">modelo1<-lm(y~x1+x2+x1*x2)</span><br>
</div>
<span style="color: rgb(237, 92, 87);">summary(modelo1)</span><br>
</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="color: rgb(32, 31, 30); font-family: "Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 14.6667px; background-color: rgb(255, 255, 255); display: inline !important"><span><br>
</span></span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="color: rgb(32, 31, 30); font-family: "Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 14.6667px; background-color: rgb(255, 255, 255); display: inline !important"><span><span>anova(modelo)<br>
</span>
<div>library(alr3)<br>
</div>
<span>pureErrorAnova(modelo)</span><br>
</span></span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-size: 14.6667px;">Gera um R^2 ajustado de 0.987 e falta de ajuste não significativa (p>0,05) pela ANOVA.</span><br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-size: 14.6667px;"><br>
</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
mod <- lm(y ~ x1+x2++x1*x2+I(x1^2)+I(x2^2))<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-size: 14.6667px;">
<div>summary(mod)<br>
</div>
<div>anova(mod)<br>
</div>
<span>pureErrorAnova(mod)</span><br>
</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-size: 14.6667px;"><span><b># EM <span style="font-family: Calibri, Helvetica, sans-serif; background-color: rgb(255, 255, 255); display: inline !important">pureErrorAnova(mod) DA ERRO E NÃO GERA OS DADOS . VOCÊS SABEM COMO CONSIGO ARRUMAR
O SCRIPT? PRECISO DOS DADOS DA FALTA DE AJUSTE PARA COMPARAR, EMBORA O R2 AJUSTADO TENHA DADO 0,991 PARA ESSE AJUSTE. PARECE SER MELHOR QUE O ANTERIOR. </span></b></span></span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="color: rgb(237, 92, 87);">Trabalhando com o pacote RSM para superfície de resposta</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>x1 <- c(-1,+1,-1,+1,0,0,0)<br>
</span>
<div>x2 <- c(-1,-1,+1,+1,0,0,0)<br>
</div>
<div>y <- c(6.07,6.67,6.25,6.7,6.43,6.47,6.44)<br>
</div>
<div>dados1<-data.frame(x1,x2,y)<br>
</div>
<span>dados1</span><br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><br>
</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><b>Não consegui gerar o planejamento pelo RStudio. Se eu usar ChemReact, por exemplo, os dados são diferentes e não consigo inserir os meus. Consigo apenas da forma demonstrada acima.</b></span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>library(rsm)<br>
</span>
<div>ajuste <- rsm(y ~ FO(x1, x2), data = dados1)<br>
</div>
<span>summary(ajuste)</span><br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><span>ajuste2 <- rsm(y ~ FO(x1, x2)+ TWI(x1, x2), data = dados1)<br>
</span><span>summary(ajuste2)</span><br>
</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><span><br>
</span></span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><span><b>o ajuste abaixo gera estruturas de confundimento, é isso? Fornece um erro ao tentar rodar.</b></span></span></div>
<div id="Signature">
<div>
<div style=""><span>ajuste3 <- rsm(y~SO(x1,x2), data = dados1)<br>
</span>
<div>summary(ajuste3)<br>
</div>
<div><br>
</div>
<div><span style="color: rgb(237, 92, 87);">graficos:</span></div>
<div><span style="color: rgb(59, 32, 77);"><b>Fiquei na dúvida porque o ajuste 2 parecia melhor (maior R2 ajustado), mas o gráfico de resíduos ficou estranho. Considerando os resíduos, o ajuste 1 parece mais adequado. O que vocês acham?</b></span></div>
<div><span>contour(ajuste2,~x1+x2, image=TRUE, img.col = terrain.colors(50),<br>
</span><span> xlabs=c("x1","x2"))</span><br>
</div>
<div><span><br>
</span></div>
<div><span><span>persp(ajuste,~x1+x2, col=rainbow(50), contours = "colors",<br>
</span>
<div> xlabs=c("x2","x3"),zlab="Resposta", theta = -30, <br>
</div>
<span> phi= 50)</span><br>
</span></div>
<div><span><span><br>
</span></span></div>
<div><span><span><b>O plot de resíduos achei estranho</b></span></span></div>
<div><span><span>plot(fitted(ajuste2),resid(ajuste2),xlab = "Valores ajustados",<br>
<div style=""> ylab = "Resíduos", main = "Resíduos x Valores ajustados")<br>
</div>
abline(h=0,lty=2,col=2)<br>
</span></span></div>
<br>
<br>
Obrigada pela ajuda
<table width="500" style="color:rgb(97,97,97); font-family:Lato,sans-serif; font-size:13px; text-align:center; border-spacing:0px; border-collapse:collapse; background-color:rgb(255,255,255)">
<tbody>
<tr>
<td>
<table style="text-align:center; border-collapse:collapse">
<tbody>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
<div style="text-align:center"><font color="#000000" face="Calibri, Helvetica, sans-serif"><br>
</font></div>
</div>
</div>
</div>
</div>
</body>
</html>