<div dir="ltr">sem saber se o plot() a q vc se refere e' um metodo especifico para o objeto de classe mcmc.list (portanto, qual o pacote q gerou mcmc.samples), nao da' pra dizer bem o que esta' acontecendo...<div>

<br></div><div>aparentemente e' um metodo sim... nesse caso, vc pode estimar a densidade (por exemplo, via density) e criar 2 matrizes (uma relativa `as coordenas X e outra Y), dai' usar matplot...</div><div><br>
</div>
<div>por exemplo:</div><div><br></div><div>set.seed(1)</div><div>obs = matrix(rnorm(1000), nc=10)</div><div>dens = apply(obs, 2, density)</div><div>matX = sapply(dens, '[[', 'x')</div><div>matY = sapply(dens, '[[', 'y')</div>

<div>matplot(matX, matY, type='l')</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 4 de abril de 2014 23:05, Daniel Marcelino <span dir="ltr"><<a href="mailto:dmsilva.br@gmail.com" target="_blank">dmsilva.br@gmail.com</a>></span> escreveu:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Caro, estou com um problema para plotar densidades de distribuição<br>
posterior e talvez alguém aqui já tenha resolvido isso de alguma<br>
maneira.<br>
Tenho um objeto de dados de classe mcmc.list com vários paramentros,<br>
entre os quais eu gostaria de sobrepor três distribuições posteriores<br>
como abaixo. Do jeito que está, contudo, não funciona porque a cada<br>
comando plot o R mantem a curva no centro e altera as escalas. Tem<br>
alguma forma de plotar essas densidades conservando as posições delas<br>
originais?<br>
Coloquei alguns dados aqui:<br>
<a href="https://dl.dropboxusercontent.com/u/1339742/mcmc.samples.Rdata" target="_blank">https://dl.dropboxusercontent.com/u/1339742/mcmc.samples.Rdata</a><br>
<br>
<br>
  quartz(width=10,height=7)<br>
 plot(mcmc.samples[, c(paste("alphaPSDB[", 2, "]", sep = ""))], trace<br>
= FALSE, density = TRUE, smooth = TRUE, lwd=3, col="blue")<br>
<br>
  par(new=TRUE)<br>
<br>
 plot(mcmc.samples[, c(paste("alphaPT[", 2, "]", sep = ""))], trace =<br>
FALSE, density = TRUE, smooth = TRUE, lwd=3, col="red")<br>
<br>
 par(new=TRUE)<br>
<br>
 plot(mcmc.samples[, c(paste("alphaPRB[", 2, "]", sep = ""))], trace =<br>
FALSE, density = TRUE, smooth = TRUE, lwd=3, col="green")<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br></div>