Elias... obrigada por responder... Então não tem como eu calcular o HPD para meus dados? Acontece que eu estou achando a distribuição empírica das minhas predições e gostaria de obter um certo limite que me dê uma alta cobertura (um limite onde tenho máxima densidade)... Como eu poderia fazer isso?? Att. Flávia Cristina ---------------------------------------------------------------------- Message: 5 Date: Tue, 21 May 2013 18:21:57 -0300 From: Elias Teixeira Krainski <eliaskrainski@yahoo.com.br> To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Ajuda: HPD em conjunto de valores preditos Message-ID: <519BE575.5040900@yahoo.com.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" Faz sentido calcular HPD de uma distribuição, não de um grid de valores  como vc fez. A função HPDinterval() aplica-se especificamente a amostras de uma distribuição. Você aplicou-a a uma sequencia de valores para os quais foi avaliada uma densidade. Ou seja, exatamente o mesmo que fazer:     HPDinterval(as.mcmc(seq(-1, 1, length=512))) On 21/05/13 14:04, Flavia Cristina M.Queiroz wrote:

Olá pessoal, Estou precisando de uma ajuda com o cálculo de intervalo de credibilidade:

Tenho um conjunto de dados (valores preditos) com os quais tenho que calcular o intervalo de credibilidade de predição. Estes valores preditos foram calculados em uma implementação de redes neurais.

Daí eu calculei o HPD usando a função HPDinterval(), do pacote coda. Fiz assim: despad <- read.table('vintemaio_pad.txt',h=F) dados <- as.matrix(despad) a=density(dados,bw=0.45,from = -1, to = 1) library(coda) arq<-as.mcmc(a$x) mcmcsamp(a$x) HPD<-HPDinterval(arq) HPD Contudo, meu intervalo é muito amplo.. tipo, meus dados vão de -1 a 1 e meu HPD deu [-0.9882583, 0.9138943]. Porque meu intervalo deu tão amplo?? O que teria que fazer de diferente?? ------------------------------------------------------------------------------- Aguardo alguma dica ou solução!! Att., Flávia Cristina