
Walmes, muito obrigado. Concordo que o problema realmente está em formular a hipóteses e não em como calcular as diferentes SQ. Seguindo suas sugestões então devo fazer as SQ Tipo II, também acho mais adequado. Mas isso não responde minha dúvida inicial: Se por algum motivo eu preciso/desejo analisar um experimento no R, usando a SQ tipo III, posso confiar que o comando Anova(objeto, type="III") do pacote car está fazendo estas reduções? E que por algum motivo (alguém sabe qual?) fornece resultados diferentes do SAS. Alguém conhece as reduções que são feitas nesta função, ou sabe onde posso encontrá-las? Atenciosamente; Tales Jesus Fernandes Doutorando em Estatística UFLA Universidade Federal de Lavras ________________________________ De: walmes . <walmeszeviani@gmail.com> Para: r-br@listas.c3sl.ufpr.br; Tales Fernandes <talesest@yahoo.com.br> Enviadas: Quarta-feira, 10 de Julho de 2013 10:03 Assunto: Re: [R-br] Soma de Quadrado Tipo III no R para analisar um fatorial desbalaceado Verifique a função drop1(). O erro é se você tem caselas pedidas então o efeito é não estimado e isso gera estrutura de confundimento etc. Você disse que a hipótese que você quer testar tem que usar a SQ tipo III, ou seja, todos os efeitos ajustados para os demais do modelo. Sejam os efeitos A e B, então você tem a interação ajustada para os efeitos principais, A*B|A e B, coisa que você testa com a SQ I. Agora, faz sentido testar o efeito de A ajustado para uma interação que contém A, A|B e A*B? Existe um modelo onde B interage com A mas sem o efeito principal de A? Não existe, então porque essa hipótese é relevante? Para mim existe mais mistério na insistências das pessoas em porque usar/ensinar essa SQ III do que como calculá-la. Recomendo que leia http://www.stats.ox.ac.uk/pub/MASS3/Exegeses.pdf À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================