создать матрицу образцов в R

У меня есть распределение вероятностей X, и я хотел бы создать выборки из 100 наблюдений: я использую sample(X,size=100,replace=TRUE). Я хотел бы построить среднее значение PDF выборки для 100 1000 10000 выборок, поэтому я попытался создать матрицы наблюдений, используя matrix(sample(X,size=100,replace=TRUE),nrow=100,ncol=100), но это сгенерировало бы один и тот же образец во всех столбцах. Любые идеи о том, как создать новый образец для каждого столбца?


person KJohn    schedule 10.10.2014    source источник
comment
Или просто сэмплируйте 1e4 и используйте dim, например temp <- sample(X,size=1e4,replace=TRUE) ; dim(temp) <- c(100, 100)   -  person David Arenburg    schedule 10.10.2014


Ответы (1)


как насчет этого? Замените 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
person DMT    schedule 10.10.2014