Python, SimPy: ¿Cómo generar un valor a partir de una distribución de probabilidad triangular?

Question

Quiero ejecutar una simulación que use como parámetro un valor generado a partir de una distribución de probabilidad triangular con un límite inferior A, modo B y un límite superior C. ¿Cómo puedo generar este valor en Python? ¿Hay algo tan simple como expovariate (lambda) (del azar) para esta distribución o tengo que codificar esto?

Answer 1

Si descarga el paquete NumPy, tiene una función numpy.random.triangular (izquierda, modo, derecha [, tamaño]) que hace exactamente lo que está buscando.

Answer 2

Puesto que, yo era la documentación de azar de Python 2.4 echaba de menos esta:

random.triangular (baja, alta, de modo) ¶ devolver un número aleatorio de punto flotante N tal que bajo < = N < = alto y con el modo especificado entre esos límites. Los límites bajo y alto por defecto son cero y uno. El argumento de modo se establece de manera predeterminada en el punto medio entre los límites, dando una distribución simétrica. Nuevo en la versión 2.6.

Answer 3

Digamos que su distribución no fue manejada por NumPy o la biblioteca estándar de Python.

En situaciones donde el rendimiento no es muy importante, el muestreo de rechazo es un truco útil para obtener sorteos de una distribución que no tiene utilizando una que sí tiene.

Para su distribución triangular, que podría hacer algo como

from random import random, uniform 

def random_triangular(low, high, mode): 
    while True: 
     proposal = uniform(low, high) 
     if proposal < mode: 
      acceptance_prob = (proposal - low)/(mode - low) 
     else: 
      acceptance_prob = (high - proposal)/(high - mode) 
     if random() < acceptance_prob: break 
    return proposal 

Puede representar algunas muestras

pylab.hist([random_triangular(1, 6, 5) for t in range(10000)]) 

para asegurarse de que todo se ve bien.