Marginsplot with hist option

[Nihat Mugurtay]

Dear friends,

I need to create a marginsplot with a histogram on the horizontal ax, reflecting the distribution of X variable.

Here is my code:

logit secondops v2x_polyarchy c.pve##c.agreeus c.loggdp, vce(cluster COWcode), if year>2017 & year<2022
margins, dydx(c.pve) at (agreeus=(0(.1)1))
marginsplot, recastci(rarea)

How can I add a histogram of X variable to the plot? There should be a command like this:

marginsplot, recastci(rarea) hist(c.agreeus)

* I am using STATA 18.