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


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/ 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.