02:31:45	Christian Barreto:	# RP con gumbel-r

fitbest=getattr(scipy.stats,"gumbel_r") # obtenienen atributos

P = fitbest.fit(data)

yrp= np.linspace(0,1000,100001)
rpd=[]

for i in yrp:
    rpd.append(fitbest.ppf([1-(1/i)], *P)[0] )
    
results_rpd=pd.DataFrame()
results_rpd["RP"] = yrp
results_rpd["PP"] = np.array(rpd)
02:34:17	Christian Barreto:	results_rpd[results_rpd["RP"]==50]
02:34:20	Christian Barreto:	results_rpd[results_rpd["RP"]==100]
02:34:24	Christian Barreto:	results_rpd[results_rpd["RP"]==200]
03:07:58	Irving:	gamma
03:49:59	Jorge Poma  Dz11:	gracias
03:50:00	Stefany Amado Menauth DZ6:	Gracias