binomialYnormal=function(tamano=30,ic=F){ binomial=1;poisson=0;normalbinomial=1;normalpoisson=0 altura=dbinom(0,tamano,1/tamano) for(i in seq(1,tamano-1,by=1)){ #altura=dbinom(i,tamano,i/tamano) posicion=tamano;lado=2 if(i>tamano/2) {posicion=0;lado=4} .x=0:tamano .xcont=seq(0,tamano,by=.25) .ybinom=dbinom(.x,tamano,i/tamano) .ypoi=dpois(.x,i) plot(.x,.ybinom,type="h",lwd=3,frame.plot=F,ylab="Probabilidad",xlab="",col="black",ylim=c(0,altura)) if(poisson==1) points(.x,.ypoi,type="h",col="red",lwd=2,lty=1) if(normalbinomial==1) lines(.xcont,dnorm(.xcont,mean=i,sd=sqrt(i*(1-i/tamano))),col="gray",lwd=3,lty=1) if(normalpoisson==1) lines(.xcont+.15,dnorm(.xcont,mean=i,sd=sqrt(i)),col="blue") text(posicion,altura*.9,bquote(paste("Binomial(n=",.(signif(tamano,3)),", p=",.(signif(i/tamano,2)),")")),pos=lado) if(normalbinomial==1) text(posicion,altura*.7,bquote(paste("Normal(",mu,"=np=",.(signif(i,3)),", ",sigma^2,"=npq=",.(signif(i*(1-i/tamano),2)),")")),pos=lado) if(poisson==1) text(posicion,altura*.7,bquote(paste("Poisson(",mu,"=np=",.(signif(i,2)),")")),pos=lado,col="red") if(normalbinomial==1) points(.x,.ybinom,type="h",lwd=1,frame.plot=F,ylab="Probabilidad",xlab="",col="black") readline()} } binomialYpoisson=function(tamano=30){ altura=dbinom(tamano,2*tamano,1/2) for(i in c(2,4,5,10,20,40,50,100,200,500,1000)){ n=tamano*i p=1/i posicion=3*tamano;lado=2 #{posicion=0;lado=4} .x=0:(3*tamano) .ybinom=dbinom(.x,n,p) .ypoi=dpois(.x,tamano) #plot(.x,altura*.x/length(.x),col="white") plot(.x,.ybinom,type="h",lwd=5,frame.plot=F,ylab="Probabilidad",xlab="",col="black",ylim=c(0,altura)) points(.x,.ypoi,type="h",col="red",lwd=3,lty=1) text(posicion,altura*.9,bquote(paste("Binomial(n=",.(signif(n,3)),", p=",.(signif(tamano,3)),"/",.(signif(n,2)),"=",.(signif(p,3)),")")),pos=lado) text(posicion,altura*.7,bquote(paste("Poisson(",mu,"=",sigma^2,"=np=",.(signif(tamano,3)),")")),pos=lado,col="red") readline()} }