Trabalhando com o arquivo de dados

Nesse exercício será analisado o experimento fatorial $2^3$ conduzido em um delineamento em blocos completos casualizados apresentado nas notas de aula. Esse experimento refere-se a adubação de pés de café, utilizando-se os nutrientes N (nitrogênio), P (fósforo) e K (potássio) com dois níveis cada.

O arquivo de dados pode ser obtido aqui..

ex05<-read.table("exemplo05.txt",header=T)

Lembre-se que cada coluna indicadora dos fatores e blocos deve ser um fator no R. Por essa razão, os vetores n, p, k e bloc precisam ser transformados para fatores.

> is.factor(ex05$n)
[1] FALSE
> is.factor(ex05$p)
[1] FALSE
> is.factor(ex05$k)
[1] FALSE
> is.factor(ex05$bloc)
[1] FALSE

Para facilitar a digitação, pode-se incluir no caminho de procura o objeto ex05 com o comando attach(ex05).

Agora, deve-se transformar as colunas em fatores

 ex05$n<-as.factor(ex05$n)
 ex05$p<-as.factor(ex05$p)
 ex05$k<-as.factor(ex05$k)
 ex05$bloc<-as.factor(ex05$bloc)

Com o comando summary obtém-se um resumo dos dados do objeto ex05.

> summary(ex05)
 n      p      k      bloc       resp     
 0:24   0:24   0:24   1:8   Min.   :1876  
 1:24   1:24   1:24   2:8   1st Qu.:3171  
                      3:8   Median :3743  
                      4:8   Mean   :3863  
                      5:8   3rd Qu.:4521  
                      6:8   Max.   :6829

Explorando a variável resposta com o comando tapply

> ex05.mn <- tapply(resp, n, mean)
> ex05.mn

> ex05.mp <- tapply(resp, p, mean)
> ex05.mp

> ex05.mk <- tapply(resp, k, mean)
> ex05.mk