R/hei_2010_PerDay_ssum.R

In yikeshu0611/nhanesR: Tools for "NHANES" Database

hei_2010_PerDay_ssum <- function(years,
                                 day='2',
                                 dietary='tot',
                                 varLabel=FALSE,
                                 energy=TRUE,
                                 component=TRUE,
                                 density=FALSE,
                                 version=2010){
  if (dietary=='iff') join <- c('seqn','line') else join <- 'seqn'

  fped <- fped_read(years = years,day = day,dietary = dietary,version=2010)

  tsv <- nhs_tsv(ifelse(day=='1',
                        sprintf('drx%s|dr1%s',dietary,dietary),
                        sprintf('drx%s|dr2%s',dietary,dietary)),years = years,cat = FALSE)
  dt <- nhs_read(tsv,
                 "drxtkcal,drxikcal,dr1tkcal,dr2tkcal,dr1ikcal,dr2ikcal:kcal",
                 "drxtsfat,drxisfat,dr1tsfat,dr2tsfat,dr1isfat,dr2isfat:sfat",
                 "drxtalco,drxialco,dr1talco,dr2talco,dr1ialco,dr2ialco:alco",
                 "drdtsodi,drdisodi,dr1tsodi,dr1isodi,dr2tsodi,dr2isodi:sodi",
                 "drxtmfat,drximfat,dr1tmfat,dr2tmfat,dr1imfat,dr2imfat:mfat",
                 "drxtpfat,drxipfat,dr1tpfat,dr2tpfat,dr1ipfat,dr2ipfat:pfat",
                 codebook = FALSE,varLabel = FALSE,cat = FALSE)
  colnames(dt) <- rename_line(colnames(dt))
  colnames(dt) <- rename_fdcd(colnames(dt))
  dt <- drop_col(dt,'fdcd')

  dt$monopoly <- dt$mfat + dt$pfat
  dt$maxalcgr <- 13*(dt$kcal/1000)
  dt$exalccal <- ifelse(dt$alco <= dt$maxalcgr,0,7*( dt$alco - dt$maxalcgr))


  indat <- dplyr::inner_join(dt,fped,join)

  indat$emptycal10 <- indat$addsugc+indat$solfatc+indat$exalccal


# leg2010a
  afterleg <- leg2010a(indat = indat,
                       kcal = 'kcal',
                       allmeat ='pf_total',
                       seaplant = 'seaplant',
                       v_total = 'v_total',
                       v_drkgr = 'v_drkgr',
                       legumes = 'v_legumes')

# hei_2010
  hei_2010(
    indat=afterleg,
    kcal = 'kcal',
    lv_total = 'legume_added_v_total',
    lbeangrn = 'legume_added_beangrn',
    f_total = 'f_total',
    wholefrt = 'f_whole',
    g_whl = 'g_whole',
    d_total = 'd_total',
    lallmeat = 'legume_added_allmeat',
    lseaplant = 'legume_added_seaplant',
    monopoly = 'monopoly',
    sfat = 'sfat',
    sodi = 'sodi',
    g_nwhl = 'g_refined',
    emptycal10 = 'emptycal10',
    varLabel = varLabel,
    energy = energy,
    component = component,
    density = density,
    join=join
  )
}