генерировать мультинормальные данные в R

Question

Я пытался генерировать мультинормальные данные в R, используя mvrnorm, однако я получил ошибки как:

> epsiloni <- mvrnorm(n = 1, rep(0,8), diag(1), tol = 1e-6, empirical = FALSE, EISPACK = FALSE)
Error in mvrnorm(n = 1, rep(0, 8), diag(1), tol = 1e-06, empirical = FALSE,  
: incompatible arguments

Но он отлично работает для других ожиданий и отклонений, таких как

betai <- mvrnorm(n = 1, mu, D, tol = 1e-6, empirical = FALSE, EISPACK = FALSE)

, где mu и D указаны иначе.

Единственная разница между этими двумя - только ожидания и дисперсия, но я не видел ничего плохого в среднем иДисперсия в epsiloni, они просто мультистандартные нормальные.

Спасибо!

Answer 1

Установка Sigma на diag(1), которая является матрицей 1x1, является проблемой.Если вам нужна единичная ковариационная матрица для 8-мерного вектора, вам нужно diag(8).

mvrnorm(n = 1, rep(0, 8), diag(8))
# [1] -0.3554192  0.6051595  0.3926595  0.2752819  0.8610572  0.2679094 -1.3581420 -0.2814057