到中的实际死亡数D与定一个标准的死亡人数ED与E之比称为死亡比standarized

Mortality

今天我们介绍一下怎么通过手动计算校正年龄、性别后的标准化死亡率

我们先导入R包和数据也不能算完全手动还是需要用到survival包。

library(survival)

bc-read.csv(E:/r/test/smr1.csv,sep,,headerTRUE)

head(bc,6)

1821.558这是一个很简单的数据sex为性别age是年龄entry_date为诊断也就是进入这个队列的时间status为结局变量futime为生存时间。

公众号回复SMR1可以获得该数据

先处理一下生存时间

bc$entry_date-as.Date(bc$entry_date)我们首先要算出它的haz

每人年的平均人口死亡率d/(pyrs)其中d是死亡人数pyrs是人年可以通过survival包的survexp函数计算。

sex这个指标不能少如果全是女的也要设置.

bc$risk1

TRUE))survexp函数中选项还有一个费率表就是用来调整生成的预测概率的可以理解为校正功能。

survexp.usr就是1940年至2014年按年龄、性别和种族划分的美国人口。

这里有黑人和白人的数据集我这里只取白人的

survexp.uswhite

如果你想算某个地方的生存率可以使用当地人口普查数据这里我们把费率表加进去。

bc$risk2

0.007464535可以看到加入费率表的概率和不加的是不一样的。

算出了预测概率后我们就可以进一步计算了。

sum(bc$status)

9.101355这样全部结果就计算出来啦计算结果我们使用survexp.fr包来验证一下

library(survexp.fr)

bc$entry_date-as.Date(bc$entry_date)

SMR(futime,