我是新来的静力学和收到以下问题,需要回答在R语言:
模拟一个i.i.d过程{Xt}t=1,···,n,样本大小n = 1000,模拟时间N = 500,服从标准正态分布Xt ∝ Normal(0,1),计算样本均值X(1),···,X(N),其中X(i)是第i次模拟的样本均值,绘制X(1),···,X(N)的直方图。
我的想法是:
1.样本量n=1000,则我应该
set.seed(1)#设置种子X1〈- rnorm(1000)#模拟X1
1.计算X1-XN的样本平均值
结果平均值〈-平均值(x1)
1.绘制平均值X1-XN的直方图
图(结果,平均值,类型=“h”)
但是,我不确定如何处理模拟时间N = 500?我生成的图只是1条直方图,所以我非常确定应该使用模拟时间。
这里模拟的目的是什么?如果我的想法在iid的情况下是正确的?谢谢
1条答案
按热度按时间7gs2gvoe1#
使用正态分布中的随机数,基本(
stats
)r代码为rnorm,默认值的平均值为0,标准差为1。我们从中得到500个样本。然后,取这1000个数字的向量的平均值。我们对replicate
重复该操作1000次,并将结果放入直方图中。