Boa noite a todos.
Estou trabalhando com a função mvrnorm, passando como parâmetros mvrnorm(n=1, matriz1[1:5202],matriz2[5202:5202]), sendo
as matrizes (matriz1 e matriz2) matrizes hipotéticas, onde estão mostrando suas dimensões.
Em experimentos realizados, estou dividindo as matrizes realizando quatro chamas ao invés de uma como a seguir:
1) mvrnorm(n=1, matriz1[1:2601], matriz2[1:2601, 1:2601])
2) mvrnorm(n=1, matriz1[1:2601], matriz2[2602:5202, 1:2601])
3) mvrnorm(n=1, matriz1[2602:5202], matriz2[1:2601, 2602:5202])
4) mvrnorm(n=1, matriz1[2602:5202], matriz2[2602:5202, 2602:5202])
O problema está nas chamadas 2 e 3, nas quais estou recebendo seguinte mensagem : " 'Sigma' is not positive definite "
Alguém pode me auxiliar quanto a essa mensagem?
Obrigado, Carlos Beleti.