# See https://academic.oup.com/ije/article/46/1/348/2622842

require(Hmisc)
d <- csv.get("sicily.csv")
d <- upData(d,
            labels=c(year    = 'Year',
                     month   = 'Month',
                     aces    = 'Acute Coronary Episodes per Month',
                     smokban = 'Smoking Ban',
                     pop     = 'Sicily Population in 10000s',
                     stdpop  = 'Age Standardized Population',
                     time    = 'Time Since Start of Study',
                     rate    = 'Standardized Rate'),
            units=c(time='months'),
            rate=aces / stdpop * 10^5)
head(d)
contents(d)
describe(d)
sicily <- d
Save(sicily)
cat(html(contents(d)), file='sicily.html')