Data 605 -Week 9 Discussion

trials_p <- runif(10000)
head(trials_p)

## [1] 0.05845939 0.08356704 0.49671653 0.61431857 0.04468008 0.28490501

length(trials_p)

## [1] 10000

min(trials_p)

## [1] 0.000385677

max(trials_p)

## [1] 0.9999889

trials_p_100 <- function(){
  p_vec = vector()
  for(i in 1:100){
    trials_p <- runif(10000)
    success = 0
    for(trial in trials_p){
      if(trial <= 0.3){
      success = success + 1
      }
    }
    p_vec[i] = success/length(trials_p)
  }
return(p_vec)
}
hist(trials_p_100(),border="red")