> a=-0.25;
> b=1.5;
> i=1;
> q=0.2;
> p<-array(0,dim=1);
> p[0]=0.3298969;
> EY=0;
> EY2=0;
> n=10;
> for(i in 1:10){
+ p<-p[(i-1)]*(a+(b/i));
+ EY=EY+p*i*q
+ EY2=EY2+p*i*i*q
+ i<-i+1
+ }
Error in p <- p[(i - 1)] * (a + (b/i)) : replacement has length zero
这个该如何解决呢?谢谢!!!