if (FALSE) {
K <- 50
df <- data.frame(
estimator = c("Estimator 1", "Estimator 2"),
parameter = rep(c("mu", "sigma"), each = K),
truth = 1:(2*K),
estimate = 1:(2*K) + rnorm(4*K)
)
estimator_labels <- c("Estimator 1" = expression(hat(theta)[1]("·")),
"Estimator 2" = expression(hat(theta)[2]("·")))
parameter_labels <- c("mu" = expression(mu), "sigma" = expression(sigma))
plotestimates(df, parameter_labels = parameter_labels, estimator_labels)}
Run the code above in your browser using DataLab