Estou testando a independência de duas variáveis, A e B, estratificadas por C. A e B são variáveis binárias e C é categórica (5 valores). Executando o teste exato de Fisher para A e B (todos os estratos combinados), recebo:
## (B)
## (A) FALSE TRUE
## FALSE 1841 85
## TRUE 915 74
OR: 1.75 (1.25 -- 2.44), p = 0.0007 *
onde OR é a razão de chances (estimativa e intervalo de confiança de 95%) e *
significa que p <0,05.
Executando o mesmo teste para cada estrato (C), recebo:
C=1, OR: 2.31 (0.78 -- 6.13), p = 0.0815
C=2, OR: 2.75 (1.21 -- 6.15), p = 0.0088 *
C=3, OR: 0.94 (0.50 -- 1.74), p = 0.8839
C=4, OR: 1.48 (0.77 -- 2.89), p = 0.2196
C=5, OR: 3.38 (0.62 -- 34.11), p = 0.1731
Por fim, executando o teste Cochran-Mantel-Haenszel (CMH), usando A, B e C, recebo:
OR: 1.56 (1.12 -- 2.18), p = 0.0089 *
O resultado do teste CMH sugere que A e B não são independentes em cada estrato (p <0,05); no entanto, a maioria dos testes dentro do estrato não é significativa, o que sugere que não temos evidências suficientes para descartar que A e B são independentes em cada estrato.
Então, que conclusão é certa? Como relatar a conclusão dados esses resultados? C pode ser considerado uma variável confusa ou não?
EDIT: Realizei o teste de Breslow-Day para a hipótese nula de que a razão de chances é a mesma entre os estratos e o valor de p foi de 0,1424.