Plot change in a single coordinate between two histories

plot_coordinate_change(
  mu_full_history,
  u_final,
  coord,
  history_i,
  history_j,
  type = c("both", "overlay", "diff"),
  highlight = c("both", "i", "j", "none"),
  highlight_u = NULL,
  res = NULL,
  main = NULL,
  xlab = "u (grid position)",
  ylab_overlay = expression(mu[j](u)),
  ylab_diff = expression(Delta * mu[j](u)),
  lty_i = 1,
  lty_j = 2,
  lty_diff = 1,
  pch_i = 16,
  pch_j = 17,
  pch_diff = 16,
  cex_pts = 0.9,
  add_zero_line = TRUE,
  legend_loc = "topright",
  bty = "n",
  ylim_overlay = NULL,
  pad_ylim = 0.04,
  include_highlight_in_ylim = TRUE,
  label_points = FALSE
)