 
  I have a program which make calculate them of 4 estimator, B, BE, wB and wBE.
   
  Now I want to varied n example: "n<-seq(1,100,1)", and to make a graph of "B" and "Be" according to "n", to make a comparison.
   
  Tank you very much.
   

n<-50

a<-0.2 

 theta<-rexp(n,a)

 x<-rpois(length(theta),theta)

  t<-rexp(1,a)

   xn<-rpois(1,t)

b<-mean(x)

 

B<-(xn+1)/(1+a)   

BE<-(1+xn)*b/(1+b)

r<-b/(1+b)

wB<-1/(a*(1+a))

wBE<-wB+(2/(a^2)+5/a+4)*(b/(1+b)-1/(1+a))^2

 

 

B

BE

wB

wBE

 		
---------------------------------

	[[alternative HTML version deleted]]

