data(kfm, package = "ISwR")
require(sjlabelled, quietly = TRUE)
data(energy, package = "ISwR")
energy <- energy |>
var_labels(
expend = "Energy expenditure (MJ/day)",
stature = "Stature"
)
p <- energy |>
strip_error(x = stature, y = expend)
gg_star(p, x1 = 1, y1 = 13, x2 = 2, y2 = 13.2, y3 = 13.4, "**") +
ylim(0, 13.6)
Run the code above in your browser using DataLab