Question

我正在寻找以小图方式绘制两组（升高的温度和控制温度）之间值的差异的方法。我希望能够以与原始图形相同的方式在x轴上绘制月份。

以下是我用来绘制当前绘制的1的以下脚本。这显示y轴=植物生长，x轴=月。

脚本：当前条形图

Tempmean<- data %>% group_by (Treatment, Month) %>% summarize (TTmean = mean(Amean, na.rm=TRUE), TTsd=sd(Amean,na.rm=TRUE))


p<-ggplot(data=Tempmean, aes(x=factor(Month), y=TTmean, fill=Treatment)) +
  geom_bar(stat="identity", position="dodge", colour="black" , size = 0.25, width=0.5) + geom_errorbar(aes(ymin=TTmean-TTsd, ymax=TTmean+TTsd), width=.1,
                position=position_dodge(.5)) + scale_fill_manual(values=c("darkgray","darkolivegreen")) + scale_x_discrete(breaks=6:8,labels=c("June","July","August")) + scale_y_continuous(limits=c(0,20), breaks=seq(0,20,2))
p

这是我正在使用的数据2。我希望从aCO2的TTmean中获取eCO2的TTmean。

数据：

    structure(list(Treatment = c("aCO2", "aCO2", "aCO2", "eCO2","eCO2", "eCO2"), Month = c(6L, 7L, 8L, 6L, 7L, 8L), TTmean = c(10.1922587348143,10.1061784054575, 8.27148533916994, 12.0261355594138,10.8954781586458, 10.9468200269188), TTsd =c(7.04936647397141,4.18653008350561, 1.50026716071241, 3.25471492346035, 0.742036555955107, 2.00464198948226)), row.names = c(NA, -6L), class = c("grouped_df", "tbl_df", "tbl", "data.frame"), vars = "Treatment", drop = TRUE, indices = list(0:2, 3:5), group_sizes = c(3L, 3L), biggest_group_size = 3L, labels = structure(list(Treatment = c("aCO2", "eCO2")), row.names = c(NA, -2L),class = "data.frame", vars = "Treatment", drop = TRUE))

Answer 1

这应该可以解决问题。

我放弃了TTsd，因为您似乎不需要它。诀窍是spread()数据，以便您可以轻松计算值的差异。我计算为aCO2减去eCO2；但您可以在mutate()

中进行更改

library(tidyverse)
Tempmean %>%
    select(-TTsd) %>%
    # either group_by Month, or just ungroup entirely
    group_by(Month) %>%
    spread(Treatment, TTmean) %>%
    mutate(T_diff = aCO2 - eCO2) %>%
    ggplot(aes(factor(Month), T_diff)) %+%
    geom_bar(position = "dodge", stat = "identity", size = 0.25, width=0.5) %+%
    scale_x_discrete(breaks=6:8,labels=c("June","July","August"))

如何在多个采样日内绘制两组之间的差异？

1 个答案: