望求教积分相关问题。我想要求函数G的积分,积分结果报错是什么原因?
(其中n=15;y(j)表示y按照大小排序,即代码中的YO;yita=1或0是一个示性函数)
G<-function(y){
f=1
for(i in 1:15){
if(YO[i,2]<=y && YO[i,3]==1){
f=f*((15-i)/(15-i+1))
return(1-f) }
}
}
> integrate (G,0,0.1)
Error in integrate(G, 0, 0.1) :
evaluation of function gave a result of wrong length