请问，r语言怎么做100次系统抽样，想要求系统抽样总体均值简单估计量的循环平均方差，

请问，r语言怎么做100次系统抽样，想要求系统抽样总体均值简单估计量的循环平均方差

# 生成一个长度为100的随机样本
sample <- rnorm(100)

# 定义一个函数，用于计算系统抽样总体均值简单估计量的循环平均方差
systematic_sampling_var <- function(n, k) {
  # n为总体大小，k为每隔k个元素抽取一个样本
  # 生成一个长度为n的数列
  population <- 1:n
  # 计算总体均值
  mu <- mean(population)
  # 初始化循环平均方差
  cmv <- 0
  # 进行100次系统抽样
  for (i in 1:100) {
    # 计算第i次抽样的起始位置
    start <- ((i-1)*k) %% n + 1
    # 从起始位置开始每隔k个元素抽取一个样本
    sample <- population[start:(start+k-1) %% n + 1]
    # 计算样本均值
    xbar <- mean(sample)
    # 更新循环平均方差
    cmv <- cmv + (xbar - mu)^2
  }
  # 返回循环平均方差
  return(cmv/100)
}

# 进行100次系统抽样，每隔10个元素抽取一个样本，并计算循环平均方差
cmv <- replicate(100, systematic_sampling_var(1000, 10))

# 输出循环平均方差的平均值和标准差
cat("循环平均方差的平均值为", mean(cmv), "\n")
cat("循环平均方差的标准差为", sd(cmv), "\n")