выборка Python из разных дистрибутивов

Question

Я мог бы задать это ранее, но вопрос может быть неясным.Я пытаюсь создать выборочное распределение медианы выборки, проводя 300 повторных случайных выборок, рисуя размер выборки 50 из 3 различных нормальных распределений:

subpop1: mean = 100, std dev = 40 (14 of the 50 sample from subpop1)
subpop2: mean = 200, std dev = 70 (20 of the 50 sample from subpop2)
subpop3: mean = 300, std dev = 80 (16 of the 50 sample from subpop3)

как тогда я могу подойти к этой проблеме?это то, что я сделал до сих пор:

repeat = 300
samplesize_list = [14, 20, 16] ]
std_list = [40, 70, 80]
mean_list = [100, 200, 300]
repeat_median = np.empty(repeat, dtype = float)
for j in range(len(samplesize_list)):
    size = samplesize_list[j]
    for m in range(len(mean_list)):
        mean = mean_list[m]
        for z in range(len(std_list)):
            std = std_list[m]
            for i in range(repeat): 
                sample_data = np.random.normal(mean, std, size)
                repeat_median[i] = np.median(sample_data)
sns.distplot(repeat_median, color = 'blue')
plt.show()

я не уверен, где я ошибся, как и во вводном курсе по питону, мне нужна помощь с моим кодированием!

Answer 1

Я не знаком с графикой, но по математике:

import random
import numpy as np

groups = [
    {'label': 'sub_one', 'mean': 100, 'std_dev': 40, 'size': 14},
    {'label': 'sub_two', 'mean': 200, 'std_dev': 70, 'size': 20},
    {'label': 'sub_three', 'mean': 300, 'std_dev': 80, 'size': 16} 
]

def median(mean, std_dev):
    data = np.random.normal(mean, std_dev)
    get_median = np.median(data)
    return get_median


group_all = []

for i in range(300):
    for i in range(groups[0]['size']):
        group_all.append(median(groups[0]['mean'], groups[0]['std_dev']))

    for i in range(groups[1]['size']):
        group_all.append(median(groups[1]['mean'], groups[1]['std_dev']))

    for i in range(groups[2]['size']):
        group_all.append(median(groups[2]['mean'], groups[2]['std_dev']))

print(len(group_all))

(xenial)vash@localhost:~/python/stack_overflow$ python3.7 median.py 
15000