1. a fun;çãO É mvrnorm() (note o r antes do norm) 2. os mecanismos de procura te mostrariam que está no pacote MASS e que éoutras funcoes parecidas em outros pacotes Em Qui, 2012-02-16 às 10:20 -0200, Mauro Sznelwar escreveu:
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. . _______________________________________________ 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.