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 "saturado"</div><div>
<div>quine.hi <- 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 <- update(quine.hi, . ~ . - Eth:Sex:Age:Lrn)</div><div>dropterm(quine.nxt, test = 'F')</div>
<div><br></div><div># funcao stepAIC -- o critério por default é o AIC!</div><div>quine.stp <- 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 'F'!</div><div>dropterm(quine.stp, test = 'F')</div><div><br></div><div>quine.3 <- update(quine.stp, . ~ . - Eth:Age:Lrn)</div>
<div>dropterm(quine.3, test = 'F')</div><div><br></div><div>quine.4 <- update(quine.3, . ~ . - Eth:Age)</div><div>dropterm(quine.4, test = 'F')</div><div><br></div><div>quine.5 <- update(quine.4, . ~ . - Age:Lrn)</div>
<div>dropterm(quine.5, test = 'F')</div><div><br></div><div>No final o que me deixou confuso, vou transcrever em inglês mesmo: "...        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. ..."</div>
<div><br></div><div>Aí que veio meu questionamento: Ele usa a função stepAIC mas depois vai em busca do 'F' e faz essa argumentação da "escolha" 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"><<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>></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 "na mão", 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>