@@ -132,18 +132,21 @@ tibble(
132132) |>
133133 rowwise() |>
134134 mutate(
135- ridge = brlavaan:::pen_ridge(x),
136- ridge_bounded = brlavaan:::pen_ridge_bound(x, lb = lb, ub = ub),
137- huber_bounded = brlavaan:::pen_huber(x, lb = lb, ub = ub),
138- huber_left_bounded = brlavaan:::pen_huber(x, lb = lb, ub = Inf)
135+ Ridge = brlavaan:::pen_ridge(x, (lb+ub)/2 ),
136+ Ridge_bounded = brlavaan:::pen_ridge_bound(x, lb = lb, ub = ub),
137+ Huber_bounded = brlavaan:::pen_huber(x, lb = lb, ub = ub),
138+ Huber_left_bounded = brlavaan:::pen_huber(x, lb = lb, ub = Inf)
139139 ) |>
140140 pivot_longer(-x) |>
141141 ggplot(aes(x, value, color = name)) +
142142 annotate("rect", xmin = lb, xmax = ub, ymin = -Inf, ymax = Inf, alpha = 0.2) +
143143 geom_line() +
144144 geom_vline(xintercept = c(lb, ub), linetype = "dashed") +
145145 labs(y = "f(x)") +
146- theme_minimal()
146+ theme_minimal() +
147+ labs(color = "Penalty") +
148+ scale_color_discrete(labels = c("Huber bounded", "Huber left bounded", "Ridge", "Ridge bounded"))
149+
147150```
148151
149152### Comparing Huber penalty thresholds
0 commit comments