> 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
这个该如何解决呢？谢谢！！！