# NOT RUN {
# Retrieve public hourly atmospheric pressure data for the last 30 days at
# Reefton Ews station
# Subtract 30 days from today's date to get the start date
last_month = paste(as.character(Sys.Date() - 30), 0)
reefton_pressure = cf_query(cf_user(), cf_datatype(7, 1, 1), cf_station(),
start_date = last_month)
class(reefton_pressure) # cfPressure object
# Plot the atmospheric pressure data using the defaults
plot(reefton_pressure)
# Enlarge the text and add the observations as points
library(ggplot2) # for element_text() and geom_point()
plot(reefton_pressure, ggtheme = "bw", text = element_text(size = 16)) +
geom_point(size = 3, shape = 1)
# Save the plot as a png to the current working directory
library(ggplot2) # for ggsave()
ggsave("my_pressure_plot.png")
# }
Run the code above in your browser using DataLab