Dear All,

I estimate a multilevel logit model using melogit. The result is reported below:

Array

I would like to save var(centeraf), var(_cons) and cov(centeraf,_cons) such that I can use them to plot the total variance of the model using:

Code:
twoway function var(_cons) + covar(centeraf, _cons)*centeraf + var(centeraf)*centeraf^2, range (-10 10)

I think I have seen sometimes ago how I can save them (maybe using scalar) but I cannot remember where I saw it.

Since I am estimating more models, I would like to save those element after each estimation rather than imputing them manually every time.

Any help would be appreciated.

Thanks in advance,

Dario