---

title: "Figure 7"

author: Tobias Roider

date: "Last compiled on `r format(Sys.time(), '%d %B, %Y, %X')`"

output: 

  rmdformats::readthedown:

editor_options: 

  chunk_output_type: console

---

```{r options, include=FALSE, warning = FALSE}

library(knitr)

opts_chunk$set(echo=TRUE, tidy=FALSE, include=TRUE, message=FALSE, 

               dpi = 100, cache = FALSE, warning = FALSE)

opts_knit$set(root.dir = "../")

options(bitmapType = "cairo")

```

# Read data, functions and packages

```{r read data}

source("R/ReadPackages.R")

source("R/Functions.R")

source("R/ReadData.R")

source("R/ThemesColors.R")

source("R/Helpers.R")

```

# Representative rLN image

```{r create rLN image, fig.height=3}

plots_codex <- list()

dpi <- 400

rLN <- codex_annotation %>% filter(unique_region== "191_4reg001") %>% 

  filter(x>500, x<7500) %>% 

  filter(y>500, y<7500) %>% 

  mutate(Merged_all_simple=ifelse(Merged_final %in% c("Granulo", "Macro", "DC"), "Myeloid", Merged_final)) %>% 

  mutate(Merged_all_simple=ifelse(Merged_all_simple %in% c("MC", "NKT", "PC", "NK"), "Other", Merged_all_simple)) %>% 

  filter(((x-mean(.$x))^2+(y-mean(.$y))^2)<2500^2)

image_rln <- ggplot()+

  geom_point_rast(data=rLN %>% filter(Merged_all_simple=="B"), aes(x=x,y=y), shape=21, size=0.25, stroke=0, alpha=1, raster.dpi =dpi, 

                  color=colors_codex[["B"]], fill=colors_codex[["B"]])+

  geom_point_rast(data=rLN %>% filter(Merged_all_simple!="B"), aes(x=x,y=y, fill=Merged_all_simple, color=Merged_all_simple), 

                  shape=21, size=0.25, stroke=0, alpha=1, raster.dpi=dpi)+

  scale_color_manual(values = colors_codex, limits=limits_codex, labels=labels_codex, name=NULL)+

  scale_fill_manual(values = colors_codex, limits=limits_codex, labels=labels_codex, name=NULL)+

  guides(fill=guide_legend(nrow = 6, override.aes = list(size=1.5, color="white", stroke=0.1)))+

  ggtitle(unique(rLN$Entity))+

  coord_fixed()+

  theme_void()+

  theme(legend.position = "right",

        legend.box.background = element_rect(fill = "black"),

        legend.box.margin = unit(units = "cm", c(0, 0, 0, -0.25)),

        legend.spacing.x = unit("cm", x = 0.1),

        legend.key.height = unit("cm", x = 0.34),

        legend.key.width = unit("cm", x = 0.2),

        legend.text = element_text(color="white", size=7),

        plot.title = element_text(color="white", hjust=0.1, size=8, 

                                  margin = unit(units = "cm", c(0,0,-0.6,0)), face = "bold"),

        plot.margin = unit(units = "cm", c(0.35, 0.1, 0.1, 0.1)),

        plot.background = element_rect(fill = "black", color="black"),

        panel.background = element_rect(fill = "black", color="black"))

image_rln

#ggsave(width = 7, height = 4.75, units = "cm", filename = "Figure7_p1.pdf")

```

## Mini B-cell plot

```{r mini B cell plot, fig.height=1, eval=FALSE, include=FALSE}

ggplot(rLN %>% filter(Merged_all_simple=="B"), aes(x=x,y=y))+

  geom_point_rast(raster.dpi = dpi, alpha=0.5, shape=".", color="grey75")+

  guides(color=guide_legend(override.aes = list(size=2,alpha=0.75)))+

  scale_y_continuous(expand = c(0,0))+#, limits = c(min(rLN$y), max(rLN$y)+edgeFt))+

  scale_x_continuous(expand = c(0,0))+#, limits = c(min(rLN$x), max(rLN$x)+edgeFt))+

  coord_fixed()+

  mytheme_codex+

  theme(plot.margin = unit(units = "cm", c(0,0,0,0)))

#ggsave(width = 5, height = 5, units = "cm", filename = "Figure7_mini.pdf")

```

# Neighborhood (NH) plots

## Load NH analysis and PCA

