Olá!
Quando utilizei o código abaixo houve uma descrepância no resultado.
from scipy.stats import binom
p = 1 / 6 # Probabilidade de sair numero 5 n = 10 # Total de lançamentos k = [3,4,5,6]
print("{0:.2%}".format(binom.pmf(k, n, p).sum()))
Resultado = 22,45%
Olá Alex, tudo bem? Espero que sim!
Realizei os testes e as funções de probabilidades acumuladas binom.sf() e binom.cdf() retornaram 22.47732...% sendo arredondado para 22.48%. De fato o caso feito de forma manual mostrou uma leve diferença e retornou 22.45056...%.
Mas essas diferenças se devem por questões de arredondamentos e de somas e não vão impactar tanto na resolução.
Bons estudos!