horvitz thompson

Análisis de la varianza HT

En este código los pesos no son proporcionales a la variable respuesta y se tiene tamaño alto de población siendo el pesos mayor cerca de 0.1.

library(MASS)
N=560000;n=300
xi=rgamma(N,shape=5,rate=1/20)
fi=c(56000,50000,45000,40000,35000,30000,25000,20000,15000,10000,5000,1000,500,300,rgamma(N-14,shape=10,rate=1/10))
pi=fi/sum(fi)
pi_i=n*pi
ri=xi/pi
vtxf=(N-n)/(N)*var(ri)/n
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx1=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx2=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx3=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx4=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx5=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx6=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx7=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx8=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx9=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c1=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c2=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c3=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c4=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c5=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c6=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c7=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c8=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c9=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
y1=sample(ri,size=n,replace=FALSE,prob=pi)
y2=sample(ri,size=n,replace=FALSE,prob=pi)
y3=sample(ri,size=n,replace=FALSE,prob=pi)
y4=sample(ri,size=n,replace=FALSE,prob=pi)
y5=sample(ri,size=n,replace=FALSE,prob=pi)
y6=sample(ri,size=n,replace=FALSE,prob=pi)
y7=sample(ri,size=n,replace=FALSE,prob=pi)
y8=sample(ri,size=n,replace=FALSE,prob=pi)
y9=sample(ri,size=n,replace=FALSE,prob=pi)
y10=sample(ri,size=n,replace=FALSE,prob=pi)
c10=rbind(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10)
cx10=rbind(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10)
cxx=rbind(cx1,cx2,cx3,cx4,cx5,cx6,cx7,cx8,cx9,cx10)
T=sum(xi)
l1=apply(cxx,1,mean)-1.88*sqrt(apply(cxx,1,var)/n*(1-n/N))
l2=apply(cxx,1,mean)+1.88*sqrt(apply(cxx,1,var)/n*(1-n/N))
e0=(-abs((l1-T)*(l2-T))/((l1-T)*(l2-T))+1)/2
l1=apply(cxx,1,mean)-2.38*sqrt(apply(cxx,1,var)/n*(1-n/N))
l2=apply(cxx,1,mean)+2.38*sqrt(apply(cxx,1,var)/n*(1-n/N))
e1=(-abs((l1-T)*(l2-T))/((l1-T)*(l2-T))+1)/2
l1=apply(cxx,1,mean)-3.02*sqrt(apply(cxx,1,var)/n*(1-n/N))
l2=apply(cxx,1,mean)+3.02*sqrt(apply(cxx,1,var)/n*(1-n/N))
e2=(-abs((l1-T)*(l2-T))/((l1-T)*(l2-T))+1)/2
l1=apply(cxx,1,mean)-3.49*sqrt(apply(cxx,1,var)/n*(1-n/N))
l2=apply(cxx,1,mean)+3.49*sqrt(apply(cxx,1,var)/n*(1-n/N))
e3=(-abs((l1-T)*(l2-T))/((l1-T)*(l2-T))+1)/2
rbind(c(mean(e0),0.90),c(mean(e1),0.95),c(mean(e2),0.98),c(mean(e3),0.99))

##       [,1] [,2]
## [1,] 0.945 0.90
## [2,] 0.979 0.95
## [3,] 1.000 0.98
## [4,] 1.000 0.99