stata_challenge

stata_challenge.R

# How to make this plot in Stata?

library(tidyverse)
library(estimatr)
library(randomizr)

# make some fake data
dat <-  tibble(Y = rnorm(100), 
               Z = complete_ra(100, conditions = c("Control", "Treament")))

# obtain group means and confidence intervals
summary_df <-
  dat %>% group_by(Z) %>%
  do(tidy(lm_robust(Y ~ 1, data = .))) %>%
  mutate(Y = estimate)
  
# plot
ggplot(summary_df, aes(Z, Y)) +
  geom_point(size = 3) +
  geom_errorbar(aes(ymin = conf.low, ymax = conf.high), width = 0) +
  geom_point(data = dat, position = position_jitter(width = 0.1), alpha = 0.3) +
  theme_bw()