# 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)))
# overall summary
summarise_specs(results)
# Summary of specific analytical choices
summarise_specs(results, # data frame
x, y) # analytical choices
# Summary of other parameters across several analytical choices
summarise_specs(results,
subsets, controls,
var = p.value,
stats = list(median = median,
min = min,
max = max))
# Unnamed vector instead of named list passed to `stats`
summarise_specs(results,
controls,
stats = c(mean = mean,
median = median))
Run the code above in your browser using DataLab