In the Seurat wrapper-package "Nebulosa," is there a way to split a plot_density() visualization by e.g. treatment groups, like, for example, the "split.by" option in Seurat's FeaturePlot()? I have been looking through the documentation but haven't found anything explicit. If not, how might I (a very basic R user) be able to achieve this? For reference, I have linked the plot_density() documentation here (https://rdrr.io/bioc/Nebulosa/man/plot_density.html).