R heatmap script from Bart 10-24-13