
Junior além do comentário do Benilton, eu na sequencia investigaria com cuidado como esta matriz é construída (segundo qual modelo) e se há oa menos garantias de que seja teoricamente válida (positiva definida) . Note que teoricamente válida nao implica em numericamente computável, ainda mais se tratando de uma matriz de dimensao bastante elevada estando tudo ok voce poderia considerar outras alternativas para a geracao o que precisa, por exemplo tentando a decomposição por auto valores/vetores Em Qui, 2012-02-16 às 00:46 -0200, Junior Beleti escreveu:
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. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.