#unequal probability

Steven Thompson 是我这学期的professor, 需要资料的私我
#Example 1
y_i<-c(60,60,14,1)
p_i<-c(0.05,0.05,0.02,0.01)
n=4
tau_p<-1/n*sum(y_i/p_i)
est_var_p<-1/((n-1)*n)*sum((y_i/p_i-tau_p)^2)
#Example 2
N<-100
y_i<-c(50,14,50,25)
p_i<-c(0.05,0.01,0.05,0.02)
n=4
tau_p2<-1/n*sum(y_i/p_i)
est_var_p2<-1/((n-1)*n)*sum((y_i/p_i-tau_p2)^2)
mu_p<-1/N*tau_p2
est_var_mu<-1/N^2*est_var_p2