应用gmm方法时的无效类型(列表)消息

时间:2017-03-19 15:58:52

标签: r gmm

时刻条件函数只是exp(-g / r)-1,其中g是AAA类债券指数的已知系列日收益率,r是通过gmm得出的rikiness度量。我的代码如下:

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中的错误(formula = gmat~1,drop.unused.levels = TRUE):     变量'gmat'****

的类型(列表)无效

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

对于那些想要复制结果的人,请查看如上所述的源数据。 正确的r是1.590,可以通过excel中的目标搜索来解决,目标函数为:(平均值(exp(-g / r)-1))^ 2和目标值:0(容差:1e-12)

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

0 个答案:

没有答案