Python中的蒙特卡洛模擬

Question

我只是為蒙特卡洛模擬編寫了一個簡單的代碼：

def loss(r, loc, arg, scale, lam):
    X = []
    for x in range(27000): 
        if(r < poisson.cdf(x, lam)):  
            out = 0
        else:
            out = lognorm.rvs(s=arg,loc=loc, scale=scale)
        X.append(out)
    return np.sum(X)  

losses = []
for _ in range(2000):
    r = np.random.random()
    losses.append(loss(r, loc, arg, scale, lam))
E = np.sum(losses)/len(losses)
print(E)

plt.hist(losses, bins='auto')

但是現在，總和僅由對數正態分布的隨機變量組成-是否有可能進行合並，例如2個蒙特卡洛模擬（1個對數正態和1個伽馬）並將其繪制在一個直方圖中？

在此先感謝並致以誠摯的問候

Answer 1

檢查允許您同時作為線程運行函數的線程模塊： https : //www.tutorialspoint.com/python/python_multithreading.htm

要顯示多個圖，可以使用matplotlib的子圖功能： https ://matplotlib.org/gallery/subplots_axes_and_figures/subplot.html