[R-br] Teste de homocedasticidade e ANOVA

[Eder David Borges da Silva]

Flávia,
trat <- gl(5,10)
resp <- rnorm(trat)

m.aov <- aov(resp~trat)
summary(m.aov)

### Conforme sugerido pelo Andre,
require(lawstat)
?levene.test
levene.test(m.aov$residuals,trat)
### Normalmente, pessoalmente uso bartlett
bartlett.test(m.aov$residuals~trat)
### Para normalidade, Shapiro
shapiro.test(m.aov$residuals)
Att
Éder

Em 4 de janeiro de 2012 16:47, Flávia Machado <flavia.dfm em gmail.com> escreveu:
> Galera, alguém sabe me dizer como fazer teste de homocedasticidade
> (levene's test) no R?
>
> Abç,
>
> 2012/1/4 Flávia Machado <flavia.dfm em gmail.com>:
>> Pessoal, trabalho com dieta de muriquis (primata da Mata Atlântica) e
>> quero, com os dados que obtive, verificar qual foi o alimento mais
>> consumido em um determinado período de tempo (todo período de estudo,
>> estação seca e e. chuvosa). São quatro categorias alimentícias (Fruto,
>> Folha, Broto e Flor) e calculei, para cada uma delas, as frequências
>> média de consumo para os três períodos de tempo. Então, cada linha do
>> meu arquivo corresponde a frequência de consumo de cada um desses em
>> um determinado dia, no determinado período.
>>
>> Fiz uma ANOVA, mas não soube interpretar o valor do teste e o
>> bloxplot. De qualquer forma, pode ter dado errado porque eu não sabia
>> que tinha que fazer o teste de homocedasticidade.
>> Não entendi o que preciso colocar para fazer o "levene's test". Na
>> ajuda do R diz:
>>
>> levene.test(y, group, option = c("mean", "median","trim.mean"), trim.alpha = 1)
>>
>> Arguments
>>
>> your data
>> groupfactor of the data
>> optionThe option must be "mean" corresponding to classical Levene's
>> test (default), "median" corresponding to the robust Brown-Forsythe
>> Levene-type test or "trim.mean" corresponding to the robust
>> Levene-type test using the trimmed mean.
>> trim.alphathe fraction (0 to 0.5) of observations to be trimmed from
>> each end of 'x' before the mean is computed.
>>
>> Porém, continuo sem entender... Alguém pode me dar um help?
>>
>> Abraços!
>> Att.
>> --
>>
>> Flávia Machado
>
>
>
> --
>
> Flávia Machado
> _______________________________________________
> R-br mailing list
> R-br em 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.