How to add ICC values in esttab after mixed command

[Luigi Rombi]

Hi all, I am facing some trouble since I am not able to display icc values of mixed models with esttab command.
Unfortunately icc values are not automatically imported in scalars of mixed command (available with ereturn list) therefore I can "see" them only after having launched the estat icc command and I can also note that they are stored in r memory. In essence, this is where I am stuck:

*Store model estimates

 quietly mixed DEPVAR || L3id: ||L2id:

eststo NULLMODEL

*Compute icc's and note that they are in r

 quietly estat icc

return list

scalars:
r(icc3) = .441199262185503
r(se3) = .0953898080584049
r(icc2) = .8234534948356549
r(se2) = .0302703105433652
r(level) = 95

macros:
r(label3) : "L3id"
r(label2) : "L2id|L3id"

matrices:
r(ci3) : 1 x 2
r(ci2) : 1 x 2

*Launch esttab command

 esttab NULLMODEL, label star(* 0.05 ** 0.001 *** 0.0001) b(%9.4f) mtitles aic bic scalars(ll chi2_c p_c k_f k_r) 

------------------------------------
                              (1)   
                        NULLMODEL   
------------------------------------
DEPVAR                        
Constant                  11.5751***
                          (97.67)   
------------------------------------
lns1_1_1                            
Constant                  -0.7544***
                          (-3.93)   
------------------------------------
lns2_1_1                            
Constant                  -0.8261***
                         (-29.77)   
------------------------------------
lnsig_e                             
Constant                  -1.2123***
                        (-116.92)   
------------------------------------
Observations                 5433   
AIC                     4409.5147   
BIC                     4435.9157   
ll                      -2.20e+03   
chi2_c                  7456.0281   
p_c                        0.0000   
k_f                        1.0000   
k_r                        3.0000   
------------------------------------
t statistics in parentheses
* p<0.05, ** p<0.001, *** p<0.0001

Hence, my question is: How do I add ICC values among the scalars in the lowest part of the table?