如何在R中更有效地计算积分？

Question

我在R中有一个功能，需要将近1秒才能运行。所以我觉得这不重。它看起来像这样：

DataList <- cons1mpdf(tradedays[i])
  tempQ <- DataList[[1]]
  Qpdf <- function(X) pmax(tempQ(X),0)
  norINTEG <- function(X)  Qpdf(X) / X^(-RRA)

然后我想从这个函数中取一个积分，并将结果用在这样的混合化问题中：

a <-  integral( norINTEG ,  DataList[[2]] ,  DataList[[3]] )

  realPDF <- function(X)  norINTEG(X) / a

  zBerkowitz[i] <-  qnorm(  pmax( 0,  integral(realPDF ,  DataList[[2]] ,   head(nextmonth.close[date == tradedays[i]] , 1)    ) ) )

但是a的计算大约花了差不多1分钟RRA大于3.所有这些代码都在for循环中，这个积分导致我的循环运行在很多时间。问题变得更糟，因为我想在最大化问题中使用此结果，这将需要几天才能运行。 这是我想解决的问题的一部分：

你认为我做错了什么？我该如何改进呢？