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)
- rlanguageserver: Disable specific linter and avoid line-breaking
- Replace truncated word in a string
- Add label to right y axis in R ggplot
- Popup offset without using markers
- Multi marker depending a variable in a data frame
- Drawing numerous polygons
- Select icons or CircleMarkers based on factor variable
- Why is a map not displayed when I call leaflet in a function
- Display SpatialPolygonsDataFrame on leaflet map
- Leaflet in R not showing variable circle colors
- Overlay two population pyramid plots in one ggplot graph
- How to most efficiently filter a dataframe conditionally of values in another one, in the tidyverse framework?
- R group_by multiple conditionnal aggregation functions dplyr
- Conditionally split strings cells in a dataframe R
- How to have Shiny app render multiple plots
- Disable Shiny button conditionally to a rHandsontable format
- Is there a way to render dynamic headings in Quarto using R?
- What is the time complexity of sample?
- Add geom_text with x bar and round values
- Color-fill ggplot2 maps with a single color
- lapply a ggplot function on multiple inputs with consistent groups
- Is there a way to save text in R without it actively implementing \n?
- How to get number of unique observations using tbl_summary
- Using if_else in a mutate command yields an error about false having the wrong size
- Rotate x-axis value in dendogram using fviz_dend
- Export data frames in list to xlsx with named sheets
- Query all primary roads from PostGis Database filled with OSM Data within a boundary
- Write a list in r to an excel file
- Calculating net assets with missing values: discrepancy in results
- Error in curl::curl_fetch_memory(url, handle = handle) : Couldn't connect to server in R: oauth2.0_token())