在R內進行直接年齡標準化 (Direct Age Standardization in R)

    年齡標準化在流行病學的應用很廣泛. 

    它主要用於: 在各地之間進行率的比較時, 而這些率又受到年齡因素影響很大, 各地之間年齡段的比例又有差異時所用.

    年齡的標準化又分為直接標準化和間接標準化兩種. 直接標準化常用在"外在比", 如: 各國之間的比較; 間接標準化則常用在"內在比", 如某國內各省之間的比較... 又以直接標準化為多用!

    由於各地人口之年齡段比例有差異, 所以要先找出一個公認的標準作為加權數. 在世界各地的比較上, 常用WHO的年齡段比例. 中國常用1982年的人口構成.

年龄组（岁） 1964年人口构1 1982年中国人口构成2 0 14.52 2.07 1～4 7.36 5～9 13.65 11.03 10～14 12.52 13.13 15～19 9.01 12.49 20～24 7.37 7.41 25～29 7.37 9.22 30～34 6.77 7.27 35～39 5.97 5.4 40～44 5.17 4.82 45～49 4.47 4.72 50～54 3.84 4.07 55～59 3.27 3.38 60～64 2.56 2.73 65～69 1.69 2.12 70～74 1.07 1.43 75～79 0.55 0.86 80～84 0.26 0.37 85及以上 0.13 合计 100 100

    在R有”epitools”軟件包, 內有功能是專門計算直接年齡標準化率的!

library(readxl)
mor<-read_excel("ageadj.xlsx") #讀入excel文件
library(epitools) #載入epitools的統計包
whostd.prec<-mor$who/100 #將WHO的標準年齡段百分率化回小數
mor.age<-ageadjust.direct(mor$death,mor$pop,stdpop = whostd.prec) #計算直接年齡標準化率
mor.aa<-round(mor.age*10^5,1) #將率化回10萬

mor文件的內容, who: 世衛的年齡段百分率, death: 死亡人數, pop: 總人數

結果

crude.rate   adj.rate        lci        uci 
     887.3      496.2      494.1      498.3 

R的計算結果與Excel的是一樣呵~

 參考資料:

1.http://www.who.int/healthinfo/paper31.pdf