ma_results.Rmd

---
title: "Meta analysis plot organization"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

## Lets look at what MA outputs look like

```{r explore the directory}
library(PEcAn.settings)
# pecan.xml <- '/home/kristinariemer/biocro_darpa_me034_results/pecan.xml'
# settings <- PEcAn.settings::read.settings(pecan.xml)
# pft_outdir <- settings$pfts$pft$outdir
pft_outdir <- '/home/kristinariemer/biocro_darpa_me034_results_april2020/pft/SetariaWT_ME034'
dir(pft_outdir)


```


```{r explore the directory}
library(ggplot)
library(ggmcmc)

setwd(pft_outdir)
load('trait.data.Rdata', 
     'prior.distns.Rdata', 
     'post.distns.MA.Rdata', 
     'trait.mcmc.Rdata')

traits <- rownames(prior.distns)
traits_with_data <- unique(trait_data$vname)

for(trait %in% traits){
  ggmcmc::ggs_caterpillar(trait.mcmc[[trait]])
  
}

```