1

我有一个概率分布 X，我想创建 100 个观察的样本：我sample(X,size=100,replace=TRUE)想在 100,1000,10000 个样本上绘制样本均值 PDF，所以我尝试使用创建观察矩阵，matrix(sample(X,size=100,replace=TRUE),nrow=100,ncol=100)但它会生成相同的所有列中的样本。关于如何为每列创建新样本的任何想法？

4

1 回答 1

3

这个怎么样？用您的示例调用替换 rnorm。这将为每一列取一个新样本

replicate(3,rnorm(10))

#                          [,1]                     [,2]                    [,3]
# [1,] -0.439366440511456290974  0.349113310500896667499  2.10467702915785226381
# [2,]  0.788892611945899879800  0.572377925929974273878  0.92566383997665424577
# [3,]  0.098359807623723205516 -0.642162545019581476602  0.28636140673186011307
# [4,] -3.063133170307587249681  1.322694510750672014510  0.66340500173312999532
# [5,]  0.255018412772398617161  1.492588176987205361712  1.11444057062233659039
# [6,] -1.069621910039232570711 -1.460604130070508821504 -0.81534768620081377044
# [7,] -1.036421328330551894226  1.525817374339748067058  0.47070620500783272311
# [8,] -0.139135286049327872027 -0.065015174557339946992  0.21483758566831215320
# [9,] -0.370005496738202488416  1.573987068922320986530 -1.21431499328084857581
#[10,] -0.070508137614489943545  1.657541962601124518883  0.45886687983031809734
于 2014-10-10T14:21:54.963 回答