(其中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


雷达卡





京公网安备 11010802022788号







