{"id":29212,"date":"2024-02-21T17:43:57","date_gmt":"2024-02-21T09:43:57","guid":{"rendered":"http:\/\/www.biocloudservice.com\/wordpress\/?p=29212"},"modified":"2024-02-21T17:43:59","modified_gmt":"2024-02-21T09:43:59","slug":"%e4%b8%80%e6%96%87%e5%ad%a6%e4%bc%9a%ef%bc%8canova%e5%8f%8c%e5%9b%a0%e7%b4%a0%e6%96%b9%e5%b7%ae%e5%88%86%e6%9e%90","status":"publish","type":"post","link":"http:\/\/www.biocloudservice.com\/wordpress\/?p=29212","title":{"rendered":"\u4e00\u6587\u5b66\u4f1a\uff0cANOVA\u53cc\u56e0\u7d20\u65b9\u5dee\u5206\u6790"},"content":{"rendered":"

\u4eca\u5929\u5c0f\u679c\u901a\u8fc7\u8fd9\u4e00\u7bc7\u6587\u7ae0\u5e26\u5927\u5bb6\u5b66\u4f1a\u548c\u4e86\u89e3\u53cc\u56e0\u7d20\u65b9\u5dee\u5206\u6790(ANOVA),\u8bdd\u4e0d\u591a\u8bf4\uff0c\u8bf7\u770b\u4e0b\u9762\u4ee3\u7801\uff1a<\/p>\n

    \n
  1. \u5b89\u88c5\u9700\u8981\u7684R\u5305<\/li>\n<\/ol>\n

    install.packages(\u201cggplot2\u201d)<\/p>\n

    install.packages(\u201ctidyverse\u201d)<\/p>\n

    install.packages(“multcompView”)<\/p>\n

    install.packages(\u201cggsci\u201d)<\/p>\n

      \n
    1. \u8f7d\u5165\u9700\u8981\u7684R\u5305<\/li>\n<\/ol>\n

      library(ggplot2)<\/p>\n

      library(tidyverse)<\/p>\n

      library(multcompView)<\/p>\n

      library(ggsci)<\/p>\n

        \n
      1. \u4ee3\u7801\u5c55\u793a<\/li>\n<\/ol>\n

        soil <- read.table(file=”soil.txt”, sep = ‘\\t’, header = TRUE)<\/p>\n

        \"Dingtalk_20230302154627\"<\/p>\n

        group <- read.table(file=”group.txt”, sep = ‘\\t’, header = TRUE)<\/p>\n

        \"Dingtalk_20230302154741\"<\/p>\n

        #\u5de6\u8fde\u63a5\u5408\u5e76\u4e24\u4e2a\u6587\u4ef6<\/p>\n

        soil <- left_join(soil, group, by = ‘sample’)<\/p>\n

        #\u7b5b\u9009\u9700\u8981\u7684\u5217<\/p>\n

        chao1 <- select(soil ,c(‘sample’, ‘treat’, ‘times’, ‘chao1’))<\/p>\n

        #\u5c06\u8fd9\u4e24\u4e2a\u56e0\u7d20\u7c7b\u578b(trate\u548ctimes)\u6539\u4e3a\u56e0\u5b50<\/p>\n

        chao1$treat <- factor(chao1$treat)<\/p>\n

        chao1$times <- factor(chao1$times)<\/p>\n

        #\u53cc\u56e0\u7d20\u65b9\u5dee\u5206\u6790<\/p>\n

        anova<- aov(chao1~treat*times, data = chao1)<\/p>\n

        #\u8ba1\u7b97mean\u548csd(\u5206\u7ec4\u7edf\u8ba1mean\u548ctimes\uff09<\/p>\n

        dt <- chao1 %>% group_by(treat,times) %>%<\/p>\n

        summarise(w=mean(chao1), sd = sd(chao1)) %>%<\/p>\n

        arrange(desc(w)) %>%<\/p>\n

        ungroup()<\/p>\n

        #\u591a\u91cd\u914d\u5bf9\u6bd4\u8f83<\/p>\n

        cld <- TukeyHSD(anova) %>% multcompLetters4(anova,.)<\/p>\n

        le <- as.data.frame.list(cld$`treat:times`)<\/p>\n

        #\u83b7\u53d6\u5b57\u6bcd\u6807\u8bb0<\/p>\n

        dt$Tukey <- le$Letters<\/p>\n

        #\u5f00\u59cb\u7ed8\u5236\u5e26\u6709\u663e\u8457\u6027\u5b57\u6bcd\u6807\u8bb0\u7684\u67f1\u72b6\u56fe<\/p>\n

        ggplot(dt, aes(times, w)) +<\/p>\n

        geom_bar(stat = “identity”,aes(fill=times),show.legend = FALSE,width=0.5) +<\/p>\n

        geom_errorbar(aes(ymin = w-sd, ymax=w+sd), width = 0.1) +<\/p>\n

        geom_text(aes(label = Tukey, y = w + sd), vjust = -0.5)+<\/p>\n

        facet_wrap(~treat)+<\/p>\n

        scale_fill_npg()+<\/p>\n

        labs(x=””,y=””)+<\/p>\n

        theme_test()+<\/p>\n

        theme(panel.background = element_blank(),<\/p>\n

        axis.line = element_line(color = “black”),<\/p>\n

        axis.text = element_text(size = 10,color = “black”),<\/p>\n

        axis.text.y = element_text(size = 10),<\/p>\n

        axis.ticks.x = element_blank())<\/p>\n

        \"Dingtalk_20230302155410\"<\/p>\n

        \u6700\u7ec8\u901a\u8fc7\u4e00\u4e2a\u793a\u4f8b\u6570\u636e\u8fdb\u884c\u4e86\u53cc\u56e0\u7d20\u65b9\u5dee\u5206\u6790\uff0c\u4eca\u5929\u5c0f\u679c\u7684\u5206\u4eab\u5c31\u5230\u8fd9\u91cc\u5965\uff0c\u4e0b\u671f\u518d\u89c1\u3002<\/p>\n","protected":false},"excerpt":{"rendered":"

        \u4eca\u5929\u5c0f\u679c\u901a\u8fc7\u8fd9\u4e00\u7bc7\u6587\u7ae0\u5e26\u5927\u5bb6\u5b66\u4f1a\u548c\u4e86\u89e3\u53cc\u56e0\u7d20\u65b9\u5dee\u5206\u6790(ANOVA),\u8bdd\u4e0d\u591a\u8bf4\uff0c\u8bf7\u770b\u4e0b\u9762\u4ee3\u7801\uff1a \u5b89\u88c5\u9700\u8981\u7684R\u5305 […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":""},"categories":[1],"tags":[],"jetpack_featured_media_url":"","_links":{"self":[{"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts\/29212"}],"collection":[{"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=29212"}],"version-history":[{"count":1,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts\/29212\/revisions"}],"predecessor-version":[{"id":29216,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts\/29212\/revisions\/29216"}],"wp:attachment":[{"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=29212"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=29212"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=29212"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}