R/leg2010a.R

In yikeshu0611/nhanesR: Tools for "NHANES" Database

leg2010a <- function(indat,kcal,allmeat,seaplant,v_total,v_drkgr,legumes){
    indat$mbmax=2.5*(indat[,kcal]/1000)

    ck1 <- indat[,allmeat] < indat[,"mbmax"]
    indat[ck1,"meatleg"] <- indat[ck1,legumes]*4
    indat[ck1,"needmeat"] <- indat[ck1,"mbmax"]-indat[ck1,allmeat]
    ck2 <- indat[,"meatleg"] <= indat[,"needmeat"]
    ck12 <- ck1 & ck2
    indat$legtype[ck12] <- 'allmeat'
    indat$legume_added_allmeat[ck12] <- indat[ck12,allmeat]+indat[ck12,"meatleg"]
    indat$legume_added_seaplant[ck12] <- indat[ck12,seaplant]+indat[ck12,"meatleg"]
    indat$legume_added_v_total[ck12]<-indat[ck12,v_total]
    indat$legume_added_beangrn[ck12]<-indat[ck12,v_drkgr]
    ck1n2 <- ck1 & !ck2
    indat$legtype[ck1n2] <- 'meat/veg'
    indat$extrmeat[ck1n2]<-indat[ck1n2,"meatleg"]-indat[ck1n2,"needmeat"]
    indat$extrleg[ck1n2]<-indat[ck1n2,"extrmeat"]/4
    indat$legume_added_allmeat[ck1n2]<-indat[ck1n2,allmeat]+indat[ck1n2,"needmeat"]
    indat$legume_added_seaplant[ck1n2] <- indat[ck1n2,"seaplant"]+indat[ck1n2,"needmeat"]
    indat$legume_added_v_total[ck1n2]<-indat[ck1n2,v_total]+indat[ck1n2,"extrleg"]
    indat$legume_added_beangrn[ck1n2]<-indat[ck1n2,v_drkgr]+indat[ck1n2,"extrleg"]
    ckn1 <- !ck1
    indat$legtype[ckn1] <- 'allveg'
    indat$legume_added_allmeat[ckn1] <- indat[ckn1,allmeat]
    indat$legume_added_seaplant[ckn1] <- indat[ckn1,"seaplant"]
    indat$legume_added_v_total[ckn1] <- indat[ckn1,v_total]+indat[ckn1,legumes]
    indat$legume_added_beangrn[ckn1] <- indat[ckn1,v_drkgr]+indat[ckn1,legumes]

    indat
}