simulation_study/plots/ground_truth_CDS.R

In JacobHelwig/covdepGE: Covariate Dependent Graph Estimation

# visualizing the true graphs
library(covdepGE)
library(extrafont)
library(ggplot2)
library(ggpubr)
library(latex2exp)
loadfonts(device = "win")
dat <- generateData(100, 1, 1, 1)
true_graphs <- lapply(lapply(lapply(dat$true_precision, `!=`, 0), `*`, 1), `-`, diag(100))
graph_inds <- c("$1$,...,$n_0$", "$n_0+1$,...,$2n_0$", "$2n_0+1$,...,$3n_0$")
titles <- paste0("CDS ", 1:3, ", obs. ", graph_inds)
titles <- lapply(titles, TeX)
true_list <- lapply(1:length(true_graphs), function(k)
  matViz(true_graphs[[k]]) +
    ggtitle(titles[[k]]) +
    coord_cartesian(xlim = c(0, 4), ylim = c(0, 4)) +
    theme(plot.title = element_text(size = 18),
          axis.text = element_text(size = 18),
          text = element_text(family = "Times")
          ))
(true_graph <- ggarrange(plotlist = true_list, nrow = 1, common.legend = T))
ggsave("ground_truth_CDS.pdf", true_graph, height = 4, width = 11)