[R-br] Plotar várias curvas de densidades num único objeto gráfico

[Daniel Marcelino]

Interessante Benilton, sua solução ajudou bastante resolver. Converti o objeto para matrix e daí foi mais fácil seguir a sua sugestão.

Valeu,
Daniel

On Apr 5, 2014, at 12:24 AM, Benilton Carvalho <beniltoncarvalho em gmail.com> wrote:

> 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...
> 
> 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...
> 
> por exemplo:
> 
> set.seed(1)
> obs = matrix(rnorm(1000), nc=10)
> dens = apply(obs, 2, density)
> matX = sapply(dens, '[[', 'x')
> matY = sapply(dens, '[[', 'y')
> matplot(matX, matY, type='l')
> 
> 
> 
> Em 4 de abril de 2014 23:05, Daniel Marcelino <dmsilva.br em gmail.com> escreveu:
> Caro, estou com um problema para plotar densidades de distribuição
> posterior e talvez alguém aqui já tenha resolvido isso de alguma
> maneira.
> Tenho um objeto de dados de classe mcmc.list com vários paramentros,
> entre os quais eu gostaria de sobrepor três distribuições posteriores
> como abaixo. Do jeito que está, contudo, não funciona porque a cada
> comando plot o R mantem a curva no centro e altera as escalas. Tem
> alguma forma de plotar essas densidades conservando as posições delas
> originais?
> Coloquei alguns dados aqui:
> https://dl.dropboxusercontent.com/u/1339742/mcmc.samples.Rdata
> 
> 
>   quartz(width=10,height=7)
>  plot(mcmc.samples[, c(paste("alphaPSDB[", 2, "]", sep = ""))], trace
> = FALSE, density = TRUE, smooth = TRUE, lwd=3, col="blue")
> 
>   par(new=TRUE)
> 
>  plot(mcmc.samples[, c(paste("alphaPT[", 2, "]", sep = ""))], trace =
> FALSE, density = TRUE, smooth = TRUE, lwd=3, col="red")
> 
>  par(new=TRUE)
> 
>  plot(mcmc.samples[, c(paste("alphaPRB[", 2, "]", sep = ""))], trace =
> FALSE, density = TRUE, smooth = TRUE, lwd=3, col="green")
> _______________________________________________
> R-br mailing list
> R-br em 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.
> 
> _______________________________________________
> R-br mailing list
> R-br em 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140405/9a50e24f/attachment.html>