I need to generate a unit curve that is going to look like a right skewed gaussian and I have the following constraints: The X axis is Days (variable but usually 45+) All values on the Y axis sum to 1 The peak will always occur around day 4 or 5 Example: Is there a way to do this programmatically in python?

Reputation: 2112

Python generate random right skewed gaussian with constraints

I need to generate a unit curve that is going to look like a right skewed gaussian and I have the following constraints:

Example:

Is there a way to do this programmatically in python?

Upvotes: 0

Reputation: 16174

as noted by @Severin, a gamma looks to be a reasonable fit. e.g:

import matplotlib.pyplot as plt
import numpy as np
import scipy.stats as sps

x = np.linspace(75)
plt.plot(x, sps.gamma.pdf(x, 4) '.-')
plt.show()

if they really need to sum to 1, rather than integrate, I'd use the cdf and then use np.diff on the result

Upvotes: 1