plotAbundanceDensity options

[antagomir]

Can we consider case where plotAbundanceDensitycould be used to visualize a single feature (e.g. Prevotella), as a function of external covariate, such as BMI group (discrete) or BMI (continuous)? Or perhaps this requires its own function, or is too simple to implement separately.

As in

library(miaTime)
data(hitchip1006, package="miaTime")
tse <- hitchip1006
df <- as.data.frame(colData(tse))
df$Prevotella <- assay(tse)["Prevotella oralis et rel.",]
ggplot(df, aes(x = bmi_group, y = Prevotella)) + geom_jitter(width=0.2)