模拟过程中特征中的非方矩阵

Question

我使用mvn包进行了MASS模拟。

模拟应具有X和Y的均值和方差以及所选的相关性。到目前为止这是有效的：

    K=0.1
    Sim <- data.frame(X,Y)
    CovM <- matrix(data=c(var(X),(-k)*sd(X)*sd(Y),(-k)*sd(X)*sd(Y),var(Y),ncol=2,nrow=2)
    MVSim <- mvrnorm(n, mu=colMeans(Sim), Sigma=CovM, empirical=TRUE)

现在我想用for循环复制该过程：

    for(i in seq(1,m, by=2)){
      MVSimOB[, i:(i+1)] <- mvrnorm(n, mu=colMeans(Sim), Sigma=CovM, empirical=TRUE)

但是我收到以下错误消息：

  Sigma = matrix(c(var(SimYOB),  : 
  Indices out of borders

解决方案可能非常简单。但我不明白为什么我首先没有得到错误，但复制过程不起作用。 非常感谢您的帮助！