Considere o experimento em faixas de dados de produção de milho descrito na apostila do curso. Clique aqui para ver e copiar o arquivo com conjunto de dados.
O comando read.table é utilizado para ler os dados no R.
ex10 <- read.table("exemplo10.txt", header=T)
ex10
dim(ex10)
[1] 48 4
names(ex10)
[1] "espA" "denB" "bloco" "resp"
ex10$espA <- as.factor(ex10$espA)
ex10$denB <- as.factor(ex10$denB)
ex10$bloco <- as.factor(ex10$bloco)
summary(ex10)
espA denB bloco resp
1:12 1:16 1:12 Min. :16.38
2:12 2:16 2:12 1st Qu.:17.20
3:12 3:16 3:12 Median :17.50
4:12 4:12 Mean :17.44
3rd Qu.:17.68
Max. :18.27
# Interação
attach(ex10)
par(mfrow=c(1,2))
interaction.plot(denB, espA, resp)
interaction.plot(espA, denB, resp)
ADILSON DOS ANJOS 2005-11-07