-
Notifications
You must be signed in to change notification settings - Fork 0
/
gse123902_cytokines_checkpoint.R
142 lines (128 loc) · 7.25 KB
/
gse123902_cytokines_checkpoint.R
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
library(Seurat)
library(ggplot2)
library(cowplot)
library(scater)
library(scran)
library(BiocParallel)
library(BiocNeighbors)
library(data.table)
library(clustree)
library(dplyr)
setwd("~/gse123902/CytoTRACE")
hms_cluster<-readRDS("n4_cluster_id.rds")
DimPlot(hms_cluster, reduction = "umap", label = TRUE, pt.size = 0.5)
DimPlot(hms_cluster, reduction = "umap", label = FALSE, pt.size = 0.5)
features = c("CXCL15","CXCL14","CXCL12","CXCL10","CXCL9","CXCL5","CXCL3","CXCL2","CXCL1",
"CCL1","CCL5","CCL4","CCL3","CCL2","CCL7",
"IL17A", "IL12A", "IL10", "IL6", "IL4", "IL2","IL1B",
"TNF","IFNG","TGFB1","CSF2","CSF1","CD40","NOS2","ARG1","S100A8","S100A9")
DotPlot(hms_cluster,features=features)+RotatedAxis()
features = c("CXCR4","CXCR3","CXCR2","CCR1", "CCR5", "CCR4", "CCR2",
"IL17RC", "IL17RA", "IL12RB2","IL12RB1", "IL10RB","IL10RA",
"IL6RA","IL4RA","IL2RB","IL2RG","IL2RA",
"TNFRSF1B", "TNFRSF1A","IFNGR2","IFNGR1",
"TGFBR3","TGFBR2" ,"TGFBR1","CSF3R","CSF2RA","CSF1R","IL1R1")
DotPlot(hms_cluster,features=features)+RotatedAxis()
features = c("CD27","CD28","CD40","TNFRSF4","TNFRSF9","TNFRSF18","ICOS","CD266",
"ADORA2A","BTLA","CTLA4","LAG3","PDCD1","HAVCR2","VSIR","TIGIT","CD70",
"CD80","CD86","CD4oLG","TNFSF9","TNFSF4","ICOSL","NECTIN2",
"PVR","LGALS3","CD274","PDCD1IG2","CEACAM1","LGALS9")
DotPlot(hms_cluster,features=features)+RotatedAxis()
setwd("~/gse123902/CytoTRACE")
hms_cluster<-readRDS("p8_cluster_id.rds")
DimPlot(hms_cluster, reduction = "umap", label = TRUE, pt.size = 0.5)
DimPlot(hms_cluster, reduction = "umap", label = FALSE, pt.size = 0.5)
features = c("CXCL15","CXCL14","CXCL12","CXCL10","CXCL9","CXCL5","CXCL3","CXCL2","CXCL1",
"CCL1","CCL5","CCL4","CCL3","CCL2","CCL7",
"IL17A", "IL12A", "IL10", "IL6", "IL4", "IL2","IL1B",
"TNF","IFNG","TGFB1","CSF2","CSF1","CD40","NOS2","ARG1","S100A8","S100A9")
DotPlot(hms_cluster,features=features)+RotatedAxis()
features = c("CXCR4","CXCR3","CXCR2","CCR1", "CCR5", "CCR4", "CCR2",
"IL17RC", "IL17RA", "IL12RB2","IL12RB1", "IL10RB","IL10RA",
"IL6RA","IL4RA","IL2RB","IL2RG","IL2RA",
"TNFRSF1B", "TNFRSF1A","IFNGR2","IFNGR1",
"TGFBR3","TGFBR2" ,"TGFBR1","CSF3R","CSF2RA","CSF1R","IL1R1")
DotPlot(hms_cluster,features=features)+RotatedAxis()
features = c("CD27","CD28","CD40","TNFRSF4","TNFRSF9","TNFRSF18","ICOS","CD266",
"ADORA2A","BTLA","CTLA4","LAG3","PDCD1","HAVCR2","VSIR","TIGIT","CD70",
"CD80","CD86","CD4oLG","TNFSF9","TNFSF4","ICOSL","NECTIN2",
"PVR","LGALS3","CD274","PDCD1IG2","CEACAM1","LGALS9")
DotPlot(hms_cluster,features=features)+RotatedAxis()
setwd("~/gse123902/CytoTRACE")
hms_cluster<-readRDS("m5_cluster_id.rds")
DimPlot(hms_cluster, reduction = "umap", label = TRUE, pt.size = 0.5)
DimPlot(hms_cluster, reduction = "umap", label = FALSE, pt.size = 0.5)
features = c("CXCL15","CXCL14","CXCL12","CXCL10","CXCL9","CXCL5","CXCL3","CXCL2","CXCL1",
"CCL1","CCL5","CCL4","CCL3","CCL2","CCL7",
"IL17A", "IL12A", "IL10", "IL6", "IL4", "IL2","IL1B",
"TNF","IFNG","TGFB1","CSF2","CSF1","CD40","NOS2","ARG1","S100A8","S100A9")
DotPlot(hms_cluster,features=features)+RotatedAxis()
features = c("CXCR4","CXCR3","CXCR2","CCR1", "CCR5", "CCR4", "CCR2",
"IL17RC", "IL17RA", "IL12RB2","IL12RB1", "IL10RB","IL10RA",
"IL6RA","IL4RA","IL2RB","IL2RG","IL2RA",
"TNFRSF1B", "TNFRSF1A","IFNGR2","IFNGR1",
"TGFBR3","TGFBR2" ,"TGFBR1","CSF3R","CSF2RA","CSF1R","IL1R1")
DotPlot(hms_cluster,features=features)+RotatedAxis()
features = c("CD27","CD28","CD40","TNFRSF4","TNFRSF9","TNFRSF18","ICOS","CD266",
"ADORA2A","BTLA","CTLA4","LAG3","PDCD1","HAVCR2","VSIR","TIGIT","CD70",
"CD80","CD86","CD4oLG","TNFSF9","TNFSF4","ICOSL","NECTIN2",
"PVR","LGALS3","CD274","PDCD1IG2","CEACAM1","LGALS9")
DotPlot(hms_cluster,features=features)+RotatedAxis()
#average value
hms_cluster<-readRDS("n4_cluster_id.rds")
DimPlot(hms_cluster, reduction = "umap", label = TRUE, pt.size = 0.5)
DimPlot(hms_cluster, reduction = "umap", label = FALSE, pt.size = 0.5)
features=c("CD274","PDCD1","PDCD1LG2","TIGIT","NECTIN2","NECTIN3","PVR","IL7","IL7R",
"CXCR4","CXCL12","FTL","SCARA5","CCL5","ACKR1","CCR1","CCR5","ACKR2","GPR75",
"LAG3","FGL1","HAVCR2","LGALS9","CCL1","CCL2","CCL3","CCL4","CCL7","CCR2","CCR3",
"CCR8","CD226","CD27","CD28","CD40","CD40LG","CD70","CD80","CD86","CD96","CSF1",
"CSF1R","CSF2","CTLA4","CXCL1","CXCL10","CXCL14","CXCL2","CXCL3","CXCL5","CXCL9",
"CXCR1","CXCR2","CXCR3","DPP4","GMCSFR","ICOS","ICOSLG","IFNG","IFNGR2","KIR2DL5A",
"LGALS3","MERTK","NGFR","P4HB","TNFRSF18","TNFRSF4","TNFRSF9","TNFSF18","TNFSF4","TNFSF9")
a<-DoHeatmap(subset(hms_cluster, downsample = 100), features = features, size = 3)
b<-a$data
write.table(b,"b")
sce <- hms_cluster
table(Idents(sce))
DimPlot(sce,label = T)
av <-AverageExpression(sce , assays = "RNA")
av=av[[1]]
write.table(av,"av")
hms_cluster<-readRDS("p8_cluster_id.rds")
DimPlot(hms_cluster, reduction = "umap", label = TRUE, pt.size = 0.5)
DimPlot(hms_cluster, reduction = "umap", label = FALSE, pt.size = 0.5)
features=c("CD274","PDCD1","PDCD1LG2","TIGIT","NECTIN2","NECTIN3","PVR","IL7","IL7R",
"CXCR4","CXCL12","FTL","SCARA5","CCL5","ACKR1","CCR1","CCR5","ACKR2","GPR75",
"LAG3","FGL1","HAVCR2","LGALS9","CCL1","CCL2","CCL3","CCL4","CCL7","CCR2","CCR3",
"CCR8","CD226","CD27","CD28","CD40","CD40LG","CD70","CD80","CD86","CD96","CSF1",
"CSF1R","CSF2","CTLA4","CXCL1","CXCL10","CXCL14","CXCL2","CXCL3","CXCL5","CXCL9",
"CXCR1","CXCR2","CXCR3","DPP4","GMCSFR","ICOS","ICOSLG","IFNG","IFNGR2","KIR2DL5A",
"LGALS3","MERTK","NGFR","P4HB","TNFRSF18","TNFRSF4","TNFRSF9","TNFSF18","TNFSF4","TNFSF9")
a<-DoHeatmap(subset(hms_cluster, downsample = 100), features = features, size = 3)
b<-a$data
write.table(b,"b")
sce <- hms_cluster
table(Idents(sce))
DimPlot(sce,label = T)
av <-AverageExpression(sce , assays = "RNA")
av=av[[1]]
write.table(av,"av")
hms_cluster<-readRDS("m5_cluster_id.rds")
DimPlot(hms_cluster, reduction = "umap", label = TRUE, pt.size = 0.5)
DimPlot(hms_cluster, reduction = "umap", label = FALSE, pt.size = 0.5)
features=c("CD274","PDCD1","PDCD1LG2","TIGIT","NECTIN2","NECTIN3","PVR","IL7","IL7R",
"CXCR4","CXCL12","FTL","SCARA5","CCL5","ACKR1","CCR1","CCR5","ACKR2","GPR75",
"LAG3","FGL1","HAVCR2","LGALS9","CCL1","CCL2","CCL3","CCL4","CCL7","CCR2","CCR3",
"CCR8","CD226","CD27","CD28","CD40","CD40LG","CD70","CD80","CD86","CD96","CSF1",
"CSF1R","CSF2","CTLA4","CXCL1","CXCL10","CXCL14","CXCL2","CXCL3","CXCL5","CXCL9",
"CXCR1","CXCR2","CXCR3","DPP4","GMCSFR","ICOS","ICOSLG","IFNG","IFNGR2","KIR2DL5A",
"LGALS3","MERTK","NGFR","P4HB","TNFRSF18","TNFRSF4","TNFRSF9","TNFSF18","TNFSF4","TNFSF9")
a<-DoHeatmap(subset(hms_cluster, downsample = 100), features = features, size = 3)
b<-a$data
write.table(b,"b")
sce <- hms_cluster
table(Idents(sce))
DimPlot(sce,label = T)
av <-AverageExpression(sce , assays = "RNA")
av=av[[1]]
write.table(av,"av")