如何遍历第三个变量？

debugcn 发表于 Dev

起诉

我已经修改了以下代码，从https://statsandr.com/blog/how-to-do-at-test-or-anova-for-many-variables-at-once-in-r-and-以更好的方式/＃将结果进一步传播到我的数据集中：

Day<-c(1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2)
Group<-c("A","A","A","B","B","B","C","C","C","A","A","A","A","B","B","B","C","C","C")
Rain<-c(4,4,6,5,3,4,5,5,3,6,6,6,5,3,3,3,2,5,2)
UV<-c(6,6,7,8,5,6,5,6,6,6,7,7,8,8,5,6,8,5,7)

dat<-data.frame(Day,Group,Rain,UV)
x <- which(names(dat) == "Group") 
y <- which(names(dat) == "Rain"
           | names(dat) == "UV") 
method1 <- "kruskal.test" 
method2 <- "wilcox.test" 
my_comparisons <- list(c("A", "B"), c("A", "C"), c("B", "C")) # comparisons for post-hoc test

library(ggpubr)
for (i in y) {
  for (j in x) {
      p <- ggboxplot(dat,
                   x = colnames(dat[j]), y = colnames(dat[i]),
                   color = colnames(dat[j]),
                   legend = "none",
                   palette = "npg",
                   add = "jitter"
    )
    print(
      p + stat_compare_means(aes(label = paste0(..method.., ", p-value = ", ..p.format..)),
                             method = method1, label.y = max(dat[, i], na.rm = TRUE)
      )
      + stat_compare_means(comparisons = my_comparisons, method = method2, label = "p.format") 
    )
  }
}

如何通过不同的“日”进一步重复此功能？谢谢。

LC数据科学家

我想您想每天查看结果，对吗？您可以添加第三个循环，如下所示：

for (h in unique(dat$Day)) {
  for (i in y) {
    for (j in x) {
        dat_tmp <- dat[dat$Day == h,] # create a subset of the data for each day
        p <- ggboxplot(dat_tmp,
                     x = colnames(dat_tmp[j]), y = colnames(dat_tmp[i]),
                     color = colnames(dat_tmp[j]),
                     legend = "none",
                     palette = "npg",
                     add = "jitter"
      )
      print(
        p + stat_compare_means(aes(label = paste0(..method.., ", p-value = ", ..p.format..)),
                               method = method1, label.y = max(dat_tmp[, i], na.rm = TRUE)
        )
        + stat_compare_means(comparisons = my_comparisons, method = method2, label = "p.format") 
      )
    }
  }
}

我在代码中添加了第三个循环，并在dat_tmp循环内部创建了循环，该循环成为您每天进行分析所使用的数据集。

本文收集自互联网，转载请注明来源。

如有侵权，请联系[email protected] 删除。