Situación problema (diapositiva 16)

Los siguientes datos corresponden con un estudio que se llevó a ca- bo en el Virginia Polytechnic Institute and State University acerca del desarrollo de una relación entre las raíces de los árboles y la acción de un hongo. Los minerales de los hongos se transfieren a los árboles, y los azúcares de los árboles a los hongos. Se plantaron dos mues- tras de 10 plantones de roble rojo norteño en un invernadero, una de ellas contenía plantones tratados con nitrógeno y la otra plantones sin tratamiento. Todas las demás condiciones ambientales se mantu- vieron constantes. Todos los plantones contenían el hongo Pisolithus tinctorus. Los pesos en gramos de los tallos se registraron después de 140 días y los datos se presentan en la siguiente tabla.

Marco de datos

data1 <- c(1,2,3,4,5,6,7,8,9,10,0.32,0.53,0.28,0.37,0.47,0.43,0.36,
                            0.42,0.38,0.43,0.26,0.43,0.47,0.49,0.52,0.75,0.79,0.86,0.62,0.46)

dtaf <- data.frame(matrix(data1, ncol = 3, nrow = 10))
dtaf
##    X1   X2   X3
## 1   1 0.32 0.26
## 2   2 0.53 0.43
## 3   3 0.28 0.47
## 4   4 0.37 0.49
## 5   5 0.47 0.52
## 6   6 0.43 0.75
## 7   7 0.36 0.79
## 8   8 0.42 0.86
## 9   9 0.38 0.62
## 10 10 0.43 0.46

Diagrama de puntos

r1 <-ggplot(dtaf, aes(x = "", y = X2)) +
  geom_point(color=2)+ggtitle("Diagrama de puntos, Var =0.00529") + 
  ylab("Sin Nitrogeno") +   
  stat_summary(fun = "mean")


r2 <-ggplot(dtaf, aes(x = "", y = X3)) +
  geom_point(color=6)+
  ggtitle("Diagrama de puntos, Var =0.0348") + 
  ylab("Con Nitrogeno") +   
  stat_summary(fun = "mean")


plot_grid(r1, r2,
          labels = c('A', 'B'), # O "AUTO" o "auto"
          label_fontfamily = "serif",
          label_fontface = "bold",
          label_colour = "dodgerblue2")
## Warning: Removed 1 row containing missing values or values outside the scale range
## (`geom_segment()`).
## Removed 1 row containing missing values or values outside the scale range
## (`geom_segment()`).

Boxplots

p1 <- ggplot(dtaf, aes(x="", y=X2))+geom_boxplot()+geom_jitter()+
  geom_boxplot(fill = 2,           # Color caja
               alpha = 0.5,        # Transparencia
               color = 1,          # Color del borde
               outlier.colour = 2)+
  ggtitle("Boxplot") + 
  ylab("Sin Nitrogeno") +   
  stat_summary(fun = "mean") 

p2 <- ggplot(dtaf, aes(x="", y=X3))+geom_boxplot()+geom_jitter()+
  geom_boxplot(fill = 3,           # Color caja
               alpha = 0.5,        # Transparencia
               color = 1,          # Color del borde
               outlier.colour = 2)+
  ggtitle("Boxplot") + 
  ylab("Con Nitrogeno") +   
  stat_summary(fun = "mean")
  
  
plot_grid(p1, p2,
          labels = c('A', 'B'), # O "AUTO" o "auto"
          label_fontfamily = "serif",
          label_fontface = "bold",
          label_colour = "dodgerblue2")
## Warning: Removed 1 row containing missing values or values outside the scale range
## (`geom_segment()`).
## Removed 1 row containing missing values or values outside the scale range
## (`geom_segment()`).

Medidas descriptivas

summary(dtaf)
##        X1              X2               X3        
##  Min.   : 1.00   Min.   :0.2800   Min.   :0.2600  
##  1st Qu.: 3.25   1st Qu.:0.3625   1st Qu.:0.4625  
##  Median : 5.50   Median :0.4000   Median :0.5050  
##  Mean   : 5.50   Mean   :0.3990   Mean   :0.5650  
##  3rd Qu.: 7.75   3rd Qu.:0.4300   3rd Qu.:0.7175  
##  Max.   :10.00   Max.   :0.5300   Max.   :0.8600
var(dtaf$X2)
## [1] 0.005298889
var(dtaf$X3)
## [1] 0.03487222