El siguiente documento muestra el procesamiento de los datos del experimento en conejos de la base rabbit en la libreria faraway.
require(faraway)
data("rabbit")
head(rabbit)
## treat gain block
## 1 f 42.2 b1
## 2 b 32.6 b1
## 3 c 35.2 b1
## 4 c 40.9 b2
## 5 a 40.1 b2
## 6 b 38.1 b2
Como se observa en la tabla, la base cuenta con 3 variables: el tratamiento (treat) que corresponde a la dieta del conejo, la ganancia de peso (gain) y la camada del conejo (block).
require(ggplot2)
g1=ggplot(rabbit,aes(x=gain))+geom_histogram()+theme_bw()
g1
summary(rabbit$gain)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 20.60 35.73 39.30 38.47 42.08 51.90
Se observa en la figura del histograma de la ganancia de peso que los conejos presentan un rango amplio de ganancia de peso que va desde los 20 hasta los 51, y un promedio de peso de 38.
g2=ggplot(rabbit,aes(y=gain,x=treat,fill=treat))+geom_boxplot()+theme_bw()
g2
En la figura se observa que el tratamiento a, presenta el mejor desempeño en terminos de ganancia de peso, seguido del f. por el contrario los tratamientos b y e presental el menor desempeño en ganancia. Adicionalmente el tratamiento e presenta una alta variabilidad es decir conejos con ganancias muy diferentes pero en general bajas comparadas con el a.