Code for Probability Mass Function in R Studio Code