0

当我试图在其中一个参数中进行一些计算时,我发生了 mvrnorm() 函数的一个奇怪行为。

这确实可以完成工作:

mvrnorm(
  n=200,
  mu = c('x'=0,'y'=0), 
  Sigma = matrix(c(1,0.5,0.5,1),nrow = 2), 
  empirical = TRUE)

然而,用 n = (1-0.8)*1000 替换 n = 200 并不能完成这项工作:

mvrnorm(
  n=(1-0.8)*1000,
  mu = c('x'=0,'y'=0), 
  Sigma = matrix(c(1,0.5,0.5,1),nrow = 2), 
  empirical = TRUE)

Error in eS$vectors %*% diag(sqrt(pmax(ev, 0)), p) %*% t(X) : 
  non-conformable arguments
In addition: Warning message:
In matrix(rnorm(p * n), n) :
  data length [399] is not a sub-multiple or multiple of the number of rows [199]

为什么会这样?

4

0 回答 0