Teste exato de Fisher

Testes A/B em R

Lauryn Burleigh

Data Scientist

Especificações

Não paramétrico para dados categóricos
- Comeriam de novo: pizza de Cheese ou Pepperoni
Amostra < 1000
Frequência < 5 em pelo menos 20%
Use o teste exato de Fisher para p-valor exato
Tabela de frequência 2x2

Distribuição hipergeométrica

Verossimilhança da tabela de frequências

table(Pizza$Topping, Pizza$EatAgain)

              No Yes
  Cheese     140 360
  Pepperoni   80 420

Hipóteses

Hipótese nula: não há relação entre comer de novo e o tipo de cobertura

Fotografia de fatias de pizza sendo pegadas por várias mãos.

Amostra

library(pwr)
pwr.chisq.test(w = 0.1, df = 2, 
               power = 0.80, 
               sig.level = 0.05)

     Cálculo de poder do qui-quadrado 
              w = 0.1
              N = 963.4689
             df = 2
      sig.level = 0.05
          power = 0.8
NOTA: N é o número de observações

table(Pizza$Topping, Pizza$EatAgain)

              No Yes
  Cheese     140 360
  Pepperoni   80 420

Teste

freqtbl <- table(Pizza$Topping, 
                 Pizza$Time)
fisher.test(freqtbl)

    Fisher's Exact Test for Count Data
dados:  freqtbl
p-value = 0.001809
hipótese alternativa: a odds ratio verdadeira é diferente de 1
intervalo de confiança de 95 por cento:
 0.4569013 0.8417208
estimativas da amostra:
odds ratio 
 0.6209987

Tamanho do efeito

freqtbl <- table(Pizza$Topping, 
                 Pizza$Time)
fisher.test(freqtbl)

Pequeno: 1,5
Médio: 2,5
Grande: 4
Pizza de Pepperoni x 0,62 em relação à de Cheese

Fisher's Exact Test for Count Data
dados:  freqtbl
p-value = 0.001809
hipótese alternativa: true odds 
ratio is not equal to 1
intervalo de confiança de 95 por cento:
 0.4569013 0.8417208
estimativas da amostra:
odds ratio 
 0.6209987

Poder

library(pwr)
pwr.2p.test(h = 0.62, n = 500, 
            sig.level = 0.002)

    Poder para diferença de proporções 
 cálculo para distribuição binomial

              h = 0.62
              n = 500
      sig.level = 0.002
          power = 1
    alternativa = bicaudal

NOTA: mesmos tamanhos de amostra

Vamos praticar!

Testes A/B em R