[R] Error bars using data from bugs()

[toby909 at gmail.com]

only one out of many many

dt = cbind(sim$mean$u2+sim$mean$beta1,sim$sd$u2)
dt = dt[order(dt[,1]),]
bounds = 
cbind(c(dt[,1]-1.96*dt[,2],dt[,1]+1.96*dt[,2],dt[,1]),rep(1:length(sim$mean$u2),3))
bounds = bounds[order(bounds[,2]),]
plot(bounds)

T



Matthew Krachey wrote:
> I'm trying to compare the results of several models using output from bugs(). I want to summarize the models via boxplots or lineplots of the parameters. I know how to use errorbars from a regression, but how do I implement the credible intervals as error bars into these plots
> 
> Any help is much appreciated
> 
> Matthew Krachey
>  
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>