(Não consegui postar no tópico existente)
Várias vezes foi falado que o zscore (ou stats) é quantas vezes o desvio padrão que a média amostral está da média esperada, mas o correto é o ERRO PADRÃO, o que fica óbvio considerando a própria fórmula do zscore -- o desvio padrão dividido pela raíz quadrada do tamanho da amostra é exatamente o erro padrão. Persistindo a confusão gerada logo no começo sobre erro e desvio padrão.
Uma condição para usar o teste Z é conhecer o valor de desvio padrão da população. Mas mesmo que informado, a função ztest do statsmodels.stats.weightstats não pede esse valor, então não tem a saída correta. Fazendo o cálculo manual conseguimos o zscore correto.