Ingresando la base de datos de excel

library(readxl)
datos2022 <- read_excel("../baseDatos/datos2022.xls", 
    sheet = "Hoja1", col_types = c("date", 
        "text", "numeric", "numeric", "numeric", 
        "numeric", "numeric", "numeric"))
datos2022$fecha = as.Date(datos2022$fecha)

#View(datos2022)
library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.0.5
library(dplyr)
## Warning: package 'dplyr' was built under R version 4.0.5
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(knitr)
## Warning: package 'knitr' was built under R version 4.0.5
library(reshape2)

adultosCu <- ggplot(datos2022, aes(x=fecha, y=adultoCuadros)) +
   ggtitle ("Incidencia de Thrips Adultos por Cuadro") +
  geom_line(size=0.5) +
  geom_line( color="#69b3a2") +
  xlab("Semanas por año ") +
  ylab(" % Adultos por Cuadro")+
  scale_x_date(date_minor_breaks = "1 week",
               date_breaks = '5 week',
               date_labels = '%Y/%W')+
  theme_bw()+
  theme(axis.text.x = element_text(angle=60, hjust=1),
        panel.grid.minor.x = element_line(color=gray(0.7)),
        panel.grid.major.y = element_line(color = gray(0.8)))
adultosCu

Porcentaje de incidencia de Thips adultos en camas

adultosCa <- ggplot(datos2022, aes(x=fecha, y=adultoCamas)) +
   ggtitle ("Incidencia de Thrips Adultos por Cama") +
  geom_line(size=0.5) +
  geom_line( color="#006633") +
  xlab("Semanas por año ") +
  ylab("% Adultos por Cama")+
  scale_x_date(date_minor_breaks = "1 week",
               date_breaks = '5 week',
               date_labels = '%Y/%W')+
   theme_bw()+
  theme(axis.text.x = element_text(angle=60, hjust=1),
        panel.grid.minor.x = element_line(color=gray(0.7)),
        panel.grid.major.y = element_line(color = gray(0.8)))
        

adultosCa

Porcentaje de incidencia de Thips larvas en cuadros

LarvasCu <- ggplot(datos2022, aes(x=fecha, y=larvaCuadros)) +
   ggtitle ("Incidencia de Thrips Larvas por Cuadros") +
  geom_line(size=0.5) +
  geom_line( color="#FD6467") +
  xlab("Semanas por año ") +
  ylab(" % Larvas por cuadro")+
  scale_x_date(date_minor_breaks = "1 week",
               date_breaks = '5 week',
               date_labels = '%Y/%W')+
   theme_bw()+
  theme(axis.text.x = element_text(angle=60, hjust=1),
        panel.grid.minor.x = element_line(color=gray(0.7)),
        panel.grid.major.y = element_line(color = gray(0.8)))
        

LarvasCu

Porcentaje de incidencia de Thips Larvas en Camas

larvasCa <- ggplot(datos2022, aes(x=fecha, y= larvaCamas)) +
   ggtitle ("Incidencia de Thrips Larvas por Cama") +
  geom_line(size=0.5) +
  geom_line( color="#F1BB7B") +
  xlab("Semanas por año ") +
  ylab("% Larvas por Cama")+
  scale_x_date(date_minor_breaks = "1 week",
               date_breaks = '5 week',
               date_labels = '%Y/%W')+
   theme_bw()+
  theme(axis.text.x = element_text(angle=60, hjust=1),
        panel.grid.minor.x = element_line(color=gray(0.7)),
        panel.grid.major.y = element_line(color = gray(0.8)))
        

larvasCa

Porcentaje de incidencia de Virus en Cuadros

