Liping prevalences - Comparisons of Oesity in WIC and PEDs

Overall Plot

qplot(year,value, data=dmoverall, group=variable, color=variable) + 
  geom_line() + theme_df() + scale_x_continuous(breaks=seq(2000,2010,2)) +
   # scale_colour_manual(values=cbbPalette) + 
  scale_color_colorblind() +
  ylab("Obesity Prevalence (%)\n") + xlab('\nYear') + 
  theme(axis.text.x = element_text(angle=45))

Plot by age

This shows the prevalence by study (2 lines) by age (3 panels). Note that the WIC excess tends to increase over the time period, and that the difference is most evident among 2-year-olds.

Since the early 1970s, the prevalence of obesity in the U.S. has increased by more than a factor of 3 among 6- to 19-year-olds and by about 70% (1,2) among 2- to 5-year-olds. In response to these secular trends, much data on the weight and height of children have been collected in schools (3–7) and clinics (8–11). Because the quality of these measured weight and height data can vary substantially, investigators must decide how to address extreme values that could represent erroneous data. This process would influence estimates of the prevalence of obesity (12), trends over time, and associations with other characteristics.

qplot(year,value, data=dage, group=variable, color=variable) + geom_line() +
  facet_wrap(~age) + theme_df() + 
  scale_x_continuous(breaks=seq(2000,2010,2)) +
  scale_colour_manual(values=cbbPalette) + 
  ylab("Obesity Prevalence (%)\n") + xlab('\nYear') + 
  theme(axis.text.x = element_text(angle=45))

Plot by sex

dx=set(dm, sex%in% cs(boys,girls)); #dx
qplot(year,value, data=dx, group=variable, color=variable) + geom_line() +
  facet_wrap(~sex) + theme_df() + 
  scale_x_continuous(breaks=seq(2000,2010,2)) + 
  scale_colour_manual(values=cbbPalette) +
  theme(axis.text.x = element_text(angle=45))

Plot by race

dx=set(dm, race %in% cs(amer_in,asian,blacks,hisp,whites)); #dx
qplot(year,value, data=dx, group=variable, color=variable) + geom_line() +
  facet_wrap(~race) + theme_df() + 
  scale_x_continuous(breaks=seq(2000,2010,2)) +
  scale_color_colorblind() +
  theme(axis.text.x = element_text(angle=45, size = 11))