DotPlot(subset_data,features =
c(
#b cells c0 c16
'Cd79a',
'Cd79b',
#macrophage c1
'Chil3',
'Lpl',
'Ccl6',
'Abcg1',
'Ear2',
#t cell c2
'Cd3d',
'Cd8b1',
'Cd3g',
#fibro c3
'Mgp',
'Nbl1',
'Sparcl1',
'Cxcl14',
'Igfbp7',
'Bgn',
"Gsn",
#fibro c4
'Npnt',
'Gpx3',
'Inmt',
'Mfap4',
'G0s2',
"Ogn",
#fibro c5
'Hhip',
'Mustn1',
'Aspn',
'Enpp2',
'Igfbp3',
'Grem2',
'Igfbp5',
'Lum',
'Tagln',
#fibro c6
'Dcn',
'Clec3b',
'Igfbp4',
'Col1a2',
'Mmp3',
'Cygb',
'Aebp1',
'Gsn',
'Cxcl12',
'Pi16',
'Ccl11',
'Serpinf1',
'Fbln1',
'Serping1',
'Timp3',
'Col14a1',
#endo c7
'Ptprb',
'Tmem100',
'Epas1',
'Calcrl',
'Egfl7',
'Cldn5',
'Cdh5',
#macro/mono c8
'F13a1',
'S100a4',
'Ms4a6c',
'Ccr2',
#macrohpage 12
'C1qa',
'C1qb',
'C1qc',
'Pf4',
'Ms4a7',
#t cell c9
'Cd3g',
'Cd3e',
#nk c10
'Gzma',
'Ccl5',
'Nkg7',
# at c17
'Sftpa1',
'Sftpb',
'Sftpd'
) %>%unique()
) +RotatedAxis()