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]] > > ______________________________________________ > [EMAIL PROTECTED] 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. > ______________________________________________ R-help@r-project.org 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.
