使用Pearson相关和线性回归的Bonferroni校正


9

我正在针对3个DV的5个IV(5个人格特质,性格外向,和agree,尽责,神经质,开放)运行统计数据,包括PCT态度,CBT态度,PCT与CBT态度。我还添加了年龄和性别,以查看还有其他影响。

我正在测试以查看人格特征是否可以预测DV的态度。

最初,我对所有变量都使用了Pearson相关性(45个测试)。

主要发现是外向性与PCT态度在p = 0.05相关。但是,当我进行45次测试时,我对Bonferroni进行了alpha = 0.05 / 45 = 0.001的校正,因此这一发现无关紧要。

然后,我对所有变量进行了简单的线性回归,对于PCT态度,外向性再次很重要。如果我进行Bonferroni校正,那么它再次显得微不足道。

问题:

  1. 我需要Bonferroni纠正Pearson的相关性吗?
  2. 如果我这样做了,因此对PCT的态度无忧无虑,那么进行线性回归还有意义吗?
  3. 如果我进行线性回归,是否还需要对此进行Bonferroni校正?
  4. 我只报告校正后的值还是未校正和校正后的值?

4
问题1与外观相关,您将发现(相关性),问题3与(问题3 )相关。为多次比较在多元回归中调整p值是一个好主意吗?。更一般而言,此查询的结果可能令人感兴趣。
chl 2012年

Answers:


4

我认为Chl在没有直接回答问题的情况下为您指出了很多好的材料和参考资料。我给出的答案可能有点争议,因为我知道一些统计学家不相信多重性调整,而许多贝叶斯主义者也不相信p值。实际上,我曾经听过Don Berry说过,尤其是在控制I型错误的自适应设计中使用贝叶斯方法不是问题。在看到确保有害药物不会进入市场对FDA实际上有多么重要之后,他后来又收回了这一点。

我的回答是是,不是。如果您进行45次测试,您肯定需要针对多重性进行调整,但不必针对Bonferroni进行调整,因为它可能过于保守。当您为关联进行数据挖掘时,I型错误的膨胀显然是一个引起关注的问题,该主题已被引用的文章“外观,您将找到关联”。这三个链接都提供了很好的信息。我认为缺少的是Westfall和Young很好地开发的p值调整的重采样方法。您可以在我的引导书中找到示例,也可以在其重采样书中找到完整的详细信息。我的建议是考虑使用引导程序或置换方法进行p值调整,并可能考虑在严格的家庭错误率之上的错误发现率。

链接到Westfall和Young:http://www.amazon.com/Resampling-Based-Multiple-Testing-Adjustment-Probability/dp/0471557617/ref=sr_1_1?s = books&ie = UTF8&qid = 1343398751&sr = 1-1&keywords = peter + 西部荒野

Bretz等人最近在多重比较上的新书:http ://www.amazon.com/Multiple-Comparisons-Using-Frank-Bretz/dp/1584885742/ref=sr_1_2?s=books&ie=UTF8&qid=1343398796&sr=1-2&keywords =彼得+韦斯特洛斯

我的书包含第8.5节的内容和大量的引导程序参考:http ://www.amazon.com/Bootstrap-Methods-Practitioners-Researchers-Probability/dp/0471756210/ref=sr_1_2?s=books&ie=UTF8&qid=1343398953&sr =1 -2&keywords = Michael + chernick


+1 Westfall&Young末尾格雷厄姆·马丁(Graham Martin)的蒙克豪森(Munchausen)统计表格的复制品以一种非常引人入胜的方式说道。您可以在Amazon“向内看”功能中阅读此内容。(看到亚马逊为这本150 美元的书开出7 美元的折价价,真是很有趣。)
胡扯

@whuber我想我曾经看过一幅动画片,展示了男爵被他的靴子从湖中拉出自己的样子。埃夫隆(Efron)称其为引导程序是明智的,因为许多人怀疑它是否可以在统计数据中完成,就像许多人怀疑男爵的传说一样!
Michael R. Chernick

4

在我看来,这是探索性研究/数据分析,而不是确定性的。就是说,这听起来并不像您从某种理论开始,即说出于某种原因外向性应该与PCT有关。因此,我不必担心alpha调整,因为我认为它与CDA更为相关,也不会认为您的发现一定是正确的。相反,我认为它可能会是真实的,并根据我对当前主题的了解来尝试这些想法/可能性。看到此发现后,它是正确的还是您对此表示怀疑?如果这是真的,对当前的理论意味着什么?会很有趣吗?会重要吗?考虑到可能需要的时间,精力和费用,是否值得进行一项新的(确定性)研究来确定其是否正确?请记住,进行Bonferroni校正的原因是,我们希望当变量太多时会出现一些东西。因此,我认为可以进行启发式研究,即使事实证明事实并非如此,这项研究是否也能提供足够的信息'?如果您认为这不值得,则此关系将保持在“可能”类别中,然后继续进行,但如果值得这样做,请对其进行测试。


如果他真的了解探索性数据分析是什么,并且他不太重视大的相关性,那么我会同意你的观点。但是人们会承认,他们只是在进行探索性分析以过滤出弱关联,但是当看到有希望的事物时却会感到非常兴奋。那是人性的一部分。我认为以FDR为标准进行调整是控制兴奋的明智方法。
Michael R. Chernick

1
@MichaelChernick,我不一定与您不同意。我只是想提出另一种观点,而且我经常喜欢给出一个全景图,半哲学性的“这是什么”的观点。许多从业者可能陷入对他们而言似乎不可思议的细节上的困惑,而对它们的理解却根深蒂固。
gung-恢复莫妮卡

1
这里没有分歧,我理解你的意思。我只是想补充一点,如果我们能够适应并接受统计原理,并且不对结果产生既定的兴趣而亲自投入研究,那么我们可以按照您所说的去做。但这很难做到。想象一下,在一家制药公司工作过,该公司已经花费了数百万美元用于某种特定药物的临床研究,但失败了。医务主任将要您找出20个不同的亚组,然后找到一个能起作用的亚组。
Michael R. Chernick

1
亚组分析是临床研究中最具争议的方面之一。如果没有多重性调整,就无法将其合法化,事后进行就很难将其出售给FDA。这只是我最近几年的经历中的一个例子,这使我对忽略多重性的建议很敏感。
Michael R. Chernick

-1

根据以下内容:http : //birnlab.psychiatry.wisc.edu/resources/fMRI_TestRetest_Documentation.pdf

确定显着性时,请为多个比较校正p值。例如,经Bonferroni校正的p值是p值除以比较总数,在这种情况下,比较数是m(m-1)/ 2个唯一连接。

例如,您的相关性的截止p值为0.05,并假定您的相关性表为100 * 100。然后,您的p值应调整为0.05 /(100 * 99/2)。

线性回归与上面类似,应用Bonferroni校正。

我知道答案似乎与您的要求无关。在这种情况下,请让我知道,我将尽力澄清。希望对您有所帮助。

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.