Datos tomados de: https://www.ers.usda.gov/data-products/adoption-of-genetically-engineered-crops-in-the-us/

Código R:

library(tidyverse)
 # data from: https://www.ers.usda.gov/data-products/adoption-of-genetically-engineered-crops-in-the-u-s/
GEcrops <- read_csv("https://www.ers.usda.gov/webdocs/DataFiles/47649/BiotechCropsAllTables2022.csv?v=2361")
GEcrops$Value <- as.numeric(GEcrops$Value)
GEcropssummary <- GEcrops %>% group_by(Attribute,Table,Year) %>%
   summarize(mean=mean(Value, na.rm = TRUE),
             sum=sum(Value,na.rm=TRUE))
#https://stackoverflow.com/questions/21187603/replace-characters-from-a-column-of-a-data-frame-r
GEcropssummary$Attribute <- gsub("All GE varieties \\(percent of all soybeans planted\\)","soya",GEcropssummary$Attribute)
GEcropssummary$Attribute <- gsub("All GE varieties \\(percent of all corn planted\\)","maíz",GEcropssummary$Attribute)
GEcropssummary$Attribute <- gsub("All GE varieties \\(percent of all upland cotton planted\\)","algodón",GEcropssummary$Attribute)
GEcropssummary <- GEcropssummary[GEcropssummary$Attribute %in% c("maíz","soya","algodón"),]
library(hrbrthemes)
produccionOGM <- GEcropssummary %>%
  ggplot(aes(x=Year, y=mean, group=Attribute, 
             color=Attribute)) + 
  geom_line()+geom_point() +
  labs(x="Año",y="% de producción de OGMs") + 
  labs(color='Especie') + ggtitle("Producción de OGMs en Estados Unidos") +
  theme_ipsum() + theme(text = element_text(size = 20), 
                        axis.title.x = element_text(size = 20),
                        axis.title.y = element_text(size = 20), 
                        plot.title = element_text(size = 25))
ggsave("OGMs en EEUU2.svg",width = 6.5, height = 5)