-
Notifications
You must be signed in to change notification settings - Fork 1
/
Example_Code.Rmd
51 lines (41 loc) · 1.62 KB
/
Example_Code.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
---
title: "BM HCA MPAL transfer label"
output: html_notebook
---
This is an [R Markdown](http://rmarkdown.rstudio.com) Notebook. When you execute code within the notebook, the results appear beneath the code.
Try executing this chunk by clicking the *Run* button within the chunk or by placing your cursor inside it and pressing *Cmd+Shift+Enter*.
```{r}
rm(list=ls())
#Load neccessary packages
suppressPackageStartupMessages({
library(monocle3)
library(reticulate)
library(dplyr)
library(Matrix)
library(ggplot2)
library(xfun)
library(pals)
library(RColorBrewer)
library(Seurat)
library(ggplot2)
library(parallel)
library(scCustomize)
library(SeuratObject)
library(SeuratDisk)
library(tibble)
library(pbmcapply)
library(viridis)
})
Sys.setenv('R_MAX_VSIZE'=32000000000)
#Sys.getenv()
```
#load HCA bone marrow data set
```{r}
hca.immune <- Connect(filename = "/Volumes/fh/scratch/delete90/furlan_s/owalt/BM_HCA/1M-immune-human-immune-10XV2.loom", mode = "r")
hca.seurat <- as.Seurat(hca.immune)
Idents(l6.seurat) <- "ClusterName"
VlnPlot(l6.seurat, features = c("Sparc", "Ftl1", "Junb", "Ccl4"), ncol = 2)
```
Add a new chunk by clicking the *Insert Chunk* button on the toolbar or by pressing *Cmd+Option+I*.
When you save the notebook, an HTML file containing the code and output will be saved alongside it (click the *Preview* button or press *Cmd+Shift+K* to preview the HTML file).
The preview shows you a rendered HTML copy of the contents of the editor. Consequently, unlike *Knit*, *Preview* does not run any R code chunks. Instead, the output of the chunk when it was last run in the editor is displayed.