График среднего и SD плотности за один год в виде гистограммы в R

Question

Я пытаюсь создать гистограмму в R, которая будет показывать среднюю плотность в 20 году для разных размеров области уменьшения плотности (1, 5, 10, 15, 20, 25 квадратных миль) при разных коэффициентах сбора урожая (30%, 35%, 40%, 45%, 50%).Ниже приведен пример моих данных (называемых "deerharvest"):

    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", я могу получить некоторую описательную статистику (среднее, sd, 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 года, для каждого размера и нормы урожая.Есть ли способ сделать это без предварительного подбора данных?Мне нужно было бы сделать этот же график для нескольких наборов данных.Любое понимание очень ценится!