[R-br] metodologia para validação de tratamentos

[Felipe]

Para validar a ANOVA, análise de resíduos:
Considerando um delineamento inteiramente casualizado.

#### HOMOCEDASTICIDADE

boxplot(modelo$res ~ d$trat, ylab='Residuals')
### Gráfico de dispersão dos resíduos vs tratamento
plot.default(d$trat, mod1$res, ylab='Residuals', 
xlab='Tratamento',col='darkblue',cex=.75,pch=16)
#### Testando Homocedasticidade
# H0: As variâncias são comuns
# H1: As variâncias não são comuns
bartlett.test(modelo1$res, d$trat)

#### NORMALIDADE DOS RESÍDUOS
hist(modelo$res,prob=T,main='',ylab='',xlab='',bg='white',border='seagreen',breaks='FD',axes=F)
axis(1)
rug(modelo$res,col='seagreen')
lines(density(modelo$res),col='red')

stem(modelo$res)

qqnorm(modelo$res,ylab='Residuals',col='darkblue',pch=16)
qqline(modelo$res,col='red')
## Teste de Normalidade
# H0:Os resíduos seguem distribuição Normal
# H1:Os resíduos NÃO seguem distribuição Normal
shapiro.test(modelo$res)


###### INDEPENDÊNCIA DOS RESÍDUOS
par(mfrow=c(1,2))
plot(modelo$fit, modelo$res, ylab='Residuals', xlab='Fitted 
Values',col='darkblue',cex=.75,pch=16)
title('Residual vs Fitted Values')

plot(modelo$fit, order(modelo$res), ylab='Residuals', xlab='Fitted 
Values',col='darkblue',cex=.75,pch=16)
title('Residual vs Fitted Values')

####### VERIFICAÇÃO DE OUTLIERS
par(mfrow=c(2,2))
plot(modelo)

names(anova(modelo))
var <- anova(modelo)$Mean[2]
var
res <- modelo$res
resd <- (res/sqrt(var))
boxplot(resd,col='darkblue',pch=16)
plot.default(d$trat,resd,xlab='Tratamento',pch=16,type='p',col='red')
title('Standard Residuals')



On 18-03-2015 16:36, Fernando Antonio de souza wrote:
> Caros amigos,
>
> Estou analisando dados que comparam diferentes metodologias. Os dados 
> referem a medidas de áreas de olho de lombo e estou comparando 4 
> metodologias utilizadas para medila (controle, A, B, C). Eu realizei a 
> anova e utilizei o teste de dunnet para fazer as comparações dos 
> tratamento com o grupo controle.
>
> Gostaria de saber quais avaliações estatísticas a mais eu posso fazer 
> para melhorar esta análise? Quais análises de resíduos eu posso utilizar?
>
> abçs
>
> -- 
> =======================================================================
> Fernando Souza
> Zootecnista, DSc. Produção Animal
> celular: (+55) 82 8113-8781
> e-mail:nandodesouza em gmail.com <mailto:e-mail%3Anandodesouza em gmail.com>
> https://producaoanimalcomr.wordpress.com/
> ========================================================================
>
>
> _______________________________________________
> 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150319/47f755df/attachment.html>