仅在 R 中的 ANOVA 循环上打印显着结果答案

【问题标题】：Printing Significant Result only on an ANOVA loop in R仅在 R 中的 ANOVA 循环上打印显着结果
【发布时间】：2014-11-17 04:13:05
【问题描述】：

我的数据如下所示：

     KAFFE   SAFT    FERSKVARE   HERMETIKK ..... K
      35    23.91     17.3         70.2          1
      53    30.13     27.4         75.7          1
      35    34.54     17.3         81.1          2
      32    31.93     47.3         92.7          2
      34    54.69     17.2         88.4          2 
      35    34.69     19.1         81.5          3 
      54    25.63     18.3         92.1          3 
      34    46.65     18.6         76.6          3
      36    38.21     23.3         79.9          3

我使用循环为每个数据创建方差分析和事后测试

rcluso is my data frame
K is the categorical variable located at the end column



sink("AnovaTest.doc")
      for(i in 1:(ncol(rcluso)-1)){
        columns <- names(rcluso[i])
        anovaresult<- summary(aov(rcluso[,i]~K,data=rcluso))
        posthocresult <- TukeyHSD(aov(rcluso[,i]~K,data=rcluso))

    print(columns)
    print(anovaresult)
    print(posthocresult)
    }
    sink()

问题是，我只想打印出方差分析的显着结果，并且我还想为每个显着结果创建一个条形图。

你们能帮我解决这个问题吗？提前致谢。

【问题讨论】：

标签： r loops plot anova

【解决方案1】：

如果你得到 anova() 输出而不是 summary() 你可以访问 p 值。

更换你的

    anovaresult<- summary(aov(rcluso[,i]~K,data=rcluso))

与

    anovaresult<- anova(aov(rcluso[,i]~K,data=rcluso))

然后您可以使用 if() 语句有条件地打印结果，例如

  if(anovaresult$Pr[1] < 0.05){ print(anovaresult) }

【讨论】：