Importar Datos y Libreria

En la tabla se observa la base de datos del experimento de germinación de semillas (germ) en función de la humedad del suelo (moisture). Adicionalmente se tiene la variable cobertura (covered) que indica si la caja de las semillas esta cubierta o no. El experimento busca evaluar la relación entre humedad del suelo y germinación.

require(faraway)
data("seeds")
head(seeds,6)
germ moisture covered
22 1 no
41 3 no
66 5 no
82 7 no
79 9 no
0 11 no

A. Exploratorio

En la tabla se observa que el promedio de germinación de las semillas es de 47% con una desviación de 28%, lo cual nos muestra una variación muy alta en la germinación que se supone puede ser explicada por la humedad del suelo. La Figura 1 nos muestra que en experimento 8 de los ensayos no germinaron y algunos lograron un % de germinación alto (mas de 80%). Finalmente la Figura 2 nos muestra que existe una relación no lineal entre humedad del suelo y la germinación indicando que las humedades mas favorables para la germinación se encuentran entre 3 y 7%. Por encima de 7% se observa que la semilla no germina.

promedio=mean(seeds$germ,na.rm=TRUE)
desviacion=sd(seeds$germ,na.rm=TRUE)

data.frame(promedio,desviacion)
promedio desviacion
47.80851 28.04316
require(ggplot2)

p=ggplot(data = seeds,mapping = aes(x=germ))
p+geom_histogram()+ggtitle("Figura 1. Histograma de Germinación")+theme_bw()

p=ggplot(data = seeds,mapping = aes(x=moisture,y=germ))
p+geom_point()+geom_smooth()+ggtitle("Figura 2. Relación Germinación VS Humedad del Suelo")+theme_bw()

La Figura 3 nos muestra que la relación entre germinación y humedad del suelo esta mediada por la cobertura indicando que en ambos casos se logra un tope de germinación pero en diferentes valores de humedad del suelo. Por ejemplo para las cajas cubiertas se recomiendan humedades entre 3 y 5% mientras que para las descubiertas entre 7 y 9% de humedad.

p=ggplot(data = seeds,mapping = aes(x=moisture,y=germ,colour=covered))
p+geom_point()+geom_smooth()+ggtitle("Figura 3. Relación Germinación VS Humedad del Suelo por Cobertura")+theme_bw()