# basic checks, and code to share with Yimeng
# devtools::install_github("donboyd5/pp.prototypes")
library("pp.prototypes")
library("plyr")
library("dplyr")
library("tidyr")
library("ggplot2")
library("knitr")
#data(package="pp.prototypes")
?actives
?retirees
?salgrowth.hist
?salgrowth.assume
glimpse(actives)
count(actives, planname)
count(actives, age)
count(actives, ea) %>% data.frame
actives %>% filter(planname=="average") %>% select(age, ea, nactives) %>% spread(ea, nactives) %>% kable(digits=2)
actives %>% filter(planname=="average") %>% select(age, ea, salary) %>% spread(ea, salary) %>% kable(digits=0)
actives %>% filter(planname=="underfunded") %>% select(age, ea, nactives) %>% spread(ea, nactives) %>% kable(digits=2)
actives %>% filter(planname=="underfunded") %>% select(age, ea, salary) %>% spread(ea, salary) %>% kable(digits=0)
actives %>% group_by(planname) %>%
summarise(age=sum(age*nactives) / sum(nactives),
nactives.sum=sum(nactives), # rename so that it does not mess up next calculation
salary=sum(salary*nactives) / sum(nactives))
tmp <- actives %>% group_by(planname, age) %>%
summarize(nactives.sum=sum(nactives), salary=sum(nactives*salary) / nactives.sum)
qplot(age, nactives.sum, data=tmp, colour=planname, geom=c("point", "line"))
qplot(age, salary, data=tmp, colour=planname, geom=c("point", "line")) # a little odd
glimpse(retirees)
count(retirees, planname)
count(retirees, age)
retirees %>% group_by(planname) %>%
summarise(age=sum(age*nretirees) / sum(nretirees),
nretirees.sum=sum(nretirees), # rename so that it does not mess up next calculation
benefit=sum(benefit*nretirees) / sum(nretirees))
tmp <- retirees %>% group_by(planname, age) %>%
summarize(nretirees.sum=sum(nretirees), benefit=sum(nretirees*benefit) / nretirees.sum)
qplot(age, nretirees.sum, data=tmp, colour=planname, geom=c("point", "line"))
qplot(age, benefit, data=tmp, colour=planname, geom=c("point", "line")) # a little odd
glimpse(salgrowth.hist)
qplot(age, sscale.hist.rate, data=salgrowth.hist, colour=planname, geom=c("point", "line"))
glimpse(salgrowth.assume)
qplot(age, sscale.assume.rate, data=salgrowth.assume, colour=planname, geom=c("point", "line"))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.