问题如下:
我正在尝试使用R来计算上面的集成。因为我是R的新手,我不确定我的代码是否正确地计算了这个集成。如果有人能帮我解决这个R代码,我将不胜感激。我的代码如下。
x0<- c(1,3,0,0,0,3,1,6)
x1<- c(2,3,6,0,0,5,6,6)
n0<- 6
n1<- 6
N<-6
k<-8
f<- function(x0,x1,n0,n1,m,p){
u<- (p+(1-p)*(m+1)*((x0[1]+x1[1])/(n0+m+1))*(sum(sapply(0:m,function(s) {sum(sapply(0:m,function(t){((factorial(m)^2)*factorial(t+s)*factorial(x0[1]+s)/(factorial(2*m)*factorial(m+n0)*factorial(n1-t)))}))}))))
for (j in 2:k)
u<- u*(p+(1-p)*(m+1)*((x0[1]+x1[1])/(n0+m+1))*(sum(sapply(0:m, function(s) {sum(sapply(0:m,function(t){((factorial(m)^2)*factorial(t+s)*factorial(x0[1]+s)/(factorial(2*m)*factorial(m+n0)*factorial(n1-t)))}))}))))
return(u)}
int1<- function (p){sum(sapply(1:N,function(m)(1/N)*f(x0,x1,n0,n1,m,p)))}
A<- integrate(Vectorize(int1),0,1)$val
A