#install.packages("tidyverse")
#install.packages("summarytools")
#install.packages("leaflet")
# Instalamos y/o cargamos paquetes
library(readxl)
library(tidyverse)
library(summarytools)
library(leaflet)
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
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
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.
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.
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()`).
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.