以R的条形图绘制一年的密度平均值和标准偏差

Question

我正在尝试在R中创建一个条形图，以显示在不同收获率下不同密度降低区域（1、5、10、15、20、25平方英里）不同尺寸的20年的平均密度。 30％，35％，40％，45％，50％）。以下是我的数据示例（称为“收割”）：

    total.sink.density    size    year   harvest.rate
    26.4                  1       1      0.30
    25.7                  1       1      0.35
    22.9                  1       1      0.40
    20.1                  1       1      0.45
    18.7                  1       1      0.50
     .                    .       .        .
     .                    .       .        .
     .                    .       .        .
    6.5                   25      20     0.30
    7.2                   25      20     0.35
    6.9                   25      20     0.40
    6.5                   25      20     0.45
    7.0                   25      20     0.50

使用在线上找到的“ summarySE”代码，我可以获得一些描述性统计信息（平均值，标准偏差，se和ci）：

 newdeerharvest <- summarySE(deerharvest, 
 measurevar="total.sink.density", groupvars=c("harvest.rate","size"))

我认为ggplot将完成我想做的事情。我尝试过使用以下代码：

ggplot(newdeerharvest, aes(x=size, y=total.sink.density, fill=harvest.rate)) + 
geom_bar(position=position_dodge(), stat="identity") +
geom_errorbar(aes(ymin=total.sink.density-sd, ymax=total.sink.density+sd)

但是，这将在整个20年中提供收益。我想针对每种尺寸和收成率绘制仅20年的平均值。有什么方法可以完成而无需先对数据进行子集设置？我需要为多个数据集制作相同的图形。任何见解都将不胜感激！