selected_genes <- c (
"Insr" , "Meis1" , "Igf1r" , "Nrf1" , "Prlr" , "Camk1d" , "Lars2" , "Cacna2d3" , # 0 ARC
"Dcn" , "Ptgds" , "Nupr1" , "Igfbp2" , "Nrarp" , "Ctnna2" , "Ldlr" , "Mmp14" , # 1 LHA
"Nkx6-2" , "Cirbp" , "Selenop" , # 2 MnPO
"Foxg1" , "Crym" , "Sema3c" , "Meis2" , "Dbp" , "Drd1" , "Drd2" , "Drd3" , "Drd4" , "Drd5" , # 3 POA
"Egr1" , "Ttr" , "Btg2" , "Mbnl3" , "Pgf" , "Irs4" , "Gpr101" , "Nr3c2" , "Agtr1" , # 4 PVN
"Rfx4" , "Dbx2" , "Prokr2" , "Cebpb" , "Zic1" , "Zic2" , "Ccn1" , "Gata4" , "Klf10" , # 5 SCN
"Tbx3" , "Fis1" , "Ndn" , "Psmc5" , "Drap1" , "Pcsk1n" , "Rtn1" , # 6 VMH
"Emx2" , "Sgip1" , "Myoc" , "Hivep3" , "Dcc" , "Ltbp1" , "Ralyl" , "Egfem1" , # 7 VPH
"Klf4" , "Atf3" , "Nrg1" , "Cdk8" , "Pitx2" , "Otx1" , "Otx2" , "Peg3" , "Isl1" , "Myt1l" , "Nr5a1" , "Otp" , "Six6" , "Nkx2-1" , "Qrfpr" , "Hcrtr1" , "Hcrtr2" , "Tacr1" , "Trhr" , "Tshr" , "Grpr" , "Npy1r" , "Npy2r" , "Npy5r" ,
"Gfap" , "Fgf1" , "Fgfr3" , "Hepacam" , "Hif1" , "Htra1" , "Lxn" , "Ndrg2" , "Ntn1" , "Nfia" , "Slit2" , "Aqp4" , "S100a1" , "S100a6" , "S100b" , "Slc1a2" , "Slc1a3" , "Slc38a1" , "Vegfa" , "Fos" , "Fosb" , "Jun" , "Junb" , "Jund" , "Ier2" , "Socs3" , "Pde10a" , "Fbln5"
)
# Iterate_FeaturePlot_scCustom(
# seurat_object = srt,
# gene_list = selected_genes,
# reduction = "densmap_xgboost",
# na_cutoff = NA,
# pt.size = 4,
# order = TRUE,
# raster = TRUE,
# raster.dpi = c(960, 960),
# alpha_na_exp = 0.1,
# alpha_exp = 0.75,
# single_pdf = FALSE,
# file_path = plots_dir,
# file_type = ".pdf",
# file_name = sprintf("%s_selected", bioproject),
# colors_use = srt@misc$expr_Colour_Palette_gnuplot_r
# )
# Iterate_FeaturePlot_scCustom(
# seurat_object = srt,
# gene_list = gene_int[gene_int %in% all.genes],
# reduction = "densmap_xgboost",
# na_cutoff = NA,
# pt.size = 4,
# order = TRUE,
# raster = TRUE,
# raster.dpi = c(960, 960),
# alpha_na_exp = 0.1,
# alpha_exp = 0.75,
# single_pdf = FALSE,
# file_path = plots_dir,
# file_type = ".pdf",
# file_name = sprintf("%s_selected", bioproject),
# colors_use = srt@misc$expr_Colour_Palette_gnuplot_r
# )