积分最大值

Question

我在参数a，b，p，q上具有一个积分。 我对数据x的概率密度函数是

GB2PdfCalc=function(x){
  out=(abs(a)*x^(a*p-1))/(b^(a*p)*beta(p,q)*(1+(x/b)^(a))^(p+q))
  return(out)
}

我想使用GB2PdfCalc对此进行积分：

    pdfIntegration<-function(x){
      [...]
      integral[i]=integrate(GB2PdfCalc,x[k],x[j])[1]
      integral=as.vector(integral,mode="numeric")                       
      return(integral)
    }

pdfIntegration(x)

但是，我想根据a，b，p，q最大化此积分。最终，我想计算一个最大似然函数。使用mle（pdfIntegration）或optim（）对我没有帮助。 有人有主意吗？