set.seed并复制100次

Question

我use set.seed=(7)并使用replicate(n=100, rnorm(4))从标准正态分布中找到大小为4的随机样本，现在我有4行100列的数据。 如何使用复制函数运行每列的平均值（此处为100列）？ 另外，如何使用复制函数运行以下公式？

(1/(n-1)(sum((x-mean(x))^2))

Answer 1

只需将mean()添加到replicate()。

set.seed(7)

replicate(n=100, mean(rnorm(4)))

对于第二个问题，请在前面定义一个函数。

fun <- function(x, n) (1 / (n - 1) * (sum((x - mean(x))^2)))

set.seed(7)
replicate(n=100, fun(x=rnorm(4), n=100))

Answer 2

您可以这样做：

Test Success