def METHOD_NAME():
    sampling_rate = 100
    ppg = nk.ppg_simulate(
        duration=500,
        sampling_rate=sampling_rate,
        heart_rate=70,
        frequency_modulation=0.025,
        ibi_randomness=0.15,
        drift=0.5,
        motion_amplitude=0.25,
        powerline_amplitude=0.25,
        burst_amplitude=0.5,
        burst_number=3,
        random_state=0,
        show=True,
    )
    # Process the data
    df, info = nk.ppg_process(ppg, sampling_rate=sampling_rate)
    epochs = nk.epochs_create(
        df, events=[0, 15000], sampling_rate=sampling_rate, epochs_end=150
    )
    epochs_ppg_intervals = nk.ppg_intervalrelated(epochs)
    assert "PPG_Rate_Mean" in epochs_ppg_intervals.columns
    ppg_intervals = nk.ppg_intervalrelated(df)
    assert "PPG_Rate_Mean" in ppg_intervals.columns