```{r neighbourhood analysis}

# Read results from neighborhood analysis

# Please run file: analysis/NeighborhoodAnalysis.Rmd

load("output/Neighborhood_results.RData")

# Add codex annotation

codex_annotation <- left_join(codex_annotation, nn_classes, by="unique_cell_id")

codex_annotation

```

## In situ NH plot of rLN

```{r overview nn, fig.height=3.5}

plot_overview <- codex_annotation %>% 

  filter(unique_region=="191_4reg001") %>% 

  filter(x>500, x<7500) %>% 

  filter(y>500, y<7500) %>% 

  filter(((x-mean(.$x))^2+(y-mean(.$y))^2)<2500^2) %>% 

  ggplot()+

  ggrastr::geom_point_rast(aes(x=x,y=y,color=Region, fill=Region), shape=21, size=0.25, stroke=0, alpha=1, raster.dpi =300)+

  scale_color_manual(values = colors_nn)+

  scale_fill_manual(values = colors_nn)+

  guides(color=guide_legend(override.aes = list(size=3)))+

  coord_fixed()+

  theme_void()+

  theme(legend.position = "none",

        legend.title = element_blank(),

        legend.text = element_text(size=6),

        legend.spacing.x = unit("cm", x = 0.1),

        legend.key.height = unit("cm", x = 0.4),

        legend.key.width = unit("cm", x = 0.2))

plot_overview

#ggsave(width = 5.3, height = 5.3, units = "cm", filename = "Figure7_p2.pdf")

```

## PCA

```{r pca nn, fig.height=3.5, fig.width=3.5}

df_loadings <- pca_codex$rotation[c("Stromal cells", "Macro", "B", "TFH", "FDC", "TTOX", "CD4T", "Treg"), c("PC1", "PC2")] %>% 

  data.frame() %>% mutate(x=0, y=0) %>% 

  rownames_to_column("Ident")

scaling <- 7.5

plot_pca <- pca_codex$x %>% 

  data.frame() %>% 

  rownames_to_column("unique_cell_id") %>% 

  left_join(., nn_classes) %>% 

  filter(unique_cell_id %in% rLN$unique_cell_id) %>%

  sample_frac(0.3) %>% 

  ggplot(aes(x=PC1, y=PC2))+

  ggrastr::geom_point_rast(size=0.5, alpha=1, shape=21, stroke=0, aes(color=Region, fill=Region), raster.dpi = 400)+

  geom_segment(data=df_loadings, aes(x=0, xend=7.5*PC1, y=0, yend=7.5*PC2), 

               arrow = arrow(type = "closed",  length = unit(units = "cm", 0.1)), size=0.25)+

  #ggrepel::geom_text_repel(data=df_loadings, aes(x=7.5*PC1, y=7.5*PC2, label=Ident), size=2.5, segment.size=0.25)+

  guides(color=guide_legend(override.aes = list(size=3, alpha=1)))+

  scale_color_manual(values = colors_nn)+

  scale_fill_manual(values = colors_nn)+

  ylim(-7.5,5)+

  mytheme_1+

  #coord_fixed(clip = "off")+

  theme(legend.position = "none")

plot_pca

#ggsave(width = 5, height = 5.1, units = "cm", filename = "Figure7_p3.pdf")

```

## NH composition

