Walmes,<div><br></div><div>Do princípio:</div><div><br></div><div>Estou seguindo o livro MASS, segue o CMR (transcição do livro) * dados quine!</div><div><br></div><div># aqui ajusta o modelo &quot;saturado&quot;</div><div>
<div>quine.hi &lt;- aov(log(Days + 2.5) ~ . ^4, quine)</div><div># update do modelo sem a interação quádrupla - não significativa</div><div>quine.nxt &lt;- update(quine.hi, . ~ . - Eth:Sex:Age:Lrn)</div><div>dropterm(quine.nxt, test = &#39;F&#39;)</div>
<div><br></div><div># funcao stepAIC -- o critério por default é o AIC!</div><div>quine.stp &lt;- stepAIC(quine.nxt,</div><div>                     scope = list(upper = ~ Eth * Sex * Age * Lrn, lower = ~ 1),</div><div>                       trace = FALSE)</div>
<div>quine.stp$anova</div><div><br></div><div># A partir daqui ele vai observando os efeitos e tirando conforme o &#39;F&#39;!</div><div>dropterm(quine.stp, test = &#39;F&#39;)</div><div><br></div><div>quine.3 &lt;- update(quine.stp, . ~ . - Eth:Age:Lrn)</div>
<div>dropterm(quine.3, test = &#39;F&#39;)</div><div><br></div><div>quine.4 &lt;- update(quine.3, . ~ . - Eth:Age)</div><div>dropterm(quine.4, test = &#39;F&#39;)</div><div><br></div><div>quine.5 &lt;- update(quine.4, . ~ . - Age:Lrn)</div>
<div>dropterm(quine.5, test = &#39;F&#39;)</div><div><br></div><div>No final o que me deixou confuso, vou transcrever em inglês mesmo: &quot;...        This suggests, correctly, that selecting terms on the basis of AIC can be somewhat permissive in its choice of terms, being roughly equivalent to choosing an F -cutoff of 2. We can proceed manually   or by setting k = 4 in stepAIC. We obtain a model equivalent to Sex/(Age + Eth*Lrn) which is the same as that found by Aitkin (1978), apart from his choice of α = 1 for the displacement constant. ...&quot;</div>
<div><br></div><div>Aí que veio meu questionamento: Ele usa a função stepAIC mas depois vai em busca do &#39;F&#39; e faz essa argumentação da &quot;escolha&quot; de um k tal que iguale o modelo do Aitkin (1978)!</div><div>
<br></div><div>Vi seu chunk number 7, para mim, agora fica claro k = 2 o critério é AIC, k = log(nrow(dados) o critério é o BIC, o que me deixa confuso são esses update() e dropterm()!</div><div><br></div><div>abraço,</div>
<div>FH</div><div><br></div><br><div class="gmail_quote">2012/1/25 Walmes Zeviani <span dir="ltr">&lt;<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>&gt;</span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<font face="trebuchet ms,sans-serif">Fernando,<br><br>Então existem duas funções para stepwise que funcionam de forma diferente, a step() e a stepAIC()? Não entendi. Pelo que sei no R não há stepwise pelo teste F, eles são baseados em critérios de informação. No passado fiz uma busca e parece-me que o step com F está implementado em outro pacote. Um ponto que existe é: tanto F quanto AIC/BIC são funções da soma de quadrados dos resíduos, então, os modelos alcançados não serão, partindo disso, drasticamente diferentes. Sobre essa de tirar variáveis &quot;na mão&quot;, você poderia passar um CMR?<br>

<br>À disposição.<span class="HOEnZb"><font color="#888888"><br>Walmes.<br><br clear="all"></font></span></font><div class="HOEnZb"><div class="h5"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">

<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">

<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: <a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573" target="_blank">(+55) 41 3361 3573</a></span><br style="font-family:trebuchet ms,sans-serif">

<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family:trebuchet ms,sans-serif">

<span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">

<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>

</div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>