<div dir="ltr">Você precisa estudar os manuais do R (ñ o RStudio) para entender a assim denominada "<i>formula interface</i>" que é praticamente uma sublinguagem dentro do R.<br><br>De fato:<div><font face="monospace">> </font><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:16px">modelo<-lm(y~x1*x2)</span></div><div><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:16px">> </span><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:16px">modelo1<-lm(y~x1*x2)</span></div><div>Produzem regressões idênticas.<span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:16px"><br></span></div><div>Como seria idêntico:</div><div><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:16px">> </span><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:16px">modelo2<-lm(y~x1+x2+x1:x2)</span><br></div><div>Talvez você estivesse pensando em(?):</div><div><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:16px">> </span><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:16px">modelo3<-lm(y~x1:x2)</span><br></div><div><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:16px">> modelo3</span></div><div><font face="monospace">Call:<br>lm(formula = y ~ x1:x2)<br><br>Residuals:<br> 1 2 3 4 5 6 7 <br>-0.325357 0.199643 -0.220357 0.304643 -0.002857 0.037143 0.007143 <br><br>Coefficients:<br> Estimate Std. Error t value Pr(>|t|) <br>(Intercept) 6.43286 0.09079 70.852 1.06e-08 ***<br>x1:x2 -0.03750 0.12011 -0.312 0.767 <br>---<br>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<br><br>Residual standard error: 0.2402 on 5 degrees of freedom<br>Multiple R-squared: 0.01912, Adjusted R-squared: -0.1771 <br>F-statistic: 0.09748 on 1 and 5 DF, p-value: 0.7675</font><br></div><div><br></div><div></div><div>HTH</div><div>--</div><div>Cesar Rabak</div><div> </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Jun 7, 2020 at 7:31 PM Katieli Química por (R-br) <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div><font color="#201f1e"><span style="font-size:14.6667px">Caros amigos usuários do RStudio</span></font></div>
<div><font color="#201f1e"><span style="font-size:14.6667px"><br>
</span></font></div>
<div><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><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><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><font color="#201f1e"><span style="font-size:14.6667px"><br>
</span></font></div>
<div><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><font color="#201f1e"><span style="font-size:14.6667px"><br>
</span></font></div>
<div><font color="#201f1e"><span style="font-size:14.6667px">Dados:</span></font></div>
<div><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><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-size:14.6667px;background-color:rgb(255,255,255);display:inline"><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-size:14.6667px;background-color:rgb(255,255,255);display:inline">
<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-size:14.6667px;background-color:rgb(255,255,255);display:inline"><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-size:14.6667px;background-color:rgb(255,255,255);display:inline"><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">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="gmail-m_854140627346616142Signature">
<div>
<div><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> 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>
</div>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div>