```{r, include=FALSE}

df_nh <- 

  codex_annotation %>% 

  add_prop(vars = c("Region", "Merged_final"), group.vars = 1) %>% 

  group_by(Merged_final) %>% 

  dplyr::mutate(Prop=scale(Prop)[,1]) 

pheat_nh <- df_nh %>% 

  pivot_wider(names_from = "Merged_final", values_from = "Prop") %>% 

  column_to_rownames("Region") %>% 

  pheatmap::pheatmap(silent = T)

plot_nn_rLN <- ggplot(df_nh, aes(x=Merged_final, y=Region, fill=Prop))+

  geom_tile()+

  scale_fill_gradientn(colours = colorRampPalette(colors = c("#762a83", "#f7f7f7", "#1b7837"))(100), limits=c(-3, 3),

                       name="Scaled\nAbundance", breaks=c(-3,-1.5,0,1.5,3))+

  scale_x_discrete(limits=pheat_nh$tree_col$labels[pheat_nh$tree_col$order], expand = c(0,0),

                   labels=c("Plasma cells", "Mast cells", "Stromal cells", "Granulocytes", "NK cells", expression('T'[Pr]), 

                            expression('T'[REG]), expression('Memory T'[TOX]), expression('CD8'^'+'~'naive'), "DC", 

                            expression('CD4'^'+'~'naive'), expression('Memory T'[H]), "NK T-cells", "Macrophages",

                            expression('Exh. T'[TOX]), "B cells", "FDC", expression('T'[FH])))+

  scale_y_discrete(limits=rev(pheat_nh$tree_row$labels[pheat_nh$tree_row$order]), expand = c(0,0), name="Neighborhoods")+

  geom_vline(xintercept = seq(1.5, 17.5, 1), color="white", size=0.25)+

  geom_hline(yintercept = seq(1.5, 14.5, 1), color="white", size=0.25)+

  mytheme_1+

  coord_cartesian(clip = "off")+

  theme_bw()+

  mytheme_1+

  theme(legend.position = "right",

    axis.title.x = element_blank(),

        axis.title.y = element_text(vjust = 11),

        axis.text.x = element_text(angle=45, hjust=1, size=7),

        legend.text = element_text(size = 7, color="black"),

        legend.title = element_text(size = 7, color="black", vjust = 0.8),

        legend.key.height = unit(0.3, "cm"),

        legend.key.width = unit(0.3, "cm"),

        legend.box.spacing = unit(0.1, "cm"),

        axis.text.y = element_blank(),

        axis.ticks.y = element_blank(),

        plot.tag = element_text(margin = unit(c(0,0.45,0,0), units = "cm")),

        plot.margin = unit(c(0,0.25,0,0.25), "cm"))+

  labs(tag = "E")

order_y <- rev(pheat_nh$tree_row$labels[pheat_nh$tree_row$order])

for(i in 1:length(pheat_nh$tree_row$order)) {

  plot_nn_rLN <- plot_nn_rLN+

    annotation_custom(grob = rectGrob(gp = gpar( fill=colors_nn[order_y][i], lex=1, col="white")), 

                      ymin = seq(0.5, length(colors_nn)-0.5, 1)[i], 

                      ymax = seq(1.5, length(colors_nn)+0.5, 1)[i],

                      xmin = 0, xmax = -1.1)+

    annotation_custom(grob = textGrob(label = paste0("N", order_y)[i], gp = gpar(cex=0.6)), 

                      ymin = seq(0.5, length(colors_nn)-0.5, 1)[i], 

                      ymax = seq(1.5, length(colors_nn)+0.5, 1)[i],

                      xmin = 0, xmax = -1.1)

}

```

## NH proportions

