Untitled

from statsmodels.stats.proportion import proportions_ztest

# Dados do problema
p_esperada = 0.8
n_amostra = 1000
sobreviventes = 720
alpha = 0.01

# Teste de proporções
count = sobreviventes
nobs = n_amostra
stat, pval = proportions_ztest(count, nobs, p_esperada)

# Resultado
if pval < alpha:
    resultado = f"Ao nível de 1%, a afirmação das autoridades sanitárias não é verdadeira (p-valor = {pval:.5f})"
else:
    resultado = f"Ao nível de 1%, a afirmação das autoridades sanitárias é verdadeira (p-valor = {pval:.5f})"

print(resultado)