Code from my R videos available at https://vimeo.com/user19035707/videos.