###############################################################################
# CARGAMOS LAS LIBRERIAS DE FUNCIONES HECHAS QUE VAMOS A USAR
###############################################################################
import numpy as np #PARA HACER CUENTAS
import xlrd #PARA CARGAR ARCHIVOS DE EXCEL
import matplotlib.pyplot as plt #PARA GRAFICAR
import scipy as sc #PARA ANALISIS DE DATOS
plt.close('all') #CIERRA CUALQUIER GRAFICO QUE HAYA ABIERTO

###############################################################################
# DEFINIMOS NOSOTROS OTRAS DE LAS FUNCIONES QUE VAYAMOS A USAR
###############################################################################
# FUNCION PARA ABRIR UNA COLUMNA DE UN ARCHIVO DE EXCEL
def abrir_columna_excel(nombre_archivo,columna):
  wb = xlrd.open_workbook(nombre_archivo) 
  Hoja1 = wb.sheet_by_index(0) 
  N_datos=len(Hoja1.col_values(columna))
  datos=[]
  for i in range(1,N_datos):
    if type(Hoja1.cell_value(i, columna))==float:
      datos.append(Hoja1.cell_value(i, columna))
  return datos

###############################################################################
# ESCRIBIMOS LA SECUENCIA DE PROCESOS DE NUESTRO PROGRAMA
###############################################################################

N_datos=1000
datos=np.random.normal(10,1,N_datos)
unos=np.ones(len(datos))

mediana=np.median(datos)

plt.plot(datos,unos,'o')
plt.plot([mediana,mediana],[0.9,1.1])

Q1=np.percentile(datos,25)
Q2=np.percentile(datos,50)
Q3=np.percentile(datos,75)

RIC=Q3-Q1
lim_inf=Q1-1.5*RIC
lim_sup=Q3+1.5*RIC

plt.figure()
plt.plot([Q1,Q1],[0.9,1.1])
plt.plot([Q2,Q2],[0.9,1.1])
plt.plot([Q3,Q3],[0.9,1.1])
plt.plot([lim_inf,lim_inf],[0.9,1.1])
plt.plot([lim_sup,lim_sup],[0.9,1.1])

plt.figure()
plt.boxplot(datos,vert=False)