我正在使用spdep软件包在R中进行一些探索性空间分析。
我遇到了一个选项,用于调整使用函数计算的空间关联性本地指标(LISA)的plocalmoran
值。根据文档,它的目标是:
...用于多个测试的概率值调整。
在p.adjustSP
我进一步阅读的文档中,可用选项包括:
调整方法包括Bonferroni校正(“ bonferroni””),其中p值乘以比较次数。Holm(1979)(“” holm“”),Hochberg(1988)('“” hochberg“'),Hommel(1988)('” hommel“')和Benjamini&Hochberg(1995)也进行了四次较不保守的校正。 ('“ fdr”')。还包括传递选项(“ none”)。
前四种方法旨在严格控制家庭错误率。似乎没有理由使用未经修改的Bonferroni校正,因为它受到Holm方法的支配,该方法在任意假设下也有效。
当假设检验是独立的或非负相关时,Hochberg和Hommel的方法是有效的(Sarkar,1998; Sarkar和Chang,1997)。Hommel的方法比Hochberg的方法更强大,但是差异通常很小,并且Hochberg p值的计算速度更快。
Benjamini,Hochberg和Yekutieli的“ BH”(又名“ fdr”)和“ BY”方法控制错误发现率,即错误发现的预期比例在被拒绝的假设中。误发现率的条件不如家庭式错误率严格,因此这些方法比其他方法更有效。
出现的几个问题:
- 简而言之-这项调整的目的是什么?
- 是否有必要使用此类更正?
- 如果是,如何从可用选项中进行选择?