R code for smith magenis data