Add beginnings of script to prepare the Zhao data

code/prep-zhao.R

+# This scripts prepares and saves the `zhao` data.
+
+args <- R.utils::commandArgs(
+  trailingOnly = TRUE,
+  asValues = TRUE
+)
+
+suppressMessages({
+  library(scRNAseq)
+  library(scater)
+  library(scran)
+})
+
+set.seed(07122021)
+
+zhao <- readRDS(args$raw_data)
+
+# Filter to only B cells
+
+# ...