revise README.md
create stan code
create avg_estimate.R
run the avg_estimate.R
modidy avg_estimate
data2013_score.R