R/myhyper.R

In z99jacc/MATH4753ouchavv0029: Math 4753

#' Hypergeometric
#'
#' @param iter
#' @param N
#' @param r
#' @param n
#'
#' @return
#' @export
#'
#' @examples
myhyper=function(iter=100,N=20,r=12,n=5){

  sam.mat=matrix(NA,nr=n,nc=iter, byrow=TRUE)
  succ=c()
  for( i in 1:iter){
    sam.mat[,i]=sample(rep(c(1,0),c(r,N-r)),n,replace=FALSE)
    succ[i]=sum(sam.mat[,i])
  }
  succ.tab=table(factor(succ,levels=0:n))
  barplot(succ.tab/(iter), col=rainbow(n+1), main="HYPERGEOMETRIC simulation", xlab="Number of successes")
  succ.tab/iter
}
myhyper(iter=1000,n=19, N=20,r=12)

dhyper(x=0:19, m=12, n=8, k=19)