pik=inclusionprobabilities(Tot,n)#泊松抽样函数包计算入样概率
m=UPpoisson(pik)  
C = getdata(l,m)
g=getdata(l$ID,m)
M=C[,-(1:3)]
M=as.matrix(M)
C1=t(M)
C1=as.matrix(C1)
R=data.frame(ID,C1)
W=getdata(R,m)
W1=W[,-(1:2)]
W2=as.matrix(W1)
library(MASS)
V=ginv(W2)
#计算抽样误差
K=C1%*%V%*%M #重构后的矩阵
E=rr1-K
Y=sum(E*E)/sum(rr1*rr1)
大家好，这是一段抽样代码，Y表示抽样误差，我现在想加上一个循坏，比如连续抽20次，最后可以得到抽20次误差的平均值，请问可以怎么加循环。