updating tourism data

Report_book_files/Performance_indicators.qmd

@@ -110,7 +110,7 @@ Indicator 30
 
 ![diving trips](/indicator_plots/prop_diving_trips_plot_final.png){fig-alt="a plot" fig-align="center"}
 
-![Non-selective gears](/indicator_plots/Proportion of trips using non-selective gears.png){fig-alt="a plot" fig-align="center"}
+![Non-selective gears](/indicator_plots/Proportion%20of%20trips%20using%20non-selective%20gears.png){fig-alt="a plot" fig-align="center"}
 
 ## 5.5 Governance
 

indicator_processing/automated_download/cruise_visitors.R → indicator_processing/automated_download/cruise_air_visitors.R

@@ -1,11 +1,13 @@
-# Code to plot number of cruise ship visitors Puerto Rico & USVI 
+# Code to plot number of cruise ship visitors and air passengers Puerto Rico & USVI 
 
-# Last updated 6/19/2024 by Carissa Gervasi
+# Last updated 10/17/2024 by Carissa Gervasi
 
 rm(list = ls())
 dev.off()
 
 
+# Puerto Rico - cruise passengers
+
 url1<-'https://www.bde.pr.gov/BDE/PREDDOCS/I_CRUISE.XLS'
 
 library(readxl)
@@ -31,6 +33,34 @@ cruise_PR = as.numeric(df2t$V2)
 PR_cruise <- data.frame(Year = yrs_PR, Cruise_passengers = cruise_PR)
 
 
+
+# Puerto Rico - air passengers (source = PR ports authority)
+
+url1<-'https://www.bde.pr.gov/BDE/PREDDOCS/I_CARGO.XLS'
+
+library(readxl)
+library(httr)
+packageVersion("readxl")
+
+GET(url1, write_disk(tf <- tempfile(fileext = ".xls")))
+df <- read_excel(tf, 2L)
+str(df)
+
+
+df2 = df[c(52,63),]
+
+df2t = as.data.frame(t(df2))
+
+df2t = df2t[-c(1,2,12),]
+
+
+yrs_PR = as.integer(df2t$V1)
+air_PR = as.numeric(df2t$V2)
+
+# Create a data frame
+PR_air <- data.frame(Year = yrs_PR, Air_passengers = air_PR)
+
+
 ############################
 
 #USVI
@@ -78,14 +108,20 @@ table_data <- strsplit(table_text, " +")[[1]]
 years <- c(2008,2009,2010,2011,2012,2013,2014,2015,2016,2017,2018,2019,2020,2021,2022)
 
 # Extract and clean the total visitors data by removing commas
-passengers <- gsub(",", "", table_data[87:101])
-passengers <- as.numeric(passengers)
+cruise_passengers <- gsub(",", "", table_data[87:101])
+cruise_passengers <- as.numeric(cruise_passengers)
+
+# Extract and clean the total visitors data by removing commas
+air_passengers <- gsub(",", "", table_data[66:80])
+air_passengers <- as.numeric(air_passengers)
 
 # Create a data frame
-USVI_cruise <- data.frame(Year = years, Cruise_passengers = passengers)
+USVI_cruise <- data.frame(Year = years, Cruise_passengers = cruise_passengers)
+USVI_air <- data.frame(Year = years, Air_passengers = air_passengers)
 
 # Print the data frame
 print(USVI_cruise)
+print(USVI_air)
 
 
 

indicator_processing/non_automated/oceanNAICS.R

@@ -73,7 +73,7 @@ save(ind, file = "indicator_objects/oceanNAICS.RData")
 ########################################################################################
 
 
-#### ************ As of 7/5/24, can't get the API for this working. For now, will just plot the data Seann uploaded to Github. At the end of this script is some code to work with the API but it isn't yet working.
+#### ************ As of 7/5/24, can't get the API for this working. For now, will just plot the data downloaded from the data portal. At the end of this script is some code to work with the API but it isn't yet working.
 
 rm(list = ls())
 dev.off()
@@ -199,7 +199,7 @@ print(naics_df)
 ## Access the API
 
 library(devtools)
-#install_github('mikeasilva/blsAPI')
+install_github('mikeasilva/blsAPI')
 library(blsAPI)
 library(jsonlite)