Mostrar código-fonte de Teste exato de Fisher

{{Estatística sidebar}}
O '''teste exato de Fisher''' é um teste de [[significância estatística]] utilizado na análise de [[Tabela de contingência|tabelas de contingência]].<ref>{{Citar periódico|ultimo=Fisher|primeiro=R. A.|ano=1922|titulo=On the interpretation of χ<sup>2</sup> from contingency tables, and the calculation of P|periodico=Journal of the Royal Statistical Society|volume=85|numero=1|paginas=87–94|doi=10.2307/2340521|jstor=2340521}}</ref><ref>{{Citar livro|ultimo1=Fisher|primeiro1=R.A.|ano=1954|titulo=Statistical Methods for Research Workers|editora=Oliver and Boyd|isbn=0-05-002170-2}}</ref><ref>{{Citar periódico|ultimo=Agresti|primeiro=Alan|ano=1992|titulo=A Survey of Exact Inference for Contingency Tables|periodico=Statistical Science|volume=7|numero=1|paginas=131–153|doi=10.1214/ss/1177011454|jstor=2246001}}</ref> Embora na prática ele seja empregado quando os tamanhos das [[Amostra (estatística)|amostras]] são pequenos, é válido para todos os tamanhos de amostra. É nomeado em homenagem a seu inventor, [[Ronald Fisher]], e é um de uma classe de [[testes exatos]], assim chamados por conta da significância do desvio de uma [[hipótese nula]] (e.g., [[Valor-p|p-valor]]) que pode ser calculada exatamente, ao invés de depender de uma aproximação que se torna exata no limite conforme o tamanho da amostra cresce para o infinito, como em muitos testes estatísticos.

