data, analysis R script, and output added