#------------------------------------------- # Texto sem acentos! by Cristian #------------------------------------------- ## Planejamento de um experimento ## #------------------------------------------- # Delineamento inteiramente # casualizado #------------------------------------------- set.seed(1234) sample(rep(c("A", "B", "C", "D"), 5)) #------------------------------------------- ## Fazer o croqui correspondente # colocando seu numero USP como # semente #------------------------------------------- ## Entrada dos dados para analise #------------------------------------------- # Copiar o arquivo aula2.csv para seu # computador #------------------------------------------- dados=read.csv2("aula2.csv") summary(dados) attach(dados) y dados$y #------------------------------------------- ## Analise Exploratoria ## #------------------------------------------- #install.packages("lattice") library(lattice) dotplot(y~trat) boxplot(y~trat) #------------------------------------------- # Estatisticas descritivas #------------------------------------------- n = tapply(y,trat, length); n soma = tapply(y,trat,sum); soma media = tapply(y,trat,mean); media variancia = tapply(y,trat,var); variancia desv.padr = tapply(y,trat,sd); desv.padr dist.int = tapply(y,trat,IQR); dist.int #------------------------------------------- # Criando uma funcao que calcula a # amplitude #------------------------------------------- f1=function(x) max(x)-min(x) amplitude = tapply(y,trat,f1); amplitude #------------------------------------------- resumo=rbind(n, soma, media, variancia, desv.padr, amplitude,dist.int) rownames(resumo)=c("n", "Soma", "Media", "Variancia", "Desvio-padrao", "Amplitude", "Amplitude Interquartilica"); round(resumo,2) #------------------------------------------- ## Analise de variancia ## #------------------------------------------- modelo=aov(y ~ trat) anova(modelo) qf(0.95,3,16, lower.tail = TRUE) qf(0.05,3,16, lower.tail = FALSE) detach(dados) #-------------------------------------------