Greetings!
I have a longitudinal dataset consisting of 10,051 observations nested within 413 subjects. I am using melogit to estimate a multilevel logistic regression model with random intercept and random slopes. If my model looks like this:

Code:
melogit dv iv1 iv2||id: iv1 iv2, cov(unst)
How can I use the bootstrap command to generate a 95% confidence interval for the random effect coefficients [var(iv1), var(iv2), cov(iv1, iv2)]?

Any suggestions of resources or example code is greatly appreciated.