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)
## Loading required package: faraway
## Warning: package 'faraway' was built under R version 4.1.2
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 |
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 indica una alta variación en la germinación, que se supone puede ser explicadaa por la humedad del suelo. La figura 1 nos muestra que en el experimento 8 de los ensayos no germinaron y algunos lograron un alto porcentaje de germinación (más de 80%). Finalmente la figura 2 nos muestra que existe una relacion 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 del 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)
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 4.1.2
## Warning in register(): Can't find generic `scale_type` in package ggplot2 to
## register S3 method.
p=ggplot(data = seeds,mapping = aes(x=germ))
p+geom_histogram()+ggtitle("Figura 1. Histograma de Germinación")
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## Warning: Removed 1 rows containing non-finite values (stat_bin).
p=ggplot(data = seeds , aes(x=moisture,y=germ))
p+geom_point()+geom_smooth()+ggtitle("Figura 2. Relación Germinación vs Humedad del suelo")
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
## Warning: Removed 1 rows containing non-finite values (stat_smooth).
## Warning: Removed 1 rows containing missing values (geom_point).
La Figura 3 nos 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()
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
## Warning: Removed 1 rows containing non-finite values (stat_smooth).
## Warning: Removed 1 rows containing missing values (geom_point).