n <- 10 A <- matrix(rnorm(n^2), n, n) A <- A %*% t(A) condMVN(mean=rep(1,n), sigma=A, dependent=c(2,3,5), given=integer(),X.given=integer())

Question

产生此错误：

Error in solve.default(D) : 'a' is 0-diml 此错误是什么意思？

condMVNorm程序包使在R内部轻松实现条件法线。这是一个愚蠢且可复制的示例：n

Answer 1

结果：

$condMean
[1] 1 1 1

$condVar
         [,1]      [,2]      [,3]
[1,] 9.489041  2.157854  6.143564
[2,] 2.157854  6.314409 -1.636732
[3,] 6.143564 -1.636732 10.131375