[R-br] Erro: 'Sigma' is not positive definite

[Paulo J Ribeiro Jr]

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 em 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.