R/monthly_bicycle_rental.R

In zeee1/Tashu_Rpackage: Analysis and Prediction of Bicycle Rental Amount

#' Visualize the change of bicycle rental amount by temperature and each month.
#'
#' A function drawing a plot that shows change of temperature and bicycle rental ratio in each month.
#'
#' @export
#' @importFrom lubridate month
#' @importFrom dplyr summarise group_by %>%
#' @importFrom graphics par plot axis mtext legend
#' @examples
#' \dontrun{monthly_bicycle_rental()}

monthly_bicycle_rental <- function() {
  check_data()

    # Compute monthly rental ratio
  tashu_record <- tashudata::tashu
  tashu_record$month <- month(tashu_record$RENT_DATE)
  rental_by_month <- tashu_record %>% dplyr::count(month) %>% dplyr::rename(rentcount = n)
  rental_by_month$ratio <- rental_by_month$rentcount/sum(rental_by_month$rentcount) * 100

    # Compute monthly average temperature
    tashu_weather <- tashudata::weather
    tashu_weather$month <- month(tashu_weather$Datetime)

    temperature_by_month <- tashu_weather %>% group_by(month) %>% dplyr::summarise(temperature = mean(Temperature),.groups = 'drop')


    option <- par(mar = c(5, 5, 5, 5))
    on.exit(par(option))

    plot(temperature_by_month$month,
         temperature_by_month$temperature,
         axes = FALSE,
         ylim = c(0, 30),
         pch = 16,
         type = "b",
         xlab = "",
         ylab = "",
         col = "blue",
        main = "The monthly average temperature and ratio of rental")
    # Draw month axis
    axis(1, temperature_by_month$month)
    mtext("Month", side = 1, col = "black", line = 2.5)

    # Draw temperature axis
    axis(2, ylim = c(0, 30), las = 1)
    mtext("Average temperature(Celsius)", side = 2, line = 2.5)

    # Allow a second plot on the same graph
    par(new = TRUE)

    plot(rental_by_month$month,
         rental_by_month$ratio,
         xlab = "", ylab = "",
         ylim = c(0, 15),
         pch = 18,
         axes = FALSE,
         type = "b",
         col = "red"
         )

    # Draw rental ratio axis
    axis(4, ylim = c(0, 15), col = "black", las = 1)
    mtext("Ratio of bicycle Rental(%)", side = 4, line = 2.5)

    # Draw legend information
    legend("topright", cex = 0.8,
           legend = c("Average Temperature", "Ratio of bicycle Rental"),
           text.col = c("blue", "red"),
           pch = c(16, 15),
           col = c("blue","red"))


}