#####Histogramas

##### Ejemplo con faithful

par(mfrow=c(3,1)) 

hist(faithful$waiting,col= "green",probability=T,
     main = "Histograma de tiempo entre erupciones segun Sturges",ylab="Frecuencias Relativas", xlab=" ",nclass=nclass.Sturges)

hist(faithful$waiting,col= "blue",probability=T,,
     main = "Histograma de tiempo entre erupciones segun Scott",ylab="Frecuencias Relativas", xlab=" ",nclass=nclass.scott)

hist(faithful$waiting,col= "grey",probability=T,,
     main = "Histograma de tiempo entre erupciones segun Friedman Diaconis",ylab="Frecuencias Relativas", xlab=" ",nclass=nclass.FD)

####Los 3 dan lo mismo!!!

####Otro ejemplo:

set.seed(34)
x<-rexp(30,1)
 hist(x,col= "green",probability=T, main = "Histograma segun Sturges",ylab="Frecuencias Relativas", breaks="Sturges")
  hist(x,col= "blue",probability=T,main = "Histograma segun Scott",ylab="Frecuencias Relativas", breaks="scott")
 hist(x,col= "orange",probability=T,main = "Histograma segun FD",ylab="Frecuencias Relativas",breaks="fd")



#####Pongo distintos puntos de corte

eru<-faithful$eruptions

min(eru)
max(eru)


#Tomo 1.6 y 5.1
 5.1-1.6

corte1=seq(1.6,5.1,0.7)
corte1

hist(eru,col= "green",probability=T, main = "Histograma con 5 intervalos",ylab="Frecuencias Relativas", breaks=corte1)
 
corte2=seq(1.6,5.1,0.5)

 hist(eru,col= "blue",probability=T,main = "Histograma con 7 intervalos",ylab="Frecuencias Relativas", breaks=corte2)

corte3=seq(1.6,5.1,0.35)
 
hist(eru,col= "orange",probability=T,main = "Histograma con 10 intervalos",ylab="Frecuencias Relativas",breaks=corte3)


### cuando los intervalos son de distinta longitud hace el histograma correcto

 par(mfrow=c(1,1))
corte4=c(1.6,2.6,4,4.5,5.1)
hist(eru,col= "orange",main = "Histograma de erupciones",breaks=corte4)

####Para ver como toma los puntos de corte 

hist(eru,plot=F)


######Densidades

hist(eru,col= "light blue",probability=T, ylim=c(0,0.7),main = "Histograma y curva estimada de densidad de erupciones",ylab="Frecuencias Relativas",xlab="")
 
lines(density (eru))

density (eru)

#Cambiamos el ancho de ventana

lines(density (eru,bw=0.15),col="blue")

lines(density (eru,bw=0.5),col="red")

lines(density (eru,bw=0.07),col="brown")

lines(density (eru,bw=1),col="green")



####Boxplot

boxplot(eru,col=3,main="Boxplot de erupciones")

boxplot(eru,col=3,main="Boxplot de erupciones",horizontal=T)



### Otros ejemplos 

par(mfrow=c(3,2))


hist(x,col= "blue",probability=T,main = "",ylab="Frecuencias Relativas", breaks="scott")
title("Ejemplo de asimetría a derecha")
boxplot(x,col= "blue")


set.seed(64)
z<-rnorm(150)title("Ejemplo de simetría")
hist(z,col= "light blue",probability=T,main = "",ylab="Frecuencias Relativas", breaks="scott")
title("Ejemplo de simetría")
boxplot(z,col= "light blue")


y<--x+5
hist(y,col= "gray",probability=T,main = "",ylab="Frecuencias Relativas", breaks="scott")
title("Ejemplo de asimetría a izquierda")
boxplot(y,col= "gray")



#######Tallo hoja
stem(eru)

##Comparamos con el histograma
par(mfrow=c(1,1))
hist(eru,main="",col="grey")



########Funciones


Mi.Func<-function(x)
        {
          par(mfrow=c(2,1))
         hist(x, breaks = "Sturges",
      probability = T,
      col ="blue",
     main ="Histograma y boxplot",
     xlim = c(min(x)-0.5,max(x)+0.5), ylim = NULL,
     xlab = names(x), ylab="Frecuencias Relativas")
boxplot(x,col="blue")
#Algunas medidas 
Prom<-mean(x,na.rm=T)
Desvio<-sd(x,na.rm=T)
list(Promedio=Prom,DesvíoEstandar=Desvio)
             }

Mi.Func(eru)