```{r neighborhood composition, fig.height=3}

roi <- c(4,2,1,7)

df_freq_nh <- codex_annotation %>% 

  left_join(., nn_classes) %>% 

  add_prop(vars = c("Entity", "Region", "unique_region"), group.vars = 3) %>% 

  fill_zeros(names_from = "Region", values_from = "Prop") 

pvalues <- df_freq_nh %>% 

  compare_means(data=., formula = Prop ~ Entity, ref.group = "rLN", 

                group.by = "Region", p.adjust.method = "BH") %>% 

  filter(p.adj<0.07) %>% 

  mutate(p.adj_s=format(p.adj, scientific = TRUE, digits=1)) %>% 

  mutate(p.adj_f=case_when(p.adj > 0.01 ~ as.character(round(p.adj, 2)),

                           p.adj==0.01 ~ "0.01",

                           p.adj < 0.01 ~ p.adj_s),

         Entity=group2) %>% 

  filter(!is.na(p.adj)) %>% 

  mutate(Entity=factor(Entity, levels = c("rLN", "DLBCL", "MCL", "FL", "MZL"))) %>% 

  filter(Region %in% roi) %>% 

  mutate(Region=factor(Region, levels = roi)) %>% 

  arrange(Region,Entity) %>% 

  mutate(Y=c(0.68,0.1,0.35,0.55,0.77,0.09,0.20,0.09)) %>% 

  mutate(Region_nn=paste0("N", Region)) %>% 

  mutate(Region_nn=factor(Region_nn, levels = paste0("N", roi))) 

# Points to modify facet scales

d <- data.frame(Entity="rLN", Region=roi, Y=c(0.7, 0.55, 0.79, 0.35)) %>% 

  mutate(Region=factor(Region, levels = roi)) %>% 

  mutate(Region_nn=paste0("N", Region)) %>% 

  mutate(Region_nn=factor(Region_nn, levels = paste0("N", roi))) 

df_medianLines <- df_freq_nh %>%

  filter(Entity=="rLN") %>% 

  group_by(Region) %>% 

  dplyr::summarise(MedianProp=median(Prop)) %>% 

  filter(Region %in% roi) %>% 

  mutate(Region=factor(Region, levels = roi)) %>% 

  mutate(Region_nn=paste0("N", Region)) %>% 

  mutate(Region_nn=factor(Region_nn, levels = paste0("N", roi))) 

plot_freq_nn <- df_freq_nh %>% 

  filter(Region %in% roi) %>% 

  mutate(Region_nn=paste0("N", Region)) %>% 

  mutate(Region_nn=factor(Region_nn, levels = paste0("N", roi))) %>% 

  mutate(Region=factor(Region, levels = roi)) %>% 

  ggplot(aes(x=Entity, y=Prop)) +

  geom_hline(data=df_medianLines, aes(yintercept=MedianProp),

             size=0.25, linetype="dashed", color="grey60")+

  geom_boxplot(width=0.5, outlier.alpha = 0, size=0.25)+

  ggbeeswarm::geom_beeswarm(size=0.8, shape=21, stroke=0.1, cex = 1.75, aes(fill=Region))+

  geom_text(data=pvalues, inherit.aes = F, aes(y=Y, x=Entity, label=p.adj_f), size=2.5)+

  geom_point(data = d, alpha=0, aes(x=Entity, y=Y))+

  scale_fill_manual(values = colors_nn)+

  scale_x_discrete(limits=c("rLN", "DLBCL", "MCL", "FL", "MZL"))+

  facet_wrap(~Region_nn, strip.position = "right", scales = "free_y")+

  ylab("% of total area")+

  mytheme_1+

  theme(strip.text.y = element_text(angle = 0, size=7, margin = unit(units = "cm", c(0.075,0.075,0.075,0.075))),

        axis.text.x = element_text(angle=45, hjust=1),

        axis.title.x = element_blank(),

        plot.margin = unit(c(0,0,0,0.1), "cm"))+

  labs(tag = "F")

g <- ggplot_gtable(ggplot_build(plot_freq_nn))

g$grobs[[22]]$grobs[[1]]$children[[1]]$gp$fill <- colors_nn["2"]

g$grobs[[23]]$grobs[[1]]$children[[1]]$gp$fill <- colors_nn["7"]

g$grobs[[24]]$grobs[[1]]$children[[1]]$gp$fill <- colors_nn["4"]

g$grobs[[25]]$grobs[[1]]$children[[1]]$gp$fill <- colors_nn["1"]

plot_nn_rLN+wrap_ggplot_grob(g)+plot_layout(widths = c(1.6,1.2))

ggsave(width = 18.3, height = 6, units = "cm", filename = "Figure7_p4.pdf")

```

# B-NHL examples

## Cells colored by NH

```{r, fig.height=3}

regions_nn <- c("191_2reg006", "191_3reg003", "191_4reg002")

names(regions_nn) <- c("MCL", "FL", "DLBCL")

plot_nn <- list()

for(r in regions_nn){

  df_tmp <- codex_annotation %>% 

    filter(unique_region==r) %>% 

    filter(x>500, x<7500) %>% 

    filter(y>500, y<7500) %>% 

    filter(((x-mean(.$x))^2+(y-mean(.$y))^2)<2500^2)

plot_nn[[r]] <- df_tmp %>% 

    ggplot()+

    geom_point(data = data.frame(x=min(df_tmp$x)+2500, y=min(df_tmp$y)+2500), stroke=2,

                             aes(x=x,y=y), shape=21, color="white", fill="white", size=72.5)+

    ggrastr::geom_point_rast(aes(x=x,y=y,color=Region, fill=Region), shape=21, size=0.25, stroke=0, alpha=1, raster.dpi = 400)+

    scale_color_manual(values = colors_nn)+

    scale_fill_manual(values = colors_nn)+

    guides(color=guide_legend(override.aes = list(size=3)))+

    coord_fixed(clip = "off")+

    theme_void()+

    theme(legend.position = "none",

          legend.title = element_blank(),

          plot.background = element_rect(fill=NA, color=NA),

          legend.text = element_text(size=6),

          legend.spacing.x = unit("cm", x = 0.1),

          legend.key.height = unit("cm", x = 0.4),

          legend.key.width = unit("cm", x = 0.2))

#ggsave(plot_nn[[r]], width = 6, height = 6, units = "cm", filename = paste("Figure7_", r, ".pdf"))

}

plot_nn

```

