library(tidyverse)
library(Matrix)
library(cowplot)
DefaultAssay(Macro) <- 'RNA'#使用的数据是RNA中的，不是整合后的数据
M1_features <- list(c('Il23',
'Tnf',
'Cxcl9',
'Cxcl10',
'Cxcl11',
'Cd86',
'Il1a',
'Il1b',
'Il6',
'Irf1',
'Irf5',
'Ccl5',
'Cd40',
'Ido1',
'Kymu',
'Ccr7'))
M2_feature <- list(c('Il4r',
'Ccl4',
'Ccl13',
'Ccl20',
'Ccl17',
'Ccl18',
'Ccl22',
'Ccl24',
'Lyve1',
'Vegfa',
'Vegfb',
'Vegfc',
'Vegfd',
'Egf',
'Ctsa',
'Ctsb',
'Ctsc',
'Ctsd',
'Tgfb1',
'Tgfb2',
'Tgfb3',
'Mmp14',
'Mmp19',
'Mmp9',
'Clec7a',
'Wnt7b',
'Fasl',
'Tnfsf12',
'Tnfsf8',
'Cd276',
'Vtcn1',
'Msr1',
'Fn1',
'Irf4'))
Angiogenesis <- list(c('Ccnd2','Ccne1','Cd44','Cxcr4','E2f3','Edn1','Ezh2','Fgf18','Fgfr1',
'Fyn','Hey1','Itgav','Jag1','Jag2','Mmp9','Notch1','Pdgfa','Ptk2',
'Spp1','Stc1','Tnfaip6','Tymp','Vav2','Vcan','Vegfa'))
Phagocytosis <- list(c('Mrc1','Cd163','Mertk','C1qb'))
Inscore <- AddModuleScore(Macro,
features = Phagocytosis,
ctrl = 100,
name = 'Phagocytosis_score')
colnames(Inscore@meta.data)
VlnPlot(Inscore, features = 'Phagocytosis_score1',pt.size=0,cols =c('#E64B35FF', '#4DBBD5FF', '#00A087FF', '#3C5488FF' ,'#F39B7FFF'))+
#width控制箱体宽度，col控制边框颜色，fill控制填充颜色
geom_boxplot(width=.2,col='black',fill='white')+
NoLegend()
##查看ggsci中的颜色
pal_npg(palette = c('nrc'), alpha = 1)(8)
[1] '#E64B35FF' '#4DBBD5FF' '#00A087FF' '#3C5488FF' '#F39B7FFF' '#8491B4FF' '#91D1C2FF' '#DC0000FF'
