Coded live in my talk at the New York Open Statistical Programming Meetup - View it on GitHub