a simple function to simulate geometric distribution for a given p of success and number of trials, producing a ggplot bar plot of the distribution

bernoulli_function.R

require(ggplot2)
#bernoulli
bernoulli_trial <- function(p_succ,n_trials){ # probability of success (1) and number of trials
  p_insucc     <-  1-p_succ # rpobability of failure
  trials       <<- c() 
  p_cumulative <<- c()
  for(i in 1:n_trials){
    trials       <<- rbind(trials,i)
    p_cumulative <<- rbind(p_cumulative,(((p_insucc)^(i -1))*p_succ)) # (p of success within n trials) = ((1-p)^n-1)*(p)
  }
  results <<- data.frame(trials,p_cumulative)
  g <<- ggplot(results,aes(x=trials,y=p_cumulative))+
    geom_bar(stat ='identity')
  g
}