Генерация реплик К-С

Теперь вам нужна функция для получения реплик статистики Колмогорова–Смирнова из целевого распределения f. Реализуйте функцию с сигнатурой draw_ks_reps(n, f, args=(), size=10000, n_reps=10000). Здесь n — количество точек данных, а f — функция для генерации выборок из целевой функции распределения (CDF). Например, чтобы проверить соответствие экспоненциальному распределению, передайте np.random.exponential в качестве f. Эта функция, как правило, принимает аргументы, которые передаются в виде кортежа. Так, чтобы получить выборки из экспоненциального распределения со средним x_mean, используйте ключевое слово args=(x_mean,). Ключевые аргументы size и n_reps задают соответственно количество выборок из целевого распределения и количество реплик.

Это упражнение является частью курса

Примеры применения статистического мышления

Посмотреть курс

Инструкции к упражнению

Напишите функцию с сигнатурой draw_ks_reps(n, f, args=(), size=10000, n_reps=10000), которая выполняет следующее.
- Генерирует size выборок из целевого распределения f. Помните: чтобы передать args в функцию выборки, используйте конструкцию f(*args, size=size). Сохраните результат в переменной x_f.
- Инициализирует массив реплик reps как пустой массив с n_reps элементами.
- Выполняет цикл for, повторяя следующие действия n_reps раз.
  - Извлекает n выборок из f. Снова используйте *args в вызове функции. Сохраните результат в переменной x_samp.
  - Вычисляет статистику К-С с помощью dcst.ks_stat() — функции, написанной в предыдущем упражнении и удобно сохранённой в модуле dcst. Сохраните результат в массиве reps.
- Возвращает массив reps.

Интерактивное практическое упражнение

Попробуйте выполнить это упражнение, дополнив этот пример кода.

____ ____:
    # Generate samples from target distribution
    x_f = ____
    
    # Initialize K-S replicates
    reps = ____
    
    # Draw replicates
    for i in range(n_reps):
        # Draw samples for comparison
        x_samp = ____
        
        # Compute K-S statistic
        ____[i] = ____

    return reps

Редактировать и запускать код