-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathREADME.Rmd
309 lines (185 loc) · 8.32 KB
/
README.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
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
---
output: github_document
---
<!-- README.md is generated from README.Rmd. Please edit that file -->
```{r, include = FALSE}
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.path = "man/figures/README-",
out.width = "100%"
)
```
# SeuratExplorer
<!-- badges: start -->
[![Lifecycle: experimental](https://img.shields.io/badge/lifecycle-experimental-orange.svg)](https://lifecycle.r-lib.org/articles/stages.html#experimental)
<!-- badges: end -->
## 为什么做这个R包
> 目前还没有很好的用于可视化Seurat分析结果的工具,当生物信息分析员将结果交给用户后,用户如果没有任何R语言基础,还是比较难去自行进行结果检索和再分析,这个R包可以帮助此类用户进行交互,实现分析结果的可视化,方便用户自行出图。用户仅需要在自己电脑上配置好R和Rstudio,然后安装运行此软件即可,无需其他操作。
> **本质上就是把R包``Seurat``里的部分绘图工具和差异分析工具进行可视化!**
> 此R包参考了开源程序:
1. [Hla-Lab/SeuratExplorer](https://github.com/rwcrocker/SeuratExplorer/): An interactive R shiny application for exploring scRNAseq data processed in Seurat.
2. [junjunlab/scRNAtoolVis](https://github.com/junjunlab/scRNAtoolVis): Some useful function to make your scRNA-seq plot more beautiful.
## A Live Demo
Support Rds file size no more than 5GB. [Open](http://www.nibs.ac.cn:666/SeuratExplorer/)
## Installation
You can install the development version of ``SeuratExplorer`` like so:
```{r eval=FALSE}
# install dependency
if (!require("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("ComplexHeatmap")
# install SeuratExplorer
if(!require(devtools)){install.packages("devtools")}
install_github("fentouxungui/SeuratExplorer")
```
**建议将各个软件包升级到最新版本!否则可能出现兼容问题**
Run App:
```{r eval=FALSE}
library(SeuratExplorer)
launchSeuratExplorer()
```
## Functions
### Load data
- support data processed by Seurat V5 and older versions. it may takes a while for loading data.
```{r echo = FALSE, message=FALSE, warning=FALSE}
knitr::include_graphics("inst/extdata/www/loading-before.png")
knitr::include_graphics("inst/extdata/www/loading-after.png")
```
### Cell Meta
- 支持下载cell meta信息,可用于统计特定分辨率下的细胞比例。
```{r echo = FALSE, message=FALSE, warning=FALSE}
knitr::include_graphics("inst/extdata/www/cell-metadata.png")
```
### Dimensional Reduction Plot
- 支持选择Dimension Reductions
- 支持选择不同的细胞分群方案
- 支持split
- 支持调整图像长宽比
- 支持是否显示cluster label
- 支持调整label大小
- 支持调整点的大小
- 支持pdf下载,所见即所得
Example plots:
```{r echo = FALSE, message=FALSE, warning=FALSE}
knitr::include_graphics("inst/extdata/www/dimplot.png")
knitr::include_graphics("inst/extdata/www/Dimplot-splited.png")
```
### Feature Plot
- 支持输入多个基因名,以空格分割,基因名可不区分大小写。支持从excel拷贝多个基因!
- 支持选择Dimension Reductions
- 支持split
- 支持自定义最大和最小表达值对应的颜色
- 支持调整图像长宽比
- 支持调整点的大小
- 支持pdf下载,所见即所得
```{r echo = FALSE, message=FALSE, warning=FALSE}
knitr::include_graphics("inst/extdata/www/featureplot-1.png")
knitr::include_graphics("inst/extdata/www/featureplot-1-with-label.png")
knitr::include_graphics("inst/extdata/www/Featureplot-split.png")
```
### Violin Plot
- 支持输入多个基因名,以空格分割,基因名可不区分大小写。支持从excel拷贝多个基因!
- 支持选择不同的细胞分群方案
- 支持split
- 支持stack和flip图像, 和指定颜色的映射
- 支持自定义最大和最小表达值对应的颜色
- 支持调整点的大小
- 支持调整x轴和y轴label的字体大小
- 支持调整图像长宽比
- 支持pdf下载,所见即所得
```{r echo = FALSE}
knitr::include_graphics("inst/extdata/www/Vlnplot-1.png")
knitr::include_graphics("inst/extdata/www/Vlnplot-2.png")
knitr::include_graphics("inst/extdata/www/vlnplot-splited-1.png")
knitr::include_graphics("inst/extdata/www/vlnplot-splited-2.png")
```
### Dot Plot
- 支持输入多个基因名,以空格分割,基因名可不区分大小写。支持从excel拷贝多个基因!
- 支持选择不同的细胞分群方案和指定所使用的clusters
- 支持split
- 支持对cluster/idents进行聚类
- 支持旋转坐标轴lable和旋转整个图像
- 支持调整点的大小和透明度
- 支持调整x轴和y轴label的字体大小
- 支持调整图像长宽比
- 支持pdf下载,所见即所得
```{r echo = FALSE}
knitr::include_graphics("inst/extdata/www/Dotplot-1.png")
knitr::include_graphics("inst/extdata/www/dotplot-splited-1.png")
```
### Heatmap Cell Level
- 支持输入多个基因名,以空格分割,基因名可不区分大小写。支持从excel拷贝多个基因!
- 支持选择不同的细胞分群方案
- 支持调整cluster label的大小、xy轴的位置和旋转角度
- 支持调整Group bar的高度
- 支持调整Group间的间隙大小
- 支持调整基因名的字体大小
- 支持调整图像长宽比
- 支持pdf下载,所见即所得
```{r echo = FALSE}
knitr::include_graphics("inst/extdata/www/Heatmap-1.png")
```
### Heatmap Group Averaged
- 支持输入多个基因名,以空格分割,基因名可不区分大小写。支持从excel拷贝多个基因!
- 支持选择不同的细胞分群方案
- 支持调整cluster label的大小、xy轴的位置和旋转角度
- 支持调整Group bar的高度
- 支持调整Group间的间隙大小
- 支持调整基因名的字体大小
- 支持调整图像长宽比
- 支持pdf下载,所见即所得
```{r echo = FALSE}
knitr::include_graphics("inst/extdata/www/averagedheatmap.png")
```
### Ridge Plot
- 支持输入多个基因名,以空格分割,基因名可不区分大小写。支持从excel拷贝多个基因!
- 支持选择不同的细胞分群方案
- 支持调整列的数目
- 支持Stack和指定颜色的映射
- 支持调整x轴和y轴label的字体大小
- 支持调整图像长宽比
- 支持pdf下载,所见即所得
```{r echo = FALSE}
knitr::include_graphics("inst/extdata/www/ridgeplot-1.png")
knitr::include_graphics("inst/extdata/www/Ridgeplot-2.png")
knitr::include_graphics("inst/extdata/www/Ridgeplot-3.png")
```
### Plot Cell Percentage
- 支持调整图像长宽比
- 支持分面
- 支持pdf下载,所见即所得
```{r echo = FALSE}
knitr::include_graphics("inst/extdata/www/cellratio-1.png")
knitr::include_graphics("inst/extdata/www/cellratio-2.png")
```
### Cluster Markers and DEGs Analysis
- 支持同时计算所有群的marker基因
- 支持计算指定群之间的差异基因
- 支持先做subset,再计算指定群之间的差异基因
- 支持自定义部分计算参数
- 支持结果下载
```{r echo = FALSE}
knitr::include_graphics("inst/extdata/www/DEGs-1.png")
knitr::include_graphics("inst/extdata/www/DEGs-2.png")
knitr::include_graphics("inst/extdata/www/DEGs-3.png")
knitr::include_graphics("inst/extdata/www/DEGs-4.png")
```
### Top Expressed Features
> 首先逐个细胞计算得出每个细胞里的高表达基因[即该基因的UMI比例大于所设定的阈值]; 然后分群(celltype列)汇总高表达基因,得到每个高表达基因在每群里发生高表达的细胞数目,以及表达比例的均值和和中位值。
### Feature Summary
> 对于感兴趣的基因list分组做简要的统计分析,得到该基因在每群里的表达阳性比例,表达水平的均值和中位值
### Feature Correlation Analysis
- ``Find Top Correlated Gene Pairs`` 用于找出top 1000 correlated gene pairs
- ``Find Correlated Genes for A Gene`` 用于找出对于感兴趣的某个基因,找出与其表达相关性强的基因
- ``Calculate Correlation for A Gene List`` 用于对一组感兴趣的基因list,计算所有组合的表达相关性
- 如结果为空,或未找到对应的Gene Pairs,是由于输入基因的表达值太低,被去除了
- 支持分区计算correlation值
- 支持pearson和spearman算法
```{r echo = FALSE}
knitr::include_graphics("inst/extdata/www/featurecorrelation.png")
```
## Session Info
```{r echo = FALSE}
sessionInfo()
```