Integrating the conditional density produces the conditional distribution.
Here, for each ordered year, the surface gives the estimated share
of Italian regions with per-capita GDP at or below each value.
One npcdist() call estimates the whole conditional distribution.
Invert the conditional distribution at a chosen probability and you obtain a conditional quantile.
The full script is on Gallery Video Demos.