How to plot aregImpute

How can I plot kernel density plots to test the similarity between the original dataset and the imputed datasets? similar to the mice package densityplot()

In addition, how do I plot add the ecdf for each variables x,y,z below for example:
mi<-aregImpute(~x+y+z)
plot(mi)

This Q would be better for stackoverflow.com or stats.stackexchange.com. But to almost answer your question see the chapter on the Titanic in RMS.