Welcome to Axel Michaud's Site

# Data Visualization of the probability mass function 
# (pmf) of a Binomial Random Variable

# range Rx={0,1,2,...,n}
# p is probability of a success for a given trial
n <- 500
p <- 0.3

# dbinom : function for binomial pmf 
# X~binomial(n,p)
# function of n and p
x.vec <- 0:n
prob.vec <- dbinom(x.vec, n, p)

# type h creates vertical lines
plot(x.vec, prob.vec, type = 'h')