Dear Statalist,

As a robustness check for my main model, I attempt to run a set of quantile regressions for an instrumental variable model with the fixed effect. The instrumental variable model includes both birth year and birth state fixed effect. The standard errors are clustered at the birth state level.

In my main analysis, the model is defined by the following code:
Code:
ivreg2 cognition (infl_pneumonia_rate1000=inter) female agesurvey black postsulfa `postAtheen' i.birth_year i.birth_state if  eduqt == 5, cluster(birth_state) partial(i.birth_year i.birth_state)
Currently, I am trying to use IVQREG2 to do the IV quantile regression with the fixed effect. Clustered SE was done using bootstrap. The code I used to estimate the IV quantile regression with FE is:
Code:
bootstrap, cluster(birth_state) idcluster(birth_state2) seed(22): ivqreg2 cognition infl_pneumonia_rate1000 female agesurvey black postsulfa `postAtheen' by* bs*, inst(inter female agesurvey black postsulfa `postAtheen' by* bs*) q(0.1(0.1).9)
The problem here is that running IVQREG2 under bootstrap is very time-consuming. For example, my STATA has run two hours to execute the code above but is still processing. Are there any better ways to estimate an IV quantile regression with FE but are more time-efficient?

Thanks in advance!