Estou trabalhando com algumas análises espaciais exploratórias em R usando o pacote spdep.
Encontrei uma opção para ajustar os valores- p dos indicadores locais de associação espacial (LISA) calculados usando a localmoran
função. De acordo com os documentos, ele visa:
... ajuste do valor de probabilidade para vários testes.
Além disso, nos documentos de p.adjustSP
eu li que as opções disponíveis são:
Os métodos de ajuste incluem a correção de Bonferroni ('"bonferroni"'), na qual os valores de p são multiplicados pelo número de comparações. Quatro correções menos conservadoras também estão incluídas por Holm (1979) ('holm'), Hochberg (1988) ('' hochberg ''), Hommel (1988) ('' hommel '') e Benjamini e Hochberg (1995) ('"fdr"'), respectivamente. Uma opção de passagem ('"none"') também está incluída.
Os quatro primeiros métodos são projetados para fornecer um forte controle da taxa de erro familiar. Parece não haver razão para usar a correção de Bonferroni não modificada, porque é dominada pelo método de Holm, que também é válido sob premissas arbitrárias.
Os métodos de Hochberg e Hommel são válidos quando os testes de hipóteses são independentes ou quando não estão associados negativamente (Sarkar, 1998; Sarkar e Chang, 1997). O método de Hommel é mais poderoso que o de Hochberg, mas a diferença é geralmente pequena e os valores de p de Hochberg são mais rápidos de calcular.
O método "BH" (aka "fdr") e "BY" de Benjamini, Hochberg e Yekutieli controlam a taxa de falsas descobertas, a proporção esperada de descobertas falsas entre as hipóteses rejeitadas. A taxa de descoberta falsa é uma condição menos rigorosa que a taxa de erro familiar, portanto, esses métodos são mais poderosos que os outros.
Algumas perguntas que apareceram:
- Em palavras simples - qual é o objetivo desse ajuste?
- É necessário usar essas correções?
- Se sim - como escolher entre as opções disponíveis?