Olá pessoal! Para entender melhor o teste qui-quadrado decidi utilizá-lo para verificar a honestidade de um dado de 6 lados. Gostaria de saber se o fiz corretamente, dos pontos de vista da formulação das hipóteses, aplicação do teste e conclusão.
#H0: a probabilidade de cair cada número é a mesma ou p[1] = p[2] = p[3] = p[4] = p[5] = p[6]
#H1: a probabilidade de cair cada número é diferente
import pandas as pd
from scipy.stats import chi
F_Observada_1 = [12, 15, 10, 10, 8, 5]
F_Esperada_1 = [10, 10, 10, 10, 10, 10]
significancia = 0.05
confianca = 1- significancia
n = 60
k = 6
df = k-1
chi_2 = 0
for i in range(k):
chi_2 += (F_Observada_1[i] - F_Esperada_1[i]) ** 2 / F_Esperada_1[i]
chi_2_alpha = chi.ppf(confianca, df) ** 2
chi_2_alpha
print('Chi² = ', chi_2)
print('Chi²_alpha =', chi_2_alpha)
print(chi_2 > chi_2_alpha)
if chi_2 > chi_2_alpha:
print('Rejeitar H0')
else: print('Aceitar H0')
Com a saída:
Chi² = 5.8
Chi²_alpha = 11.07049769351635
False
Aceitar H0
Obrigado!