#install.packages("tidyverse")
#install.packages("summarytools")
#install.packages("leaflet")
# Instalamos y/o cargamos paquetes
library(readxl)      
library(tidyverse)   
library(summarytools) 
library(leaflet)

1. LA BASE DE DATOS

Base de datos CAPIA

Función read_excel()

capiass <-read_excel("capiass.xlsx")
capiass 
## # A tibble: 30 × 8
##        N fertilizante altura_de_planta N_hojas N_mazorcas Largo_marzorca_1
##    <dbl> <chr>                   <dbl>   <dbl>      <dbl>            <dbl>
##  1     1 orgánico                 2.45      14          1               30
##  2     2 orgánico                 2.35      15          2               23
##  3     3 orgánico                 2.49      19          1               32
##  4     4 orgánico                 2.25      15          0               30
##  5     5 orgánico                 2.5       16          1               30
##  6     6 orgánico                 2.6       15          1               30
##  7     7 orgánico                 1.8       17          0                0
##  8     8 orgánico                 1.95      17          2               30
##  9     9 orgánico                 1.96      19          2               29
## 10    10 orgánico                 1.8       17          0               23
## # ℹ 20 more rows
## # ℹ 2 more variables: Largo_marzorca_2 <dbl>, Largo_marzorca_3 <dbl>
glimpse(capiass)
## Rows: 30
## Columns: 8
## $ N                <dbl> 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16…
## $ fertilizante     <chr> "orgánico", "orgánico", "orgánico", "orgánico", "orgá…
## $ altura_de_planta <dbl> 2.45, 2.35, 2.49, 2.25, 2.50, 2.60, 1.80, 1.95, 1.96,…
## $ N_hojas          <dbl> 14, 15, 19, 15, 16, 15, 17, 17, 19, 17, 15, 19, 17, 1…
## $ N_mazorcas       <dbl> 1, 2, 1, 0, 1, 1, 0, 2, 2, 0, 0, 2, 3, 2, 2, 1, 2, 0,…
## $ Largo_marzorca_1 <dbl> 30, 23, 32, 30, 30, 30, 0, 30, 29, 23, 31, 32, 32, 29…
## $ Largo_marzorca_2 <dbl> NA, 18, 0, 0, 0, 0, 0, 15, 20, 0, 0, 21, 24, 26, 21, …
## $ Largo_marzorca_3 <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10, 0, 0, 0, 0, 0…

#La base de datos CAPIA tiene 30 casos y 7 variables ( fertilizante, altura de planta, número de hojas, número de mazorca, largo de mazorca 1, largo de mazorca 2 y largo de mazorca 3).

##2. Datos estadísticos, media, valor mínimo y valor máximo para las variables estudiadas

summary(capiass)
##        N         fertilizante       altura_de_planta    N_hojas     
##  Min.   : 1.00   Length:30          Min.   :1.280    Min.   :13.00  
##  1st Qu.: 8.25   Class :character   1st Qu.:1.820    1st Qu.:15.00  
##  Median :15.50   Mode  :character   Median :2.100    Median :16.50  
##  Mean   :15.50                      Mean   :2.133    Mean   :16.40  
##  3rd Qu.:22.75                      3rd Qu.:2.450    3rd Qu.:17.75  
##  Max.   :30.00                      Max.   :2.900    Max.   :19.00  
##                                                                     
##    N_mazorcas   Largo_marzorca_1 Largo_marzorca_2 Largo_marzorca_3 
##  Min.   :0.00   Min.   : 0.00    Min.   : 0.000   Min.   : 0.0000  
##  1st Qu.:0.25   1st Qu.:24.25    1st Qu.: 0.000   1st Qu.: 0.0000  
##  Median :1.00   Median :30.00    Median : 0.000   Median : 0.0000  
##  Mean   :1.10   Mean   :25.10    Mean   : 7.552   Mean   : 0.3333  
##  3rd Qu.:2.00   3rd Qu.:31.75    3rd Qu.:20.000   3rd Qu.: 0.0000  
##  Max.   :3.00   Max.   :39.00    Max.   :28.000   Max.   :10.0000  
##                                  NA's   :1

Para la variable estudiada altura de planta, la media es de 2,10 metros, el valor mínimo 1,28 metro y el valor máximo 2,90 metros.

En el número de hojas, la media es de 16, 5, con un valor mínimo de 13 y un máximo de 19 y para el número de marzorcas la media es de 1 mazrca por planta.

