popEpi包裹按年龄分层的SIR

Question

我试图使用popEpi软件包按已达到的年龄来计算SIRs，我想按达到的年龄(0-5；5-10;10-15 . 60-65；65+)将结果分层。我想不出该怎么做。为了根据后续时间对结果进行分层，我使用了以下语法：

c <- lexpand( DatasetSMN, status = status, birth = bi_date, exit = ex_date, entry = dg_date,
               breaks = list(per = 1960:2013, age = 1:100, fot = c(0,5,10,15,20,25,30,35,Inf)), 
               aggre = list(fot, agegroup = age, year = per, sex) )

se <- sir( coh.data = c, coh.obs = 'from0to1', coh.pyrs = 'pyrs',
           ref.data = incdata, ref.rate = 'value', conf.type = "univariate", 
           adjust = c('agegroup','year','sex'), print ='fot')

这是可行的，但我不知道如何适应这一点，以获得按年龄分层的结果。

提前感谢！

Answer 1

可以将产出按已达到的年龄划分如下：

update(se, print = list(agegroup))

通过执行以下操作，也可以跨多个变量做到这一点：

update(se, print = list(sex, agegroup))

或者，我们甚至可以做一些不同的事情：

update(se, print = list(age.cat = cut(agegroup, c(0,50,75,100), dig.lab=4)))