我正在尝试使用重复循环来生成100个样本大小为n = 100的泊松分布的数据集,我想按行和列排列结果,但这只是向我重复以显示最后一个数据集,而不是所有数据集。同时,我还将尝试找出获取100个数据集的均值,方差和MSE的方法。
set.seed(124)
a <- 1
repeat{
b = rpois(100, lambda = 3)
Storage100 <- matrix(data=b,nrow=100,ncol=1)
a = a+1
print(b)
if (a>100){break
}
}
Storage100
我希望我的100个数据集可以像第一列中的第一组数据,第二列中的第二组数据一样显示。...
答案 0 :(得分:1)
使用replicate
和simplify
作为TRUE
来获得尺寸为100 X 100的矩阵,其中每一列代表分布。
set.seed(124)
m1 <- replicate(100, matrix(data=rpois(100, lambda = 3),ncol = 1), simplify = TRUE)
要获取每一列的mean
,我们可以使用colMeans
(感谢@ jay.sf)
colMeans(m1)