以正态分布复制样本

Question

将种子设置为7。执行以下计算100次。 从标准正态分布中找到大小为4的随机样本。从样本数据的平均值中找到平均平方偏差（使用公式σ^ 2作为enter image description here 然后找到这100个差异估计的平均值。

对于σ^ 2，我的代码是(1/n) * sum((x-mean(xi))^2)。 首先，我不知道我的σ^ 2代码是否有意义。

然后我将代码运行到这个问题

set.seed(7)
xi <- rnorm(4)
mean(xi)
sd <- (1 / 4) * sum( (xi - mean(xi) )^2 )
q2 <- replicate(n=100, expr = {x=rnorm(4, mean(xi), sd)})
mean(q2)

但是我找不到正确的答案。