var1<-rnorm(1000,-1,1) var2<-rexp(150,5) var3<-runif(1500) part<-50 h<-(max(var1)-min(var1))/part hist(var1, seq(min(var1),max(var1), by=h), probability=TRUE) #hacemos una función histojuli<-function(x, part) { tam<-length(x) h<-(sort(x)[tam]-sort(x)[1])/part plot(hist(x, seq(sort(x)[1],sort(x)[tam], by=h))) } histojuli(var1, 8) histojuli(var3,30) histojuli(var3,15) par(mfrow=c(1,2))