Consider returning median output type from linear_pool()

[eahowerton]

In the hubEnsembles manuscript, we have a plot of quantile average and linear pool ensembles from the incident example data (Figure 5), where the median values are used as the point estimate. However, linear_pool() doesn't return results for median output type (I think), so I ended up having to add this clunky few lines of code, where I pull output_type_id == 0.5, rename the output_type to "median" and append it to the end of the linear_pool results.

model_output_plot <- linear_pool_ens |>
  dplyr::filter(output_type_id == 0.5, target == "wk inc flu hosp") |>
  dplyr::mutate(output_type = "median", output_type_id = NA)
model_output_plot <- dplyr::bind_rows(linear_pool_ens, model_output_plot)

[eahowerton]

Noting this is related to hubverse-org/hubData#26

[eahowerton]

New functions in hubUtils will be useful: hubverse-org/hubUtils#150