plot_fitted: Plot the Bayesian curve fitted values
In drkowal/fosr: Function-on-Scalars Regression for High Dimensional Data

Plot the curve posterior means with posterior credible intervals (pointwise and joint), the observed data, and true curves (if known)

plot_fitted(
  y,
  mu,
  postY,
  y_true = NULL,
  t01 = NULL,
  include_joint_bands = FALSE
)

`y`	the `n x 1` vector of time series observations
`mu`	the `n x 1` vector of fitted values, i.e., posterior expectation of the mean
`postY`	the `nsims x n` matrix of posterior draws from which to compute intervals
`y_true`	the `n x 1` vector of points along the true curve
`t01`	the observation points; if NULL, assume `n` equally spaced points from 0 to 1
`include_joint_bands`	logical; if TRUE, compute simultaneous credible bands