<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:8px"><div id="yiv6536683256"><div id="yui_3_16_0_1_1416654216114_7711"><div style="color:#000;background-color:#fff;font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:8px;" id="yui_3_16_0_1_1416654216114_7710"><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2539" style="">Bom dia, estou rodando um experimento do tipo Box Behnken e usei a função reg e comparei com o obtido pela função rsm, o problema é que o que os dois consideram significativo difere, aonde posso estar errando? Testei também a função AIC para o modelo que ele fornece e dá diferente pela função AIC().</i></div><div class="yiv6536683256" style="" id="yiv6536683256yui_3_16_0_1_1416654216114_2562"><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div></div><div id="yiv6536683256yui_3_16_0_1_1416654216114_2538"> </div><div class="yiv6536683256signature" id="yiv6536683256yui_3_16_0_1_1416654216114_2341"><span style="font-weight:bold;font-style:italic;" id="yiv6536683256yui_3_16_0_1_1416654216114_2340">Ari Clecius Alves de Lima</span><br style="font-weight:bold;font-style:italic;"><span style="font-weight:bold;font-style:italic;" id="yiv6536683256yui_3_16_0_1_1416654216114_2563">Engenheiro Químico</span><br style="font-weight:bold;font-style:italic;"><span style="font-weight:bold;font-style:italic;" id="yiv6536683256yui_3_16_0_1_1416654216114_2543">Me. Engenharia Civil</span><br style="font-weight:bold;font-style:italic;"><span style="font-weight:bold;font-style:italic;" id="yiv6536683256yui_3_16_0_1_1416654216114_2342">(085)88412345</span><br style="font-weight:bold;font-style:italic;"><span style="font-weight:bold;font-style:italic;" id="yiv6536683256yui_3_16_0_1_1416654216114_2542">(085)33669042</span></div><div class="yiv6536683256signature" id="yiv6536683256yui_3_16_0_1_1416654216114_2341"><span style="font-weight:bold;font-style:italic;"><br></span></div><div class="yiv6536683256signature" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" dir="ltr"><span style="font-weight:bold;font-style:italic;" id="yiv6536683256yui_3_16_0_1_1416654216114_2565">#CMR</span></div><div class="yiv6536683256signature" id="yiv6536683256yui_3_16_0_1_1416654216114_2341"><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><br></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">library(DoE.wrapper)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">Design.1_2<- bbd.design(nfactors= 3 , ncenter= 4 ,randomize= FALSE,seed=7480, </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""> factor.names=list( A=c(0,400),B=c(600,2200),C=c(1,10) ) )</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">res=c(73.450,75.725,63.980,67.670,76.070,77.345,66.005,67.335,84.740,81.595,74.830,70.325,72.460,73.170,73.925,73.800)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">length(res)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#Medidas resumo</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">Design.1_2<- add.response(Design.1_2,res, replace=FALSE)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">plancoded<- code.design(Design.1_2)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">plancoded</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="" id="yui_3_16_0_1_1416654216114_7714">LinearModel.1 <- aov(res ~(x1+x2+x3)^2+ I(x1^2) + I(x2^2) + I(x3^2), </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""> data=plancoded)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">sreg = step(LinearModel.1,direction="both")</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#Melhor modelo</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#Step: AIC=14.07</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="" id="yui_3_16_0_1_1416654216114_7715"># Df Sum of Sq RSS AIC</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#<none> 16.074 14.074</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#+ x1:x2 1 0.501 15.574 15.568</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#+ x2:x3 1 0.462 15.612 15.607</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#+ x1:x3 1 0.001 16.073 16.073</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#- x1 1 9.181 25.255 19.303</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#- I(x2^2) 1 9.310 25.384 19.385</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="" id="yui_3_16_0_1_1416654216114_7716">#- I(x3^2) 1 36.195 52.269 30.941</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#- x2 1 79.223 95.297 40.551</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#- I(x1^2) 1 86.793 102.867 41.774</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#- x3 1 212.747 228.821 54.566</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">summary(sreg)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""># Df Sum Sq Mean Sq F value Pr(>F) </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x1 1 9.18 9.18 5.140 0.04959 * </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x2 1 79.22 79.22 44.357 9.27e-05 ***</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x3 1 212.75 212.75 119.118 1.72e-06 ***</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#I(x1^2) 1 86.79 86.79 48.596 6.53e-05 ***</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#I(x2^2) 1 9.31 9.31 5.213 0.04831 * </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="" id="yiv6536683256yui_3_16_0_1_1416654216114_2527">#I(x3^2) 1 36.20 36.20 20.266 0.00149 ** </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#Residuals 9 16.07 1.79 </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#---</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="" id="yiv6536683256yui_3_16_0_1_1416654216114_2528">#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">Design.1_2.coded<-plancoded</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">rsmModel.2 <- rsm(res ~ FO(x1, x2, x3) + TWI(x1, x2, x3) + PQ(x1, x2, x3), </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""> data=Design.1_2.coded)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">summary(rsmModel.2)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""># Estimate Std. Error t value Pr(>|t|) </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#(Intercept) 73.33875 0.79347 92.4275 1.081e-10 ***</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x1 1.07125 0.56107 1.9093 0.104810 </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x2 -3.14688 0.56107 -5.6087 0.001370 ** </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x3 -5.15688 0.56107 -9.1911 9.350e-05 ***</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x1:x2 0.35375 0.79347 0.4458 0.671360 </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x1:x3 0.01375 0.79347 0.0173 0.986736 </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x2:x3 -0.34000 0.79347 -0.4285 0.683246 </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x1^2 -4.65813 0.79347 -5.8705 0.001081 ** </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x2^2 1.52563 0.79347 1.9227 0.102878 </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#x3^2 3.00813 0.79347 3.7911 0.009062 ** </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#---</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#Multiple R-squared: 0.9664, Adjusted R-squared: 0.916 </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#F-statistic: 19.17 on 9 and 6 DF, p-value: 0.0009315</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">LinearModel.2 <- aov(res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2), </i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""> data=plancoded)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">summary(LinearModel.2 )</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">AIC(LinearModel.2)</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#[1] 61.47997</i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style=""><br class="yiv6536683256" style=""></i></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><i class="yiv6536683256" style="">#valor obtido anteriormente Step: AIC=14.07</i></div></div><div class="yiv6536683256" id="yiv6536683256yui_3_16_0_1_1416654216114_2341" style=""><span style="font-weight:bold;font-style:italic;" class="yiv6536683256"><br class="yiv6536683256" style=""></span></div></div></div></div></div></body></html>