added script for plotting rate het data