analysis.r: add sample R script