-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathQC.Rmd
48 lines (33 loc) · 1 KB
/
QC.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
---
title: "Quality Control"
author: "D. Ford Hannum Jr."
date: "7/16/2020"
output:
html_document:
toc: true
toc_depth: 3
number_sections: false
theme: united
highlight: tango
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = FALSE, message = FALSE)
library(Seurat)
library(ggplot2)
library(data.table)
library(MAST)
```
```{r reading in the data}
wbm <- readRDS('./data/wbm_clustered_filtered_named.rds')
```
# Plotting Different QC Metrics by Clusters
```{r plots on plots}
VlnPlot(wbm, features = 'nFeature_RNA', split.by = 'cluster_IDs')
VlnPlot(wbm, features = 'nCount_RNA', split.by = 'cluster_IDs')
VlnPlot(wbm, features = 'percent.mt', split.by = 'cluster_IDs')
```
```{r feature plots}
FeatureScatter(wbm, feature1 = 'nCount_RNA', feature2 = 'nFeature_RNA')
FeatureScatter(wbm, feature1 = 'nCount_RNA', feature2 = 'percent.mt')
FeatureScatter(wbm, feature1 = 'percent.mt', feature2 = 'nFeature_RNA')
```