计算R中的积分

我计算累积分布函数，其结果应位于[0,1]。计算CDF的方程式为：

\开始{align}

F=\int{\hat{a}}}{x}\frac{2}{\hat{b}-\hat{a}}-\sum\nolimits{k=0}{N-1}C{k}-\text{cos}\bigg（\big（y-\hat{a}\big）\frac{k\pi}{b}-\hat}\bigg

\结束{align}

在哪里

• Ck

  是一个向量

• cos术语

  是一个向量
• 长度（

  ck

  ）=长度（

  cos-term

  ）=N

我确信这个等式是正确的，但恐怕我的代码是错误的

这是我的密码：

integrand<-function(x,myCk)
{
 (2/(b-a))*(t(myCk)%*%as.matrix(cos((x-hat.a)*uk)))
}

f <- function(x){integrand(x,myCk)}

# define a vectorized version of this function 
fv <- Vectorize(f,"x")

res<-integrate(fv,upper = r,lower = hat.a, subdivisions = 2000)$value

---

integrated不需要ˋtˋ和ˋas.matrixˋ。如何确保cos项是非负的？cos项可以是负的@Karsten W。那么这个被积函数怎么可能是一个密度呢？这个公式是由本文提出的，它是方程（14）@Karsten W.也许被积函数中的
a
应该是
hat.a
。否则，代码在我看来是正常的。你为什么认为计算结果是错误的？你画了被积函数的图吗？不需要ˋtˋ和ˋas.matrixˋ。如何确保cos项是非负的？cos项可以是负的@Karsten W。那么这个被积函数怎么可能是一个密度呢？这个公式是由本文提出的，它是方程（14）@Karsten W.也许被积函数中的
a
应该是
hat.a
。否则，代码在我看来是正常的。你为什么认为计算结果是错误的？你画了被积函数的图吗？