Prévia do material em texto
EXERCÍCIO PARA ENTREGAR – DBC 14/05/2018
EXERCÍCIO PARA ENTREGAR: Com a finalidade de aumentar a produção de lã de suas ovelhas, por meio de uma alimentação mais
apropriada, um criador separou 28 ovelhas de sua criação. Como eram de idades diferentes elas foram divididas em sete grupos, de modo que
dentro de cada grupo existiam quatro ovelhas com idades similares e homogêneas para as demais características. Em cada grupo foi
realizado um sorteio para distribuir inteiramente ao acaso quatro tipos de alimentação. O experimento iniciou-se no momento de se realizar
uma nova tosquia, obtendo os seguintes resultados expressos em unidades de medidas de lã por animal.
Produção de lã segundo a alimentação ingerida e os grupos homogêneos
## Entrada de dados
dbc <- read.table("XXXXXXXX.txt", h=T, d=","); dbc
TR <- dbc[,1]; TR
BL <- dbc [,2]; BL
Y <- dbc [,3]; Y
FTR <- as.factor(TR) #Toda fonte de variação deve ser um fator
FBL <- as.factor(BL) #Toda fonte de variação deve ser um fator
mod <- aov(Y~FTR+FBL) #Modelo aditivo DBC
## TESTES DE HOMOCEDASTICIDADE.
𝐻0: Os resíduos dos tratamentos ___________________________ homocedásticos (homogeneidade das variâncias).
𝐻1: Os resíduos dos tratamentos ___________________________ homocedásticos (homogeneidade das variâncias).
# Teste de Bartlett
bartlett.test(Y~FTR)
Conclusão. Os resíduos dos tratamentos ___________ homocedásticos, pois 𝑝 −valor = _________ no Teste de ______________
## TESTES DE NORMALIDADE.
𝐻0: Os resíduos do experimento ___________________________ distribuição normal.
𝐻1: Os resíduos do experimento ___________________________ distribuição normal.
# Teste de Shapiro-Wilk
rs <- rstudent(mod); shapiro.test(rs)
Conclusão. Os resíduos do experimento _____________ distribuição normal, pois 𝑝 −valor = _________ no Teste de _________
## TESTE F
# Para Tratamento:
𝐻0: Os alimentos testados possuem efeitos ________________________na a produção de lã de ovelhas.
𝐻0: Os alimentos testados possuem efeitos ________________________na a produção de lã de ovelhas.
# Para Bloco:
𝐻0: Os grupos (blocos) testados possuem efeitos ________________________na a produção de lã de ovelhas.
𝐻0: Os grupos (blocos) testados possuem efeitos ________________________na a produção de lã de ovelhas.
# Tabela de Análise de Variância
anv <- anova(mod); anv
Conclusão.
o Para tratamento
Os alimentos testados possuem efeitos ____________________na a produção de lã de ovelhas., pois 𝑝 −valor = _______________.
Assim, o teste é _________________________ ao nível de significância de _______%. Deve-se ______________________ a hipótese
nula ___________________________________ e concluir que os efeitos dos tratamentos ___________________ entre si ao nível
de significância _______%. Essas diferenças __________________ ser atribuídas ao acaso e _________________ ao efeito dos
tratamentos, com um grau de confiança de ___________________%.
o Para bloco
Os grupos testados possuem efeitos _________________________na a produção de lã de ovelhas, pois 𝑝 −valor = _____________.
Assim, o teste é _________________________ ao nível de significância de _______%. Deve-se ______________________ a hipótese
nula ___________________________________ e concluir que os efeitos dos grupos ___________________ entre si ao nível de
significância _______%. Essas diferenças __________________ ser atribuídas ao acaso e _________________ ao efeito dos
tratamentos, com um grau de confiança de ___________________%.
## TESTE DE TUKEY
#install.packages("agricolae");
require(agricolae)
t_tukey <- HSD.test(mod, "FTR", group=T,alpha=0.05)
t_tukey
𝑚____ = ___________
𝑚____ = ___________
𝑚____ = ___________
𝑚____ = ___________
Médias seguidas de pelo menos uma letra em comum não diferem entre si teste de Tukey, ao nível de significância de 5%
Conclusão:
Verifica-se, em média, uma inferioridade do(s) tratamento(s) _____________ em relação ao(s) tratamento(s) _______________.
Verifica-se, em média, uma superioridade do(s) tratamento(s) ____________ em relação ao(s) tratamento(s) _______________.
## Coeficiente de Variação do Experimento
cv <- 100*sqrt(anv$"Mean Sq"[length(anv$"Mean Sq")])/mean(Y,na.rm=T); cv
Conclusão:
O coeficiente de variação do experimento é de _________%, dessa maneira, conclui-se que o experimento é _______________.