楼主您没给一些程序指令什么的，我不好回。
总之， r(S) 是以向量矩阵储存的，所以您要用矩阵的指令调用，
然后，由于aic与bic都是一个纯量，它是一个值，请从已存向量矩阵调用，
如果您须要在变量看到这个值，您可以利用指令 g 变量名=纯量名
但这个的缺点是，您的模型变化的个数多，变量也跟著多，占用记忆体也较多。

以下是利用Stata的例子，假设三个不同lag期的arch模型的aic与bic调用。
webuse wpi1
forv i=1/3{
arch D.ln_wpi, arch(1/`i')
set mo off
estat ic
mat m`i'=r(S)
sca aic`i'=m`i'[1,5]
sca bic`i'=m`i'[1,6]
g aic`i'=aic`i'
g bic`i'=bic`i'  
}
  

哦哦！！非常受用！非常感谢！