---
title: "corr interro 2"
output: pdf_document
---

```{r, message = FALSE}
library(SummarizedExperiment)
library(tidyverse)
```

```{r}
se <- readRDS("data/se.rds")
```

```{r}
rowData(se)
cbind(assay(se), rowData(se)) %>% 
  as_tibble() %>% 
  pivot_longer(names_to = "samples", values_to = "gene_expression",-c(gene,location)) %>% 
full_join(as_tibble(colData(se))) %>%
  ggplot(aes(x = group,y = gene_expression)) +
  geom_boxplot() +
  facet_wrap(~ gene)
  
```

