#=============================================================================== #. Name : Análise de experimento em parcela subdividida #. Author : J.C.Faria # Date : 25/10/2022 00:48:23 # Aim : Ilustrar análise em PSD # Mail : <<>> #=============================================================================== # Banzatto e Kronka (1992), apresentam o ensaio citado por Steel e Torrie (1980), no # qual são comparadas 4 variedades de aveia e 4 tratamentos de sementes (3 produtos # químicos + testemunha não tratada) quanto aos efeitos sobre a produção. # # As variedades utilizadas (fator A) foram: # A1 – Vicland 1 – infectada com o fungo Helminthosporium victoriae # A2 – Vicland 2 – não infectada # A3 – Clinton – resistente a H. victoriae # A4 – Branch – resistente a H. victoriae # # As variedades foram distribuídas aleatoriamente nas parcelas de cada um dos quatro # blocos do ensaio. # # Os tratamentos de sementes utilizados (fator B) foram: # B1 – Testemunha (não tratado) # B2 – Ceresan M # B3 – Panogen # B4 – Agros # # Os tratamentos de sementes foram aleatoriamente distribuídos nas quatro # subparcelas de cada parcela. #. Leitura dos dados dad <- read.table('https://lec.pro.br/download/R/dados/steel_torrie_1980.txt', dec=',', head=T, stringsAsFactors=TRUE) str(dad) head(dad, 10) #. Análise av <- aov(pro ~ blk + A + B + A*B + Error(blk:A), data=dad) summary(av) #.. Gráfico dos efeitos par(mfrow=c(2, 1)) with(dad, interaction.plot(A, B, pro, xlab='A', ylab='Pro', col='blue3', lwd=2, ylim=c(min(pro), max(pro)))) with(dad, interaction.plot(B, A, pro, xlab='B', ylab='pro', col='blue3', lwd=2, ylim=c(min(pro), max(pro)))) #. Teste de Tukey library(TukeyC) #.. A/B - Supparcela dentro de cada nível da parcela #...A1/B A1_B <- TukeyC(av, which='A:B', fl1=1) summary(A1_B) #...A2/B A2_B <- TukeyC(av, which='A:B', fl1=2) summary(A2_B) #...A3/B A3_B <- TukeyC(av, which='A:B', fl1=3) summary(A3_B) #...A4/B A4_B <- TukeyC(av, which='A:B', fl1=4) summary(A4_B) #.. B/A - Parcela dentro de cada nível da subparcela #...B1/A B1_A <- TukeyC(av, error='Within/blk:A', which='B:A', fl1=1) summary(B1_A) #...B2/A B2_A <- TukeyC(av, error='Within/blk:A', which='B:A', fl1=2) summary(B2_A) #...B3/A B3_A <- TukeyC(av, error='Within/blk:A', which='B:A', fl1=3) summary(B3_A) #...B4/A B4_A <- TukeyC(av, error='Within/blk:A', which='B:A', fl1=4) summary(B4_A) #. Comparando as MSD #... Coparando subparcela dentro de parcela A1_B$out$MSD # 8.58 OK #... Coparando parcela dentro de subparcela (média ponderada de: res(a) e res(b)) B1_A$out$MSD # 11.02 OK