Fernando Bastos
25 de novembro de 2018
dados <- read.table("Exer_9_1.txt", header=T)
head(dados)## ï..a b bloco resp
## 1 1 1 1 42.9
## 2 1 1 2 41.6
## 3 1 1 3 28.9
## 4 1 1 4 30.8
## 5 1 2 1 53.8
## 6 1 2 2 58.5dim(dados)## [1] 64 4names(dados) <- c("Parcelas","Sub","blocos","Resp")
str(dados)## 'data.frame': 64 obs. of 4 variables:
## $ Parcelas: int 1 1 1 1 1 1 1 1 1 1 ...
## $ Sub : int 1 1 1 1 2 2 2 2 3 3 ...
## $ blocos : int 1 2 3 4 1 2 3 4 1 2 ...
## $ Resp : num 42.9 41.6 28.9 30.8 53.8 58.5 43.9 46.3 49.5 53.8 ...dados$Parcelas <- as.factor(dados$Parcelas)
dados$Sub <- as.factor(dados$Sub)
dados$blocos <- as.factor(dados$blocos)
str(dados)## 'data.frame': 64 obs. of 4 variables:
## $ Parcelas: Factor w/ 4 levels "1","2","3","4": 1 1 1 1 1 1 1 1 1 1 ...
## $ Sub : Factor w/ 4 levels "1","2","3","4": 1 1 1 1 2 2 2 2 3 3 ...
## $ blocos : Factor w/ 4 levels "1","2","3","4": 1 2 3 4 1 2 3 4 1 2 ...
## $ Resp : num 42.9 41.6 28.9 30.8 53.8 58.5 43.9 46.3 49.5 53.8 ...summary(dados)## Parcelas Sub blocos Resp
## 1:16 1:16 1:16 Min. :28.30
## 2:16 2:16 2:16 1st Qu.:44.90
## 3:16 3:16 3:16 Median :52.30
## 4:16 4:16 4:16 Mean :52.81
## 3rd Qu.:62.38
## Max. :75.40dados.av <- aov(Resp ~ blocos + Parcelas*Sub + Error(blocos/Parcelas), data=dados)
summary(dados.av)##
## Error: blocos
## Df Sum Sq Mean Sq
## blocos 3 2843 947.6
##
## Error: blocos:Parcelas
## Df Sum Sq Mean Sq F value Pr(>F)
## Parcelas 3 2848.0 949.3 13.82 0.00102 **
## Residuals 9 618.3 68.7
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## Sub 3 170.5 56.85 2.799 0.05386 .
## Parcelas:Sub 9 586.5 65.16 3.208 0.00595 **
## Residuals 36 731.2 20.31
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1library("ExpDes.pt")##
## Attaching package: 'ExpDes.pt'## The following object is masked from 'package:stats':
##
## ccfpsub2.dbc(dados$Parcelas, dados$Sub, dados$blocos, dados$Resp, quali = c(TRUE, TRUE), mcomp = "tukey", fac.names = c("Parcelas", "Sub"), sigT = 0.05, sigF = 0.05)## ------------------------------------------------------------------------
## Legenda:
## FATOR 1 (parcela): Parcelas
## FATOR 2 (subparcela): Sub
## ------------------------------------------------------------------------
##
## ------------------------------------------------------------------------
## $`Quadro da analise de variancia\n------------------------------------------------------------------------\n`
## GL SQ QM Fc Pr(>Fc)
## Parcelas 3 2848.0 949.34 13.8188 0.001022 **
## Bloco 3 2842.9 947.62 13.7938 0.001029 **
## Erro a 9 618.3 68.70
## Sub 3 170.5 56.85 2.7987 0.053859 .
## Parcelas*Sub 9 586.5 65.16 3.2082 0.005945 **
## Erro b 36 731.2 20.31
## Total 63 7797.4
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## ------------------------------------------------------------------------
## CV 1 = 15.69515 %
## CV 2 = 8.534077 %
##
##
##
## Interacao significativa: desdobrando a interacao
## ------------------------------------------------------------------------
##
## Desdobrando Parcelas dentro de cada nivel de Sub
## ------------------------------------------------------------------------
## GL SQ QM Fc valor.p
## Parcelas : Sub 1 3.00000 1404.1825 468.06083 14.442654 9e-06
## Parcelas : Sub 2 3.00000 412.9700 137.65667 4.247584 0.014042
## Parcelas : Sub 3 3.00000 324.7650 108.25500 3.340355 0.034095
## Parcelas : Sub 4 3.00000 1292.5700 430.85667 13.294668 1.7e-05
## Erro combinado 26.77883 867.8545 32.40823
## ------------------------------------------------------------------------
##
##
## Parcelas dentro de Sub 1
## ------------------------------------------------------------------------
## Teste de Tukey
## ------------------------------------------------------------------------
## Grupos Tratamentos Medias
## a 4 61.925
## ab 3 53.925
## b 2 50.85
## c 1 36.05
## ------------------------------------------------------------------------
##
## Parcelas dentro de Sub 2
## ------------------------------------------------------------------------
## Teste de Tukey
## ------------------------------------------------------------------------
## Grupos Tratamentos Medias
## a 4 63.425
## ab 2 55.375
## b 3 51.375
## b 1 50.625
## ------------------------------------------------------------------------
##
## Parcelas dentro de Sub 3
## ------------------------------------------------------------------------
## Teste de Tukey
## ------------------------------------------------------------------------
## Grupos Tratamentos Medias
## a 4 57.675
## ab 3 55.875
## ab 2 53.1
## b 1 45.85
## ------------------------------------------------------------------------
##
## Parcelas dentro de Sub 4
## ------------------------------------------------------------------------
## Teste de Tukey
## ------------------------------------------------------------------------
## Grupos Tratamentos Medias
## a 4 61.25
## a 3 56.05
## a 2 54.3
## b 1 37.3
## ------------------------------------------------------------------------
##
##
## Desdobrando Sub dentro de cada nivel de Parcelas
## ------------------------------------------------------------------------
## GL SQ QM Fc valor.p
## Sub : Parcelas 1 3 583.49188 194.49729 9.575873 8.7e-05
## Sub : Parcelas 2 3 45.21187 15.07062 0.741987 0.534047
## Sub : Parcelas 3 3 56.95687 18.98562 0.934738 0.433943
## Sub : Parcelas 4 3 71.34188 23.78063 1.170814 0.33435
## Erro b 36 731.20250 20.31118
## ------------------------------------------------------------------------
##
##
## Sub dentro de Parcelas 1
## ------------------------------------------------------------------------
## Teste de Tukey
## ------------------------------------------------------------------------
## Grupos Tratamentos Medias
## a 2 50.625
## ab 3 45.85
## bc 4 37.3
## c 1 36.05
## ------------------------------------------------------------------------
## ------------------------------------------------------------------------
##
##
## Sub dentro de Parcelas 2
## ------------------------------------------------------------------------
## De acordo com o teste F, as medias desse fator sao estatisticamente iguais.
## ------------------------------------------------------------------------
## Niveis Medias
## 1 1 50.850
## 2 2 55.375
## 3 3 53.100
## 4 4 54.300
## ------------------------------------------------------------------------
##
## Sub dentro de Parcelas 3
## ------------------------------------------------------------------------
## De acordo com o teste F, as medias desse fator sao estatisticamente iguais.
## ------------------------------------------------------------------------
## Niveis Medias
## 1 1 53.925
## 2 2 51.375
## 3 3 55.875
## 4 4 56.050
## ------------------------------------------------------------------------
##
## Sub dentro de Parcelas 4
## ------------------------------------------------------------------------
## De acordo com o teste F, as medias desse fator sao estatisticamente iguais.
## ------------------------------------------------------------------------
## Niveis Medias
## 1 1 61.925
## 2 2 63.425
## 3 3 57.675
## 4 4 61.250
## ------------------------------------------------------------------------