#Instalación de paquetes
{r} install.packages("tidyverse") install.packages("readxl") install.packages("summarytools") install.packages("leaflet")
{r} library("tidyverse") library("readxl") library("summarytools") library("leaflet") library("dplyr")
#2. Análisis Descriptivo #Responder a las siguientes preguntas
utilizando todo lo visto en el módulo:
#A. ¿Cuántos casos (o filas) tiene en total la base de datos?
#Función read_excel()
{r} TRIGO <-read_excel("TRIGO.xlsx")
{r} glimpse(TRIGO) #Visualización de la tabla de
datos
{r} TRIGO
##Respuesta:La base de datos tiene 4,170 observaciones
#B.Cúantos casos (o filas) tiene la base de datos de su localidad?
{r} RAF <-TRIGO %>% filter (Localidad == "RAF") RAF
##Respuesta:La Localidad Rafaela tiene 366 observaciones
#C.¿Cómo se distribuye el rendimiento de trigo en su localidad? Elabore un gráfico y saque conclusiones.
{r} ggplot(RAF, aes(Rendimiento)) + geom_histogram(color = "yellow", fill = "darkgreen")+ geom_vline(xintercept = 4361, color = "red")+ labs(title = "Rendimiento de trigo en la localidad de Rafaela, Buenos Aires", x = "Rendimiento (Kg/ha)", y = "Frecuencia")
###Para una mejor visualización del comporamiento de los datos se agrupo
en 13 columnas segun el siguiente gráfico
{r} ggplot(RAF, aes(Rendimiento)) + geom_histogram(bins = 13,color = "gray", fill = "lightblue")+ geom_vline(xintercept = 4361, color = "red" )+ geom_freqpoly(color = "red", linewidth = 0.2)+ labs(title = "Rendimiento de trigo en la localidad de Rafaela, Buenos Aires", x = "Rendimiento (Kg/ha)", y = "Frecuencia")
##Respuesta: Según lo visualizado en el Histograma, la distribución del rendimiento en la localidad Rafaela es bimodal; Un conjunto de datos se encuentra entre los 3.300 y 3.600 kg/ha y el otro conjunto de datos entre los 4.800 y 5.300 kg/ha, aproximadamente. Esta distribución de los datos puede estar indicando que existe más de un proceso en la recolecciónn de los mismo, posiblemente se manifieste el efecto del fertilizante en el rendimiento. La media del rendimiento (4361Kg) se ubica a la derecha del conjunto de datos modal. la distribución de los datos no es normal
#D.¿Cómo varía el rendimiento promedio en su localidad para los diferentes tratamientos? ¿Qué tratamiento recomendaría (Con fungicida o Sin fungicida). En todos los gráficos realizar las modificaciones que consideren adecuadas (ejes, tamaño de letra, color, etc.) para mejorar la presentación.
{r} ggplot(RAF, aes(Tratamiento, Rendimiento, color = Tratamiento)) + geom_boxplot() + stat_summary(fun = mean, color = "black", size = 0.8, shape = 4) + labs(title = "Distribución del rendimiento por tratamiento", x = "Tratamiento", y = "Rendimiento (Kg/ha)") + theme_classic() + theme(plot.title = element_text(size = 20, hjust = 0.5), axis.title.x = element_text(size = 14, color = "black"), axis.title.y = element_text(size = 16, color = "black"), axis.text.x = element_text(size = 12, color = "black"), axis.text.y = element_text(size = 12, color = "black"), legend.position = ("right"), legend.text = element_text(size = 14), legend.title = element_text(size = 16)) + scale_y_continuous(limits = c(200, 8000),breaks = seq(200, 8000, by = 1000)) + scale_color_manual(values = c("SinFung" = "blue", "ConFung" = "orange"))
##Respuesta:
##El boxplot muestra la distribución del rendimiento de trigo, en Kg/ha, para los tratamientos “Con fungicida” (ConFung) y “Sin fungicida” (SinFung) realizados en la localidad de Rafaela. ##El Rendimiento promedio para ambos tratamientos se encuentran entre los 4300 kg/ha y 4500 kg/ha, observandose una mejor performance del rendimiento en el tratamiento SinFung. Los promedios de los rendimientos, en ambos casos se encuentran muy próximos a la mediana ##El Rendimeinto en las parcelas con fungicida presenta un rango algo mas amplio con valores mínimos bajos y valores máximos respecto al tratamiento Sin Fungicida. Por otra parte, en el tratamiento sin fungicida hay una menor dispersión en los rendimientos; la mitad de los casos tiene un valor superior a la mediana y la otra mitad tiene un valor inferior. En ambos tratamientos no se observan valores atípico.
##Recomendación ##En la localidad de Rafaela, No se recomienda realizar aplicación de fertilizante para emjorar los rendimientos,expresados en Kg/ha, en trigo
#E. Agregar una ubicación ficticia de un campo dentro de su localidad
usando el paquete leaflet.
{r} leaflet() %>% addProviderTiles(providers$Esri.WorldImagery) %>% setView(lng = -61.460535, lat = -31.285141, zoom = 15) %>% addMarkers(lng = -61.460535, lat = -31.285141)
#F. Crear una tabla de frecuencias para la variable Rendimiento y analizar los resultados.
{r} TF_Rendimiento_RAF <- cut(RAF$Rendimiento, breaks = seq(min(RAF$Rendimiento), max(RAF$Rendimiento), by = 1000), right = FALSE, include.lowest = TRUE) freq(TF_Rendimiento_RAF, report.nas = FALSE, justify = "center" )
##Respuesta
##Los rendimientos de trigo en la localidad Rafaela varían entre 1350kg/ha y 7350 kg/ha, lo que muestra una amplia dispersión en los datos, aunque el 75,6% de los datos relevados del rendiemnto varia entre 3350y6350 kg/ha. Estos intervalos comprenden la mayoría de los datos, lo que sugiere que el rendimiento promedio de los genotipos de trigo se encuentra en este rango. ##Los intervalos con mayor frecuencia son: 3350-4350 kg/ha, con 119 observaciones concentrando el 32% de los datos; con 87 observaciones, el intervalo 4350-5500 kg/ha (23.85% de los datos).El 20,55% de los datos se encuentra por debajo de los 3350 kg/ha y Tan solo el 3,01% supera el rendimeito de 6350kg/ha lo que muestra que los valores más altos son poco frecuentes y pueden estar influenciado por condiciones optimas en cuanto a clima o de manejo. La distribución muestra una tendencia central en torno a 3350-4350 kg/ha, con una disminución gradual en las frecuencias hacia los extremos.
#G. Obtener medidas de resumen para la variable Rendimiento según el tratamiento aplicado.
###Cambiamos el nombre de los casos
{r} RAF <- RAF %>% mutate(Tratamiento_RAF = case_when( Tratamiento == "SinFung" ~ "Sin Funguicida", Tratamiento == "ConFung" ~ "Con Funguicida")) RAF
##Medidas de resumen por tratameinto para la localidad Rafaela
{r} RAF %>% group_by(Tratamiento_RAF) %>% descr(Rendimiento, headings = FALSE, justify = "center")
#H. ¿Qué genotipo recomendaría sembrar en la localidad? justifique la respuesta
{r} Rendimiento_promedio <- RAF %>% group_by(Genotipo, Tratamiento) %>% summarise(Rendimiento_promedio = mean(Rendimiento, na.rm = TRUE)) %>% arrange(desc(Rendimiento_promedio)) # Ordenar de mayor a menor rendimiento
{r} Rendimiento_promedio
#Seleccionar el genotipo de mayor rendimiento
{r} Mejor_genotipo <- Rendimiento_promedio %>% group_by(Tratamiento) %>% slice(1) # Tomar el genotipo con mayor rendimiento por tratamiento Mejor_genotipo
##Respuesta ##El genotipo que recomendaría es ACA905 ya que tiene un rendimeinto del 5543 kg/ha. Esta recomenación se sustenta en que el genotipo ONIX, el segundo genotipo con mayor rendimiento, tratado con fungincida tan solo supera a ACA905 en 300kg/ha. Si se considera los costos que implica la plicación del fungicida, debería ser muy superior el rendimiento como para justificar su utilización en la localidad de Rafaela.
{r setup, include=FALSE} knitr::opts_chunk$set(echo = TRUE)
This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.
When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:
{r cars} summary(cars)
You can also embed plots, for example:
{r pressure, echo=FALSE} plot(pressure)
Note that the echo = FALSE parameter was added to the
code chunk to prevent printing of the R code that generated the
plot.