用R语言实现以下代码 在R中编写总体为25个单元，排成五行五列的形式用简单随机抽样方法计算总体单元的均值估计与方差估计。

以下是R语言代码实现：

# 创建总体为25个单元的数据框
pop <- data.frame(matrix(rnorm(25), nrow=5))

# 计算每个单元的均值
pop_means <- apply(pop, 1, mean)

# 简单随机抽样n个样本
n <- 5
sample_index <- sample(1:nrow(pop), n, replace=FALSE)
sample <- pop[sample_index, ]

# 计算简单随机抽样样本的均值
sample_means <- apply(sample, 1, mean)

# 计算简单随机抽样样本的方差
sample_vars <- apply(sample, 1, var)

# 计算总体单元的均值估计
mean_est <- mean(sample_means)

# 计算总体单元的方差估计
var_est <- ((n-1)/n) * sum(sample_vars) + (n^2-n)/(n*(n-1)) * var(sample_means)

# 输出均值估计与方差估计
cat("总体单元的均值估计为：", mean_est, "\n")
cat("总体单元的方差估计为：", var_est, "\n")

其中，`pop`是总体数据框，`pop_means`是每个单元的均值，`n`是简单随机抽样的样本大小，`sample_index`是抽样得到的样本的索引，`sample`是简单随机抽样得到的样本，`sample_means`是样本均值，`sample_vars`是样本方差，`mean_est`是总体单元的均值估计，`var_est`是总体单元的方差估计。