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 . <[email protected]>
Para: [email protected]; Tales Fernandes <[email protected]>
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: [email protected]
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
