R语言NHANES数据分析(1)

#NHANES数据分析
library(maditr)
library(nhanesdata)
data_a<-nhanesdata::demo %>%
  left_join(.,nhanesdata::bmx,by="SEQN") %>%
  left_join(.,nhanesdata::bpq,by="SEQN")
colnames(data_a)
data_b=data_a[,c("RIDAGEYR","DMDMARTL","SDMVSTRA","SDMVPSU","WTINT2YR","DMDEDUC2","DMDCITZN","BMXBMI","BPQ020")] %>% na.omit()
#install.packages("finalfit")
data_c<-subset(data_b,data_b$BPQ020 %in% c(1,2))
data_c$BPQ020<-data_c$BPQ020-1

#基本分析
library(survey)
svy_design <- survey::svydesign(strata= ~SDMVSTRA, id = ~SDMVPSU, weights = ~WTINT2YR, nest = TRUE,data = data_c)
result<-survey::svyglm(BPQ020~BMXBMI, family = binomial(), design = svy_design)

#单因素和多因素分析
library(finalfit)
dep="BMXBMI"
exp=c("DMDEDUC2","RIDAGEYR")
uni=svy_design%>%
  svyglmuni(dep,exp)%>%
  fit2df(estimate_suffix="(univariable)")
multi=svy_design%>%
  svyglmmulti(dep,exp)%>%
  fit2df(estimate_suffix="(multivariable)")
data_c%>%
  summary_factorlist(dep,exp,fit_id=TRUE)%>%
  ff_merge(uni)%>%
  ff_merge(multi)%>%
  select(-fit_id,-index)%>%
  dependent_label(data_c,dep)->table1
table1