Health Clusters Regression

#Create new variables based on clusters identified in EGA
acute.clean <- transform(acute.clean,
                         HEALTH.1 = HCCARVAS + HCDIAB + HCHIBP + HCHICHL,
                         HEALTH.2 = HCALZHEIM + HCPRESULC + HCORALDENTAL + HCDYSPH,
                         MENTAL.HEALTH = DXANXIETY + DXMOOD + DXPSYCHOTIC + DXMIOTHER,
                         BEHAVIOR = DISBEH + SELFINJ + UNCPBEH + DXBEHAVIOR)

#Poisson regression with groups as predictor variables
model.groups <- glm(ER_Total~GENDER + AGE+ RACE_BLACK + RACE_LATINO + LEVELID + EXPRESS +  HEALTH.1 + HEALTH.2 + MENTAL.HEALTH + BEHAVIOR + TYPEHOME, family = "poisson", data = acute.clean)
summary(model.groups)
summ(model.groups)