R/cp_GQage.R

Defines functions cp_GQage

Documented in cp_GQage

#'Creates a bar chart showing the Group Quarters population by age
#'
#'This function creates a ggpplot2 chart that shows the GQ population
#'by age at the 2010 Census.  The chart
#'gives an idea of how important GQ is for an age group in a place.
#'
#' @param fips is the fips code for the county to be charted
#' @param base is the base text size for the ggplot2 object and codemog_theme()


cp_GQage=function(fips, state="08", base=12){
require(codemog, quietly = TRUE)
require(ggplot2, quietly = TRUE)
require(grid, quietly = TRUE)
require(gridExtra, quietly = TRUE)
require(scales, quietly = TRUE)
require(dplyr, quietly = TRUE)
require(tidyr, quietly = TRUE)



d10=codemog_api(data="pco1",db="c2010",geonum=paste("1",state , fips,sep=""),meta="no")
d10[,7:ncol(d10)]=as.numeric(as.character(d10[,7:ncol(d10)]))
d10t=codemog_api(data="p12",db="c2010",geonum=paste("1",state , fips,sep=""),meta="no")
d10t[,7:56]=as.numeric(as.character(d10t[,7:56]))

gq=d10%>%
  mutate(age1=pco1003+pco1022,
         age2=pco1004+pco1023,
         age3=pco1005+pco1024,
         age4=pco1006+pco1025,
         age5=pco1007+pco1026,
         age6=pco1008+pco1027,
         age7=pco1009+pco1028,
         age8=pco1010+pco1029,
         age9=pco1011+pco1030,
         age10=pco1012+pco1031,
         age11=pco1013+pco1032,
         age12=pco1014+pco1033,
         age13=pco1015+pco1034,
         age14=pco1016+pco1035+pco1017+pco1036+pco1018+pco1037+pco1019+pco1038+pco1020+pco1039)%>%
  select(geoname:geonum, age1:age14)%>%
  gather(variable, value, -geoname:-geonum)%>%
  mutate(agecat=ordered(as.factor(variable), levels=c("age1", "age2", "age3", "age4",
                                                      "age5", "age6", "age7", "age8", "age9",
                                                      "age10","age11","age12", "age13", "age14"),
                        labels=c("Under 5\nyears", "5 to 9", "10 to 14","15 to 19" ,"20 to 24", "25 to 29","30 to 34",
                                 "35 to 39", "40 to 44", "45 to 49", "50 to 54", "55 to 59",
                                 "60 to 64", "65\nand Over")),
         population="GQ")
tp=d10t%>%
  mutate(age1=p12003+p12027,
         age2=p12004+p12028,
         age3=p12005+p12029,
         age4=p12006+p12030+p12007+p12031,
         age5=p12008+p12009+p12010+p12032+p12033+p12034,
         age6=p12011+p12035,
         age7=p12012+p12036,
         age8=p12013+p12037,
         age9=p12014+p12038,
         age10=p12015+p12039,
         age11=p12016+p12040,
         age12=p12017+p12041,
         age13=p12018+p12019+p12042+p12043,
         age14=p12020+p12021+p12022+p12023+p12024+p12025+
           p12044+p12045+p12046+p12047+p12048+p12049)%>%
  select(geoname:geonum,age1:age14)%>%
  gather(variable, value, -geoname:-geonum)%>%
  mutate(agecat=ordered(as.factor(variable), levels=c("age1", "age2", "age3", "age4",
                                                      "age5", "age6", "age7", "age8", "age9",
                                                      "age10","age11","age12", "age13", "age14"),
                        labels=c("Under 5\nyears", "5 to 9", "10 to 14","15 to 19" ,"20 to 24", "25 to 29","30 to 34",
                                 "35 to 39", "40 to 44", "45 to 49", "50 to 54", "55 to 59",
                                 "60 to 64", "65\nand Over")),
    population="Total")
hh=tp%>%
  spread(population, value)%>%
  inner_join(spread(gq, population, value))%>%
  mutate(hh=Total-GQ)


p=hh%>%
  ggplot()+
  geom_bar(aes(x=agecat, y=Total, fill=rgb(31,74,126, max=255)),stat="identity")+
  geom_bar(aes(x=agecat, y=GQ, fill=rgb(191,32,38, max=255)),stat="identity")+
  scale_y_continuous(label=comma)+
  scale_fill_identity(name="Population Type", guide="legend", labels=c("Total Population","Group Quarters Population"))+
  theme_codemog(base_size=base)+
  theme(legend.text=element_text(size=rel(1.05)), axis.text=element_text(size=rel(1.10)))+
  labs(x="Age", y="Population", title=paste(d10$geoname, "Population by Age and Type, 2010 \nSource: U.S. Census Bureau"))
return(p)


}
ColoradoDemography/codemogProfile documentation built on Aug. 30, 2017, 5 a.m.