I am looking for a way to compute average marginal effects with clustered standard errors which i seem to be having a few problems with. My model is as follows: