Dear STATA pros,

For my research, using the India NSS data, I'm observing the gap in consumption of major food groups by religion and employing RIF QR counterfactual decomposition methods.

A study by (Srinivasan, Chittur S.; Zanello, Giacomo; Shankar, Bhavani, 2013) https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3729423/ assesses the rural-urban differentials in HAZ scores by first estimating the distributions of HAZ scores separately for rural and urban children in each country using kernel smoothing techniques. I would like to generate figures like Figures 1 and 2 (shown on page 9) which show the cumulative distribution functions for urban and rural HAZ scores in Bangladesh and Nepal as well as aggregate the results of QR-CD analysis.

My question is how do you implement this on STATA? What is the command?

I have plotted the cumulative distribution functions for hindu and non-hindu milk consumption but I am struggling with how to plot the "counterfactual" curve on the same graph.

**plotting cumulative distribution function
cumul milk if hindu==1, gen(x)
cumul milk if hindu==0, gen(y)
stack x milk y milk, into(c milk) wide clear
line x y milk, sort

Thank you for your help.
Samira.