Producing box whisker plots from a model with two fixed factors

[leahbecu]

Hi Brian,

I ran a model with two fixed factors, region and depth. I have the posterior plots, and in the output options I have used "return_obj=TRUE".

I can adjust the posterior plots using for example: NorthDeep<-gpostD5.6$both[[1]][[1]]+ ggplot2::theme_classic()

But I'm wondering if it's possible to create box whisker plots from the data instead. I have looked at the code and vignette for adjusting plots (the mantis shrimp and wolves examples). That code uses p.fac1 which doesn't exist for my data.

I have also tried: both<-jags.D5.6$BUGSoutput$sims.list$both but that just ends up with null data.

Do you have any insights into how I can extract this data for boxplots.

Thanks, Leah