R中的GMM与正常CDF的积分

时间:2017-08-06 23:00:08

标签: r integral gmm

我正在尝试运行GMM估算,其中包括正常CDF的集成,其中参数存在于函数和积分区间中。代码的要点如下:

g1 <- function(b,x){
  e <- b[1] +b[2]*x$x1
  r <- e + b[3]*x$x2
  n <- b[4]+b[5]*x$x3
  ncdf <- function(z){
    return((pnorm((log(z)-b[6]*log(e))/(b[7]/x$x4)))^n)
  }
  m1 <- x$y-integrate(ncdf,-Inf,r)/ncdf(r)
  f <- cbind(m1)
  return(f)
}

init = rep(0,7)

res<-gmm(g1,bids,init)

对于轻微的变化,我要么

Error in integrate(ncdf, -Inf, r) : 'bound' must be of length one 

Error in integrate(ncdf, -Inf, r) : a limit is NA or NaN

帮助将不胜感激。提前谢谢。

1 个答案:

答案 0 :(得分:1)

所以我可以使用Vectorize()

来解决它
g1 <- function(b,x){
  e <- b[1] +b[2]*x$x1
  r <- e + b[3]*x$x2
  n <- b[4]+b[5]*x$x3
  ncdf <- function(a,b,c,d){
    return((pnorm((a-b)/c))^d)
  }

  intn <- function(z){
    return(integrate(ncdf,-Inf,z,b=m,c=s,d=n)$value)
  }

  vint <- Vectorize(intn)

  d <- b[6]*x$x6
  s <- b[7]*x$x7

  m1 <- x$y-vint(r)/ncdf(r,m,s,n)
  f <- cbind(m1)
  return(f)
}