0

矩条件函数就是 exp(-g/r)-1,其中 g 是 AAA 级债券指数的已知日收益率序列,r 是通过 gmm 导出的风险度量。我的代码如下:

View(Source)
library(gmm)
data(Source)
x <- Source[1:5200,"AAA"]
m <- function(r,x)
 {m.1 <- exp(-x[,"AAA"]/r)-1}
summary(gmm(m,x,t0=1,method="BFGS",control=1e-12))

这会产生以下错误消息:

****model.frame.default 中的错误(公式 = gmat ~ 1,drop.unused.levels = TRUE):变量“gmat”的类型(列表)无效****

谁能帮我弄清楚出了什么问题?非常感谢!

对于那些想要复制结果的好心人,请找到上面提到的源数据。正确的 r 为 1.590 ,可以通过 excel 中的目标搜索来解决,目标函数为 :(average(exp(-g/r)-1) )^2 ,目标值:0(公差:1e-12)

https://docs.google.com/spreadsheets/d/1AnTERQd2jm9ttKDZa7On3DLzEZUWaz5Km3nKaB7K18o/edit?usp=sharing

4

0 回答 0