<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">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: <br>
HPDinterval(as.mcmc(seq(-1, 1, length=512)))<br>
<br>
On 21/05/13 14:04, Flavia Cristina M.Queiroz wrote:<br>
</div>
<blockquote
cite="mid:1369155882.21299.YahooMailNeo@web126204.mail.ne1.yahoo.com"
type="cite">
<div style="color:#000; background-color:#fff; font-family:times
new roman, new york, times, serif;font-size:18pt">
<div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.800000190734863px; line-height:
13.600000381469727px;">Olá pessoal, </div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.727272033691406px; line-height:
13.600000381469727px; background-color: transparent;
font-style: normal;"><span style="background-color:
transparent; line-height: 13.600000381469727px;">Estou
precisando de uma ajuda com o cálculo de intervalo de
credibilidade:</span></div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.727272033691406px; line-height:
13.600000381469727px; background-color: transparent;
font-style: normal;"><span style="background-color:
transparent; line-height: 13.600000381469727px;"><br>
</span></div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.727272033691406px; line-height:
13.600000381469727px; background-color: transparent;
font-style: normal;"><span style="background-color:
transparent; line-height: 13.600000381469727px;">Tenho um
conjunto de dados (valores preditos) com os quais tenho
que calcular o intervalo de credibilidade de predição.</span><br>
</div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.727272033691406px; line-height:
13.600000381469727px; background-color: transparent;
font-style: normal;"><span style="background-color:
transparent; line-height: 13.600000381469727px;">Estes
valores preditos foram calculados em uma implementação de
redes neurais.</span></div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.727272033691406px; line-height:
13.600000381469727px; background-color: transparent;
font-style: normal;"><span style="background-color:
transparent; line-height: 13.600000381469727px;"><br>
</span></div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.727272033691406px; line-height:
13.600000381469727px; background-color: transparent;
font-style: normal;">Daí eu calculei o HPD usando a função
HPDinterval(), do pacote coda.<br>
Fiz assim:</div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.800000190734863px; line-height:
13.600000381469727px;">despad <-
read.table('vintemaio_pad.txt',h=F)<br>
dados <- as.matrix(despad)<br>
a=density(dados,bw=0.45,from = -1, to = 1)</div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.800000190734863px; line-height:
13.600000381469727px;">library(coda)<br>
arq<-as.mcmc(a$x) <br>
mcmcsamp(a$x)<br>
HPD<-HPDinterval(arq)<br>
HPD</div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.800000190734863px; line-height:
13.600000381469727px;">Contudo, meu intervalo é muito
amplo.. tipo, meus dados vão de -1 a 1 e meu HPD deu
[-0.9882583, 0.9138943].</div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.800000190734863px; line-height:
13.600000381469727px;">Porque meu intervalo deu tão amplo??
O que teria que fazer de diferente??<br>
-------------------------------------------------------------------------------</div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.800000190734863px; line-height:
13.600000381469727px;">Aguardo alguma dica ou solução!!</div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.800000190734863px; line-height:
13.600000381469727px;">Att.,</div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.800000190734863px; line-height:
13.600000381469727px;">Flávia Cristina</div>
<div style="margin-top: 10px; color: rgb(51, 51, 51);
font-family: 'lucida grande', tahoma, verdana, arial,
sans-serif; font-size: 12.727272033691406px; line-height:
13.600000381469727px; background-color: transparent;
font-style: normal;"><br>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>