R语言NHANES数据分析（2）

#交互分析
library(jstable)
help(package="nhanesA")
help(package="jstable")
str(data_c)
svy_design$variables$DMDCITZN<-as.factor(svy_design$variables$DMDCITZN)
TableSubgroupMultiGLM(BPQ020~BMXBMI,data=svy_design,var_subgroups = c("DMDCITZN"))

#中介分析
library(mediation)
#install.packages("mma")
help(package="mma")
svy_design <- survey::svydesign(strata= ~SDMVSTRA, id = ~SDMVPSU, weights = ~WTINT2YR, nest = TRUE,data = data_c)
result<-survey::svyglm(BPQ020~BMXBMI, family = gaussian(), design = svy_design)
result1<-survey::svyglm(DMDCITZN~BMXBMI, family = gaussian(), design = svy_design)
med_out <- mediate(model.m=result1, model.y=result, 
                   treat = "BMXBMI", 
                   mediator = "BPQ020",
                   boot = TRUE)
summary(med_out)

接之前代码，完成交互分析和中介分析。

结果解读：

交互分析：

中介分析：

ACME间接效应，P值大于0.05，无统计学意义。

ADE直接效应，P值小于0.05，有统计学意义，说明BMXBMI对DMDCITZN的影响主要通过直接效应。