Tem alguma bibliotéca para esta função mvnorm? 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. .

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.