#=============================================================================== #. Título: Teste de comparação de médias múltiplas - TCMM #. Curso : CET076 - Metodologia e Estatística Experimental - Agronomia # Autor : José Cláudio Faria # Data : 06/04/2022 17:49:12 #=============================================================================== # Objetivos #=============================================================================== # - Exemplificar, através do R, as análises e determinações numéricas do curso # de CET076, referentes a análise de variância e teste de comparação de médias #=============================================================================== #.. Pacotes necessários: # TukeyC (CRAN) # install.packages('TukeyC') #. DIC - balanceado # Leitura na WEB # dad_b <- read.table('https://lec.pro.br/download/R/dados/tcmm_bal.txt', # head=TRUE, # dec=',') # Leitura local # dad_b <- read.table('tcmm_bal.txt', # head=TRUE, # dec=',') # Informando os dados dad_b <- data.frame( tra=gl(4, 5, labels=LETTERS[1:4]), r=1:5, pro=c(25, 26, 20, 23, 21, 31, 25, 28, 27, 24, 22, 26, 28, 25, 29, 33, 29, 31, 34, 28) ) dad_b str(dad_b) av_b <- aov(pro ~ tra, data=dad_b) summary(av_b) #. Carregando o pacote necessário library(TukeyC) tk_b <- TukeyC(av_b, sig.lev=5/100) summary(tk_b) # Observando os todos os resultados do teste tk_b # Observando os resultados mais importantes summary(tk_b) plot(tk_b, disp='cip') #. DIC - desbalanceado # Leitura na web #dad_d <- read.table('https://lec.pro.br/download/R/dados/tcmm_des.txt', # head=TRUE, # dec=',', # na.str='.') # Leitura local #dad_d <- read.table('tcmm_bal.txt', # head=TRUE, # dec=',', # na.str='.') # Informando os dodos dad_d <- data.frame( tra = sort(rep(LETTERS[1:4], 5)), r=1:5, pro= c(NA, 26, 20, 23, 21, 31, 25, 28, 27, 24, 22, 26, NA, 25, 29, 33, 29, 31, 34, NA) ) str(dad_d) dad_d av_d <- aov(pro ~ tra, data=dad_d) summary(av_d) tk_d <- TukeyC(av_d) # Observando os todos os resultados do teste tk_d # Observando os resultados mais importantes summary(tk_d) plot(tk_d, disp='cip')