add R script for plotting
added hwe_exact call