Cara melakukan uji t dua sampel dengan python

Uji-t dua sampel digunakan untuk menguji apakah rata-rata dua populasi sama atau tidak.

Tutorial ini menjelaskan cara melakukan uji-t dua sampel dengan Python.

Contoh: dua contoh uji-t dengan Python

Peneliti ingin mengetahui apakah dua spesies tumbuhan berbeda memiliki rata-rata tinggi yang sama. Untuk mengujinya, mereka mengumpulkan sampel acak sederhana yang terdiri dari 20 tanaman dari setiap spesies.

Ikuti langkah-langkah berikut untuk melakukan uji-t dua sampel untuk menentukan apakah kedua spesies tanaman memiliki tinggi yang sama.

Langkah 1: Buat datanya.

Pertama, kita akan membuat dua tabel untuk menampung pengukuran setiap kelompok yang terdiri dari 20 tanaman:

 import numpy as np

group1 = np.array([14, 15, 15, 16, 13, 8, 14, 17, 16, 14, 19, 20, 21, 15, 15, 16, 16, 13, 14, 12])
group2 = np.array([15, 17, 14, 17, 14, 8, 12, 19, 19, 14, 17, 22, 24, 16, 13, 16, 13, 18, 15, 13])

Langkah 2: Lakukan uji-t dua sampel.

Selanjutnya, kita akan menggunakan fungsi ttest_ind() dari pustaka scipy.stats untuk melakukan uji-t dua contoh, yang menggunakan sintaks berikut:

ttest_ind(a, b, equal_var=Benar)

Emas:

• a : tabel contoh observasi kelompok 1
• b : tabel contoh hasil observasi kelompok 2
• equal_var: Jika benar, lakukan uji t 2 sampel independen standar yang mengasumsikan varian populasi yang sama. Jika salah, lakukan uji-t Welch , yang tidak mengasumsikan variansi populasi sama. Hal ini berlaku secara default.

Sebelum melakukan pengujian, kita perlu memutuskan apakah kita akan mengasumsikan bahwa kedua populasi mempunyai varian yang sama atau tidak. Secara umum, kita dapat berasumsi bahwa suatu populasi mempunyai varians yang sama jika rasio varians sampel terbesar dan varians sampel terkecil kurang dari 4:1.

 #find variance for each group
print(np.var(group1), np.var(group2))

7.73 12.26

Rasio varians sampel terbesar dan varians sampel terkecil adalah 12,26/7,73 = 1,586 , yang kurang dari 4. Artinya kita dapat mengasumsikan bahwa varians populasi adalah sama.

Jadi, kita dapat melanjutkan untuk melakukan uji t dua sampel dengan varian yang sama:

 import scipy.stats as stats

#perform two sample t-test with equal variances
stats.ttest_ind(a=group1, b=group2, equal_var=True)

(statistic=-0.6337, pvalue=0.53005)

Statistik uji-t adalah -0,6337 dan nilai p dua sisi yang sesuai adalah 0,53005 .

Langkah 3: Interpretasikan hasilnya.

Dua asumsi untuk uji t dua sampel ini adalah:

H ₀ : µ ₁ = µ ₂ (rata-rata kedua populasi adalah sama)

H _A : µ ₁ ≠µ ₂ (rata-rata dua populasi tidak sama)

Karena nilai p pengujian kami (0,53005) lebih besar dari alpha = 0,05, kami gagal menolak hipotesis nol pengujian tersebut. Kami tidak memiliki cukup bukti untuk mengatakan bahwa rata-rata tinggi tanaman antara kedua populasi tersebut berbeda.

Sumber daya tambahan

Cara Melakukan Uji-T Satu Sampel dengan Python
Cara Melakukan Uji-T Sampel Berpasangan dengan Python