Fisher disse ter concebido o teste depois de um comentário da [[Muriel Bristol|Dra. Muriel Bristol]], que afirmou ser capaz de detectar se o chá ou o leite foi adicionado primeiro em sua xícara. Ele testou seu pedido no experimento "[[dama apreciadora de chá]]".<ref name="newman">{{Citar livro|primeiro=Sir Ronald A.|ultimo=Fisher|capitulo=Mathematics of a Lady Tasting Tea|anooriginal=1935|ano=1956|titulo=The World of Mathematics, volume 3|ultimo2=Newman|primeiro2=James Roy|url=https://books.google.com/?id=oKZwtLQTmNAC&pg=PA1512&dq=%22mathematics+of+a+lady+tasting+tea%22|editora=Courier Dover Publications|isbn=978-0-486-41151-4}}</ref>

== Propósito ==
[[Ficheiro:Nice_Cup_of_Tea.jpg|thumb|256x256px|Um [[Bule|Bule de chá]], um ''creamer'' e uma [[xícara]] cheia de chá com [[leite]] - pode um degustador dizer se o leite foi colocado primeiro?]]
O teste é útil para [[Variável categórica|dados categóricos]], que resultam de classificação de objetos em duas maneiras diferentes; ele é usado para examinar a significância da associação (contingência) entre os dois tipos de classificação. Assim, no exemplo original de Fisher, um critério de classificação poderia ser se o leite ou chá foi colocado na xícara primeiro; o outro poderia ser se a Dra. Bristol pensava que o leite ou chá, fora colocado em primeiro lugar. Queremos saber se essas duas classificações são associados - isto é, se a Dra. Bristol realmente poderia dizer se o leite ou o chá fora servido em primeiro lugar. A maioria dos usos do teste exato de Fisher envolvem, como neste exemplo, uma tabela de contingência 2×2. O p-valor do teste é calculado como se as margens da tabela fossem fixas, isto é, como se no exemplo da degustação de chá, a Dra. Bristol soubesse o número de xícaras com cada tratamento (leite ou chá primeiro) que havia e poderia, portanto, fornecer estimativas com o número correto em cada categoria. Como apontado por Fisher, isso leva, sob a hipótese nula de independência a uma [[distribuição hipergeométrica]] dos números nas células da tabela.

Com grandes amostras, um [[Teste Qui-Quadrado|teste qui-quadrado]] pode ser utilizado nesta situação. No entanto, o valor de significância que ele oferece é apenas uma aproximação, pois a [[distribuição amostral]] da estatística de teste é calculada somente aproximadamente igual a teórica da distribuição qui-quadrado teórica. A aproximação é inadequada quando o tamanho das amostras é pequena, ou se os dados são muito desigualmente distribuídos entre as células da tabela, resultando na contagem de células previstas na hipótese nula (os "valores esperados") ser baixa. A regra de ouro usual para decidir se o teste de aproximação da qui-quadrado é bom o suficiente é que o teste qui-quadrado não é adequado quando os valores esperados nas células da tabela de contingência estão abaixo de 5, ou abaixo de 10 quando há apenas um grau de liberdade (esta regra é agora conhecida por ser excessivamente conservadora).<ref name="Larntz1978">{{citar periódico|doi=10.2307/2286650|ultimo=Larntz|primeiro=Kinley|ano=1978|titulo=Small-sample comparisons of exact levels for chi-squared goodness-of-fit statistics|periodico=Journal of the American Statistical Association|volume=73|numero=362|paginas=253–263|jstor=2286650}}</ref> Na verdade, para dados pequenos, esparsos, ou não balanceados, o p-valor exato e assintótico podem ser muito diferentes e podem levar a conclusões opostas sobre a hipótese de interesse.<ref name="Mehta1984">{{Citar periódico|ultimo=Mehta|primeiro=Cyrus R|ultimo2=Patel|primeiro2=Nitin R|ultimo3=Tsiatis|primeiro3=Anastasios A|ano=1984|titulo=Exact significance testing to establish treatment equivalence with ordered categorical data|periodico=Biometrics|volume=40|numero=3|paginas=819–825|doi=10.2307/2530927|pmid=6518249|jstor=2530927}}</ref><ref name="Mehta1995">Mehta, C. R. 1995. SPSS 6.1 Exact test for Windows. Englewood Cliffs, NJ: Prentice Hall.</ref> Em contraste, o teste exato de Fisher é, como seu nome indica, exato, conforme o procedimento experimental mantém os totais das linhas e colunas fixos, e pode, portanto, ser utilizado independentemente das características da amostra. Torna-se difícil calcular com amostras grandes ou tabelas bem equilibradas, mas felizmente esses são exatamente as condições em que o teste qui-quadrado é apropriado.

Para os cálculos à mão, o teste só é viável no caso de uma tabela de contingências 2×2. No entanto, o princípio do teste pode ser estendido para o caso geral de uma tabela m×n,<ref>{{citar periódico|ultimo=Mehta|primeiro=C.R.|ultimo2=Patel|primeiro2=N.R.|ano=1983|titulo=A Network Algorithm for Performing Fisher's Exact Test in ''r ''X''c'' Contingency Tables|url=|periodico=Journal of the American Statistical Association|volume=78|numero=382|paginas=427–434|doi=10.2307/2288652}}</ref><ref>[http://mathworld.wolfram.com/FishersExactTest.html mathworld.wolfram.com] Página dando a fórmula para a forma geral do teste de Fisher para tabelas de contingência ''m''&nbsp;×&nbsp;''n''</ref> e alguns pacotes estatísticos fornecem um cálculo (às vezes usando um [[método de Monte Carlo]] para obter uma aproximação) para o caso mais geral.<ref>{{citar periódico|primeiro=Cyrus R.|ultimo=Mehta|primeiro2=Nitin R.|ultimo2=Patel|titulo=ALGORITHM 643: FEXACT: a FORTRAN subroutine for Fisher's exact test on unordered r×c contingency tables|periodico=ACM Trans. Math. Softw.|volume=12|numero=2|ano=1986|paginas=154–161|doi=10.1145/6497.214326}}</ref>

== Exemplo ==
Por exemplo, uma amostra de adolescentes pode ser dividida em masculino e feminino por um lado, e aqueles que estão e não estão atualmente a estudar para um exame de estatística no outro. Vamos supor, por exemplo, que a proporção de pessoas estudando é maior entre mulheres do que entre os homens, e queremos testar se a diferença de proporções que observamos é significativa. Os dados podem parecer como a tabela abaixo:
{| class="wikitable"
!
!Homens
!Mulheres
!Total da linha
|-
!Estudiosos
|1
|9
|10
|-
!Não estudiosos
|11
|3
|14
|-
!Total da coluna
|12
|12
|24
|}
A pergunta que se faz sobre esses dados é: sabendo que 10 destes 24 adolescentes são estudiosos, e que 12 dos 24 são do sexo feminino, e supondo que a hipótese nula de que homens e mulheres têm a mesma probabilidade de estudar, qual é a probabilidade de que esses 10 estudiosos seria tão desigualmente distribuídos entre as mulheres e os homens? Se tivéssemos que escolher 10 dos adolescentes ao acaso, qual a probabilidade de que 9 ou mais deles estarem entre as 12 mulheres, e apenas 1 ou menos estarem entre os 12 homens?

Antes de prosseguir com o teste de Fisher, primeiro apresentamos alguns apontamentos. Nós representamos as células pelas letras a, b, c e d, chame os totais das linhas e colunas de totais marginais, e represente o total por n. Assim, a tabela agora tem esse aspecto:
{| class="wikitable"
!
!Homens
!Mulheres
!Total da linha
|-
!Estudiosos
|'''a'''
|'''b'''
|a+b
|-
!Não estudiosos
|'''c'''
|'''d'''
|c+d
|-
!Total da coluna
|a+c
|b+d
|a+b+c+d = n
|}
Fisher mostrou que a probabilidade de obtenção de tais valores é dada pela distribuição hipergeométrica:

<math>p = \frac{ \displaystyle{{a+b}\choose{a}} \displaystyle{{c+d}\choose{c}} }{ \displaystyle{{n}\choose{a+c}} } = \frac{(a+b)!~(c+d)!~(a+c)!~(b+d)!}{a!~~b!~~c!~~d!~~n!}</math>

onde <math> \tbinom nk </math> é o [[coeficiente binomial]] e o símbolo <math>!</math> indica o [[Fatorial|operador fatorial]]. Com os dados acima, isso nos dá:

<math>p = \frac{ \displaystyle{{10}\choose{1}} \displaystyle{{14}\choose{11}} }{ \displaystyle{{24}\choose{12}} }= \frac{10!~14!~12!~12!}{1!~9!~11!~3!~24!} \approx 0.001346076</math>

A fórmula acima dá a probabilidade hipergeométrica exata de observar esta disposição específica dos dados, assumindo os dados marginais totais, na hipótese nula de que homens e mulheres têm a mesma probabilidade de ser estudiosos. Para colocar de outra forma, se nós assumimos que a probabilidade de que um homem é um estudioso é P, a probabilidade de que uma mulher é uma estudiosa é p, e assumimos que os homens e as mulheres da nossa amostra, independentemente de se são ou não estudiosos, então esta fórmula hipergeométrica dá a [[Probabilidade condicionada|probabilidade condicional]] de se observar os valores de <math>a</math>, <math>b</math>, <math>c</math>, <math>d</math>, em quatro células, condicionalmente sobre o observado marginais (i.e., assumindo os totais das linhas e colunas mostrados nas margens da tabela). Isso permanece verdadeiro mesmo se os homens estão em nosso exemplo com diferentes probabilidades do que as mulheres. A exigência é apenas que as duas características de classificação - gênero e estudo (ou não) - não sejam associados.

Por exemplo, suponha que soubessemos as probabilidades <math>P, Q, p</math> e <math>q</math> com <math>P + Q = p + q = 1</math> tais que homens estudiosos, homens não estudiosos, mulheres estudiosas e mulheres não estudiosas tenham suas probabilidades respectivas <math>(Pp, Pq, Qp</math> e <math>Qq)</math> para cada indivíduo encontrado sob o processo de amostragem. Então se alguém for calcular as distribuições condicionais das entradas das células dadas as marginais, obteria a fórmula acima onde nem <math>p</math> ou <math>P</math> ocorram. Portanto, podemos calcular a probabilidade exata de qualquer arranjo dos 24 adolescentes nas quatro células da tabela, mas Fisher mostrou que para gerar um nível de significância, devemos considerar somente os casos onde os totais marginais são iguais aos observados na tabela, e dentre esses, somente os casos onde o arranjo é tão ou mais extremo como o arranjo observado (teste de Barnard relaxa essa restrição em um conjunto de totais marginais). Neste exemplo, existem 11 casos dessa forma. Desses, apenas um é mais extremo na mesma direção que nossos dados; Se apresenta assim:
{| class="wikitable"
!
!Homens
!Mulheres
!Total da linha
|-
!Estudiosos
|0
|10
|10
|-
!Não estudiosos
|12
|2
|14
|-
!Total da coluna
|12
|12
|24
|}
Para essa tabela (com proporções extremamente desiguais) a probabilidade é:

<math>p = \frac{ \displaystyle{{10}\choose{0}} \displaystyle{{14}\choose{12}} }{ \displaystyle{{24}\choose{12}} }= \frac{10!~14!~12!~12!}{0!~10!~12!~2!~24!} \approx 0.000033652</math>.

Para calcular a significância do dado observado, isto é, a probabilidade total de observar os dados como extremos ou mais extremos de a hipótese nula é verdadeira, temos que calcular os valores de <math>p</math> para essas duas tabelas e então somá-los. Isto nos dá um [[teste unicaudal]], com <math>p \approx 0.001346076 + 0.000033652 = 0.001379728 </math>. Por exemplo, no programa [[R (linguagem de programação)|R]], este valor pode ser obtido por <code>fisher.test(rbind(c(1,9),c(11,3)), alternative="less")$p.value</code>. Este valor pode ser interpretado como a soma das evidências providas pelos dados observados - ou qualquer tabela extrema - para a hipótese nula (de que não há diferença nas proporções de estudiosos entre homens e mulheres). Quanto menor o valor de <math>p</math>, maior a evidência para rejeitar a hipótese nula; Então aqui a evidência é forte de que homens e mulheres não são igualmente prováveis de serem estudiosos.

Para [[Teste bicaudal|testes bicaudais]] devemos considerar também tabelas que são igualmente extremas, mas na direção oposta. Infelizmente, classificação de tabelas de acordo como se são ou não "tão extremas quanto" é problemática. Uma abordagem usada pela função<code>fisher.test</code> em [[R (linguagem de programação)|R]] é computar o p valor assumindo as probabilidades para todas as tabelas menores ou iguais aquela observada. No exemplo, o p-valor bilateral é duas vezes o p-valor unilateral - mas no geral pode diferir substancialmente para tabelas com quantidades pequenas, exceto para o caso de testes estatísticos que tem uma distribuição amostral simétrica.

== Controvérsias ==
Apesar do fato de que o teste de Fisher dá os p-valores exatos, alguns autores têm argumentado que o teste é conservador, isto é, que a verdadeira taxa de rejeição é menor que o nível de significância.<ref name="Liddell-1976">{{Citar periódico|doi=10.2307/2988087|ultimo=Liddell|primeiro=Douglas|ano=1976|titulo=Practical tests of 2x2 contingency tables|periodico=The Statistician|volume=25|numero=4|paginas=295–304|jstor=2988087}}</ref><ref name="Berkson1978">{{citar periódico|ultimo=Berkson|primeiro=Joseph|ano=1978|titulo=In dispraise of the exact test|periodico=Journal of Statistic Planning and Inference|volume=2|paginas=27–42|doi=10.1016/0378-3758(78)90019-8}}</ref><ref name="DAgostino1988">{{citar periódico|doi=10.2307/2685002|ultimo=D'Agostino|primeiro=R. B.|ultimo2=Chase|primeiro2= W.|ultimo3=Belanger|primeiro3=A.|ano=1988|titulo=The Appropriateness  of  Some Common Procedures  for Testing Equality  of Two Independent  Binomial  Proportions|periodico=The  American  Statistician|volume=42|numero=3|paginas=198–202|jstor=2685002}}</ref> A aparente contradição decorre da combinação de estatística discreta com níveis de significância fixados.<ref name="Yates1984">{{citar periódico|doi=10.2307/2981577|ultimo=Yates|primeiro=F.|ano=1984|titulo=Tests of Significance  for 2  x  2 Contingency  Tables (with discussion)|periodico=Journal of the Royal Statistical Society, Series A|volume=147|numero=3|paginas=426–463|jstor=2981577}}</ref><ref name="Little1989">{{citar periódico|doi=10.2307/2685390|ultimo=Roderick|primeiro=J. A. Little|ano=1989|titulo=Testing the Equality of Two Independent Binomial Proportions|periodico=The  American  Statistician|volume=43|numero=4|paginas=283–288|jstor=2685390}}</ref> Para ser mais preciso, considere a seguinte proposta para o teste de significância no nível de 5%: rejeite a hipótese nula para cada tabela para as quais o teste de Fisher atribui um p-valor menor ou igual a a 5%. Porque o conjunto de todas as tabelas é discreto, pode não haver uma tabela para as quais igualdade é alcançada. Se <math>\alpha_e</math>é o maior p-valor menor que 5% que pode realmente ocorrer para alguma tabela, então o teste proposto efetivamente testa no nível <math>\alpha_e</math>. Para amostras de tamanho pequenos, <math>\alpha_e</math> pode ser significantemente menor que 5%.<ref name="Liddell-1976" /><ref name="Berkson1978" /><ref name="DAgostino1988" /> Enquanto este efeito ocorre para qualquer estatística discreta (não somente em tabelas de contingência ou para testes de Fisher), tem sido argumentado que o problema é composto pelas condições de teste de Fisher nas marginais.<ref>Cyrus R. Mehta and Pralay Senchaudhuri (4 September 2003). "[http://www.cytel.com/Papers/twobinomials.pdf Conditional versus Unconditional Exact Tests for Comparing Two Binomials] {{Wayback|url=http://www.cytel.com/Papers/twobinomials.pdf |date=20110708215834 }}". Retrieved 20 November 2009.</ref> Para evitar o problema, muitos autores desencorajam o uso de níveis de significância fixos quando lidar com problemas discretos.<ref name="Yates1984" /><ref name="Little1989" />

A decisão de condicionar nas margens da tabela também é controverso.<ref name="Barnard1945a">{{citar periódico|doi=10.1038/156177a0|ultimo=Barnard|primeiro=G.A|ano=1945|titulo=A New Test for 2×2 Tables|periodico=Nature|volume=156|paginas=177|numero=3954}}</ref><ref name="NatureDiscussion">{{citar periódico|ultimo=Fisher|primeiro=Ronald|ano=1945|periodico=Nature|volume=156|paginas=388|doi=10.1038/156388a0|titulo=A New Test for 2 × 2 Tables|numero=3961}}; {{citar periódico|ultimo=Barnard|primeiro=G.A|ano=1945|periodico=Nature|volume=156|paginas=783|titulo=A New Test for 2×2 Tables|doi=10.1038/156783b0|numero=3974}}</ref> Os p-valores derivados do teste de Fisher vêm da distribuição que condicionam os totais marginais. Nesse sentido, o teste é exato somente para a distribuição condicional e não para a tabela original onde os totais marginais podem mudar de experimento para experimento. É possível obter um p-valor exato para a tabela 2x2 quando as margens não são fixadas. o [[teste de Barnard]], por exemplo, permite margens aleatórias. Entretanto, alguns autores<ref name="Yates1984" /><ref name="Little1989" /><ref name="NatureDiscussion" /> (incluindo, posteriormente, o próprio Barnard)<ref name="Yates1984" /> tem criticado o teste de Barnard baseado nesta propriedade. Eles argumentam que o sucesso marginal total é (quase<ref name="Little1989" />) uma [[estatística auxiliar]], (quase) não contendo informação sobre a propriedade testada.

O ato de condicionar no sucesso marginal de uma tabela 2x2 pode ser mostrado para ignorar alguma informação nos dados sobre chances proporcionais desconhecidas.<ref name="Choi2015">{{citar periódico|ultimo=Choi|primeiro=L|ultimo2=Blume|primeiro2=J.D.|ultimo3=Dupont|primeiro3=W.D.|ano=2015|titulo=Elucidating the Foundations of Statistical Inference with 2 x 2 Tables|periodico=PLoS ONE|volume=10|numero=4|paginas=e0121263|doi=10.1371/journal.pone.0121263|url=http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0121263}}</ref> O argumento de que os totais marginais são (quase) auxiliares implica que a [[função de verossimilhança]] apropriada por fazer inferências sobre estas chances proporcionais devem ser condicionados no sucesso marginal.<ref name="Choi2015" /> Se essa informação perdida é importante para propósitos inferenciais é a essência da controvérsia.<ref name="Choi2015" />

== Alternativas ==
Uma alternativa ao teste exato, o [[teste de Barnard]], tem sido desenvolvido e os proponentes sugerem que esse método é mais poderoso, particularmente em tabelas 2 × 2. Outra alternativa é usar estimadores de [[máxima verossimilhança]] para calcular o p-valor da [[distribuição binomial]] ou [[Distribuição multinomial|multinomial]] e rejeitar ou falhar em rejeitar baseado no p-valor.{{carece de fontes|data=junho de 2017}}Para dados categóricos estratificados o [[teste de Cochran–Mantel–Haenszel]] deve ser usado ao invés do teste de Fisher. Choi et al.<ref name="Choi2015" /> propôs um p-valor derivado da razão de verossimilhança baseada na distribuição condicional da [[razão de possibilidades]] dada a razão de sucesso marginal. Este p-valor é inferencialmente consistente com testes clássicos de dados normalmente distribuídos bem como razões de verossimilhanças e intervalos de suporte baseados na função de verosimilhança condicional. Também é facilmente computável.<ref name="Choi2011">{{citar web|ultimo=Choi|primeiro=Leena|ano=2011|titulo=ProfileLikelihood: profile likelihood for a parameter in commonly used statistical models; 2011. R package version 1.1.|url=https://cran.r-project.org/web/packages/ProfileLikelihood/index.html}} See also: [http://statcomp2.vanderbilt.edu:37212/dalep/ProfileLikelihood/ Likelihood Ratio Statistics for 2 x 2 Tables] {{Wayback|url=http://statcomp2.vanderbilt.edu:37212/dalep/ProfileLikelihood/ |date=20160604023208 }} (Online calculator).</ref>

== Veja também ==
* [[Teste Qui-Quadrado]]
* [[Teste de Barnard]]

{{Referências|col=2}}

== Ligações externas ==
* [http://www.langsrud.com/stat/fisher.htm Calculate Fisher's Exact Test Online]
* [https://play.google.com/store/apps/details?id=appinventor.ai_admin.SimpleImpactCalculator Fisher's Exact Test calculator for Android devices]
* [http://statcomp2.vanderbilt.edu:37212/ProfileLikelihood/ Likelihood Ratio Statistics for 2X2 Tables]

{{Estatística}}

{{Portal3|Probabilidade e estatística}}

[[Categoria:Testes estatísticos]]