[R-br] Soma de Quadrado Tipo III no R para analisar um fatorial desbalaceado

Tales Fernandes talesest em yahoo.com.br
Quinta Julho 11 10:09:43 BRT 2013


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 em gmail.com>
Para: r-br em listas.c3sl.ufpr.br; Tales Fernandes <talesest em 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 em ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/03712892/attachment.html>


Mais detalhes sobre a lista de discussão R-br