Menerapkan koreksi Bonferroni
Setelah mengidentifikasi perbedaan yang signifikan antar kelompok terapi dengan Tukey's HSD, kita ingin menegaskan temuan tersebut menggunakan koreksi Bonferroni. Koreksi Bonferroni adalah penyesuaian statistik yang konservatif untuk mengatasi masalah perbandingan berganda. Metode ini menurunkan peluang terjadinya hasil positif palsu dengan menyesuaikan tingkat signifikansi. Dalam konteks studi Anda tentang efektivitas CBT, DBT, dan ACT, menerapkan koreksi Bonferroni akan membantu memastikan bahwa perbedaan signifikan yang Anda amati antar kelompok terapi tidak terjadi secara kebetulan.
DataFrame therapy_outcomes telah dimuat kembali bersama dengan pandas as pd, from scipy.stats import ttest_ind, dan from statsmodels.sandbox.stats.multicomp import multipletests.
Latihan ini adalah bagian dari kursus
Perancangan Eksperimen dengan Python
Petunjuk latihan
- Lakukan uji t independen untuk semua pasangan kelompok terapi dalam
therapy_pairsdan tambahkan p-value (p_val) ke daftarp_values. - Terapkan koreksi Bonferroni untuk menyesuaikan p-value dari beberapa pengujian dan cetak hasilnya.
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
p_values = []
therapy_pairs = [('CBT', 'DBT'), ('CBT', 'ACT'), ('DBT', 'ACT')]
# Conduct t-tests and collect P-values
for pair in ____:
group1 = therapy_outcomes[therapy_outcomes['Therapy_Type'] == ____]['Anxiety_Reduction']
group2 = therapy_outcomes[therapy_outcomes['Therapy_Type'] == ____]['Anxiety_Reduction']
t_stat, p_val = ____(group1, group2)
p_values.____(p_val)
# Apply Bonferroni correction
print(____(____, alpha=0.05, method='____')[1])