Bonferroni在R

Question

我正在使用图像Image中的相同方法，我需要知道如何将脚本放在R中，以便在他们提到的Duncan测试中进行Bonferroni校正，给出“（α '=α/ 3，α= 0.05）“。 我使用了以下脚本：

  library(agricolae)
  model <- aov(value~ind,data=data.frame)
  comparison <- duncan.test(model, "Ind", main= "comparison")
  duncan.test(model,"ind", console=TRUE) 

我尝试纠正p值，在你的一个答案Answer from StackOverflow之后在duncan.test中添加“alpha”，如下所示：alpha by 3 = 0.05 / 3~0.0167，但是它产生了一个错误在该计划中。我将衷心感谢您的帮助。