Z calculado negativo

Ei, Thiago! Tudo bem?

Fiz o teste por aqui, o seu código está correto aparentemente e não consegui identificar exatamente o problema nele. Diante disso, recomendo que revise o código anterior se está certinho como o da aula, para termos certeza que não afetará no resultado:

#Seleção das amostras 
homens = dados.query('Sexo == 0').sample(n = 500, random_state = 101).Renda
mulheres = dados.query('Sexo == 1').sample(n = 500, random_state = 101).Renda

#Dados do problema
media_amostra_M = mulheres.mean()
media_amostra_M

desvio_padrao_amostra_M = mulheres.std()
desvio_padrao_amostra_M

media_amostra_H = homens.mean()
media_amostra_H

desvio_padrao_amostra_H = homens.std()
desvio_padrao_amostra_H

significancia = 0.01
confianca = 1 - significancia
n_M = 500
n_H = 500
D_0 = 0


#fixação da significância do teste

probabilidade = confianca
probabilidade

z_alpha = norm.ppf(probabilidade)
z_alpha.round(2)

numerador = (media_amostra_H - media_amostra_M) - D_0

denominador = np.sqrt((desvio_padrao_amostra_H ** 2 / 500) + (desvio_padrao_amostra_M ** 2 / 500))

z = numerador / denominador

z

Após revisar, reinicie o notebook para ver se o retorno está ok. Caso persista, compartilhe o link do seu notebook para que eu possa ter a mesma visibilidade que você, assim conseguimos buscar outras soluções.

Até mais, Thiago!