View source: R/plot_all_joint_marginals.R

plot_all_joint_marginals | R Documentation |

For each sampled MCMC generation, the values estimated for a pair of parameters are logged and plotted against one another. Points are color coded by when in the analysis they were sampled, so that users can visually assess mixing. A joint marginal plot is generated for all combinations of parameters, excluding the phi parameters estimated in the beta-binomial model.

```
plot_all_joint_marginals(MCMC.output, percent.burnin = 0, thinning = 1)
```

`MCMC.output` |
The standard MCMC output file generated from a BEDASSLE run. |

`percent.burnin` |
The percent of the sampled MCMC generations to be discarded as "burn-in." If the
MCMC is run for 1,000,000 generations, and sampled every 1,000 generations, there
will be 1,000 sampled generations. A |

`thinning` |
The multiple by which the sampled MCMC generations are thinned. A |

Visualizations of the joint marginal distributions allow users to (1) assess how well the MCMC is mixing, and (2) potentially diagnose instances of non-identifiability in the model. Strong linear trends in the joint marginal, or visible "ridges" in the likelihood surface, may be indicative of parameter non-identifiability, in which multiple combinations of values of these two parameters provide equally reasonable fits to the data.

Gideon Bradburd

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.