ERGM MCMC diagnostics joint plot
How to combine MCMC diagnostic plots from two (or more) ERGM estimation in one figure as in the example below?
library(ergm)
library(latticeExtra)
data(florentine)
# Fit two ERGM models
model1 <- ergm(flomarriage ~ edges + kstar(2))
model2 <- ergm(flomarriage ~ edges + kstar(2) + triangle)
mcmc.diagnostic(model1)
mcmc.diagnostic(model2)
Solution
Obviously you can "glue" together the plots in some image-editting program. To do it in R you'll need to use coda (used by ergm internally anyway) and some lattice gymnastics, for example:
p1 <- model1$sample |>
coda::as.mcmc.list() |>
ergm_plot.mcmc.list(main = "Model 1")
p2 <- model2$sample |>
coda::as.mcmc.list() |>
ergm_plot.mcmc.list(main = "Model 2")
plot(p1, split=c(1,1,2,1) )
plot(p2, split=c(2,1,2,1), newpage = FALSE)
- Converting the row of a data.table to a vector
- ggplot multi-line log-log labels
- dplyr filter vector by position in pipe %>%
- Removing white spaces in a histogram?
- How to assign and match land cover type from Corine to a dataframe with a set of lon lat coordinates?
- Error in list2env(list, envir = .GlobalEnv) : names(x) must be a character vector of the same length as x
- How to traverse hclust internal node in R
- How are the abbreviations from the par function like mgp, mar, mai, deciphered?
- How to map over function with multiple arguments (and changing data)?
- Join/fill in data based on both row and column values
- How to plot a Log Pearson III distribution using `ggplot2` in R?
- Variable names and Easystats reports
- R plotly(): Adding regression line to a correlation scatter plot
- How to Use ifelse Statements in R
- Creating subgraphs from neighbors with max weight in R using igraph
- merge based on id and dates with missing data
- Automatically adjust the size of the text inside each bar on a bar plot using ggplot2 in R
- gghighlight and facet_wrap one always highlighted
- Mask and Crop functions in Terra R
- How to use a variable in dplyr::filter?
- Changing R Temporary Files
- How do you add a general label to facets in ggplot2?
- change colour of controls in R DT datatable
- Click Button Within Datatable to Display Popup Twice in a row
- How to display values in filter with ascending order
- Select Data From Column in Data Table
- Warning: Error in dataTableOutput: unused argument (height = "auto")
- How to have cells with rounded corners in {gt}?
- Two colours in the same box
- Change color of datatable filter dropdown box