Упорядочить другой набор данных, используя матричный код

Question

Я пытаюсь использовать повторный цикл для генерации 100 наборов данных распределения Пуассона с размером выборки n = 100, и я хотел бы упорядочить результат по строкам и столбцам, но это просто покажет, что я повторяюсь, чтобы показать мне последнийнабор данных пока не весь набор данных.В то же время я бы также попытался выяснить, как получить среднее значение, дисперсию и MSE для набора данных 100.

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 может отображаться как первыйнабор данных в первом столбце, второй набор данных во втором столбце .....

Answer 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)