R/hei_2015_PerPerson_ssum.R

In yikeshu0611/nhanesR: Tools for "NHANES" Database

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

    # day 1
    fped <- fped_read(years = years,day = '1',dietary = dietary,version = 2015)
    tsv <- nhs_tsv(sprintf('drx%s|dr1%s',dietary,dietary),years = years,cat = FALSE)
    dt <- nhs_read(tsv,
                   "drxtkcal,drxikcal,dr1tkcal,dr1ikcal,dr2tkcal,dr2ikcal:kcal",
                   "drxtsfat,drxisfat,dr1tsfat,dr1isfat,dr2tsfat,dr2isfat:sfat",
                   "drdtsodi,drdisodi,dr1tsodi,dr1isodi,dr2tsodi,dr2isodi:sodi",
                   "drxtmfat,drximfat,dr1tmfat,dr1imfat,dr2tmfat,dr2imfat:mfat",
                   "drxtpfat,drxipfat,dr1tpfat,dr1ipfat,dr2tpfat,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

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

    # day 2
    fped <- fped_read(years = years,day = '2',dietary = dietary,version = 2015)
    tsv <- nhs_tsv(sprintf('drx%s|dr2%s',dietary,dietary),years = years,cat = FALSE)
    dt <- nhs_read(tsv,
                   "drxtkcal,drxikcal,dr1tkcal,dr1ikcal,dr2tkcal,dr2ikcal:kcal",
                   "drxtsfat,drxisfat,dr1tsfat,dr1isfat,dr2tsfat,dr2isfat:sfat",
                   "drdtsodi,drdisodi,dr1tsodi,dr1isodi,dr2tsodi,dr2isodi:sodi",
                   "drxtmfat,drximfat,dr1tmfat,dr1imfat,dr2tmfat,dr2imfat:mfat",
                   "drxtpfat,drxipfat,dr1tpfat,dr1ipfat,dr2tpfat,dr2ipfat:pfat",
                   codebook = FALSE,varLabel = FALSE,cat = FALSE,Year = FALSE)
    dt$monopoly <- dt$mfat + dt$pfat
    colnames(dt) <- rename_line(colnames(dt))
    colnames(dt) <- rename_fdcd(colnames(dt))
    dt <- drop_col(dt,'fdcd')
    indat2 <- dplyr::inner_join(dt,fped,join)

    indat <- dplyr::inner_join(indat1,indat2,join)
    choice <- colnames(indat2)[-c(1:length(join))]
    for (i in 1:length(choice)) {
        which <- which(colnames(indat) %in% paste0(choice[i],c('.x','.y')))
        indat$last <- rowSums(indat[,which],na.rm = TRUE)
        indat <- indat[,-which]
        colnames(indat)[ncol(indat)] <- choice[i]
    }
    hei_2015(
        indat = indat,
        kcal = 'kcal',
        vtotalleg = 'vtotalleg',
        vdrkgrleg = 'vdrkgrleg',
        f_total = 'f_total',
        fwholefrt = 'f_whole',
        g_whole = 'g_whole',
        d_total = 'd_total',
        pfallprotleg = 'pfallprotleg',
        pfseaplantleg = 'pfseaplantleg',
        monopoly = 'monopoly',
        satfat = 'sfat',
        sodium = 'sodi',
        g_refined = 'g_refined',
        add_sugars = 'add_sugars',
        varLabel = varLabel,
        energy = energy,
        component = component,
        density = density,join=join
    )
}