库(dplyr)库(rsample)cv<-函数(x){返回(sd(x)/平均值(x))}sim<-data.frame(id=1:100000,变量=rpois(100000,5),组10=样品(1:10100000,替换=真),group100=样本(1:100100000,替换=真),group1000=样本(1:1000100000,替换=真),group10000=样本(1:10000000000,替换=真))AA<-空G<-c(10100100010000)for(k in 1:4){对于(1:G[k]中的j){如果(k==1){一个<-dplyr::滤波器(sim,group10==j)}否则,如果(k==2){A<-dplyr::过滤器(sim,group100==j)}否则,如果(k==3){A<-dplyr::过滤器(sim,group1000==j)}否则,如果(k==4){A<-dplyr::filter(sim,group10000==j)}AA<-rbind(AA,c(G[k],j、,nrow(A),平均值(A$变量),var(A$变量),sd(A$变量),cv(A$变量)))}}AA<-as.data.frame(AA)列名(AA)<-c(“组”、“组N”、“N”、”平均值“、”变异“、”SD“、”CV“)封头(AA)尾部(AA)AA10<-AA%>%过滤器(组==10)AA100<-AA%>%过滤器(组==100)AA1000<-AA%>%过滤器(组==1000)AA10000<-AA%>%过滤器(组==10000)(cv(sim$变量))par(mfrow=c(1,4))绘图(AA10$N,AA10$CV,ylim=c(0,1.2),ylab=“CV”,xlab=“样本量”,main=“10组”)abline(h=cv(sim$变量),col=“red”)绘图(AA100$N、 AA100型$CV,ylim=c(0,1.2),ylab=“CV”,xlab=“样本大小”,main=“100组”)abline(h=cv(sim$变量),col=“红色”)绘图(AA1000$N,AA1000$CV,ylim=c(0,1.2),ylab=“CV”,xlab=“样本量”,main=“1000-组”)abline(h=cv(sim$变量),col=“red”)绘图(AA10000$N、 10000澳大利亚元$CV,ylim=c(0,1.2),ylab=“CV”,xlab=“样本大小”,main=“10000-组”)abline(h=cv(sim$变量),col=“red”)开发关闭()