我正在尝试使用包使用 R 中的广义矩量法来估计模型,gmm
但出现以下错误:
ar.ols 中的错误(x,aic = aic,order.max = order.max,na.action = na.action,:'order.max' 必须 <'n.used'
我为您提供了我使用的代码以及我生成的样本数据,而不是这些数字随机数,我将拥有与财务回报相关的数量。
x = matrix(rnorm(248), nrow = 62, ncol = 4)
g0 = function(tet, x){
m1 = mean(tet[1]*(x[,1]/x[,2])^(-tet[2])*(x[,3]-1))
m2 = mean(tet[1]*(x[,1]/x[,2])^(-tet[2])*(x[,4]-x[,3]))
g = cbind(m1, m2)
return(g)
}
gmm(g0, x, c(beta = 1, gamma = 1))
你能告诉我我的代码有什么问题吗?