I'm running a multilevel mixed model with:
eststo: xtmixed y x1 x2 || country: || household:, reml
eststo: xtmixed y x1 x2 x3 x4|| country: || household:, reml

I want to obtain the estimates of variances in all the levels recorded in the esttab
This is the code I'm using now but it does not give me all the variances (from the random effects parameters), both at level 3, level2, and level 1
esttab using file.rtf,se bic aic eqlabels("" "sd(country)" "sd(_cons)" "sd(Residual)", none) title(Regression Output) nonumbers mtitles("Model 1" "Model 2") addnote("Source: Here")