## Cells colored by subset

```{r, fig.height=3}

regions_nn <- c("191_2reg006", "191_3reg003", "191_4reg002")

names(regions_nn) <- c("MCL", "FL", "DLBCL")

df_images <- list()

images <- list()

margins <- c(0.1, 0, 0.4, -0.5)

dpi <- 600

for(r in 1:3){

df_images[[r]] <- codex_annotation %>% 

  filter(unique_region==regions_nn[r]) %>% 

  filter(x>500, x<7500) %>% 

  filter(y>500, y<7500) %>% 

  mutate(Merged_all_simple=ifelse(Merged_final %in% c("Granulo", "Macro", "DC"), "Myeloid", Merged_final)) %>% 

  mutate(Merged_all_simple=ifelse(Merged_all_simple %in% c("PC", "MC", "NK", "NKT"), "Other", Merged_all_simple)) %>% 

  filter(((x-mean(.$x))^2+(y-mean(.$y))^2)<2500^2)

images[[r]] <- 

  ggplot()+

  geom_point_rast(data=df_images[[r]] %>% filter(Merged_all_simple=="B"), aes(x=x,y=y), 

                  shape=21, size=0.25, stroke=0, alpha=1, raster.dpi =dpi, 

                  color=colors_codex[["B"]], fill=colors_codex[["B"]])+

  geom_point_rast(data=df_images[[r]] %>% filter(Merged_all_simple!="B"), 

                  aes(x=x,y=y, fill=Merged_all_simple, color=Merged_all_simple), 

                  shape=21, size=0.25, stroke=0, alpha=1, raster.dpi=dpi)+

  scale_color_manual(values = colors_codex, limits=limits_codex, labels=labels_codex, name="Cell type")+

  scale_fill_manual(values = colors_codex, limits=limits_codex, labels=labels_codex, name="Cell type")+

  ggtitle(unique(df_images[[r]]$Entity))+

  coord_fixed()+

  theme_void()+

  theme(legend.position = "none",

        plot.title = element_text(color="white", hjust=0.1, size=10, 

                                  margin = unit(units = "cm", c(0,0,-1,0)), face = "bold"),

        plot.margin = unit(units = "cm", margins),

        panel.background = element_rect(fill = "black", color="black"),

        plot.background =  element_rect(fill = "black", color="black"))

}

#images

emptyplot <- ggplot()+

  geom_point_rast(data=df_images[[1]] %>% filter(Merged_final=="BC"), aes(x=x,y=y), raster.dpi = dpi, shape=".",

                  color=colors_codex[["B"]])+

  coord_fixed()+

  mytheme_codex+

  theme(panel.background = element_rect(fill = "black", color="black"),

        plot.background =  element_rect(fill = "black", color="black"))

p_full <- images[[3]]+emptyplot+images[[1]]+emptyplot+images[[2]]+emptyplot+plot_layout(widths = c(1,0.1,1,0.1,1,0.3))

p_full

#ggsave(p_full, width = 22.5, height = 6, units = "cm", filename = "Figure7_mini.pdf")

```

## Legend

```{r, fig.height=1}

plot.legend <- images[[r]]+

 guides(fill=guide_legend(nrow = 1, override.aes = list(size=1.75, color="white", stroke=0.25)))+

 guides(color=guide_legend(nrow = 1, override.aes = list(size=1.75, color="white", stroke=0.25)))+

 theme(legend.position = "bottom",

       legend.title = element_blank(),

        legend.box.background = element_rect(fill = "black"),

        legend.box.margin = unit(units = "cm", c(0, 0, 0, 0)),

        legend.spacing.x = unit("cm", x = 0.1),

        legend.key.height = unit("cm", x = 0.34),

        legend.key.width = unit("cm", x = 0.16),

        plot.margin = unit(units = "cm", c(0,0,0,0)),

        plot.background = element_rect(fill = "black", color="black"),

        panel.background = element_rect(fill = "black", color="black"),

        legend.text = element_text(color="white", size=6.5))

as_ggplot(get_legend(plot.legend))

#ggsave(width = 19, height = 1, units = "cm", filename = "Figure7_legend.pdf")

```

# Closest cells to B-cells

