workflowr::wflow_remove("analysis/zhao-data.Rmd")

analysis/zhao-data.Rmd

----
-title: "Zhao dataset"
-output: html_document
----
-
-```{r setup, include=FALSE}
-knitr::opts_chunk$set(echo = TRUE)
-```
-
-```{r libraries}
-library(tibble)
-library(dplyr)
-library(ggplot2)
-library(scran)
-library(scater)
-```
-
-```{r load-data}
-zhao <- readRDS(here::here("data", "real_data", "zhao.rds"))
-```
-
-```{r}
-seurat_res <- readRDS(here::here("results", "real_data", "zhao-seurat_wilcox.rds"))
-scanpy_res <- readRDS(here::here("results", "real_data", "zhao-scanpy_t_rankby_raw.rds"))
-```
-
-# Naive B
-
-```{r}
-scanpy_res$result %>% 
-  filter(cluster == "Naive B")
-```
-
-# MZB
-
-```{r}
-scanpy_res$result %>% 
-  filter(cluster == "MZB")
-```