r - 在 R 中使用“复制”时尝试存储估计值

Question

我是 RStudio 的新手，有一个问题希望有人能帮助我。我正在使用“复制”函数来模拟 log-GARCH 模型，我得到的是 100 次复制，我想要做的是存储估计值以便我可以计算平均值。我怎样才能做到这一点？

代码：

library(lgarch)


replicate(n=100,{ x <- lgarchSim(500, constant=0.5) mymod <- lgarch(x) }, simplify=FALSE )

这产生了以下输出：

输出

附件是复制 [99] 和 [100] 的图像，我本质上想要的是将“intercept”、“arch1”和“garch1”存储在一个列表中。

Answer 1

结果是给你一个list模型。您可以使用从每个模型lapply中提取，将它们一起放入一个矩阵并获取以获得您的平均系数：coefficientsrbindcolMeans

library(lgarch)

my_list <- replicate(n = 100, { 
  x <- lgarchSim(500, constant = 0.5)
  mymod <- lgarch(x) 
  }, simplify = FALSE )

colMeans(do.call(rbind, lapply(my_list, coefficients)))
#>   intercept       arch1      garch1       Elnz2 
#>  1.15879953  0.04965793  0.82525587 -1.26277188