```{r fig.height=2.2}

codex_freq <- codex_annotation %>% 

  add_prop(vars = c("unique_region", "Merged_final"), group.vars = 1)

nn <- run_NNanalysis(data = codex_annotation, regions = unique(codex_annotation$unique_region), 

                     plan_session = "multisession", 

                     add.prop=FALSE,

                     n_workers = 10, 

                     nn = 1)

nn_sum <- 

  nn %>% select(-name) %>% 

  left_join(codex_annotation %>% select(unique_cell_id, unique_region)) %>% 

  left_join(codex_annotation %>% select(unique_cell_id, Ident_center=Merged_final)) %>% 

  filter(Ident_center=="B", ) %>% 

  add_prop(vars = c("unique_region", "Merged_final"), group.vars = 1) %>% 

  mutate(Prop=100*Prop) %>% 

  left_join(., codex_annotation %>% select(unique_region, Entity) %>% distinct)

art_max <- c(20, 14, 11, 12)

names(art_max) <- c("rLN", "DLBCL", "MCL", "FL")

plot_nn <- list()

for(e in names(art_max)){

  selected <- 

    nn_sum %>% filter(Entity==e) %>% 

    group_by(Merged_final) %>% 

    summarise(Mean=mean(Prop), SEM=sd(Prop)/sqrt(length(Prop))) %>% 

    top_n(Mean, n = 10) %>% 

    mutate(code=Merged_final=="FDC") %>% 

    mutate(SEM=ifelse(Mean>art_max[e], NA, SEM)) %>% 

    mutate(Mean_new=ifelse(Mean>art_max[e], art_max[e], Mean)) 

  order_ <- selected %>% 

    arrange(desc(Mean_new)) %>% 

    pull(Merged_final)

  mean_pat <- nn_sum %>% filter(Entity==e) %>% 

    left_join(codex_annotation %>% select(PatientID, unique_region) %>% distinct) %>% 

    group_by(PatientID, Merged_final) %>% 

    summarise(Mean_pat=mean(Prop))

 plot_nn[[e]] <- 

    selected %>% 

    ggplot(aes(x=reorder(Merged_final, -Mean_new), y=Mean_new, fill=Merged_final, color=code))+

    geom_errorbar(aes(ymin=Mean_new, ymax=Mean_new+SEM), width=0.2, color="black", size=0.25)+

    geom_bar(stat = "identity", size=0.5, width=0.5, fill="white", color="white")+

    geom_bar(stat = "identity", size=0.25, width=0.5, alpha=0.6)+

    ggbeeswarm::geom_beeswarm(data=mean_pat, inherit.aes = F, aes(x=Merged_final, y=Mean_pat),

               color="black", stroke=0, size=0.65, alpha=0.5, cex=1.75)+

    annotation_custom(grob = textGrob(label = e, just = "right", x = 0.92, y=0.9, gp = gpar(cex=0.6)))+

    geom_segment(inherit.aes = F, 

                 aes(x=1, xend=1, y=1.02*Mean_new[1], yend=1.14*Mean_new[1]), 

                 color="black", size=0.15,

                 arrow = arrow(type = "closed", length = unit(units = "cm", 0.1)))+

    geom_text(aes(x=2.2, y=1.08*Mean_new[1], label=round(Mean[1],1)), 

              check_overlap = T, size=2.5, color="black")+

    scale_color_manual(values = c("white", "black"), limits=c(F, T))+

    scale_fill_manual(values = colors_codex_exp)+

    scale_x_discrete(limits=order_, labels=unlist(labels_codex_exp))+

    scale_y_continuous(name = "% of cells closest\nto B-cells", limits=c(0,1.15*art_max[e]))+

    mytheme_1+

    theme(legend.position = "none",

          axis.text.x = element_text(angle=45, hjust=1, size=7),

          plot.background = element_rect(fill = NA, colour = NA),

          panel.background = element_rect(fill = NA, colour = NA),

          plot.margin = unit(c(0,0.1,0,0), "cm"),

          axis.title.x = element_blank())

  if(e!="rLN") {

    plot_nn[[e]] <- plot_nn[[e]]+

      theme(axis.title.y = element_blank())}

}

plot_nn$rLN+labs(tag = "G")+plot_nn$DLBCL+

  plot_nn$MCL+plot_nn$FL+

  plot_layout(nrow = 1)

#ggsave(width = 18, height = 4.75, units = "cm", filename = "Figure7_p5.pdf")

```

# Session info

```{r}

sessionInfo()

```