Calculate Probit Estimates for samples via Loop

Question

I have a sample dataset that looks like this:

library(dplyr)
test_df <- structure(list(test = c(0L, 1L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 1L, 
0L, 1L, 0L, 1L, 0L