Create calculate_probabilities.Rscript