function dedicated to add simulation standard error bound, in development draw the bound to current plot
add_sim_bound(
data,
family = gaussian(),
M = 200,
n = 1000,
true.theta = c(1, 1, 1)/sqrt(3)
)a list of simulated data
default is gaussian()
number of simulations
sample size
the true coefficients