Blog 3

hardtoget<-haven::read_xpt("/Users/christacrumrine/Downloads/LLCP2020.XPT ")

names(hardtoget)<-tolower(gsub(pattern = "_", replacement = "",x=names(hardtoget)))

Data Preperation

#Age cut into intervals
hardtoget$age<-cut(hardtoget$ageg5yr,
                   breaks=c(0,24,39,59,79,99))


#Hours of time spent as provider 
hardtoget$hoursofcare<-Recode(hardtoget$crgvhrs1, recodes="1=1; 2=2; 3=3; 4=4; else=NA", as.factor = T)
hardtoget$hoursofcare<-relevel(hardtoget$hoursofcare, ref = "1")

#has health coverage
hardtoget$healthinsurace_coverage<- car::Recode(hardtoget$hlthpln1,
                       recodes="1='no, has coverage'; 2='yes, no coverage';else=NA",
                       as.factor=T)
hardtoget$healthinsurace_coverage<-relevel(hardtoget$healthinsurace_coverage, ref='yes, no coverage')


#sex
hardtoget$male<-as.factor(ifelse(hardtoget$colgsex==1,
                                "Male",
                                "Female"))

#marital status
hardtoget$marst<-Recode(hardtoget$marital, recodes="1='married'; 2='divorced'; 3='widowed'; 4='separated'; 5='nm';6='cohab'; else=NA", as.factor=T)


#depression
hardtoget$depression<-Recode(hardtoget$addepev3, recodes="1=1; 2=0; else=NA")



#assist personal care
hardtoget$ADLcare<-Recode(hardtoget$crgvper1, recodes="1=1; 2=0; else=NA")

#assist home tasks
hardtoget$IADLcare<-Recode(hardtoget$crgvhou1, recodes="1=1; 2=0; else=NA")

#Poor or fair self rated health
hardtoget$badhealth<-Recode(hardtoget$genhlth,
                           recodes="4:5=1; 1:3=0; else=NA")

#Drink alcohol
hardtoget$drink<-Recode(hardtoget$drnkany5, recodes="1=1; 2=0; else=NA")

#length caring for someone 
hardtoget$cargving<-Recode(hardtoget$crgvlng1,
recodes="1='0less30'; 2='1to6mnths'; 3='2sixmnthsto2yrs'; 4='3twofvyears'; 5='4fvyrsormore';7:9=NA",
as.factor=T)
hardtoget$cargving<-relevel(hardtoget$cargving, ref='0less30')


#have you been told you have diabetes
hardtoget$hvdiab<-Recode(hardtoget$prediab1,
                       recodes="1:2 ='yes'; 3 ='no'; else = NA",
                       as.factor=T)

#Ever had a stroke
hardtoget$stroke<-Recode(hardtoget$cvdstrk3, recodes="1=1; 2=0; else=NA")

#smoked
hardtoget$smoked<-Recode(hardtoget$smoke100, recodes="1=1; 2=0; else=NA")


#race/ethnicity
hardtoget$black<-Recode(hardtoget$imprace,
                       recodes="2=1; 9=NA; else=0")
hardtoget$white<-Recode(hardtoget$imprace,
                       recodes="1=1; 9=NA; else=0")
hardtoget$other<-Recode(hardtoget$imprace,
                       recodes="3:4=1; 9=NA; else=0")
hardtoget$hispanic<-Recode(hardtoget$imprace,
                          recodes="5=1; 9=NA; else=0")

hardtoget$race_eth<-Recode(hardtoget$imprace,
                          recodes="1='nhwhite'; 2='nh_black'; 3='nh_other';4='nh_multirace'; 5='hispanic'; else=NA",
                          as.factor = T)

#education level
hardtoget$educ<-Recode(hardtoget$educa,
recodes="1:2='0Prim'; 3='1somehs'; 4='2hsgrad'; 5='3somecol'; 6='4colgrad';9=NA",
as.factor=T)
hardtoget$educ<-relevel(hardtoget$educ, ref='2hsgrad')

#BMI, in the brfss_17a the bmi variable has 2 implied decimal places,
#so we must divide by 100 to get real bmi's

hardtoget$bmi<-hardtoget$bmi5/100

#Healthy days
hardtoget$healthdays<-Recode(hardtoget$menthlth, recodes = "88=0; 77=NA; 99=NA")

Research Questions:

Do people who act as a caregiver for someone report higher levels of depression?

Do people who assist with activities of daily living report higher levels than those who provide supervision or assist with only instrumental activities of daily living?

What effect does income have on the number of poor rated health days for someone who act as a caregiver for someone?