names(CO2)
## [1] "Plant" "Type" "Treatment" "conc" "uptake"
head(CO2)
## Plant Type Treatment conc uptake
## 1 Qn1 Quebec nonchilled 95 16.0
## 2 Qn1 Quebec nonchilled 175 30.4
## 3 Qn1 Quebec nonchilled 250 34.8
## 4 Qn1 Quebec nonchilled 350 37.2
## 5 Qn1 Quebec nonchilled 500 35.3
## 6 Qn1 Quebec nonchilled 675 39.2
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.1.3 ✔ readr 2.1.4
## ✔ forcats 1.0.0 ✔ stringr 1.5.0
## ✔ ggplot2 3.4.4 ✔ tibble 3.2.1
## ✔ lubridate 1.9.3 ✔ tidyr 1.3.0
## ✔ purrr 1.0.2
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
CO2%>%ggplot(aes(Treatment,uptake))+geom_boxplot(fill="blue",alpha=0.5)+stat_summary(fun=mean)
## Warning: Removed 2 rows containing missing values (`geom_segment()`).

#CO2%>%count(Plant,uptake)%>%group_by(Plant)%>%mutate(mean=mean(uptake))%>%mutate(ftc_reorder(Plant,mean))%>%ggplot(aes(Plant,mean))+geom_col()