NHANES(National Health and Nutrition Examination Survey,国家健康和营养检查调查)中的饮食炎症指数(Dietary Inflammatory Index, DII)是一种用于评估个人饮食对体内炎症水平影响的工具。DII通过量化个体饮食中各种成分的抗炎或促炎特性来计算一个综合分数。较高的DII分数意味着饮食可能促进更多的炎症,而较低的分数则表示饮食具有更强的抗炎作用。
在既往文章中,我介绍了使用dietaryindex包膳食指数,有不少粉丝私信问如何计算炎症指数,这个其实挺简单的,今天简单介绍一下把
先导入R包和数据,
library(dietaryindex)
data("DHQ3_exp")
这个是数据使用dietaryindex包自带的数据。第一行是标签,第二个是记录次数,如果没有多次记录就默认是1次,sex是性别,age是年龄,后面都是一些维生素,食物或者微量元素的摄入量
dietaryindex包计算饮食炎症指数主要是通过DII函数来计算,它的函数如下
DII(
SERV_DATA,
RESPONDENTID,
REPEATNUM = 1,
ALCOHOL_DII = NULL,
VITB12_DII = NULL,
VITB6_DII = NULL,
BCAROTENE_DII = NULL,
CAFFEINE_DII = NULL,
CARB_DII = NULL,
CHOLES_DII = NULL,
KCAL_DII = NULL,
EUGENOL_DII = NULL,
TOTALFAT_DII = NULL,
FIBER_DII = NULL,
FOLICACID_DII = NULL,
GARLIC_DII = NULL,
GINGER_DII = NULL,
IRON_DII = NULL,
MG_DII = NULL,
MUFA_DII = NULL,
NIACIN_DII = NULL,
N3FAT_DII = NULL,
N6FAT_DII = NULL,
ONION_DII = NULL,
PROTEIN_DII = NULL,
PUFA_DII = NULL,
RIBOFLAVIN_DII = NULL,
SAFFRON_DII = NULL,
SATFAT_DII = NULL,
SE_DII = NULL,
THIAMIN_DII = NULL,
TRANSFAT_DII = NULL,
TURMERIC_DII = NULL,
VITA_DII = NULL,
VITC_DII = NULL,
VITD_DII = NULL,
VITE_DII = NULL,
ZN_DII = NULL,
TEA_DII = NULL,
FLA3OL_DII = NULL,
FLAVONES_DII = NULL,
FLAVONOLS_DII = NULL,
FLAVONONES_DII = NULL,
ANTHOC_DII = NULL,
ISOFLAVONES_DII = NULL,
PEPPER_DII = NULL,
THYME_DII = NULL,
ROSEMARY_DII = NULL
)
SERV_DATA填入数据,RESPONDENTID填入ID,如果没有重复记录REPEATNUM这里填入1,之后看你收集了那个指标,对应录入就行
out<-DII(SERV_DATA=DHQ3_exp, RESPONDENTID=DHQ3_exp$`Respondent ID`, REPEATNUM = 1,
ALCOHOL_DII =DHQ3_exp$`Alcohol (g)`, VITB12_DII=DHQ3_exp$`Vitamin B12 (mcg)`, VITB6_DII =DHQ3_exp$`Vitamin B6 (mg)`)
如果你只是输入部分指标,看见警告是正常的,不必担心
一句话代码,结果就出来啦,总的指标炎症指数和单独指标的炎症指数都出来啦,非常简单。
如果想了解怎么计算的可以看这篇文章,有给出计算公式。
如果感兴趣的多,有空我也可以手动计算一下