Example drake plan

plan.R

the_plan <-
  drake_plan(
    start_date = as_datetime("2018-01-01"),
    end_date = as_datetime("2020-01-01") - seconds(1),

    permanent_limit = 12 * 60,
    auxiliary_limit = 7 * 60,

    percentile = 0.9,

    # spatial block is H3
    h3_resolution = 6,
    hexes = get_hexes(h3_resolution),

    remoteness_areas = get_remoteness_areas(.filter = list(ste_code16 = 3)),

    tagged_hexes = tag_hexes(hexes, remoteness_areas),

    stations = get_stations(),

    crews = get_crewnum_19(),

    hexes_assigned_stations = read_rds(
      file_in("//pims/dependencies/hex_assigned_stations.rds")
    ),

    hexes_assigned_crew = assign_crew(hexes_assigned_stations,
                                      crews),

    h3_geoscience_res6_source = fetch_geoscience_data(),

    h3_geoscience_res6 = read_geoscience_data(h3_geoscience_res6_source),

    qld_sa1s = get_qld_sa1s(),

    future_pop = get_population_data(),

    hex_populations = intersect_population_hexes(future_pop,
                                                 qld_sa1s,
                                                 hexes),

    oms_responses = get_oms_responses(start_date, end_date),

    # FIXME: hexify the remoteness areas / create remoteness area clusters from hexes?
    # this would enforce each incident falls into the same remoteness area and hex bin
    tagged_oms_responses = tag_responses(oms_responses, hexes),

    hex_response_times = summarise_hex_response_times(
      tagged_oms_responses,
      tagged_hexes,
      percentile
    ),

    ## Structure Fires
    hex_summary = summarise_hexes(
      hex_response_times,
      h3_geoscience_res6,
      tagged_oms_responses,
      hexes_assigned_crew,
      hex_populations
    ),

    hex_predictions = knn_prediction(hex_summary),


    report = target(
      command = {
        rmarkdown::render(knitr_in("doc/analysis.Rmd"))
        file_out("doc/analysis.html")
      }
    )
  )