feature_schex <- function(seurat_object,
                          gene_ids,
                          label = FALSE,
                          nbins = 150,
                          reduction = "UMAP",
                          assay = "SCT",
                          slot = "data",
                          action = "mean"){
                          
  require(Seurat)
  require(schex)
  
  seurat_object <- make_hexbin(seurat_object, nbins = nbins, dimension_reduction = reduction)
  
  plot_list <- list()
  
  for(gene in gene_ids){
    print(gene)
    schex_plot <- plot_hexbin_feature(seurat_object, type=slot, feature=gene, 
                        action= action, xlab="UMAP 1", ylab="UMAP 2", mod = assay, 
                        title= gene) +
      theme_cowplot() +
      theme(legend.position = "bottom")
    plot_list[[gene]] <- schex_plot
  }
  
  return(plot_list)
}