Add Stata shapefile creation code