r语言 - 如何模拟 50 个随机样本并计算每个样本的均值和方差

Question

当我模拟 50 个正态分布的随机样本并尝试计算每个模拟的均值和方差时，我得到此错误："有 50 个或更多警告(使用 warnings(( 查看前 50 个警告("。

n=100
mean=100
sd=25 
sample=NULL
meansample=NULL
sdsample=NULL
for (i in 1:50)
  {
  sample[i]=rnorm(n,mean,sd)
  meansample[i]=mean(sample[i])
  sdsample[i]=sd(sample[i])
  }
sample
meansample
sdsample`

我想问我如何正确计算平均值和标准偏差，以及为什么在执行代码时出现此错误"有 50 个或更多警告(使用 warnings(( 查看前 50 个警告("。

Answer 1

我个人会这样做。

x <- replicate(50, rnorm(100, 100, 25), simplify = FALSE)
sapply(x, mean)
sapply(x, sd)

当前代码的问题是示例将包含 50 个单独的数据集，因此您可能希望它是一个列表，将其称为 sample[[i]]。

n=100
mean=100
sd=25 
sample=list()
meansample=NULL
sdsample=NULL
for (i in 1:50)
{
  sample[[i]]=rnorm(n,mean,sd)
  meansample[i]=mean(sample[[i]])
  sdsample[i]=sd(sample[[i]])
}
sample
meansample
sdsample