Estimate the three Probit models presented in Table 1 in Chen, Y., et al. (2016). Note, do not try to match with the exact numbers presented in Table 1. Table 1 presents the maringal effects and we will calculate them in the next question. Here, use table1_reg1, table1_reg2 and table1_reg3 to store the fitted models for each column. Feel free to collect the three columns into a regression table with Stargazer. Hints: Here, use the "join" column as the outcome variable and include all treatment dummies (treatment2, treatment3, ..., treatment8) for column (1) and (2). For column 3, you need to drop treatmtent2 in your specification; Use smf.probit() to specify the probit model and .fit() operator to estimate the model; The autograder checks the point-estimates for the three fitted models