# NOT RUN {
library(dplyr)
df <- data.frame(indicatorid = rep(c(1234, 5678, 91011, 121314), each = 19 * 2 * 5),
year = rep(2006:2010, each = 19 * 2),
sex = rep(rep(c("Male", "Female"), each = 19), 5),
ageband = rep(c(0,5,10,15,20,25,30,35,40,45,
50,55,60,65,70,75,80,85,90), times = 10),
obs = sample(200, 19 * 2 * 5 * 4, replace = TRUE),
pop = sample(10000:20000, 19 * 2 * 5 * 4, replace = TRUE))
refdf <- data.frame(refcount = sample(200, 19, replace = TRUE),
refpop = sample(10000:20000, 19, replace = TRUE))
df %>%
group_by(indicatorid, year, sex) %>%
phe_smr(obs, pop, refdf$refcount, refdf$refpop, type="standard")
## OR
df %>%
group_by(indicatorid, year, sex) %>%
phe_smr(obs, pop, refdf$refcount, refdf$refpop, confidence=99.8, refvalue=100)
# }
Run the code above in your browser using DataCamp Workspace