virusCu <- ggplot(datos2022, aes(x=fecha, y= virusCuadros)) +
   ggtitle ("Incidencia de Virus por Cuadros") +
  geom_line(size=0.5) +
  geom_line( color="#C93312") +
  xlab("Semanas por año ") +
  ylab("% virus por cuadros")+
  scale_x_date(date_minor_breaks = "1 week",
               date_breaks = '5 week',
               date_labels = '%Y/%W')+
   theme_bw()+
  theme(axis.text.x = element_text(angle=60, hjust=1),
        panel.grid.minor.x = element_line(color=gray(0.7)),
        panel.grid.major.y = element_line(color = gray(0.8)))
        

virusCu

Porcentaje de incidencia de Virus en Camas

virusCa <- ggplot(datos2022, aes(x=fecha)) +
  geom_line(aes(y=virusCamas),
            size=0.5, color= "#FF0000") +
  ggtitle ("Incidencia de virus por Cama") +
  xlab("Semanas por año ") +
  ylab("% virus por Cama")+
  scale_x_date(date_minor_breaks = "1 week",
               date_breaks = '5 week',
               date_labels = '%Y/%W')+
   theme_bw()+
  theme(axis.text.x = element_text(angle=60, hjust=1),
        panel.grid.minor.x = element_line(color=gray(0.7)),
        panel.grid.major.y = element_line(color = gray(0.8)))
        

virusCa

datos2022 %>% 
  select(fecha, adultoCamas, larvaCamas, virusCamas) %>% 
  melt('fecha') %>% 
  ggplot() +
  aes(x = fecha, y = value, fill=variable)+
  geom_area(alpha=0.9, color='darkgrey')+
  ggtitle ("Incidencia de Thrips por Cama") +
  xlab("Semanas por año ") +
  ylab("% de Trhips")+
  theme_bw()

datos2022 %>% 
  select(fecha, adultoCuadros, larvaCuadros, virusCuadros) %>% 
  melt('fecha') %>% 
  ggplot() +
  aes(x = fecha, y = value, fill=variable)+
  geom_area(alpha=0.9, color='darkgrey')+
  ggtitle ("Incidencia de Thrips por Cuadros") +
  xlab("Semanas por año ") +
  ylab("% de Trhips")+
  theme_bw()

Descomposicion de datos

Adultos Cama

# Los datos se vuelven una serie de tiempo 
#start = es desde que año y mes se inicia
#frequency = es cuantos datos se tomaron en el año 
adultosCa.ts= ts(datos2022$adultoCamas, start = c(2020,01), frequency = 53)
#plot de la serie creada
plot.ts(adultosCa.ts)

#se utiliza la funcion decompose para sacar los componentes de una serie
adultosCa.componentes = decompose(adultosCa.ts)
#se grafica 
plot(adultosCa.componentes)

Adultos Cuadro

adultosCu.ts= ts(datos2022$adultoCuadros, start = c(2020,01), frequency = 53)
#plot de la serie creada
plot.ts(adultosCu.ts)

#se utiliza la funcion decompose para sacar los componentes de una serie
adultosCu.componentes = decompose(adultosCu.ts)
#se grafica 
plot(adultosCu.componentes)

Larvas Cama

larvaCa.ts= ts(datos2022$larvaCamas, start = c(2020,01), frequency = 53)
plot.ts(larvaCa.ts)

larvaCa.componentes = decompose(larvaCa.ts)
plot(larvaCa.componentes)

Larvas Cuadro

larvaCu.ts= ts(datos2022$larvaCuadros, start = c(2020,01), frequency = 53)
plot.ts(larvaCu.ts)

larvaCu.componentes = decompose(larvaCu.ts)
plot(larvaCu.componentes)

Virus Cama

virusCa.ts= ts(datos2022$virusCamas, start = c(2020,01), frequency = 53)
plot.ts(virusCa.ts)

virusCa.componentes = decompose(virusCa.ts)
plot(virusCa.componentes)

Virus Cuadro

virusCu.ts= ts(datos2022$virusCuadros, start = c(2020,01), frequency = 53)
plot.ts(virusCu.ts)

virusCu.componentes = decompose(virusCu.ts)
plot(virusCu.componentes)