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 に答える