R code for mark analysis
