Account for clustered standard errors in mediation analysis in R

Question

I have an experimental set-up that requires that I account for the clustered nature of the data. I currently do this using lm_robust in the estimatr li