ComeçarComece de graça

Teste de permutação com dados de sapos

A força média do golpe do Sapo A foi 0,71 newtons (N), e a do Sapo B foi 0,42 N, resultando em uma diferença de 0,29 N. É possível que os sapos golpeiem com a mesma força e que essa diferença observada tenha ocorrido ao acaso. Você vai calcular a probabilidade de obter pelo menos 0,29 N de diferença na força média do golpe sob a hipótese de que as distribuições das forças de golpe dos dois sapos são idênticas. Usaremos um teste de permutação com a estatística de teste sendo a diferença das médias para avaliar essa hipótese.

Para sua conveniência, os dados foram armazenados nos arrays force_a e force_b.

Este exercício faz parte do curso

Pensamento Estatístico em Python (Parte 2)

Ver curso

Instruções do exercício

  • Defina uma função com a assinatura diff_of_means(data_1, data_2) que retorne a diferença entre as médias de dois conjuntos de dados: média de data_1 menos a média de data_2.
  • Use essa função para calcular a diferença empírica das médias observada nos sapos.
  • Gere 10.000 réplicas por permutação da diferença das médias.
  • Calcule o p-valor.
  • Imprima o p-valor.

Exercício interativo prático

Experimente este exercício completando este código de exemplo.

def diff_of_means(data_1, data_2):
    """Difference in means of two arrays."""

    # The difference of means of data_1, data_2: diff
    diff = ____

    return diff

# Compute difference of mean impact force from experiment: empirical_diff_means
empirical_diff_means = ____

# Draw 10,000 permutation replicates: perm_replicates
perm_replicates = draw_perm_reps(____, ____,
                                 ____, size=10000)

# Compute p-value: p
p = np.sum(____ >= ____) / len(____)

# Print the result
print('p-value =', p)
Editar e executar o código