R/mycltb.R

In reaganrasmussen/GOODMATH4753repo: Houses functions for MATH5743

#' Histogram of a binomial clt
#'
#' @param n sample size
#' @param iter number of iterations
#' @param p probability of a success
#' @param ... take parameters from within the function
#'
#' @return a histogram
#' @export
#'
#' @examples \dontrun{mycltb(n=4,iter=10000,p=0.3)}
mycltb=function(n,iter,p,...){

  rbinom <- hist <- curve <- dnorm <- x <- NULL

  y=rbinom(n*iter,size=n,prob=p)
  data=matrix(y,nrow=n,ncol=iter,byrow=TRUE)
  w=apply(data,2,mean)
  param=hist(w,plot=FALSE)
  ymax=max(param$density)
  ymax=1.1*ymax


  hist(w,freq=FALSE,  ylim=c(0,ymax),
       main=paste("Histogram of sample mean","\n", "sample size= ",n,sep=""),
       xlab="Sample mean",...)
  curve(dnorm(x,mean=n*p,sd=sqrt(p*(1-p))),add=TRUE,col="Red",lty=2,lwd=3)

}