ЯсноПонятно24

Изначально в коде была допущена ошибка в написании знака "<" вместо "<" в операторах сравнения. Вот исправленный код:

import scipy.stats as st

def z_test_return(x1, x2, y1, y2):
    
    p1 = x1/y1
    p2 = x2/y2
    
    p_combined = (x1 + x2) / (y1 + y2)
    difference = p1 - p2 
    
    z_value = difference / (p_combined * (1 - p_combined) * (1/y1 + 1/y2)) ** 0.5
    
    distr = st.norm(0, 1)
    
    p_value = (1 - distr.cdf(abs(z_value))) * 2
    return p_value

alpha = 0.05

for i in range(0, 4):
    for j in range(2, 3):
        if j < 3:
            k = j + 1
        else:
            k = 1
    
        p_value = z_test_return(popul.iloc[i, j], popul.iloc[i, k], 
                                popul.iloc[4, j], popul.iloc[4, k])
        print('Событие:', popul.event_name.unique()[i])
        print('p-значение:', p_value)

        if p_value < alpha:
            print('Отвергаем нулевую гипотезу: между долями есть разница')
        else:
            print('Не получилось отвергнуть нулевую гипотезу, нет оснований считать доли разными')
        print(' ')

Теперь код должен правильно считать p-value и выводить результаты теста.