freq(capiass$fertilizante, justify =  "center", #Controla la alineación del texto centrada en la tabla.
       report.nas = FALSE, # no se incluiran los valores faltantes (NA) en la tabla.
       headings = FALSE, # no se muestran los encabezados o títulos de la tabla. 
       cumul = FALSE) #indica que no se mostrarán los porcentajes acumulados en la tabla.)
## 
##                    Freq     %    
## ----------------- ------ --------
##     orgánico        15    50.00  
##    tradicional      15    50.00  
##       Total         30    100.00

3. TIPOS DE FERTILIZANTES

Fertilizante orgánico y tradicional

Filtrar los datos para fertilizante

orgánico <- capiass %>% filter(fertilizante == "orgánico")
tradicional<-capiass%>% filter(fertilizante == "tradicional")

##4. TIPOS DE FERTILIZANTES SEGÚN LA VARIABLE ESTUDIADA ALTURA DE PLANTA

ggplot(capiass, aes(fertilizante,altura_de_planta)) +
geom_boxplot()

ggplot(capiass, aes(fertilizante, altura_de_planta, color = fertilizante)) +
geom_boxplot() +
stat_summary(fun = mean, color = "blue", size = 0.2, shape = 4) +
labs(x = "fertilizante",
y = "altura de planta (m)") +
theme_classic() +
theme(axis.title.x = element_text(size = 18),  # Tamaño de letra del título del eje x
axis.title.y = element_text(size = 18),  # Tamaño de letra del título del eje y
axis.text.x = element_text(size = 14),   # Tamaño de letra de las categorías del eje x
axis.text.y = element_text(size = 14),
legend.position = "none") +  # Elimina la leyenda
scale_y_continuous(limits = c(0,4)) + # Ajusta los límites del eje y
scale_color_manual(values = c("organico" = "red", "tradicional" = "orange"))
## Warning: Removed 2 rows containing missing values or values outside the scale range
## (`geom_segment()`).

Se observa una media ligeramente mayor de altura de planta cuando se aplico en el maiz el fertilizante organico.

4 - Variable Distribución de altura de planta

ggplot(orgánico, aes(altura_de_planta)) +
  geom_histogram(bins=7,color = "yellow", fill = "darkgreen")

# la mayor distribución de freciencias esta cerca a 2 y 2,5 metros de altura de planta, siendo 2,5 metros la atura mayor. Aparentemente no habría plantas con alturas cercanas al rango de 1,5 metros.

5. Valores medios de altura de planta comparando los dos fertilzantes

ggplot(capiass, aes(fertilizante,altura_de_planta, color = fertilizante)) + 
  geom_boxplot() +
  stat_summary(fun = mean, color = "black") # función para agregar la media muestral y asignarle color
## Warning: Removed 2 rows containing missing values or values outside the scale range
## (`geom_segment()`).

Con la aplicación del fertilizante orgánico se logra valores mas altos de altura de planta y un valor de media tambien mayor. Por lo que se podría recomendar esa aplicación de fertilizante para lograr mayor altura de planta.

6. Ubicación geográfica del ensayo de maiz capia.

Localidad San Jose norte, Departamento Santa María, provincia de Catamarca

library(leaflet)
leaflet() %>%
  addProviderTiles(providers$Esri.WorldImagery) %>%
  setView(lng = -66.052477, lat = -26.791282, zoom = 16) %>% 
  addMarkers(lng = -66.052477, lat = -26.791282)
capiass%>% 
  group_by(fertilizante) %>%
  descr(Largo_marzorca_1)
## Descriptive Statistics  
## Largo_marzorca_1 by fertilizante  
## Data Frame: capiass  
## N: 30  
## 
##                     orgánico   tradicional
## ----------------- ---------- -------------
##              Mean      27.53         22.67
##           Std.Dev       8.13         14.67
##               Min       0.00          0.00
##                Q1      29.00          0.00
##            Median      30.00         29.00
##                Q3      32.00         32.00
##               Max      32.00         39.00
##               MAD       1.48          7.41
##               IQR       2.50         19.50
##                CV       0.30          0.65
##          Skewness      -2.55         -0.74
##       SE.Skewness       0.58          0.58
##          Kurtosis       5.80         -1.26
##           N.Valid      15.00         15.00
##                 N      15.00         15.00
##         Pct.Valid     100.00        100.00

En cuanto al largo de mazorca, se observa que la media es mayor con fertilizante orgánico (27,53) en comparación con las plantas con fertilizante tradicional (22,67) Los valores de curtosis, superiores a 0 o menores a 0, estarían indicando que las distribuciones de ambos tipos de fertilizantes no se aproximan a una distribución normal. Lo que es correcto ya que se trata de maices andinos, que tienen una gran variabilidad genetica.