# load additional library
library(ggplot2) # for further customization of the plots
# Run specification curve analysis
results <- run_specs(df = example_data,
y = c("y1", "y2"),
x = c("x1", "x2"),
model = c("lm"),
controls = c("c1", "c2"),
subsets = list(group1 = unique(example_data$group1),
group2 = unique(example_data$group2)))
# Plot simple specification curve
plot_curve(results)
# Ribbon instead of CIs and customize further
plot_curve(results, ci = FALSE, ribbon = TRUE) +
geom_hline(yintercept = 0) +
geom_hline(yintercept = median(results$estimate),
linetype = "dashed") +
theme_linedraw()
Run the code above in your browser using DataLab