Lucas, Sem um CMR (v. o guia de postagem a respeito) este seu último e-mail vira o que os *nerds* denominaram jocosamente de « *confessional debugging* »! Por óbvio se você a variável contiver dados advindos dos resíduos de uma ANOVA o teste testará apenas eles. Por outro lado se a variável resposta diferisse apenas no valor médio e as variâncias fossem similares os testes não deveriam divergir do ponto de vista estatístico, ou seja, o grande resultado deveria bater, mesmo se na minúcia os números não fossem exatamente iguais. Testes como esse tendem a ser "sensíveis demais" se a amostra for muito grande, outro dado que nos permitiria ajudá-lo mais seria você no CMR postar também a tabela ANOVA decorrente da análise de variância nos dados, e rodar os diagnósticos gráficos sobre a ANOVA, que podem lhe dar *insights* sobre outros pressupostos da ANOVA. HTH -- Cesar Rabak 2018-03-08 15:17 GMT-03:00 Lucas Stempkowski via R-br < r-br@listas.c3sl.ufpr.br>:
Obrigado Marcus,
No primeiro comando que usei #1 considerei a anova (ex.av) onde res seriam os resíduos.
bartlett.test(ex.av$res~trat) #1
No segundo comando #2 não considerei a anova.
bartlett.test(ex$resp, ex$trat) #2
No primeiro comando, pelo que entendi agora, o teste considerou os resíduos..... no segundo os dados de forma geral, não sei se estou certo.
*--*
*-----------------------------------------------------*
Lucas Antonio Stempkowski
Eng. Agronômo - UFFS
Mestrando em Produção Vegetal CAV/UDESC
+55 54 99153-0124/ 98448-4633
--------------------------------------------------
------------------------------ *De:* Marcus Nunes <marcus.nunes@gmail.com> *Enviado:* quinta-feira, 8 de março de 2018 14:27:00 *Para:* Lucas Stempkowski; a lista Brasileira oficial de discussão do programa R. *Assunto:* Re: [R-br] AJUDA: Teste de Homogeneidade
Ambas as formas são idênticas. Elas não devem dar resultados diferentes. Rode os exemplos presentes no help da função para verificar isto:
?bartlett.test bartlett.test(InsectSprays$count, InsectSprays$spray) bartlett.test(count ~ spray, data = InsectSprays)
Eu desconfio que há um erro no teu código. Em primeiro lugar, parece que há dois bancos de dados: ex e ex.av. Qual dos dois é o correto? Assumindo que ex é o correto, rode os seguintes comandos
bartlett.test(ex$resp, ex$trat) # 1 bartlett.test(res ~ trat, data=ex) # 2
e veja que os resultados foram iguais.
A diferença entre uma forma e outra está apenas na sintaxe. Pessoalmente, prefiro a forma # 2, pois acho ela mais elegante. Entretanto, esta é apenas a minha opinião pessoal.
-- Marcus Nunes Professor Adjunto Universidade Federal do Rio Grande do Norte Centro de Ciências Exatas e da Terra Departamento de Estatística Laboratório de Estatística Aplicada marcus.nunes@ccet.ufrn.br http://marcusnunes.me/
2018-03-08 14:11 GMT-03:00 Lucas Stempkowski via R-br < r-br@listas.c3sl.ufpr.br>:
Boa tarde pessoal,
Sou iniciante com o R, e estou com uma dúvida sobre o teste de bartlett. Vi alguns exemplos sobre este teste, e , alguns é utilizado o comando:
bartlett.test(ex.av$res~trat) e, em outros o comando:
bartlett.test(ex$resp, ex$trat)
testando os dois observei valores distintos, mas quando devo usar cada um deles?
Desde já agradeço Lucas
*--*
*-----------------------------------------------------*
Lucas Antonio Stempkowski
Eng. Agronômo - UFFS
Mestrando em Produção Vegetal CAV/UDESC
+55 54 99153-0124 <+55%2054%2099153-0124>/ 98448-4633
--------------------------------------------------
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.