Trabajo Practico 1 y 2. Modulo 4
2024-12-03
install.packages("tidyverse")Instalamos y/o cargamos paquetes
library("tidyverse")## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.1.4 ✔ readr 2.1.5
## ✔ forcats 1.0.0 ✔ stringr 1.5.1
## ✔ ggplot2 3.5.1 ✔ tibble 3.2.1
## ✔ lubridate 1.9.3 ✔ tidyr 1.3.1
## ✔ purrr 1.0.2
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errorslibrary("readxl")2.1 LA BASE DE DATOS
Base de datos TRIGO
Función read_excel()
TRIGO <- read_excel("TRIGO.xlsx")glimpse(TRIGO)## Rows: 4,170
## Columns: 5
## $ Año <dbl> 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007…
## $ Localidad <chr> "CHA", "CHA", "CHA", "CHA", "CHA", "CHA", "CHA", "CHA", "C…
## $ Tratamiento <chr> "ConFung", "ConFung", "ConFung", "SinFung", "SinFung", "Si…
## $ Genotipo <chr> "KleinTauro", "KleinTauro", "KLEINCASTOR", "KleinTauro", "…
## $ Rendimiento <dbl> 5604.40, 4945.05, 6106.75, 4992.15, 5086.34, 6342.23, 5949…2.2 - EQUIPO DE TRABAJO
#Equipo 2- Localidad Los Altos
# Filtrar los datos para la localidad LOS ALTOS
ALT <- TRIGO %>% filter(Localidad == "ALT")2.3 - ANALISIS EXPLORATORIO CON GGPLOT2
A - Cantidad de casos de la base de datos
glimpse(TRIGO)## Rows: 4,170
## Columns: 5
## $ Año <dbl> 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007…
## $ Localidad <chr> "CHA", "CHA", "CHA", "CHA", "CHA", "CHA", "CHA", "CHA", "C…
## $ Tratamiento <chr> "ConFung", "ConFung", "ConFung", "SinFung", "SinFung", "Si…
## $ Genotipo <chr> "KleinTauro", "KleinTauro", "KLEINCASTOR", "KleinTauro", "…
## $ Rendimiento <dbl> 5604.40, 4945.05, 6106.75, 4992.15, 5086.34, 6342.23, 5949…#La Base de Datos de trigo tiene 4170 casos y 5 variables ( Año, Localidad, Tratamiento, genotipo y Rendimiento).
B - Cantidad de casos de la localidad Los Altos
glimpse(ALT)## Rows: 360
## Columns: 5
## $ Año <dbl> 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007…
## $ Localidad <chr> "ALT", "ALT", "ALT", "ALT", "ALT", "ALT", "ALT", "ALT", "A…
## $ Tratamiento <chr> "SinFung", "SinFung", "ConFung", "SinFung", "ConFung", "Si…
## $ Genotipo <chr> "BIOINTA1001", "ONIX", "ONIX", "BIOINTA1001", "ONIX", "ONI…
## $ Rendimiento <dbl> 5821.55, 5207.25, 5778.69, 6028.69, 5807.26, 5378.68, 5571…#la localidad de Los Altos tiene 360 casos
C - Distribución del rendimiento el Los Altos
ggplot(ALT, aes(Rendimiento)) +
geom_histogram()## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
ggplot(ALT, aes(Rendimiento)) +
geom_histogram(color = "yellow", fill = "darkgreen")## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
D- Variación del rendimiento promedio en la Localidad de Los Altos para los diferentes tratamientos
#Polígono de frecuencias
ggplot(ALT, aes(Rendimiento)) +
geom_freqpoly(binwidth = 100)
ggplot(ALT, aes(Rendimiento, group = Tratamiento, colour = Tratamiento)) +
geom_freqpoly(binwidth = 200)
# Tratamiento Sin fungicida: (azul) las frecuencias mas altas se
encuentran concentradas en valores altos de rendimiento, superiores a
4.000 kg/ha #Tratamiento Sin fungicida (rojo), tambien presenta valores
de frecuencias altas cercano a 4.000 Kg/ha, se observa un pico máximo
(valor más alto del conteo) con alrededor de 30 observaciones
ggplot(ALT,aes(x = Tratamiento, y = Rendimiento, fill = Tratamiento)) +
geom_boxplot(outlier.colour = "black", outlier.shape = 19, outlier.size = 2) +
labs(title = "Rendimeinto de Trigo en la Loc. de Los Altos con diferentes Niveles de Fungicida",
x = "Tratamiento",
y = "Rendimiento") +
theme_classic()
# Ambos tratamientos presentan distribuciones asimetrícas, con mayor
concentración de valores cercanos a la media de 4.000 Kg, que es uno de
los rangos promedios mas bajos. # Comparando los dos tratamientos, aun
cuando los dos presentan valores máximos de rendimiento cercanos a 4.000
kg/ha, se recomienda el tratamiento Con fungicida, por presentar algunos
valores máximos mas altos que el tratamiento Sin Fungicidas,con
rendimientos superiores a 5.000 kg/ha. Esto podria ser importante para
tener en cuenta desde algunas practicas culturales u otro manejo, para
poder conseguir mas individuos con esos valores altos de
rendimiento.
E. Agregar una ubicación ficticia de un campo dentro de su localidad usando el paquete leaflet.
Localidad Los Altos, provincia de Catamarca
library(leaflet)
leaflet() %>%
addProviderTiles(providers$Esri.WorldImagery) %>%
setView(lng = -65.5003491, lat = -28.0477257, zoom = 16) %>%
addMarkers(lng = -65.5003491, lat = -28.0477257)##F. Crear una tabla de frecuencias para la variable Rendimiento y analizar los resultados.
Tablas para variables cuantitativas continuas
Instalamos y cargamos los paquetes
## Frequencies
## TF_Rendimiento_ALT
## Type: Factor
##
## Freq % % Cum.
## ------------------------- ------ -------- --------
## [3.13e+03,3.63e+03) 15 4.19 4.19
## [3.63e+03,4.13e+03) 86 24.02 28.21
## [4.13e+03,4.63e+03) 108 30.17 58.38
## [4.63e+03,5.13e+03) 51 14.25 72.63
## [5.13e+03,5.63e+03) 34 9.50 82.12
## [5.63e+03,6.13e+03) 45 12.57 94.69
## [6.13e+03,6.63e+03] 19 5.31 100.00
## Total 358 100.00 100.00El analisis de frecuencias muestra que los rendimientos de trigo en la localidad de Los Altos, se concentra en los tres primeros intervalos, mostrando en el % acumulado el 58% de las observaciones. El intervalo (4.130, 4.630) presenta la mayor frecuencia, con 108 observaciones (30,17%), seguido por el intrvalo (3.630, 4.130) con 86 observaciones (24,02%) Los rendimientos mas altos estan dados en el lintervalo (6.130,6.630) con solo 19 observaciones (5,31%)
G. Obtener medidas de resumen para la variable Rendimiento según el tratamiento aplicado.
## Año Localidad Tratamiento Genotipo
## Min. :2007 Length:360 Length:360 Length:360
## 1st Qu.:2009 Class :character Class :character Class :character
## Median :2010 Mode :character Mode :character Mode :character
## Mean :2010
## 3rd Qu.:2011
## Max. :2011
## Rendimiento
## Min. :3134
## 1st Qu.:4100
## Median :4479
## Mean :4694
## 3rd Qu.:5254
## Max. :6972El rendimiento de trigo en la Localidad de Los Altos, muestra un valor de Media de 4.479 Kg/ha, con una Mediana de 4479 kg/ha, lo que sugiere que hay una distribución asimétrica, con los datos ligeramente inclinados hacia valores menores. El rango intercuartílico (IQR) definido entre el primer cuartil (Q1 = 4.100) y el tercer cuartil (Q3 = 5.254). El valor mínimo es 3.134 y el máximo es 6972, lo que demuestra la presencia de una brecha considerable en los valores extremos. El rendimiento de trigo en la localidad de Los Altos, presenta una distribución equilibrada y moderadamente dispersa, con un rango amplio pero una concentración significativa de los datos entre los valores de 4.100 y 6.972.
H. ¿Qué genotipo recomendaría sembrar en la localidad? justifique la respuesta.
Paquete summarytools
## `summarise()` has grouped output by 'Genotipo'. You can override using the
## `.groups` argument.## # A tibble: 58 × 7
## # Groups: Genotipo [29]
## Genotipo Tratamiento Media_Rendimiento Desviacion Min Max n
## <chr> <chr> <dbl> <dbl> <dbl> <dbl> <int>
## 1 ONIX ConFung 5964. 297. 5779. 6307. 3
## 2 KLEINCASTOR ConFung 5924. 341. 5686. 6314. 3
## 3 BIOINTA1001 ConFung 5912. 466. 5507. 6422. 3
## 4 BIOINTA1001 SinFung 5833. 190. 5650. 6029. 3
## 5 KLEINCASTOR SinFung 5748. 159. 5564. 5843. 3
## 6 ACA801 ConFung 5522. 273. 5207. 5700. 3
## 7 ONIX SinFung 5324. 101. 5207. 5386. 3
## 8 ACA801 SinFung 5293. 96.1 5214. 5400. 3
## 9 KleinLeon ConFung 5282. 965. 4457 6593. 6
## 10 LE2331 SinFung 5200. 1031. 4100 6400. 6
## # ℹ 48 more rowsAl analizar los datos de rendimiento promedio por genotipo y tratamiento, se observa en general, que los tratamientos con fungicida (ConFung) presentan rendimientos ligeramente superiores en comparación con los tratamientos sin fungicida (SinFung). En el caso del genotipo ONIX, tiene valores mayores de rendimiento Con fungicida (5964,4) que Sin fungicida (5.323,9). Lo mismo sucede con el genotipo KLEINCASTOR, el rendimiento incrementa con la aplicación del fungicida. En cambio, para el genotipo BIOINTA1001 no hay diferencias significativas si el tratamiento es Sin fungicida o Con fungicida Se recomienda sembrar, si el tratamiento es Con fungicida, los genotipos ONIX y KLEINCASTOR. Si se prefiere el tratamiento Sin fungicida se podría utilizar el genotipo BIOINTA1001.