0

The goal is to distribute man-days effort (y axis) through the days time period (x axis) according to defined curve shapes: Bell, FrontLoad, BackLoad, U-shape

These models could be perfectly drawn with Beta distribution and playing with the 'a' and 'b' params to get desired shapes.

My (weak mathematician) problem is to convert scipy beta PDF function result density list (y1, y2, y3) into percentage loads list between 0 and 1 (y-axis loads in each period point x (day#)). Sum of these percentage loads supposed to be qual to 1 (100%). X values is the days and always > 0

The idea then to iterate through y-axis percentage values list and multiply them on total man-days value to get desired man-days effort distribution through the days.

Code snippet looks like this:

import numpy as np
from scipy.stats import beta
import matplotlib.pyplot as plt

total_demand = 100  # total man-days demand
x = np.linspace(0, 1, 100)  # 100 days range. PDF function accepts values in range 0...1
y1 = beta.pdf(x, 2, 8)  # a=2, b=8. FrontLoad distribution
y2 = beta.pdf(x, 5, 5)  # a=5, b=5. Bell-shape distribution
y3 = beta.pdf(x, 8, 2)  # a=8, b=2. BackLoad distribution

y1_demands = [total_demand*y for y in y1]   # does not work because y values are density and not percentage and 
                                            # could be greater than 1

Density result graphs:

Klever
  • 133
  • 8

1 Answers1

1

The instantaneous value of the pdf isn't bounded on [0,1]. Only the integral over any subset is. In this case, you want your subsets to be days, so you'd need to integrate over each day. That's not easy with the pdf, but with the pre-integrated cdf, it's easy - just take the running difference:

y1 = np.diff(beta.cdf(x, 2, 8))  # a=2, b=8. FrontLoad distribution
y2 = np.diff(beta.cdf(x, 5, 5))  # a=5, b=5. Bell-shape distribution
y3 = np.diff(beta.cdf(x, 8, 2))  # a=8, b=2. BackLoad distribution

I can't attach images, but that should give you your expected graphs.

Also this:

y1_demands = [total_demand*y for y in y1]

can just be:

y1_demands = total_demand*y1

numpy is intended to prevent for loops, that are slow in python. if you find yourself iterating through arrays, there's probably a much faster (and usually clearer) way to do it in numpy.

Daniel F
  • 13,620
  • 2
  • 29
  • 55
  • Hi Daniel! Thank you very much! That is exactly what I tried to get! However, would like to clarify next things: - As I got, beta.cdf function to provide reasonable result accepts x only in range 0...1, so I still need to divide 1 on total days value with linespace to gex x-axis values, am I right? - Produced with np.diff list y lists contains x-1 length values. Could you give me some background to understand why it is so? Should I and zero value at the beginning of y lists to have same length for x and y lists? Thank you in advance! – Klever May 05 '22 at 11:42
  • 1
    you should be able to call `beta.cdf(x,#,#,scale = 100)` and have the domain change to 0-100. – Daniel F May 05 '22 at 13:00
  • Hi Daniel, thank you once again! Really appreciated! – Klever May 05 '22 at 18:01