{"id":375,"date":"2022-04-09T21:07:32","date_gmt":"2022-04-09T13:07:32","guid":{"rendered":"https:\/\/linguopeng.top\/?p=375"},"modified":"2022-04-09T21:07:33","modified_gmt":"2022-04-09T13:07:33","slug":"%e5%a4%9a%e5%9b%a0%e7%b4%a0%e6%96%b9%e5%b7%ae%e5%88%86%e6%9e%90","status":"publish","type":"post","link":"https:\/\/linguopeng.top\/?p=375","title":{"rendered":"\u591a\u56e0\u7d20\u65b9\u5dee\u5206\u6790"},"content":{"rendered":"\n<pre class=\"wp-block-code\"><code>library(multcompView)\nlibrary(stats)\nlibrary(tidyverse)\nlibrary(ggsci)\n\ndata2 &lt;- ToothGrowth\ndata2$dose = as.factor(data2$dose)\n\nanova &lt;- aov(len ~ supp*dose, data = data2) #two way anova\nsummary(anova) #represents ANOVA\n\ntukey &lt;- TukeyHSD(anova)\n\ngroup_lettering &lt;- multcompLetters4(anova, tukey) #\ngroup_lettering\ngroup_lettering2 &lt;- data.frame(group_lettering$`supp:dose`$Letters)\ngroup_lettering2\n\nmean_data2 &lt;- data2 %>% \n  group_by(supp, dose) %>% \n  summarise(len_mean=mean(len), sd = sd(len)) %>%\n  arrange(desc(len_mean))\n\nmean_data2$group_lettering &lt;- group_lettering2$group_lettering..supp.dose..Letters\n\n\nggplot(mean_data2,aes(x = dose, y = len_mean,group=supp))+\n  geom_bar(position=position_dodge(0.9),stat = \"identity\", aes(fill = supp), \n           show.legend = TRUE)\n\nggplot(mean_data2, aes(x = dose, y = len_mean,group=supp))+\n  geom_bar(position=position_dodge(0.9),stat = \"identity\",\n           aes(fill = supp),show.legend = TRUE) +\n  geom_errorbar(aes(ymin = len_mean-sd, ymax=len_mean+sd),width = 0.1,\n                position=position_dodge(0.9))\n\nggplot(mean_data2, aes(x = dose, y = len_mean,group=supp))+\n  geom_bar(position=position_dodge(0.9),stat = \"identity\",\n           aes(fill = supp), show.legend = TRUE) +\n  geom_errorbar(aes(ymin = len_mean-sd, ymax=len_mean+sd),\n                width = 0.1, position=position_dodge(0.9))+\n  geom_text(aes(label = group_lettering, y = len_mean + sd),\n            vjust=-0.4,position=position_dodge(0.9))\n\n\nggplot(mean_data2, aes(x = dose, y = len_mean,group=supp))+\n  geom_bar(position=position_dodge(0.9),stat = \"identity\",\n           aes(fill = supp), show.legend = TRUE) +\n  geom_errorbar(aes(ymin = len_mean-sd, ymax=len_mean+sd),\n                width = 0.1, position=position_dodge(0.9)) + \n  geom_text(aes(label = group_lettering, y = len_mean + sd),\n            vjust=-0.4, position=position_dodge(0.9)) +\n  scale_y_continuous(expand = expansion(0),limits = c(0,35),\n                     breaks = seq(0,35,5))+\n  labs(x=NULL,y=NULL)+\n  theme_minimal()+\n  theme(\n    plot.margin = unit(c(0.2,0.2,0.2,0.2), \"cm\"),\n    panel.background = element_blank(),\n    axis.line = element_line(color = \"black\"),\n    axis.title = element_text(size = 10, color = \"black\",face = \"bold\"),\n    axis.text = element_text(size = 10,color = \"black\"),\n    axis.text.x = element_text(margin=margin(t =3)),\n    axis.text.y = element_text(size = 10),\n    axis.title.y = element_text(margin = margin(r = 10)),\n    axis.ticks.x = element_blank())+\n  # scale_fill_jco()\nscale_fill_jama()<\/code><\/pre>\n\n\n\n<figure class=\"wp-block-image size-full\"><div class='fancybox-wrapper lazyload-container-unload' data-fancybox='post-images' href='https:\/\/linguopeng.top\/wp-content\/uploads\/2022\/04\/QQ\u622a\u56fe20220409210532.png'><img class=\"lazyload lazyload-style-1\" src=\"data:image\/svg+xml;base64,PCEtLUFyZ29uTG9hZGluZy0tPgo8c3ZnIHdpZHRoPSIxIiBoZWlnaHQ9IjEiIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIgc3Ryb2tlPSIjZmZmZmZmMDAiPjxnPjwvZz4KPC9zdmc+\"  loading=\"lazy\" decoding=\"async\" width=\"577\" height=\"394\" data-original=\"https:\/\/linguopeng.top\/wp-content\/uploads\/2022\/04\/QQ\u622a\u56fe20220409210532.png\" src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAAEAAAABCAYAAAAfFcSJAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsQAAA7EAZUrDhsAAAANSURBVBhXYzh8+PB\/AAffA0nNPuCLAAAAAElFTkSuQmCC\" alt=\"\" class=\"wp-image-376\"\/><\/div><\/figure>\n\n\n\n<figure class=\"wp-block-image size-full\"><div class='fancybox-wrapper lazyload-container-unload' data-fancybox='post-images' href='https:\/\/linguopeng.top\/wp-content\/uploads\/2022\/04\/QQ\u622a\u56fe20220409210645.png'><img class=\"lazyload lazyload-style-1\" src=\"data:image\/svg+xml;base64,PCEtLUFyZ29uTG9hZGluZy0tPgo8c3ZnIHdpZHRoPSIxIiBoZWlnaHQ9IjEiIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIgc3Ryb2tlPSIjZmZmZmZmMDAiPjxnPjwvZz4KPC9zdmc+\"  loading=\"lazy\" decoding=\"async\" width=\"680\" height=\"199\" data-original=\"https:\/\/linguopeng.top\/wp-content\/uploads\/2022\/04\/QQ\u622a\u56fe20220409210645.png\" src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAAEAAAABCAYAAAAfFcSJAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsQAAA7EAZUrDhsAAAANSURBVBhXYzh8+PB\/AAffA0nNPuCLAAAAAElFTkSuQmCC\" alt=\"\" class=\"wp-image-377\"\/><\/div><\/figure>\n","protected":false},"excerpt":{"rendered":"","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[],"class_list":["post-375","post","type-post","status-publish","format-standard","hentry","category-uncategorized"],"_links":{"self":[{"href":"https:\/\/linguopeng.top\/index.php?rest_route=\/wp\/v2\/posts\/375","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/linguopeng.top\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/linguopeng.top\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/linguopeng.top\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/linguopeng.top\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=375"}],"version-history":[{"count":1,"href":"https:\/\/linguopeng.top\/index.php?rest_route=\/wp\/v2\/posts\/375\/revisions"}],"predecessor-version":[{"id":378,"href":"https:\/\/linguopeng.top\/index.php?rest_route=\/wp\/v2\/posts\/375\/revisions\/378"}],"wp:attachment":[{"href":"https:\/\/linguopeng.top\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=375"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/linguopeng.top\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=375"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/linguopeng.top\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=375"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}