1

I am very new to R, but am trying to create vertical mean lines in a barplot from the ggplot2 library.

My current graph is grouped by year, with each colored bar representing a different age. The y axis is count: mygraph

I would like to have a vertical mean line for each group by year (i.e., a mean line of age for year 1, a mean line of age for year 2, etc.)

This is what I have written thus far:

library(ggplot2)
AgeDat<-read.csv(File Name)
columnmeans<-c("MeanValue", "MeanValue", "MeanValue"..."MeanValue")
SepAgePlot<-ggplot(AgeDat, aes(factor(Year), Count, fill=(factor(Age))))+
  geom_bar(stat="identity", position="dodge")+ 
  SepAgePlot+(geom_vline(data=AgeDat, xintercept = columnmeans)+
  theme_classic()
SepAgePlot+labs (x="Year", y="count", title="my graph")

Thank you in advance for all your help!

UDPATE: This would be a sample data set: 

Year    Age Count
1964    31  15
1964    33  23
1964    34  54
1964    35  8
1964    36  44
1964    37  21
1964    38  23
1964    39  26
1964    40  23
1965    30  22
1965    31  23
1965    32  45
1965    33  55
1965    34  23
1965    35  10
1965    36  12
1965    37  16
1965    38  32
1965    39  36
1965    40  13
1966    30  27
1966    31  32
1966    32  19
1966    33  45
1966    34  35
1966    35  60
1966    36  15
1966    37  28
1966    38  56
1966    39  18
1966    40  25
1967    30  36
1967    31  32
1967    32  23
1967    33  9
1967    34  15
1967    35  0
1967    36  5
1967    37  7
1967    38  24
1967    39  31
1967    40  24

The mean age of each year is:

1964 35.4

1965 35.6

1966 35.0

1967 34.6

My goal is to create a barplot of the data grouped by year, with a line representing the mean age for each year.

user9165024
  • 69
  • 1
  • 8
  • 1
    What does your data look like? It's easier to help you if you provide a [reproducible example](https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example) that can be used to test and verify possible solutions. Is there an example somewhere of what exactly you mean by "vertical mean line"? Not sure what the desired output is. – MrFlick Jan 19 '18 at 15:07
  • Thank you for your help! To provide some clarification: the first graph in the answers section of the following link is what I am hoping to accomplish, however I would like to do this for each year of my graph and for the mean:https://stackoverflow.com/questions/37280407/ggplot2-how-to-create-vertical-line-corresponding-to-quantile-in-geom-bar-plot – user9165024 Jan 19 '18 at 15:32
  • Please either use a built-in data set or share some data reproducibly with `dput()`. See the link in MrFlick's comment above for a nice tutorial on making reproducible examples in R. – Gregor Thomas Jan 19 '18 at 15:40
  • Can you please post your sample data in an easy-to-copy-paste format? `dput()` produces this automatically. Otherwise at least remove the `>` and `.` you have at the start of each line, and the `:` in the mean data. See MrFlick's reproducible example link for more suggestions. – Gregor Thomas Jan 19 '18 at 16:26
  • I have edited the sample data per your recommendations. Thank you for your help! – user9165024 Jan 19 '18 at 17:02

1 Answers1

2

You can use the group_by and summarise functions from the dplyr package to add a column with the mean age for each year:

AgeDat <- AgeDat %>% group_by(Year) %>% summarise(MeanValue=sum(Age*Count)/sum(Count)) %>% inner_join(AgeDat, .)
AgeDat
    Year Age Count MeanValue
1   1964  30    44  34.64545
2   1964  31    44  34.64545
3   1964  32    46  34.64545
4   1964  33    35  34.64545
5   1964  34    83  34.64545
6   1964  35    70  34.64545
7   1964  36    73  34.64545
8   1964  37    85  34.64545
9   1964  38    31  34.64545
10  1964  39    39  34.64545
11  1965  30   100  34.65485
12  1965  31    19  34.65485

You can then use geom_vline to add a vertical line to each year, and facet_wrap to get the age distribution by year:

ggplot(AgeDat) + 
    geom_bar(aes(x=Age, y=Count), stat="identity", position="dodge") + 
    geom_vline(aes(xintercept=MeanValue)) + 
    facet_wrap(~Year, nrow=1) + 
    theme_classic()

enter image description here

LucyMLi
  • 657
  • 4
  • 14
  • Excellent solution! I altered some of the code to accommodate work I had already completed, however your response was integral! Thank you for your response, LucyMLi! – user9165024 Jan 22 '18 at 21:18