Setup

Setup

knitr::opts_chunk$set(echo = TRUE)
knitr::opts_knit$set(root.dir = "D:/Jose")

Load Libraries

options(scipen=20)

require(dplyr);require(doParallel);require(epitools);require(foreach)
## Loading required package: dplyr
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
## Loading required package: doParallel
## Loading required package: foreach
## Loading required package: iterators
## Loading required package: parallel
## Loading required package: epitools
require(foreign);require(kableExtra);require(magrittr);require(MASS)
## Loading required package: foreign
## Loading required package: kableExtra
## 
## Attaching package: 'kableExtra'
## The following object is masked from 'package:dplyr':
## 
##     group_rows
## Loading required package: magrittr
## Loading required package: MASS
## 
## Attaching package: 'MASS'
## The following object is masked from 'package:dplyr':
## 
##     select
require(plyr);require(psych);require(psychTools);require(SASxport)                            
## Loading required package: plyr
## ------------------------------------------------------------------------------
## You have loaded plyr after dplyr - this is likely to cause problems.
## If you need functions from both plyr and dplyr, please load plyr first, then dplyr:
## library(plyr); library(dplyr)
## ------------------------------------------------------------------------------
## 
## Attaching package: 'plyr'
## The following objects are masked from 'package:dplyr':
## 
##     arrange, count, desc, failwith, id, mutate, rename, summarise,
##     summarize
## Loading required package: psych
## Loading required package: psychTools
## Loading required package: SASxport
## 
## Attaching package: 'SASxport'
## The following objects are masked from 'package:foreign':
## 
##     lookup.xport, read.xport
require(survey)
## Loading required package: survey
## Loading required package: grid
## Loading required package: Matrix
## Loading required package: survival
## 
## Attaching package: 'survival'
## The following object is masked from 'package:epitools':
## 
##     ratetable
## 
## Attaching package: 'survey'
## The following object is masked from 'package:graphics':
## 
##     dotchart

Variable Selection

myvars=c("X.AGE.G", "X.IMPRACE","SEX" ,"MARITAL",   "INCOME2",   "EDUCA" ,"EMPLOY1", "VETERAN3" , "X.RFBMI5", "CVDCRHD4", "CVDSTRK3",  "CHCSCNCR" , "CHCOCNCR" , "CHCCOPD1" , "HAVARTH3" , "MENTHLTH", "CHCKIDNY" , "DIABETE3" , "X.STSTR" ,  "X.LLCPWT", "X.STATE", "Year");print(myvars)
##  [1] "X.AGE.G"   "X.IMPRACE" "SEX"       "MARITAL"   "INCOME2"   "EDUCA"    
##  [7] "EMPLOY1"   "VETERAN3"  "X.RFBMI5"  "CVDCRHD4"  "CVDSTRK3"  "CHCSCNCR" 
## [13] "CHCOCNCR"  "CHCCOPD1"  "HAVARTH3"  "MENTHLTH"  "CHCKIDNY"  "DIABETE3" 
## [19] "X.STSTR"   "X.LLCPWT"  "X.STATE"   "Year"

Primary Recoding

We make the following recodes.

.STATE: Not all states/territories are surveyed every year, creating gaps. We combine the FIPS codes for Census subregions to ensure that each region is represented every year.

.RFBMI5: Adults who have a body mass index greater than 25.00 Values: 1=No, 2=Yes, 9=Don’t Know/Refused/Missing Recode: 0=No, 1=Yes or Don’t Know/Refused/Missing (modal imputation, 8% missing)

.AGE.G: Six-level imputed age category Values: 1=18 to 24, 2=25 to 34, 3=35 to 44, 4=45 to 54, 5=55 to 64, 6=65 or older

.IMPRACE: Imputed race/ethnicity value Values: 1=White Non-Hispanic, 2=Black Non-Hispanic, 3=Asian Non-Hispanic, 4=American Indian/Alaskan Native Non-Hispanic, 5=Hispanic, 6=Other Race, Non-Hispanic Recode: 0=White, 1=Non-White

HISPANIC: Derived from IMPRACE Recode: 0=Non-Hispanic, 1=Hispanic

.SEX: Indicate sex of respondent Values: 1=Male, 2=Female, 9=Refused (.06%) Recode: 0=Not Known to Be Male, 1=Known to be Male (modal imputation)

.MARITAL: Are you..? Values: 1= Married, 2=Divorced, 3=Widowed, 4=Separated, 5=Never Married, 6=Unmarried Couple, 9=Refused, BLANK=Not asked or Missing Recode: 0=Did not Identify as Married, 1= Self-Identified as Married

.INCOME2: Is your annual household income from all sources: Values: 1=LT 10K, 2=LT 15K, 3=LT 20K, 4=LT 25K, 5=LT 35K, 6=LT 50K, 7=LT 75K, 8=GE 75K, 77=Don’t Know/Not Sure, 99= Refused, Not Asked or Missing Recode: 0=Not Identified as 75K or more, 1=Identified as 75K or more

EDUCA: What is the highest grade or year of school you completed? Values: 1=Never Attended School or Only Kindergarten, 2=Grades 1-8, 3=Grades 9-11, 4=Grade 12 or GED, 5=College 1-3 Years, 6=College 4 or more (graduate), 9=Refused, BLANK=Not Asked/Missing Recode: 0=Not Identified as College Grad, 1=College Graduate

EMPLOY1: Are you currently…? Values: 1=Employed for Wages, 2= Self-Employed, 3=Out of work >1 year, 4=Out of Work< 1 year, 5=A homemaker, 6= A student, 7=Retire, 8=Unable to Work, 9=Refused, BLANK=Not Asked/Missing Recode: 0=Did not Identify as Employed for Wages, 1=identified as Employed for Wages

VETERAN3: Have you ever served on active duty in the United States Armed Forces, either in the regular military or in a National Guard or military reserve unit? Values: 1=Yes, 2=No, 7=Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing Recode: 0=Not Yes, 1=Yes

The next several variables have the same recode. CVDCRHD4:(Ever told) you had angina or coronary heart disease? CVDSTRK3:(Ever told) you had a stroke. CHCSCNCR:(Ever told) you had skin cancer? CHCONCR: (Ever told) you had any other types of cancer? CHCCOPD: (Ever told) you have chronic obstructive pulmonary disease, C.O.P.D., emphysema or chronic bronchitis? HAVARTH3:(Ever told) you have some form of arthritis, rheumatoid arthritis, gout, lupus, or fibromyalgia? (Arthritis diagnoses include: rheumatism, polymyalgia rheumatica; osteoarthritis (not osteporosis); tendonitis, bursitis, bunion, tennis elbow; carpal tunnel syndrome, tarsal tunnel syndrome; joint infection, etc.) CHCKIDNY: (Ever told) you have kidney disease? Do NOT include kidney stones, bladder infection or incontinence.(Incontinence is not being able to control urine flow.) Values: 1=Yes, 2=No, 7=Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing Recode: 0=Not Yes, 1=Yes

DIABETE3: (Ever told) you have diabetes (If ´Yes´ and respondent is female, ask ´Was this only when you were pregnant?´ Values: 1=Yes, 2=Yes, Gestational Only, 3=No, 4=No, Pre-Diabetes, =Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing Recode: 0=Not Yes, 1=Yes

MENTHLTH: Now thinking about your mental health, which includes stress, depression, and problems with emotions, for how many days during the past 30 days was your mental health not good? Values: 1-30=Days, 88=None, 77=Don’Know/Not Sure, 99=Refused, BLANK=Not Asked/Missing Recode: 0=Not Positive Number between 1-30, 1=Positive Number between 1-30

#Pre-processing Function

myf=function(mydata, myvars){
  mydata=mydata[myvars]

  mydata$X.STATE[mydata$X.STATE %in% c(9, 23, 25, 33, 44, 50)==TRUE]="NEW.ENGLAND"
  mydata$X.STATE[mydata$X.STATE %in% c(17, 18, 26, 39, 55)==TRUE]="EAST.NORTH.CENTRAL"
  mydata$X.STATE[mydata$X.STATE %in% c(1, 21, 28, 47)==TRUE]="EAST.SOUTH.CENTRAL"
  mydata$X.STATE[mydata$X.STATE %in% c(34, 36, 42)==TRUE]="MIDDLE.ATLANTIC"
  mydata$X.STATE[mydata$X.STATE %in% c(8, 30, 32, 35, 49, 56, 4, 16)==TRUE]="MOUNTAIN"
  mydata$X.STATE[mydata$X.STATE %in% c(2, 6, 15, 41, 53)==TRUE]="PACIFIC"
  mydata$X.STATE[mydata$X.STATE %in% c(10, 11, 12, 13, 24, 37, 45, 51, 54)==TRUE]="SOUTH.ATLANTIC"
  mydata$X.STATE[mydata$X.STATE %in% c(19, 20, 27, 29, 31, 38, 46)==TRUE]="WEST.NORTH.CENTRAL"
  mydata$X.STATE[mydata$X.STATE %in% c(5, 22, 40, 48)==TRUE]="WEST.SOUTH.CENTRAL"
  mydata$X.STATE[mydata$X.STATE %in% c(60, 3, 81, 7, 64, 14, 66, 84, 86, 67, 89, 68, 71, 76, 69, 70, 95, 43, 72, 74, 52, 78, 79)==TRUE]="TERRITORIES"
  mydata$X.STATE=as.factor(mydata$X.STATE)
  
  mydata$X.RFBMI5[mydata$X.RFBMI5!=2]=0 #1=0, No = 0
  mydata$X.RFBMI5[mydata$X.RFBMI5==2]=1 #2=1, Yes = 1
  mydata$X.RFBMI5[is.na(mydata$X.RFBMI5)==TRUE]=0 #NA's
  
  mydata$HISPANIC=mydata$X.IMPRACE #Create Hispanic Variable
  mydata$HISPANIC[mydata$HISPANIC!=5]=0
  mydata$HISPANIC[mydata$HISPANIC==5]=1
  mydata$HISPANIC[is.na(mydata$HISPANIC)==TRUE]=0 #NA's
  
  mydata$X.IMPRACE[mydata$X.IMPRACE!=1]=0 #1=0, No = 0
  mydata$X.IMPRACE[is.na(mydata$X.IMPRACE)==TRUE]=0
  
  mydata$SEX[mydata$SEX!=1]=0 #Male=1, Others=0
  mydata$SEX[is.na(mydata$SEX)==TRUE]=0 
  
  mydata$MARITAL[mydata$MARITAL!=1]=0  #2 or higher = 0
  mydata$MARITAL[is.na(mydata$MARITAL)==TRUE]=0 
  
  mydata$INCOME2[mydata$INCOME2!=8]=0 #77 or 99 = 0
  mydata$INCOME2[mydata$INCOME2==8]=1 #GE 75K = 1
  mydata$INCOME2[is.na(mydata$INCOME2)==TRUE]=0 
  
  mydata$EDUCA[mydata$EDUCA!=6]=0 #Not College Grad
  mydata$EDUCA[mydata$EDUCA==6]=1 #College Grad
  mydata$EDUCA[is.na(mydata$EDUCA)==TRUE]=0 #
  
  mydata$EMPLOY1[mydata$EMPLOY1!=1]=0 #Not Employed for Wages
  mydata$EMPLOY1[is.na(mydata$EMPLOY1)==TRUE]=0 
  
  mydata$VETERAN3[mydata$VETERAN3!=1]=0 #Not Veteran
  mydata$VETERAN3[is.na(mydata$VETERAN3)==TRUE]=0 
  
  mydata$CVDCRHD4[mydata$CVDCRHD4!=1]=0 #Heart Disease
  mydata$CVDCRHD4[is.na(mydata$CVDCRHD4)==TRUE]=0 
  
  mydata$CVDSTRK3[mydata$CVDSTRK3!=1]=0 #Stroke
  mydata$CVDSTRK3[is.na(mydata$CVDSTRK3)==TRUE]=0 
  
  mydata$CHCSCNCR[mydata$CHCSCNCR!=1]=0 #Skin Cancer
  mydata$CHCSCNCR[is.na(mydata$CHCSCNCR)==TRUE]=0 
  
  mydata$CHCOCNCR[mydata$CHCOCNCR!=1]=0 #other Cancer
  mydata$CHCOCNCR[is.na(mydata$CHCOCNCR)==TRUE]=0 
  
  mydata$CHCCOPD1[mydata$CHCCOPD1!=1]=0 #COPD
  mydata$CHCCOPD1[is.na(mydata$CHCCOPD1)==TRUE]=0 
  
  mydata$HAVARTH3[mydata$HAVARTH3!=1]=0 #Arthritis
  mydata$HAVARTH3[is.na(mydata$HAVARTH3)==TRUE]=0 
  
  mydata$CHCKIDNY[mydata$CHCKIDNY!=1]=0 #Kidney
  mydata$CHCKIDNY[is.na(mydata$CHCKIDNY)==TRUE]=0 
  
  mydata$DIABETE3[mydata$DIABETE3!=1]=0 #Diabetes
  mydata$DIABETE3[is.na(mydata$DIABETE3)==TRUE]=0  
  
  mydata$MENTHLTH[mydata$MENTHLTH<=30]=1 #Mental Health
  mydata$MENTHLTH[mydata$MENTHLTH>30]=0
  mydata$MENTHLTH[is.na(mydata$MENTHLTH)==TRUE]=0
  mydata$MENTHLTH[is.na(mydata$MENTHLTH)==TRUE]=0  
  
  colnames(mydata)=c( "Age", "Caucasian", "Male" ,"Married", "Income.75K", 
    "College.Graduate" ,"Employed.for.Wages", "Veteran","Overweight.Obese", "Heart.Disease",
    "Stroke", "Skin.Cancer","Cancer" , "COPD" , "Arthritis", 
    "Depression", "Kidney.Disease" , "Diabetes","Stratum" ,  "Weights" , 
    "Region", "Year", "Hispanic")
  
  mydata=mydata[, c(1:2, 23, 3:22)] #Put Hispanic by Race
  for (i in 1:22){mydata[,i]=as.numeric(mydata[,i])}# for checking descriptives
  
  return(mydata)
}

print("Variable Recoding Function Loaded....")
## [1] "Variable Recoding Function Loaded...."

Prep and Print Functions

We use functions for unregistering clusters, citing, and printing.

options(survey.lonely.psu = "adjust")
myprint=function(x){x%>%kbl()%>%kable_classic(html_font = "Cambria")}
mycite=function(x){citation(x)}
unregister=function() {env <- foreach:::.foreachGlobals; rm(list=ls(name=env), pos=env)}

myprint;mycite;unregister
## function(x){x%>%kbl()%>%kable_classic(html_font = "Cambria")}
## function(x){citation(x)}
## function() {env <- foreach:::.foreachGlobals; rm(list=ls(name=env), pos=env)}

Read Data

Read Data with Cluster

Reading with \(foreach\) is significantly faster.

clust=makeCluster(18)
registerDoParallel(clust)
st=proc.time()
mytemp=foreach(i=2003:2019) %dopar% SASxport::read.xport(paste0("Y", i, ".xpt"))
proc.time()-st
##    user  system elapsed 
##   13.73   14.30   65.74
for (i in 1:17){mytemp[[i]]$Year=i+2002}
unregister()

Coding Changes

The base year for variables was 2017. Post and prior years were recoded to match the variables in this year. Changes are noted on each tab.

Year 2019

In 2019, CHCKNDY was CHCKNDY2; SEX was SEXVAR, CHCOPD1 was CHCCOPD2, HAVARTH3 was HAVARTH4, and DIABETE3 was DIABETE4.

CHCKDNY2: Not including kidney stones, bladder infection or incontinence, were you ever told you have kidney disease? 1=Y, 2=N, 7=Don’t Know /Not Sure, 9=Refused, BLANK=Not Asked or Missing CHCKDNY1: (Ever told) you have kidney disease? Do NOT include kidney stones, bladder infection or incontinence. (Incontinence is not being able to control urine flow.) 1=Y, 2=N, 7=Don’t Know /Not Sure, 9=Refused, BLANK=Not Asked or Missing. These response options were identical.

SEXVAR: Sex of Respondent, 1=Male, 2=Female (No Refused Category) SEX: Indicate Sex of Respondent, 1=Male, 2=Female, 9=Refused Sex was modeled as either known to be male or not known to be male (meaning that the 7’s and 9’s were assumed female, the modal response value)

CHCOPD2:(Ever told) (you had) chronic obstructive pulmonary disease, C.O.P.D.,emphysema or chronic bronchitis? CHCOPD1:(Ever told) you have chronic obstructive pulmonary disease, C.O.P.D., emphysema or chronic bronchitis? Identical with same responses. This applied to HAVARTH3 and 4 as well as DIABETE3 and 4.

DIABETE4: (Ever told) (you had) diabetes? (If ´Yes´ and respondent is female, ask ´Was this only when you were pregnant?´. If Respondent says pre-diabetes or borderline diabetes, use response code 4.) DIABETE3: (Ever told) you have diabetes (If ´Yes´ and respondent is female, ask ´Was this only when you were pregnant?´. If Respondent says pre-diabetes or borderline diabetes, use response code 4.)

HAVARTH4:(Ever told) you have some form of arthritis, rheumatoid arthritis, gout, lupus, or fibromyalgia? (Arthritis diagnoses include: rheumatism, polymyalgia rheumatica; osteoarthritis (not osteporosis); tendonitis, bursitis, bunion, tennis elbow; carpal tunnel syndrome, tarsal tunnel syndrome; joint infection, etc.) HAVARTH3: (Ever told) you have some form of arthritis, rheumatoid arthritis, gout, lupus, or fibromyalgia? (Arthritis diagnoses include: rheumatism, polymyalgia rheumatica; osteoarthritis (not osteporosis); tendonitis, bursitis, bunion, tennis elbow; carpal tunnel syndrome, tarsal tunnel syndrome; joint infection, etc.)

Year 2018

In 2018, there were only two changes from the base year: CHCKDNY became CHCKDNY1, SEX became SEX1. CHCKDNY1: (Ever told) you have kidney disease? (Do NOT include kidney stones, bladder infection or incontinence.) CHCKDNY: (Ever told) you have kidney disease? Do NOT include kidney stones, bladder infection or incontinence.(Incontinence is not being able to control urine flow.) The response options were identical (1=Y, 2=N, 7=Don’t Know /Not Sure, 9=Refused, BLANK=Not Asked or Missing). These variables were treated as identical.

SEX1: What is your sex? or What was your sex at birth? Was it…, 1=Male, 2=Female, 7=Don’t Know/Not Sure, 9=Refused SEX: Indicate Sex of Respondent, 1=Male, 2=Female, 9=Refused Sex was modeled as either known to be male or not known to be male (meaning that the 7’s and 9’s were assumed female, the modal response value)

Year 2017

This year was the base year for all variables.

Year 2016

In 2016, the imputed race variables was named .RACE.

.RACE=Race/ethnicity categories Responses: 1. White Non-Hispanic, 2. Black Non-Hispanic, 3. American Indian / Native Alaskan Non-Hispanic, 4. Asian Non-Hispanic, 5. Native Hawaiian / PI, 6. Other Race Non-Hispanic, 7. Multiracial Non-Hispanic, 8. Hispanic, 9. Don’t Know / Not Sure / Refused

.IMPRACE=Imputed race/ethnicity value Responses: 1. White Non-Hispanic, 2. Black Non-Hispanic, 3. Asian Non-Hispanic, 4. American Indian / Native Alaskan Non-Hispanic, 5.Hispanic, 6. Other Race Non-Hispanic

To address this, we zero out 5’s and then assign 8’s to 5’s as we are only interested in White vs. Non-White, and Hispanics vs. Non-Hispanic.

Year 2015

In 2015, the imputed race variable was named .RACE.

Year 2014

In 2014, the imputed race variable was named .RACE.

Year 2013

In 2013, the imputed race variable was named .RACE.

Year 2012

In 2012, the imputed race variable was named .RACE2. The coding was similar but switched from RACE.

.RACE=Race/ethnicity categories Responses: 1. White Non-Hispanic, 2. Black Non-Hispanic, 3. Asian Only Non-Hispanic, 4. Native Hawaiian / PI, 5. American Indian or Alaskan Native Non-Hispanic, 6. Other Race Non-Hispanic, 7. Multiracial Hispanic, 8. Hispanic, 9. Don’t know/Not sure/Refused

Versus…

.IMPRACE=Imputed race/ethnicity value Responses: 1. White Non-Hispanic, 2. Black Non-Hispanic, 3. Asian Non-Hispanic, 4. American Indian / Native Alaskan Non-Hispanic, 5.Hispanic, 6. Other Race Non-Hispanic

Further the employment variable was EMPLOY rather than EMPLOY1.

.EMPLOY Values: 1. Employed for Wages, 2. Self-Employed, 3. Out of work for > 1 year, 4. Out of work < 1 year, 5. A Homemaker, 6. A Student, 7. Retired, 8. Unable to Work, 9. Refused, BLANK Not Asked or Missing

EMPLOY1: Are you currently…? Values: 1=Employed for Wages, 2= Self-Employed, 3=Out of work >1 year, 4=Out of Work< 1 year, 5=A homemaker, 6= A student, 7=Retire, 8=Unable to Work, 9=Refused, BLANK=Not Asked/Missing

Year 2011

In 2011, the imputed race variable was named .RACE2. The employment variable was EMPLOY rather than EMPLOY1 and measured as stated above. The COPD variable was CHCCOPD rather than CHCCOPD1 and measured identically.

Year 2010

In 2010, the race, employment and COPD changes were measured identically. Note: CDC recommends against combining (joint analysis) any years before 2010 due to the addition of cellular and raking methodology in 2011. However, these are the best estimates for the population, so treating the years independently when combined solves this issue.

On top of these changes, SEX was SEXG with 1=Male and 2=Female. The recoding in the base works. VETERAN3 was VETERAN2.

VETERAN2: Have you ever served on active duty in the United States Armed Forces, either in the regular military or in a National Guard or military reserve unit? Active duty does not include training for the Reserves or National Guard, but DOES include activation, for example, for the Persian Gulf War. Values: 1=Yes Now on Active Duty, 2=Yes On Active Duty In the Past But Not the Last 12 Months, 3=Yes, on active duty in the past, but not during the last 12 months, 4=No, training for Reserves or National Guard only, 5=No, never served in the military, 7=Don’t know/Not Sure, 9=Refused, BLANK=Not Asked / Missing

Versus

VETERAN3: Have you ever served on active duty in the United States Armed Forces, either in the regular military or in a National Guard or military reserve unit? Values: 1=Yes, 2=No, 7=Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing

To address this, we recoded 1-3 for VETERAN2 as 1 and all other values as 0.

HAVARTH2 preceeded HAVARTH3. The responses were identical

HAVARTH2:Have you ever been told by a doctor or other health professional that you have some form of arthritis, rheumatoid arthritis, gout, lupus, or fibromyalgia?(Arthritis diagnoses include: rheumatism, polymyalgia rheumatica; osteoarthritis (not osteporosis); tendonitis, bursitis, bunion, tennis elbow; carpal tunnel syndrome, tarsal tunnel syndrome; joint infection, etc.) Values: 1=Yes, 2=No, 7=Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing

HAVARTH3:(Ever told) you have some form of arthritis, rheumatoid arthritis, gout, lupus, or fibromyalgia? (Arthritis diagnoses include: rheumatism, polymyalgia rheumatica; osteoarthritis (not osteporosis); tendonitis, bursitis, bunion, tennis elbow; carpal tunnel syndrome, tarsal tunnel syndrome; joint infection, etc.) Values: 1=Yes, 2=No, 7=Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing

CHCSNCR,Skin Cancer, was derivable using CNCRTYPE, a follow-on to whether an individual was told they had cancer (CNCRHAVE). What type of cancer was it? 20 and 21 were “Melanoma” and “Other Skin Cancer.”

Similarly, Other Cancer (CHCOCNCR) were derivable from CNCRHAVE by eliminating 20 and 21 (skin cancers).

COPD was not available. Kidney Disease was not available.

DIABETE2 preceded DIABETE3. Responses were identical.

DIABETE2: Have you ever been told by a doctor that you have diabetes (If “Yes” and respondent is female, ask “Was this only when you were pregnant?”. Values: 1=Yes, 2=Yes, Gestational Only, 3=No, 4=No, Pre-Diabetes, =Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing Recode: 0=Not Yes, 1=Yes

DIABETE3: (Ever told) you have diabetes (If ´Yes´ and respondent is female, ask ´Was this only when you were pregnant?´ Values: 1=Yes, 2=Yes, Gestational Only, 3=No, 4=No, Pre-Diabetes, =Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing Recode: 0=Not Yes, 1=Yes

2010 did not use raking for weights, so they are less reliable in terms of estimating the population. Thus the weighting variables were different. The stratum weights remained the same, but the individual weight variable was different.

Year 2009

Year 2009 was identical to Year 2010 in terms of variables for this study.

Year 2008

In 2008, no questions existed about arthritis, skin cancer, other cancers, COPD, or kidney disease.
VETERAN was used instead of VETERAN2. The measurement and question was identical.

VETERAN: Have you ever served on active duty in the United States Armed Forces, either in the regular military or in a National Guard or military reserve unit? Active duty does not include training for the Reserves or National Guard, but DOES include activation, for example, for the Persian Gulf War. VETERAN2: Have you ever served on active duty in the United States Armed Forces, either in the regular military or in a National Guard or military reserve unit? Active duty does not include training for the Reserves or National Guard, but DOES include activation, for example, for the Persian Gulf War.

Both responses were the same.

Aside from these changes, there were no other differences from 2008.

Year 2007

In 2007, no questions existed about skin cancer, other cancers, COPD, or kidney disease.
VETERAN wasn’t used. Instead, VETERAN1 was. These responses were identical to VETERAN3. What was old was made new again.

VETERAN1: Have you ever served on active duty in the United States Armed Forces, either in the regular military or in a National Guard or military reserve unit? Active duty does not include training for the Reserves or National Guard, but DOES include activation, for example, for the Persian Gulf War. Values: 1=Yes, 2=No, 7=Don’t Know/Not Sure, 9=Refused, Blank=Not Asked/Missing

Additionally, HAVARTH2 was in this set, but this measure was measured differently from HAVARTH3.

HAVARTH2: During the past 30 days, have you had symptoms of pain, aching, or stiffness in or around a joint? Values: 1=Yes, 2=No, 7=Don’t know/Not sure, 9=Refused HAVARTH3:(Ever told) you have some form of arthritis, rheumatoid arthritis, gout, lupus, or fibromyalgia? (Arthritis diagnoses include: rheumatism, polymyalgia rheumatica; osteoarthritis (not osteporosis); tendonitis, bursitis, bunion, tennis elbow; carpal tunnel syndrome, tarsal tunnel syndrome; joint infection, etc.) Values: 1=Yes, 2=No, 7=Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing

These variables were different and could not be compared alike.

Other variable naming changes existed, but there were no differences in the questions.

Year 2006

Another change..

CVDCRHD3/ CVDINFR3: Has a doctor, nurse, or other health professional ever told you that you had any of the following? For each, tell me “Yes”, “No”, or you’re “Not sure”: (Ever told) you had angina or coronary heart disease?

Values: 1=Yes, 2=No, 7=Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing

versus

CVDCRHD4: (Ever told) you had angina or coronary heart disease? Values: 1=Yes, 2=No, 7=Don’t Know/Not Sure, 9=Refused, BLANK=Not Asked/Missing

Year 2005

2005 looks like 2006 for the variables in question.

Year 2004

In 2004, no questions existed about skin cancer, other cancers, COPD, or kidney disease. A module about heart disease contained insufficient usable responses (256,507 blanks versus 303,822 total). Stroke was still present, but also largely unused (256,513 blanks).

Year 2003

In 2003, there was not a 6-group age variable but rather a 13 group. This was easy to collapse by recoding. No questions existed about skin cancer, other cancers, COPD, or kidney disease. DIABETES was the precursor to DIABETE2 and excluded only the pre-diabetes response. RFBMI3 was the precursor to RFBMI4 and was 1) Not at risk (_BMI3 < 2500) and 2) At Risk (2500 <= _BMI3 < 9999). This was identical to the results from RFBMI4 when recoded

Specific Year Recoding

#2019
mytemp[[17]]$CHCKIDNY=mytemp[[17]]$CHCKDNY2 #identical 
mytemp[[17]]$SEX=mytemp[[17]]$SEXVAR        #SEXVAR:  M/F, SEX:  M/F/R 
mytemp[[17]]$CHCCOPD1=mytemp[[17]]$CHCCOPD2 #identical
mytemp[[17]]$HAVARTH3=mytemp[[17]]$HAVARTH4 #identical
mytemp[[17]]$DIABETE3=mytemp[[17]]$DIABETE4 #identical

#2018
mytemp[[16]]$CHCKIDNY=mytemp[[16]]$CHCKDNY1 #identical
mytemp[[16]]$SEX=mytemp[[16]]$SEX1          #SEX1:  M/F/DK/R

#2017 is base year for coding, year 2017
#2016
mytemp[[14]]$X.IMPRACE=mytemp[[14]]$X.RACE  #different order and items
mytemp[[14]]$X.IMPRACE[mytemp[[14]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[14]]$X.IMPRACE[mytemp[[14]]$X.IMPRACE==8]=5 #assign 8 to all 5s

#2015
mytemp[[13]]$X.IMPRACE=mytemp[[13]]$X.RACE
mytemp[[13]]$X.IMPRACE[mytemp[[13]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[13]]$X.IMPRACE[mytemp[[13]]$X.IMPRACE==8]=5 #assign 8 to all 5s

#2014
mytemp[[12]]$X.IMPRACE=mytemp[[12]]$X.RACE
mytemp[[12]]$X.IMPRACE[mytemp[[12]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[12]]$X.IMPRACE[mytemp[[12]]$X.IMPRACE==8]=5 #assign 8 to all 5s

#2013
mytemp[[11]]$X.IMPRACE=mytemp[[11]]$X.RACE
mytemp[[11]]$X.IMPRACE[mytemp[[11]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[11]]$X.IMPRACE[mytemp[[11]]$X.IMPRACE==8]=5 #assign 8 to all 5s

#2012
mytemp[[10]]$X.IMPRACE=mytemp[[10]]$RACE2
mytemp[[10]]$X.IMPRACE[mytemp[[10]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[10]]$X.IMPRACE[mytemp[[10]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[10]]$EMPLOY1=mytemp[[10]]$EMPLOY

#2011
mytemp[[9]]$X.IMPRACE=mytemp[[9]]$RACE2
mytemp[[9]]$X.IMPRACE[mytemp[[9]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[9]]$X.IMPRACE[mytemp[[9]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[9]]$EMPLOY1=mytemp[[9]]$EMPLOY
mytemp[[9]]$CHCCOPD1=mytemp[[9]]$CHCCOPD

#2010
mytemp[[8]]$X.IMPRACE=mytemp[[8]]$RACE2
mytemp[[8]]$X.IMPRACE[mytemp[[8]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[8]]$X.IMPRACE[mytemp[[8]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[8]]$EMPLOY1=mytemp[[8]]$EMPLOY
mytemp[[8]]$X.RFBMI5=mytemp[[8]]$X.RFBMI4
mytemp[[8]]$SEX=mytemp[[8]]$X.SEXG.
mytemp[[8]]$VETERAN3=mytemp[[8]]$VETERAN2
mytemp[[8]]$VETERAN3[mytemp[[8]]$VETERAN3<=3]=1
mytemp[[8]]$VETERAN3[mytemp[[8]]$VETERAN3>=4]=0
mytemp[[8]]$HAVARTH3=rep('',nrow(mytemp[[8]]))
mytemp[[8]]$CHCSCNCR=rep('',nrow(mytemp[[8]]))
mytemp[[8]]$CHCOCNCR=rep('',nrow(mytemp[[8]]))
mytemp[[8]]$CHCCOPD1=rep('',nrow(mytemp[[8]]))
mytemp[[8]]$CHCCOPD1=rep('',nrow(mytemp[[8]]))
mytemp[[8]]$CHCKIDNY=rep('', nrow(mytemp[[8]]))
mytemp[[8]]$DIABETE3=mytemp[[8]]$DIABETE2
mytemp[[8]]$X.LLCPWT=mytemp[[8]]$X.FINALWT

#2009
mytemp[[7]]$X.IMPRACE=mytemp[[7]]$RACE2
mytemp[[7]]$X.IMPRACE[mytemp[[7]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[7]]$X.IMPRACE[mytemp[[7]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[7]]$EMPLOY1=mytemp[[7]]$EMPLOY
mytemp[[7]]$X.RFBMI5=mytemp[[7]]$X.RFBMI4
mytemp[[7]]$SEX=mytemp[[7]]$X.SEXG.
mytemp[[7]]$VETERAN3=mytemp[[7]]$VETERAN2
mytemp[[7]]$VETERAN3[mytemp[[7]]$VETERAN3<=3]=1
mytemp[[7]]$VETERAN3[mytemp[[7]]$VETERAN3>=4]=0
mytemp[[7]]$HAVARTH3=rep('',nrow(mytemp[[7]]))
mytemp[[7]]$CHCSCNCR=rep('',nrow(mytemp[[7]]))
mytemp[[7]]$CHCOCNCR=rep('',nrow(mytemp[[7]]))
mytemp[[7]]$CHCCOPD1=rep('',nrow(mytemp[[7]]))
mytemp[[7]]$CHCKIDNY=rep('', nrow(mytemp[[7]]))
mytemp[[7]]$DIABETE3=mytemp[[7]]$DIABETE2
mytemp[[7]]$X.LLCPWT=mytemp[[7]]$X.FINALWT

#2008
mytemp[[6]]$X.IMPRACE=mytemp[[6]]$RACE2
mytemp[[6]]$X.IMPRACE[mytemp[[6]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[6]]$X.IMPRACE[mytemp[[6]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[6]]$EMPLOY1=mytemp[[6]]$EMPLOY
mytemp[[6]]$X.RFBMI5=mytemp[[6]]$X.RFBMI4
mytemp[[6]]$SEX=mytemp[[6]]$X.SEXG.
mytemp[[6]]$VETERAN3=mytemp[[6]]$VETERAN1
mytemp[[6]]$HAVARTH3=rep('',nrow(mytemp[[6]]))
mytemp[[6]]$CHCSCNCR=rep('',nrow(mytemp[[6]]))
mytemp[[6]]$CHCOCNCR=rep('',nrow(mytemp[[6]]))
mytemp[[6]]$CHCCOPD1=rep('', nrow(mytemp[[6]]))
mytemp[[6]]$CHCKIDNY=rep('', nrow(mytemp[[6]]))
mytemp[[6]]$DIABETE3=mytemp[[6]]$DIABETE2
mytemp[[6]]$X.LLCPWT=mytemp[[6]]$X.FINALWT
mytemp[[6]]$X.STSTR=mytemp[[6]]$X.STSTR
mytemp[[6]]$X.AGE.G=mytemp[[6]]$X.AGE.G
mytemp[[6]]$X.STATE=mytemp[[6]]$X.STATE

#2007
mytemp[[5]]$X.IMPRACE=mytemp[[5]]$RACE2
mytemp[[5]]$X.IMPRACE[mytemp[[5]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[5]]$X.IMPRACE[mytemp[[5]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[5]]$EMPLOY1=mytemp[[5]]$EMPLOY
mytemp[[5]]$SEX=mytemp[[5]]$X.SEXG.
mytemp[[5]]$VETERAN3=mytemp[[5]]$VETERAN1
mytemp[[5]]$HAVARTH3=rep('',nrow(mytemp[[5]]))
mytemp[[5]]$CHCSCNCR=rep('',nrow(mytemp[[5]]))
mytemp[[5]]$CHCOCNCR=rep('',nrow(mytemp[[5]]))
mytemp[[5]]$CHCCOPD1=rep('', nrow(mytemp[[5]]))
mytemp[[5]]$CHCKIDNY=rep('', nrow(mytemp[[5]]))
mytemp[[5]]$DIABETE3=mytemp[[5]]$DIABETE2
mytemp[[5]]$X.LLCPWT=mytemp[[5]]$X.FINALWT
mytemp[[5]]$X.STSTR=mytemp[[5]]$X.STSTR
mytemp[[5]]$X.AGE.G=mytemp[[5]]$X.AGE.G
mytemp[[5]]$X.STATE=mytemp[[5]]$X.STATE
mytemp[[5]]$SEX=mytemp[[5]]$X.SEXG.
mytemp[[5]]$X.RFBMI5=mytemp[[5]]$X.RFBMI4

#2006
mytemp[[4]]$X.IMPRACE=mytemp[[4]]$RACE2
mytemp[[4]]$X.IMPRACE[mytemp[[4]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[4]]$X.IMPRACE[mytemp[[4]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[4]]$EMPLOY1=mytemp[[4]]$EMPLOY
mytemp[[4]]$X.RFBMI5=mytemp[[4]]$X.RFBMI4
mytemp[[4]]$SEX=mytemp[[4]]$X.SEXG.
mytemp[[4]]$VETERAN3=mytemp[[4]]$VETERAN
mytemp[[4]]$HAVARTH3=rep('',nrow(mytemp[[4]]))
mytemp[[4]]$CHCSCNCR=rep('',nrow(mytemp[[4]]))
mytemp[[4]]$CHCOCNCR=rep('',nrow(mytemp[[4]]))
mytemp[[4]]$CHCCOPD1=rep('', nrow(mytemp[[4]]))
mytemp[[4]]$CHCKIDNY=rep('', nrow(mytemp[[4]]))
mytemp[[4]]$CVDCRHD4=mytemp[[4]]$CVDCRHD3
mytemp[[4]]$DIABETE3=mytemp[[4]]$DIABETE2
mytemp[[4]]$X.LLCPWT=mytemp[[4]]$X.FINALWT
mytemp[[4]]$X.STSTR=mytemp[[4]]$X.STSTR
mytemp[[4]]$X.AGE.G=mytemp[[4]]$X.AGE.G
mytemp[[4]]$X.STATE=mytemp[[4]]$X.STATE
mytemp[[4]]$SEX=mytemp[[4]]$X.SEXG.
mytemp[[4]]$X.RFBMI5=mytemp[[4]]$X.RFBMI4

#2005
mytemp[[3]]$X.IMPRACE=mytemp[[3]]$RACE2
mytemp[[3]]$X.IMPRACE[mytemp[[3]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[3]]$X.IMPRACE[mytemp[[3]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[3]]$EMPLOY1=mytemp[[3]]$EMPLOY
mytemp[[3]]$X.RFBMI5=mytemp[[3]]$X.RFBMI4
mytemp[[3]]$SEX=mytemp[[3]]$X.SEXG.
mytemp[[3]]$VETERAN3=mytemp[[3]]$VETERAN
mytemp[[3]]$HAVARTH3=rep('',nrow(mytemp[[3]]))
mytemp[[3]]$CHCSCNCR=rep('',nrow(mytemp[[3]]))
mytemp[[3]]$CHCOCNCR=rep('',nrow(mytemp[[3]]))
mytemp[[3]]$CHCCOPD1=rep('', nrow(mytemp[[3]]))
mytemp[[3]]$CHCKIDNY=rep('', nrow(mytemp[[3]]))
mytemp[[3]]$CVDCRHD4=mytemp[[3]]$CVDCRHD3
mytemp[[3]]$DIABETE3=mytemp[[3]]$DIABETE2
mytemp[[3]]$X.LLCPWT=mytemp[[3]]$X.FINALWT
mytemp[[3]]$X.STSTR=mytemp[[3]]$X.STSTR
mytemp[[3]]$X.AGE.G=mytemp[[3]]$X.AGE.G
mytemp[[3]]$X.STATE=mytemp[[3]]$X.STATE
mytemp[[3]]$SEX=mytemp[[3]]$X.SEXG.
mytemp[[3]]$X.RFBMI5=mytemp[[3]]$X.RFBMI4

#2004
mytemp[[2]]$X.IMPRACE=mytemp[[2]]$RACE2
mytemp[[2]]$X.IMPRACE[mytemp[[2]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[2]]$X.IMPRACE[mytemp[[2]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[2]]$EMPLOY1=mytemp[[2]]$EMPLOY
mytemp[[2]]$X.RFBMI5=mytemp[[2]]$X.RFBMI4
mytemp[[2]]$VETERAN3=mytemp[[2]]$VETERAN
mytemp[[2]]$HAVARTH3=rep('',nrow(mytemp[[2]]))
mytemp[[2]]$CHCSCNCR=rep('',nrow(mytemp[[2]]))
mytemp[[2]]$CHCOCNCR=rep('',nrow(mytemp[[2]]))
mytemp[[2]]$CHCCOPD1=rep('', nrow(mytemp[[2]]))
mytemp[[2]]$CHCKIDNY=rep('', nrow(mytemp[[2]]))
mytemp[[2]]$CVDCRHD4=rep('', nrow(mytemp[[2]]))
mytemp[[2]]$CVDSTRK3=rep('', nrow(mytemp[[2]]))
mytemp[[2]]$DIABETE3=mytemp[[2]]$DIABETE2
mytemp[[2]]$X.LLCPWT=mytemp[[2]]$X.FINALWT
mytemp[[2]]$X.STSTR=mytemp[[2]]$X.STSTR
mytemp[[2]]$X.AGE.G=mytemp[[2]]$X.AGE.G
mytemp[[2]]$X.STATE=mytemp[[2]]$X.STATE
mytemp[[2]]$SEX=mytemp[[2]]$X.SEXG.
mytemp[[2]]$X.RFBMI5=mytemp[[2]]$X.RFBMI4

#2003
mytemp[[1]]$X.AGE.G=mytemp[[1]]$X.AGEG5YR
mytemp[[1]]$X.AGE.G[mytemp[[1]]$X.AGEG5YR>1 & mytemp[[1]]$X.AGEG5YR<=3]=2
mytemp[[1]]$X.AGE.G[mytemp[[1]]$X.AGEG5YR>3 & mytemp[[1]]$X.AGEG5YR<=5]=3
mytemp[[1]]$X.AGE.G[mytemp[[1]]$X.AGEG5YR>5 & mytemp[[1]]$X.AGEG5YR<=7]=4
mytemp[[1]]$X.AGE.G[mytemp[[1]]$X.AGEG5YR>7 & mytemp[[1]]$X.AGEG5YR<=9]=5
mytemp[[1]]$X.AGE.G[mytemp[[1]]$X.AGEG5YR>9]=6
mytemp[[1]]$X.IMPRACE=mytemp[[1]]$RACE2
mytemp[[1]]$X.IMPRACE[mytemp[[1]]$X.IMPRACE==5]=0 #zero out item 5
mytemp[[1]]$X.IMPRACE[mytemp[[1]]$X.IMPRACE==8]=5 #assign 8 to all 5s
mytemp[[1]]$EMPLOY1=mytemp[[1]]$EMPLOY
mytemp[[1]]$X.RFBMI5=mytemp[[1]]$X.RFBMI4
mytemp[[1]]$VETERAN3=mytemp[[1]]$VETERAN
mytemp[[1]]$HAVARTH3=rep('',nrow(mytemp[[1]]))
mytemp[[1]]$CHCSCNCR=rep('',nrow(mytemp[[1]]))
mytemp[[1]]$CHCOCNCR=rep('',nrow(mytemp[[1]]))
mytemp[[1]]$CHCCOPD1=rep('', nrow(mytemp[[1]]))
mytemp[[1]]$CHCKIDNY=rep('', nrow(mytemp[[1]]))
mytemp[[1]]$CVDCRHD4=rep('', nrow(mytemp[[1]]))
mytemp[[1]]$CVDSTRK3=rep('', nrow(mytemp[[1]]))
mytemp[[1]]$DIABETE3=mytemp[[1]]$DIABETES
mytemp[[1]]$X.LLCPWT=mytemp[[1]]$X.FINALWT
mytemp[[1]]$X.STSTR=mytemp[[1]]$X.STSTR
mytemp[[1]]$X.STATE=mytemp[[1]]$X.STATE
mytemp[[1]]$SEX=mytemp[[1]]$X.SEXG.
mytemp[[1]]$X.RFBMI5=mytemp[[1]]$X.RFBMI3

Reduce

st=proc.time()
mytemp=lapply(mytemp,myf,myvars)
proc.time()-st
##    user  system elapsed 
##    5.08    3.83    9.48

Combine 2011-2019

mynew=mytemp
mynew[[1]]=mynew[[2]]=mynew[[3]]=mynew[[4]]=mynew[[5]]=mynew[[6]]=mynew[[7]]=mynew[[8]]=NULL

mynew=rbind(mynew[[1]], mynew[[2]], mynew[[3]],mynew[[4]],mynew[[5]], mynew[[6]],mynew[[7]],mynew[[8]])

myjoinedsurvey=survey::svydesign(id=~1,strata =~Stratum,weights =~Weights, check.strata=FALSE,data = mynew)

Make Survey

clust=makeCluster(17)
registerDoParallel(clust)
mysurveys=foreach(i=1:17)%do%survey::svydesign(id=~1,strata = ~Stratum,weights = ~Weights,data = mytemp[[i]])
#rm(mytemp)
#rm(mynew)
unregister()

Descriptive Statistics

Age-Adjusted Tables

myf2=function(svy,dv){
  temp=survey::svytable(as.formula(paste0(dv, "~ Veteran+as.factor(Age)")), svy)
  temp1=survey::svytable(as.formula("~Veteran+as.factor(Age)"), svy)
  temp2=survey::svytable(as.formula(paste0(dv, "~as.factor(Age)")), svy)
  myvet=epitools::ageadjust.direct(count=temp[2,],pop=temp1[2,],stdpop=temp2,conf.level=.95)
  mynonvet=epitools::ageadjust.direct(count=temp[1,],pop=temp1[1,], stdpop=temp2,conf.level=.95)
  z=round(c(mynonvet[2],myvet[2]),2)
  names(z)=c(paste(dv, "Non-Vet"), paste(dv, "Vet"))
  return(z)
}

a1=c('Overweight.Obese', 'Heart.Disease','Stroke','Skin.Cancer','Cancer','COPD',
               'Arthritis','Depression','Kidney.Disease','Diabetes')

unregister()
clust=makeCluster(3)
registerDoParallel(clust)

myt1=proc.time()
  p1=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[1])),myf2(myjoinedsurvey,a1[1]))
  p2=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[2])),myf2(myjoinedsurvey,a1[2]))
  p3=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[3])),myf2(myjoinedsurvey,a1[3]))
  p4=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[4])),myf2(myjoinedsurvey,a1[4]))
  p5=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[5])),myf2(myjoinedsurvey,a1[5]))
  p6=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[6])),myf2(myjoinedsurvey,a1[6]))
  p7=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[7])),myf2(myjoinedsurvey,a1[7]))
  p8=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[8])),myf2(myjoinedsurvey,a1[8]))
  p9=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[9])),myf2(myjoinedsurvey,a1[9]))
  p10=rbind(do.call(rbind,parLapply(clust,mysurveys,myf2,a1[10])),myf2(myjoinedsurvey,a1[10]))
proc.time()-myt1
##    user  system elapsed 
##  125.31    7.20  226.30
ptot=cbind(p1,p2,p3,p4,p5,p6,p7,p8,p9,p10)
ptot=as.data.frame(ptot)
rownames(ptot)=c('Y2003','Y2004','Y2005','Y2006','Y2007','Y2008',
                 'Y2009','Y2010','Y2011','Y2012','Y2013','Y2014',
                 'Y2015','Y2016','Y2017','Y2018','Y2019','Y11-19')
myprint(ptot)
Overweight.Obese Non-Vet Overweight.Obese Vet Heart.Disease Non-Vet Heart.Disease Vet Stroke Non-Vet Stroke Vet Skin.Cancer Non-Vet Skin.Cancer Vet Cancer Non-Vet Cancer Vet COPD Non-Vet COPD Vet Arthritis Non-Vet Arthritis Vet Depression Non-Vet Depression Vet Kidney.Disease Non-Vet Kidney.Disease Vet Diabetes Non-Vet Diabetes Vet
Y2003 0.56 0.70 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN 0.37 0.31 NaN NaN 0.11 0.13
Y2004 0.57 0.70 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN 0.37 0.31 NaN NaN 0.12 0.13
Y2005 0.58 0.70 0.08 0.13 0.05 0.06 NaN NaN NaN NaN NaN NaN NaN NaN 0.36 0.30 NaN NaN 0.13 0.14
Y2006 0.58 0.71 0.08 0.14 0.05 0.06 NaN NaN NaN NaN NaN NaN NaN NaN 0.37 0.30 NaN NaN 0.13 0.15
Y2007 0.59 0.72 0.08 0.13 0.05 0.06 NaN NaN NaN NaN NaN NaN NaN NaN 0.36 0.28 NaN NaN 0.14 0.16
Y2008 0.60 0.73 0.08 0.14 0.05 0.06 NaN NaN NaN NaN NaN NaN NaN NaN 0.36 0.30 NaN NaN 0.14 0.15
Y2009 0.60 0.73 0.07 0.13 0.05 0.06 NaN NaN NaN NaN NaN NaN NaN NaN 0.36 0.30 NaN NaN 0.14 0.16
Y2010 0.61 0.73 0.08 0.14 0.05 0.06 NaN NaN NaN NaN NaN NaN NaN NaN 0.36 0.30 NaN NaN 0.14 0.16
Y2011 0.60 0.73 0.08 0.13 0.05 0.06 0.11 0.16 0.11 0.13 0.08 0.09 0.37 0.36 0.37 0.30 0.03 0.03 0.15 0.18
Y2012 0.60 0.72 0.08 0.13 0.05 0.06 0.11 0.16 0.11 0.12 0.08 0.09 0.38 0.38 0.37 0.30 0.04 0.04 0.15 0.18
Y2013 0.60 0.72 0.08 0.13 0.05 0.06 0.11 0.16 0.11 0.13 0.08 0.10 0.38 0.37 0.35 0.29 0.04 0.04 0.16 0.18
Y2014 0.59 0.72 0.08 0.13 0.05 0.06 0.11 0.17 0.11 0.12 0.09 0.10 0.39 0.37 0.35 0.29 0.04 0.04 0.16 0.18
Y2015 0.59 0.71 0.07 0.12 0.05 0.06 0.12 0.17 0.12 0.13 0.08 0.10 0.37 0.37 0.36 0.30 0.04 0.04 0.16 0.19
Y2016 0.59 0.71 0.08 0.12 0.05 0.06 0.11 0.16 0.11 0.13 0.08 0.10 0.38 0.38 0.36 0.30 0.04 0.05 0.17 0.19
Y2017 0.59 0.71 0.07 0.12 0.05 0.06 0.12 0.17 0.12 0.14 0.09 0.11 0.37 0.37 0.38 0.31 0.04 0.05 0.17 0.19
Y2018 0.60 0.72 0.07 0.12 0.05 0.07 0.12 0.17 0.12 0.14 0.09 0.11 0.38 0.39 0.39 0.32 0.05 0.05 0.17 0.20
Y2019 0.60 0.71 0.07 0.11 0.05 0.07 0.12 0.18 0.12 0.14 0.09 0.11 0.37 0.37 0.41 0.35 0.05 0.05 0.17 0.19
Y11-19 0.60 0.72 0.08 0.13 0.05 0.06 0.11 0.16 0.11 0.13 0.09 0.10 0.38 0.37 0.37 0.30 0.04 0.04 0.16 0.19
rm(p1,p2,p3,p4,p5,p6,p7,p8,p9,p10)
unregister()

Vets vs. Non-Vets

clust=makeCluster(3)
registerDoParallel(clust)
st=proc.time()
res=foreach(i=1:17, .combine="rbind")%do%survey::svytable(~Veteran,mysurveys[[i]])
proc.time()-st
##    user  system elapsed 
##    5.69    0.02    5.70
unweighted=foreach(i=1:17, .combine="rbind")%do% table(mysurveys[[i]]$variables$Veteran)
res=cbind(res,unweighted)
colnames(res)=c("Non-Veteran Weighted", "Veteran Weighted", "Non-Veteran Samples", "Veteran Samples")
rownames(res)=seq(2003,2019)
res%>%kable(format.args = list(big.mark=",",scientific = FALSE))
Non-Veteran Weighted Veteran Weighted Non-Veteran Samples Veteran Samples
2003 190,348,049 30,003,072 228,159 36,525
2004 191,637,278 29,746,086 260,982 42,840
2005 194,578,583 29,532,523 305,107 51,005
2006 198,138,945 29,118,914 304,989 50,721
2007 202,498,717 27,673,461 370,990 59,922
2008 205,615,985 27,244,684 358,433 56,076
2009 208,756,506 26,249,349 374,909 57,698
2010 211,037,577 26,048,662 390,643 60,432
2011 212,198,501 25,812,791 441,873 64,594
2012 216,959,427 26,098,283 415,817 59,870
2013 219,968,409 26,056,006 430,268 61,505
2014 220,704,167 27,778,365 402,544 62,120
2015 224,174,518 27,172,620 383,614 57,842
2016 227,144,466 27,006,670 422,384 63,919
2017 229,254,924 26,398,281 392,148 57,868
2018 230,694,063 27,379,324 381,382 56,054
2019 226,740,688 25,689,603 365,038 53,230
rm(unweighted)
rm(res)
unregister()

Models

Run Models

options(digits=3)

clust=makeCluster(18)
registerDoParallel(clust)

mymod=function(svy,dv){
  f=as.formula(paste0(dv,'~as.factor(Age)+Caucasian+Hispanic+Male+Married+Income.75K+
    College.Graduate+Employed.for.Wages+as.factor(Region)+Veteran'))
  survey::svyglm(f,design=svy,family=quasibinomial)
}

mymod1=function(svy,dv){
  f=as.formula(paste0(dv,'~as.factor(Age)+Caucasian+Hispanic+Male+Married+Income.75K+
    College.Graduate+Employed.for.Wages+as.factor(Region)+Veteran+as.factor(Year)'))
  survey::svyglm(f,design=svy,family=quasibinomial)
}

st=proc.time()
temp1=lapply(mysurveys, mymod, a1[1])
temp8=lapply(mysurveys, mymod, a1[8])
temp10=lapply(mysurveys, mymod, a1[10])

mysurveys[[1]]=mysurveys[[2]]=NULL
temp2=lapply(mysurveys, mymod, a1[2])
temp3=lapply(mysurveys, mymod, a1[3])

mysurveys[[1]]=mysurveys[[2]]=mysurveys[[3]]=mysurveys[[4]]=mysurveys[[5]]=mysurveys[[6]]=NULL
temp4=lapply(mysurveys, mymod, a1[4])
temp5=lapply(mysurveys, mymod, a1[5])
temp6=lapply(mysurveys, mymod, a1[6])
temp7=lapply(mysurveys, mymod, a1[7])
temp9=lapply(mysurveys, mymod, a1[9])

proc.time()-st
##    user  system elapsed 
##   605.5    33.3   639.1
st1=mymod1(myjoinedsurvey,a1[1])
st2=mymod1(myjoinedsurvey,a1[2])
st3=mymod1(myjoinedsurvey,a1[3])
st4=mymod1(myjoinedsurvey,a1[4])
st5=mymod1(myjoinedsurvey,a1[5])
st6=mymod1(myjoinedsurvey,a1[6])
st7=mymod1(myjoinedsurvey,a1[7])
st8=mymod1(myjoinedsurvey,a1[8])
st9=mymod1(myjoinedsurvey,a1[9])
st10=mymod1(myjoinedsurvey,a1[10])


unregister()

Odds Ratios

clust=makeCluster(3)
registerDoParallel(clust)

myodds=function(x){cbind(exp(summary(x)$coefficients[,1]), summary(x)$coefficients[,4])}
t1=parLapply(clust,temp1,myodds)
t2=parLapply(clust,temp2,myodds)
t3=parLapply(clust,temp3,myodds)
t4=parLapply(clust,temp4,myodds)
t5=parLapply(clust,temp5,myodds)
t6=parLapply(clust,temp6,myodds)
t7=parLapply(clust,temp7,myodds)
t8=parLapply(clust,temp8,myodds)
t9=parLapply(clust,temp9,myodds)
t10=parLapply(clust,temp10,myodds)

or1=cbind(myodds(st1),myodds(st2),myodds(st3),myodds(st4),myodds(st5),myodds(st6),
          myodds(st7),myodds(st8),myodds(st9),myodds(st10))


unregister()

Results

myprint2=function(x,nm,y){x%>%kbl(caption=nm, col.names=y)%>%kable_classic(html_font="Cambria")}

mynames=c(2003,'p', 2004, 'p', 2005, 'p', 2006, 'p', 2007, 'p', 2008, 'p', 2009,'p', 2010, 'p', 2011, 'p',
          2012, 'p', 2013, 'p', 2014, 'p', 2015, 'p', 2016, 'p', 2017, 'p', 2018,'p', 2019, 'p')

myprint2(do.call(cbind,t1), a1[1], mynames)
Overweight.Obese
2003 p 2004 p 2005 p 2006 p 2007 p 2008 p 2009 p 2010 p 2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 0.450 0.000 0.494 0.000 0.520 0.000 0.536 0.000 0.477 0.000 0.538 0.000 0.515 0.000 0.496 0.000 0.497 0.000 0.464 0.000 0.483 0.00 0.446 0.000 0.470 0.000 0.453 0.000 0.480 0.000 0.480 0.000 0.530 0.000
as.factor(Age)2 2.074 0.000 1.973 0.000 1.840 0.000 1.946 0.000 2.145 0.000 1.998 0.000 2.227 0.000 2.344 0.000 2.016 0.000 2.127 0.000 2.043 0.00 1.995 0.000 2.024 0.000 1.935 0.000 1.969 0.000 1.970 0.000 1.802 0.000
as.factor(Age)3 2.619 0.000 2.502 0.000 2.401 0.000 2.541 0.000 2.630 0.000 2.549 0.000 2.773 0.000 2.952 0.000 2.854 0.000 3.176 0.000 2.920 0.00 2.888 0.000 2.807 0.000 2.680 0.000 2.704 0.000 2.658 0.000 2.488 0.000
as.factor(Age)4 3.081 0.000 2.855 0.000 2.794 0.000 2.870 0.000 3.032 0.000 2.819 0.000 3.187 0.000 3.391 0.000 3.348 0.000 3.578 0.000 3.430 0.00 3.402 0.000 3.310 0.000 3.156 0.000 3.216 0.000 3.074 0.000 2.972 0.000
as.factor(Age)5 3.640 0.000 3.516 0.000 3.344 0.000 3.257 0.000 3.608 0.000 3.329 0.000 3.485 0.000 3.629 0.000 3.626 0.000 3.842 0.000 3.707 0.00 3.648 0.000 3.513 0.000 3.309 0.000 3.393 0.000 3.274 0.000 3.231 0.000
as.factor(Age)6 2.454 0.000 2.490 0.000 2.445 0.000 2.591 0.000 2.687 0.000 2.442 0.000 2.724 0.000 2.961 0.000 3.060 0.000 3.265 0.000 3.197 0.00 3.250 0.000 3.167 0.000 2.962 0.000 3.017 0.000 2.991 0.000 2.832 0.000
Caucasian 0.883 0.000 0.811 0.000 0.851 0.000 0.850 0.000 0.921 0.000 0.879 0.000 0.910 0.000 0.876 0.000 0.883 0.000 0.908 0.000 0.922 0.00 0.998 0.766 0.954 0.000 1.025 0.001 0.996 0.615 0.990 0.240 0.993 0.408
Hispanic 1.030 0.063 0.967 0.020 0.994 0.648 0.892 0.000 1.116 0.000 1.076 0.000 1.137 0.000 1.115 0.000 1.132 0.000 1.120 0.000 1.105 0.00 1.200 0.000 1.205 0.000 1.195 0.000 1.135 0.000 1.152 0.000 1.135 0.000
Male 2.134 0.000 2.123 0.000 2.070 0.000 2.037 0.000 2.082 0.000 2.036 0.000 2.054 0.000 2.080 0.000 2.008 0.000 2.011 0.000 1.938 0.00 1.866 0.000 1.919 0.000 1.888 0.000 1.790 0.000 1.735 0.000 1.611 0.000
Married 1.125 0.000 1.135 0.000 1.146 0.000 1.145 0.000 1.091 0.000 1.093 0.000 1.102 0.000 1.088 0.000 1.088 0.000 1.102 0.000 1.104 0.00 1.081 0.000 1.035 0.000 1.062 0.000 1.075 0.000 1.045 0.000 1.053 0.000
Income.75K 0.940 0.000 0.917 0.000 0.921 0.000 0.935 0.000 0.937 0.000 0.938 0.000 0.907 0.000 0.935 0.000 0.932 0.000 0.941 0.000 0.952 0.00 1.028 0.001 1.120 0.000 1.102 0.000 1.102 0.000 1.150 0.000 1.207 0.000
College.Graduate 0.684 0.000 0.673 0.000 0.703 0.000 0.696 0.000 0.706 0.000 0.723 0.000 0.712 0.000 0.700 0.000 0.696 0.000 0.687 0.000 0.682 0.00 0.695 0.000 0.678 0.000 0.680 0.000 0.685 0.000 0.709 0.000 0.702 0.000
Employed.for.Wages 1.100 0.000 1.154 0.000 1.159 0.000 1.186 0.000 1.188 0.000 1.185 0.000 1.217 0.000 1.173 0.000 1.188 0.000 1.175 0.000 1.190 0.00 1.196 0.000 1.153 0.000 1.161 0.000 1.187 0.000 1.235 0.000 1.264 0.000
as.factor(Region)2 1.041 0.043 1.127 0.000 1.077 0.000 1.074 0.000 1.140 0.000 1.126 0.000 1.062 0.000 1.070 0.000 1.049 0.001 1.044 0.005 1.085 0.00 1.055 0.001 1.080 0.000 1.142 0.000 1.062 0.000 1.090 0.000 1.089 0.000
as.factor(Region)3 0.865 0.000 0.925 0.000 0.904 0.000 0.825 0.000 0.918 0.000 0.892 0.000 0.830 0.000 0.859 0.000 0.824 0.000 0.834 0.000 0.848 0.00 0.828 0.000 0.789 0.000 0.817 0.000 0.808 0.000 0.789 0.000 0.807 0.000
as.factor(Region)4 0.787 0.000 0.828 0.000 0.805 0.000 0.798 0.000 0.836 0.000 0.805 0.000 0.794 0.000 0.817 0.000 0.773 0.000 0.767 0.000 0.786 0.00 0.799 0.000 0.768 0.000 0.780 0.000 0.801 0.000 0.822 0.000 0.792 0.000
as.factor(Region)5 0.794 0.000 0.854 0.000 0.836 0.000 0.785 0.000 0.852 0.000 0.839 0.000 0.771 0.000 0.817 0.000 0.808 0.000 0.811 0.000 0.801 0.00 0.765 0.000 0.750 0.000 0.782 0.000 0.797 0.000 0.811 0.000 0.710 0.000
as.factor(Region)6 0.915 0.000 0.930 0.000 0.893 0.000 0.810 0.000 0.821 0.000 0.858 0.000 0.774 0.000 0.806 0.000 0.784 0.000 0.788 0.000 0.729 0.00 0.699 0.000 0.718 0.000 0.759 0.000 0.718 0.000 0.761 0.000 0.740 0.000
as.factor(Region)7 0.923 0.000 0.944 0.000 0.971 0.016 0.895 0.000 0.945 0.000 0.921 0.000 0.863 0.000 0.940 0.000 0.893 0.000 0.889 0.000 0.897 0.00 0.886 0.000 0.877 0.000 0.888 0.000 0.849 0.000 0.894 0.000 0.824 0.000
as.factor(Region)8 1.044 0.273 0.988 0.737 0.954 0.167 1.038 0.272 0.918 0.006 0.906 0.002 0.865 0.000 0.920 0.007 0.922 0.006 0.915 0.004 0.957 0.15 0.957 0.158 0.945 0.079 0.991 0.765 1.053 0.114 1.153 0.000 1.133 0.000
as.factor(Region)9 1.022 0.263 1.066 0.000 1.046 0.007 1.005 0.758 1.030 0.050 1.057 0.000 1.001 0.925 0.995 0.751 0.944 0.000 0.959 0.004 0.924 0.00 0.965 0.015 0.924 0.000 0.995 0.737 0.964 0.016 1.009 0.542 0.942 0.000
as.factor(Region)10 0.960 0.013 1.012 0.429 0.979 0.135 0.921 0.000 0.947 0.000 0.988 0.352 0.962 0.003 0.969 0.015 0.959 0.000 0.939 0.000 0.954 0.00 0.979 0.096 0.953 0.000 1.022 0.079 1.018 0.167 1.014 0.286 0.957 0.001
Veteran 1.213 0.000 1.178 0.000 1.119 0.000 1.160 0.000 1.178 0.000 1.223 0.000 1.186 0.000 1.177 0.000 1.202 0.000 1.184 0.000 1.208 0.00 1.246 0.000 1.226 0.000 1.243 0.000 1.282 0.000 1.293 0.000 1.284 0.000
myprint2(do.call(cbind,t2), a1[2], mynames[5:34])
Heart.Disease
2005 p 2006 p 2007 p 2008 p 2009 p 2010 p 2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 0.006 0.000 0.005 0.000 0.004 0.000 0.004 0.000 0.003 0.000 0.002 0.000 0.003 0.000 0.003 0.000 0.003 0.000 0.005 0.000 0.003 0.000 0.003 0.000 0.004 0.000 0.004 0.000 0.004 0.000
as.factor(Age)2 1.671 0.000 1.803 0.000 2.860 0.000 2.426 0.000 2.059 0.000 2.835 0.000 2.181 0.000 2.407 0.000 2.540 0.000 1.277 0.001 2.271 0.000 2.165 0.000 1.706 0.000 2.454 0.000 1.640 0.000
as.factor(Age)3 2.880 0.000 3.611 0.000 4.810 0.000 4.977 0.000 4.698 0.000 6.527 0.000 4.992 0.000 5.636 0.000 4.555 0.000 2.917 0.000 5.360 0.000 5.845 0.000 4.257 0.000 4.712 0.000 3.911 0.000
as.factor(Age)4 6.960 0.000 8.001 0.000 11.198 0.000 11.614 0.000 11.619 0.000 18.534 0.000 13.274 0.000 13.474 0.000 13.538 0.000 8.653 0.000 13.928 0.000 12.319 0.000 9.750 0.000 9.327 0.000 10.542 0.000
as.factor(Age)5 13.714 0.000 16.916 0.000 21.206 0.000 23.802 0.000 23.740 0.000 40.789 0.000 26.200 0.000 27.464 0.000 23.469 0.000 14.990 0.000 26.501 0.000 23.806 0.000 18.725 0.000 17.931 0.000 18.754 0.000
as.factor(Age)6 19.039 0.000 25.341 0.000 30.807 0.000 36.391 0.000 35.283 0.000 63.368 0.000 38.093 0.000 38.093 0.000 37.015 0.000 23.337 0.000 39.596 0.000 35.834 0.000 27.865 0.000 26.005 0.000 28.347 0.000
Caucasian 1.041 0.100 1.129 0.000 1.196 0.000 1.126 0.000 1.146 0.000 1.028 0.202 1.066 0.002 1.093 0.000 1.084 0.000 1.219 0.000 1.092 0.000 1.184 0.000 1.203 0.000 1.094 0.000 1.068 0.004
Hispanic 1.279 0.000 1.345 0.000 1.325 0.000 1.239 0.000 1.146 0.000 0.923 0.016 0.907 0.002 0.873 0.000 0.970 0.323 0.961 0.204 0.874 0.000 0.980 0.511 0.833 0.000 0.859 0.000 0.773 0.000
Male 1.562 0.000 1.601 0.000 1.650 0.000 1.655 0.000 1.630 0.000 1.739 0.000 1.616 0.000 1.776 0.000 1.673 0.000 1.641 0.000 1.662 0.000 1.646 0.000 1.703 0.000 1.697 0.000 1.701 0.000
Married 0.940 0.001 0.939 0.000 0.876 0.000 0.921 0.000 0.894 0.000 0.849 0.000 0.892 0.000 0.858 0.000 0.865 0.000 0.859 0.000 0.862 0.000 0.870 0.000 0.896 0.000 0.917 0.000 0.854 0.000
Income.75K 0.755 0.000 0.633 0.000 0.688 0.000 0.729 0.000 0.746 0.000 0.729 0.000 0.731 0.000 0.777 0.000 0.720 0.000 0.736 0.000 0.755 0.000 0.781 0.000 0.741 0.000 0.745 0.000 0.772 0.000
College.Graduate 0.750 0.000 0.778 0.000 0.772 0.000 0.818 0.000 0.754 0.000 0.796 0.000 0.753 0.000 0.699 0.000 0.768 0.000 0.778 0.000 0.756 0.000 0.748 0.000 0.712 0.000 0.708 0.000 0.784 0.000
Employed.for.Wages 0.572 0.000 0.564 0.000 0.465 0.000 0.527 0.000 0.516 0.000 0.513 0.000 0.485 0.000 0.459 0.000 0.474 0.000 0.490 0.000 0.452 0.000 0.462 0.000 0.473 0.000 0.483 0.000 0.472 0.000
as.factor(Region)2 1.149 0.000 1.039 0.316 1.044 0.222 1.063 0.083 1.067 0.074 1.078 0.030 1.160 0.000 1.246 0.000 1.218 0.000 1.215 0.000 1.158 0.000 1.268 0.000 1.169 0.000 1.295 0.000 1.202 0.000
as.factor(Region)3 1.065 0.041 0.963 0.215 0.995 0.845 0.938 0.026 0.998 0.945 1.000 0.997 1.035 0.205 0.976 0.371 1.001 0.957 1.067 0.021 0.975 0.402 0.959 0.135 0.990 0.726 0.995 0.861 0.964 0.276
as.factor(Region)4 0.900 0.010 0.785 0.000 0.832 0.000 0.797 0.000 0.808 0.000 0.844 0.000 0.816 0.000 0.780 0.000 0.748 0.000 0.856 0.000 0.777 0.000 0.832 0.000 0.769 0.000 0.835 0.000 0.844 0.000
as.factor(Region)5 1.006 0.900 0.846 0.000 0.910 0.020 0.825 0.000 0.940 0.143 0.960 0.317 0.894 0.005 0.861 0.000 0.841 0.000 0.971 0.463 0.932 0.094 0.872 0.001 0.867 0.001 0.891 0.008 0.891 0.010
as.factor(Region)6 0.964 0.246 0.796 0.000 0.715 0.000 0.735 0.000 0.822 0.000 0.878 0.000 0.905 0.000 0.807 0.000 0.784 0.000 0.943 0.039 0.853 0.000 0.831 0.000 0.759 0.000 0.876 0.000 0.781 0.000
as.factor(Region)7 1.124 0.000 0.987 0.635 0.968 0.203 0.932 0.007 1.022 0.421 1.015 0.557 1.045 0.075 1.042 0.095 0.969 0.210 1.076 0.004 1.001 0.982 0.966 0.166 0.967 0.196 1.118 0.000 1.007 0.789
as.factor(Region)8 1.654 0.000 1.235 0.002 1.221 0.002 1.411 0.000 1.669 0.000 2.138 0.000 1.934 0.000 1.931 0.000 1.863 0.000 2.308 0.000 1.828 0.000 1.802 0.000 1.984 0.000 2.193 0.000 2.138 0.000
as.factor(Region)9 0.967 0.392 0.881 0.001 0.863 0.000 0.886 0.001 0.890 0.002 0.930 0.042 0.910 0.006 0.881 0.000 0.886 0.001 0.934 0.057 0.930 0.050 0.874 0.000 0.904 0.006 0.986 0.712 0.913 0.020
as.factor(Region)10 1.145 0.000 1.005 0.878 1.065 0.035 1.030 0.335 1.059 0.069 1.199 0.000 1.122 0.000 1.014 0.630 0.942 0.043 1.129 0.000 1.038 0.232 1.048 0.104 0.990 0.750 1.139 0.000 1.029 0.381
Veteran 1.397 0.000 1.326 0.000 1.333 0.000 1.431 0.000 1.448 0.000 1.316 0.000 1.363 0.000 1.279 0.000 1.361 0.000 1.377 0.000 1.286 0.000 1.329 0.000 1.341 0.000 1.377 0.000 1.340 0.000
myprint2(do.call(cbind,t3), a1[3], mynames[5:34])
Stroke
2005 p 2006 p 2007 p 2008 p 2009 p 2010 p 2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 0.012 0.000 0.006 0.000 0.009 0.000 0.007 0.000 0.009 0.000 0.006 0.000 0.006 0.000 0.006 0.000 0.007 0.000 0.010 0.000 0.006 0.000 0.007 0.000 0.008 0.000 0.009 0.000 0.006 0.000
as.factor(Age)2 1.443 0.000 2.612 0.000 2.397 0.000 2.710 0.000 1.188 0.062 2.550 0.000 3.130 0.000 3.334 0.000 2.763 0.000 2.320 0.000 3.129 0.000 3.203 0.000 2.338 0.000 2.320 0.000 3.470 0.000
as.factor(Age)3 2.920 0.000 4.887 0.000 3.675 0.000 5.081 0.000 3.288 0.000 5.775 0.000 6.007 0.000 6.670 0.000 5.816 0.000 4.401 0.000 7.182 0.000 6.510 0.000 4.975 0.000 5.339 0.000 6.805 0.000
as.factor(Age)4 4.980 0.000 9.425 0.000 7.980 0.000 10.077 0.000 7.281 0.000 11.346 0.000 12.124 0.000 11.493 0.000 10.975 0.000 9.129 0.000 13.744 0.000 11.522 0.000 9.089 0.000 10.451 0.000 14.730 0.000
as.factor(Age)5 8.903 0.000 15.257 0.000 12.587 0.000 15.708 0.000 11.201 0.000 15.609 0.000 17.824 0.000 17.363 0.000 15.838 0.000 12.071 0.000 20.080 0.000 16.010 0.000 13.107 0.000 14.487 0.000 21.132 0.000
as.factor(Age)6 13.661 0.000 21.576 0.000 18.340 0.000 23.838 0.000 16.729 0.000 24.451 0.000 24.326 0.000 22.580 0.000 20.751 0.000 15.701 0.000 24.455 0.000 19.177 0.000 15.655 0.000 17.441 0.000 24.709 0.000
Caucasian 0.680 0.000 0.770 0.000 0.737 0.000 0.708 0.000 0.748 0.000 0.731 0.000 0.732 0.000 0.721 0.000 0.802 0.000 0.717 0.000 0.788 0.000 0.842 0.000 0.818 0.000 0.726 0.000 0.731 0.000
Hispanic 0.730 0.000 0.778 0.000 0.810 0.000 0.733 0.000 0.766 0.000 0.673 0.000 0.720 0.000 0.623 0.000 0.675 0.000 0.622 0.000 0.627 0.000 0.666 0.000 0.603 0.000 0.703 0.000 0.590 0.000
Male 1.096 0.001 1.203 0.000 1.141 0.000 1.049 0.046 0.999 0.975 1.062 0.007 0.999 0.950 1.125 0.000 1.076 0.000 1.081 0.000 1.087 0.000 1.110 0.000 1.122 0.000 1.130 0.000 1.174 0.000
Married 0.722 0.000 0.700 0.000 0.696 0.000 0.677 0.000 0.698 0.000 0.742 0.000 0.726 0.000 0.705 0.000 0.712 0.000 0.694 0.000 0.698 0.000 0.713 0.000 0.725 0.000 0.772 0.000 0.687 0.000
Income.75K 0.662 0.000 0.470 0.000 0.495 0.000 0.553 0.000 0.565 0.000 0.578 0.000 0.623 0.000 0.667 0.000 0.531 0.000 0.545 0.000 0.549 0.000 0.604 0.000 0.579 0.000 0.550 0.000 0.564 0.000
College.Graduate 0.654 0.000 0.686 0.000 0.696 0.000 0.666 0.000 0.727 0.000 0.644 0.000 0.651 0.000 0.616 0.000 0.664 0.000 0.695 0.000 0.642 0.000 0.604 0.000 0.618 0.000 0.640 0.000 0.698 0.000
Employed.for.Wages 0.435 0.000 0.356 0.000 0.358 0.000 0.363 0.000 0.365 0.000 0.351 0.000 0.352 0.000 0.325 0.000 0.372 0.000 0.343 0.000 0.336 0.000 0.356 0.000 0.373 0.000 0.349 0.000 0.370 0.000
as.factor(Region)2 1.210 0.000 1.345 0.000 1.155 0.001 1.219 0.000 1.192 0.000 1.239 0.000 1.207 0.000 1.294 0.000 1.239 0.000 1.281 0.000 1.280 0.000 1.304 0.000 1.316 0.000 1.333 0.000 1.323 0.000
as.factor(Region)3 0.863 0.000 0.872 0.001 0.869 0.000 0.868 0.000 0.864 0.000 0.934 0.062 0.775 0.000 0.870 0.000 0.861 0.000 0.866 0.000 0.893 0.001 0.887 0.000 0.917 0.013 0.884 0.000 0.961 0.287
as.factor(Region)4 0.842 0.002 1.001 0.984 0.904 0.033 0.887 0.013 0.909 0.049 1.003 0.940 0.901 0.011 0.914 0.040 0.823 0.000 0.889 0.006 0.879 0.003 0.925 0.047 0.867 0.001 0.883 0.003 0.927 0.073
as.factor(Region)5 0.862 0.015 0.836 0.004 0.786 0.000 0.817 0.000 0.792 0.000 0.823 0.000 0.797 0.000 0.793 0.000 0.812 0.000 0.878 0.008 0.850 0.002 0.827 0.000 0.950 0.304 0.807 0.000 0.801 0.000
as.factor(Region)6 0.964 0.372 1.064 0.121 0.873 0.000 0.859 0.000 0.971 0.438 0.940 0.086 0.791 0.000 0.873 0.000 0.794 0.000 0.931 0.030 0.860 0.000 0.840 0.000 0.833 0.000 0.841 0.000 0.878 0.000
as.factor(Region)7 1.023 0.539 1.109 0.005 1.020 0.542 1.019 0.580 0.988 0.718 1.096 0.004 1.054 0.064 1.004 0.885 1.029 0.335 1.026 0.381 1.026 0.399 0.995 0.853 1.056 0.075 1.030 0.325 1.062 0.046
as.factor(Region)8 0.616 0.000 0.572 0.000 0.520 0.000 0.677 0.000 0.627 0.000 0.727 0.002 0.532 0.000 0.712 0.001 0.590 0.000 0.693 0.000 0.556 0.000 0.574 0.000 0.817 0.044 0.378 0.000 0.661 0.000
as.factor(Region)9 0.999 0.989 1.067 0.192 1.017 0.715 1.019 0.675 1.033 0.487 1.051 0.263 0.930 0.076 1.022 0.611 0.887 0.004 0.984 0.702 1.007 0.880 1.020 0.604 1.018 0.670 0.940 0.146 0.975 0.560
as.factor(Region)10 1.168 0.000 1.207 0.000 1.054 0.164 1.077 0.056 1.079 0.053 1.210 0.000 1.019 0.580 1.101 0.007 0.997 0.926 1.120 0.001 1.122 0.001 0.971 0.379 1.307 0.000 1.203 0.000 1.256 0.000
Veteran 1.202 0.000 1.208 0.000 1.300 0.000 1.268 0.000 1.393 0.000 1.239 0.000 1.334 0.000 1.268 0.000 1.287 0.000 1.260 0.000 1.243 0.000 1.271 0.000 1.265 0.000 1.295 0.000 1.218 0.000
myprint2(do.call(cbind,t4), a1[4], mynames[17:34])
Skin.Cancer
2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 0.001 0.000 0.001 0.000 0.001 0.000 0.001 0.000 0.001 0.000 0.001 0.000 0.001 0.000 0.001 0.000 0.001 0.000
as.factor(Age)2 1.219 0.005 1.510 0.000 1.756 0.000 1.620 0.000 1.429 0.000 1.924 0.000 1.479 0.000 1.375 0.000 1.160 0.085
as.factor(Age)3 2.824 0.000 3.229 0.000 4.689 0.000 4.256 0.000 3.508 0.000 4.887 0.000 2.745 0.000 3.072 0.000 3.360 0.000
as.factor(Age)4 6.632 0.000 7.237 0.000 10.899 0.000 10.238 0.000 8.926 0.000 11.166 0.000 8.254 0.000 7.677 0.000 7.510 0.000
as.factor(Age)5 12.262 0.000 13.913 0.000 20.179 0.000 18.335 0.000 16.156 0.000 20.505 0.000 14.720 0.000 14.452 0.000 14.583 0.000
as.factor(Age)6 26.185 0.000 28.817 0.000 43.563 0.000 38.750 0.000 36.712 0.000 41.614 0.000 30.895 0.000 29.971 0.000 28.971 0.000
Caucasian 5.411 0.000 5.482 0.000 6.087 0.000 5.930 0.000 6.351 0.000 5.669 0.000 7.988 0.000 6.727 0.000 8.050 0.000
Hispanic 1.312 0.000 1.205 0.000 1.017 0.740 1.185 0.001 0.995 0.932 1.072 0.143 1.316 0.000 1.310 0.000 1.708 0.000
Male 1.085 0.000 1.077 0.000 1.048 0.002 1.063 0.000 1.059 0.000 1.055 0.000 1.023 0.154 1.054 0.001 1.056 0.001
Married 1.060 0.000 1.110 0.000 1.051 0.001 1.087 0.000 1.105 0.000 1.083 0.000 1.106 0.000 1.088 0.000 1.095 0.000
Income.75K 1.262 0.000 1.185 0.000 1.164 0.000 1.159 0.000 1.195 0.000 1.132 0.000 1.172 0.000 1.182 0.000 1.154 0.000
College.Graduate 1.219 0.000 1.246 0.000 1.268 0.000 1.246 0.000 1.234 0.000 1.295 0.000 1.237 0.000 1.248 0.000 1.285 0.000
Employed.for.Wages 0.737 0.000 0.796 0.000 0.780 0.000 0.805 0.000 0.807 0.000 0.749 0.000 0.801 0.000 0.768 0.000 0.750 0.000
as.factor(Region)2 1.506 0.000 1.612 0.000 1.484 0.000 1.440 0.000 1.494 0.000 1.473 0.000 1.399 0.000 1.452 0.000 1.353 0.000
as.factor(Region)3 1.026 0.333 0.914 0.001 0.980 0.455 0.977 0.394 1.007 0.812 0.974 0.337 0.979 0.456 0.960 0.165 0.917 0.005
as.factor(Region)4 1.439 0.000 1.460 0.000 1.552 0.000 1.486 0.000 1.474 0.000 1.419 0.000 1.508 0.000 1.465 0.000 1.433 0.000
as.factor(Region)5 1.039 0.282 1.042 0.254 1.086 0.017 1.094 0.012 1.060 0.104 1.038 0.286 1.096 0.013 1.042 0.283 1.044 0.243
as.factor(Region)6 1.403 0.000 1.459 0.000 1.429 0.000 1.282 0.000 1.370 0.000 1.330 0.000 1.587 0.000 1.485 0.000 1.441 0.000
as.factor(Region)7 1.571 0.000 1.591 0.000 1.616 0.000 1.514 0.000 1.554 0.000 1.559 0.000 1.561 0.000 1.535 0.000 1.568 0.000
as.factor(Region)8 0.568 0.000 0.814 0.138 1.153 0.269 0.992 0.953 0.865 0.345 0.876 0.328 1.132 0.382 0.817 0.181 0.778 0.083
as.factor(Region)9 1.065 0.043 1.092 0.006 1.159 0.000 0.991 0.770 1.068 0.038 1.063 0.048 1.084 0.015 1.024 0.487 1.050 0.133
as.factor(Region)10 1.448 0.000 1.409 0.000 1.302 0.000 1.235 0.000 1.170 0.000 1.321 0.000 1.369 0.000 1.228 0.000 1.182 0.000
Veteran 1.314 0.000 1.288 0.000 1.295 0.000 1.360 0.000 1.289 0.000 1.276 0.000 1.334 0.000 1.266 0.000 1.400 0.000
myprint2(do.call(cbind,t5), a1[5], mynames[17:34])
Cancer
2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 0.013 0.000 0.010 0.000 0.009 0.000 0.008 0.000 0.009 0.000 0.009 0.000 0.008 0.000 0.008 0.000 0.009 0.000
as.factor(Age)2 2.130 0.000 3.043 0.000 3.009 0.000 2.720 0.000 2.393 0.000 2.771 0.000 2.701 0.000 2.797 0.000 2.576 0.000
as.factor(Age)3 3.290 0.000 3.929 0.000 4.477 0.000 4.996 0.000 3.792 0.000 5.094 0.000 4.775 0.000 5.035 0.000 4.070 0.000
as.factor(Age)4 5.507 0.000 7.282 0.000 7.739 0.000 8.156 0.000 6.935 0.000 9.012 0.000 9.141 0.000 9.574 0.000 7.584 0.000
as.factor(Age)5 9.267 0.000 12.575 0.000 12.435 0.000 13.226 0.000 11.627 0.000 13.821 0.000 14.876 0.000 14.878 0.000 12.578 0.000
as.factor(Age)6 15.729 0.000 20.339 0.000 22.566 0.000 22.688 0.000 20.144 0.000 24.052 0.000 25.651 0.000 25.875 0.000 21.441 0.000
Caucasian 1.322 0.000 1.320 0.000 1.384 0.000 1.337 0.000 1.386 0.000 1.290 0.000 1.406 0.000 1.301 0.000 1.410 0.000
Hispanic 0.760 0.000 0.854 0.000 0.763 0.000 0.795 0.000 0.719 0.000 0.705 0.000 0.960 0.126 0.935 0.010 0.918 0.001
Male 0.592 0.000 0.614 0.000 0.594 0.000 0.618 0.000 0.618 0.000 0.618 0.000 0.628 0.000 0.650 0.000 0.663 0.000
Married 0.960 0.001 0.973 0.036 0.939 0.000 0.957 0.001 0.930 0.000 0.939 0.000 0.966 0.008 0.962 0.003 0.951 0.000
Income.75K 0.957 0.008 0.887 0.000 0.968 0.049 0.951 0.002 0.945 0.001 0.977 0.142 0.994 0.699 0.941 0.000 0.972 0.081
College.Graduate 0.986 0.331 0.990 0.490 0.978 0.134 1.019 0.215 1.025 0.105 1.048 0.001 1.007 0.620 1.062 0.000 0.974 0.083
Employed.for.Wages 0.632 0.000 0.682 0.000 0.679 0.000 0.628 0.000 0.673 0.000 0.645 0.000 0.609 0.000 0.662 0.000 0.671 0.000
as.factor(Region)2 0.936 0.017 0.990 0.712 1.084 0.003 1.165 0.000 1.061 0.044 0.926 0.007 1.107 0.000 1.041 0.162 1.025 0.400
as.factor(Region)3 1.006 0.795 0.950 0.023 1.018 0.427 0.995 0.841 1.044 0.072 0.913 0.000 0.979 0.366 0.989 0.629 0.996 0.887
as.factor(Region)4 1.007 0.789 0.963 0.183 0.989 0.683 1.072 0.013 1.070 0.017 0.931 0.009 0.965 0.203 0.971 0.291 0.948 0.060
as.factor(Region)5 0.986 0.649 0.998 0.937 1.063 0.045 1.166 0.000 1.075 0.024 0.996 0.889 1.077 0.017 0.969 0.331 1.008 0.794
as.factor(Region)6 0.951 0.022 1.052 0.020 1.110 0.000 1.122 0.000 1.119 0.000 0.956 0.041 0.993 0.742 0.937 0.005 0.933 0.003
as.factor(Region)7 1.021 0.300 0.961 0.054 1.003 0.887 1.031 0.143 1.072 0.001 0.928 0.000 0.981 0.344 0.980 0.324 1.016 0.450
as.factor(Region)8 0.757 0.000 0.787 0.001 0.901 0.154 0.954 0.513 1.086 0.267 1.064 0.356 0.765 0.000 0.933 0.314 0.901 0.152
as.factor(Region)9 0.994 0.823 1.002 0.936 1.038 0.173 1.035 0.229 1.049 0.096 0.937 0.019 1.010 0.715 0.937 0.024 0.961 0.170
as.factor(Region)10 1.062 0.009 0.944 0.019 1.015 0.530 1.012 0.638 1.056 0.030 0.977 0.333 0.929 0.003 0.870 0.000 1.002 0.924
Veteran 1.633 0.000 1.551 0.000 1.549 0.000 1.453 0.000 1.483 0.000 1.448 0.000 1.525 0.000 1.475 0.000 1.468 0.000
myprint2(do.call(cbind,t6), a1[6], mynames[17:34])
COPD
2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 0.045 0.000 0.043 0.000 0.036 0.000 0.034 0.000 0.034 0.000 0.039 0.000 0.032 0.000 0.036 0.000 0.031 0.000
as.factor(Age)2 1.597 0.000 1.740 0.000 2.060 0.000 1.982 0.000 1.642 0.000 1.995 0.000 1.987 0.000 1.757 0.000 2.083 0.000
as.factor(Age)3 2.486 0.000 2.639 0.000 2.834 0.000 3.045 0.000 2.558 0.000 3.094 0.000 3.408 0.000 3.008 0.000 3.302 0.000
as.factor(Age)4 4.020 0.000 4.456 0.000 5.020 0.000 5.400 0.000 4.521 0.000 4.893 0.000 5.630 0.000 4.983 0.000 5.382 0.000
as.factor(Age)5 4.841 0.000 5.500 0.000 6.327 0.000 6.655 0.000 6.017 0.000 6.537 0.000 8.085 0.000 6.894 0.000 7.797 0.000
as.factor(Age)6 4.292 0.000 4.801 0.000 5.407 0.000 5.917 0.000 5.340 0.000 5.697 0.000 6.744 0.000 6.077 0.000 7.124 0.000
Caucasian 1.256 0.000 1.195 0.000 1.245 0.000 1.320 0.000 1.401 0.000 1.253 0.000 1.372 0.000 1.295 0.000 1.325 0.000
Hispanic 0.664 0.000 0.643 0.000 0.662 0.000 0.665 0.000 0.597 0.000 0.638 0.000 0.606 0.000 0.745 0.000 0.648 0.000
Male 0.771 0.000 0.807 0.000 0.820 0.000 0.821 0.000 0.842 0.000 0.827 0.000 0.832 0.000 0.803 0.000 0.821 0.000
Married 0.652 0.000 0.642 0.000 0.671 0.000 0.626 0.000 0.656 0.000 0.619 0.000 0.624 0.000 0.660 0.000 0.642 0.000
Income.75K 0.522 0.000 0.511 0.000 0.485 0.000 0.531 0.000 0.508 0.000 0.547 0.000 0.505 0.000 0.511 0.000 0.537 0.000
College.Graduate 0.491 0.000 0.477 0.000 0.469 0.000 0.479 0.000 0.457 0.000 0.458 0.000 0.463 0.000 0.488 0.000 0.469 0.000
Employed.for.Wages 0.467 0.000 0.449 0.000 0.478 0.000 0.459 0.000 0.453 0.000 0.462 0.000 0.446 0.000 0.468 0.000 0.503 0.000
as.factor(Region)2 1.239 0.000 1.269 0.000 1.292 0.000 1.355 0.000 1.400 0.000 1.241 0.000 1.258 0.000 1.332 0.000 1.241 0.000
as.factor(Region)3 0.879 0.000 0.929 0.002 0.922 0.000 0.879 0.000 0.956 0.075 0.855 0.000 0.832 0.000 0.856 0.000 0.879 0.000
as.factor(Region)4 0.820 0.000 0.880 0.000 0.882 0.000 0.832 0.000 0.878 0.000 0.812 0.000 0.801 0.000 0.824 0.000 0.857 0.000
as.factor(Region)5 0.905 0.002 0.879 0.000 0.919 0.010 0.889 0.000 0.908 0.006 0.813 0.000 0.865 0.000 0.803 0.000 0.792 0.000
as.factor(Region)6 0.701 0.000 0.804 0.000 0.807 0.000 0.796 0.000 0.794 0.000 0.740 0.000 0.756 0.000 0.738 0.000 0.751 0.000
as.factor(Region)7 1.032 0.115 1.016 0.439 1.017 0.409 1.046 0.029 1.053 0.018 0.950 0.010 1.032 0.131 1.033 0.129 1.014 0.505
as.factor(Region)8 0.603 0.000 0.669 0.000 0.688 0.000 0.709 0.000 0.882 0.116 0.830 0.006 0.877 0.083 1.009 0.900 0.891 0.134
as.factor(Region)9 0.838 0.000 0.872 0.000 0.841 0.000 0.853 0.000 0.903 0.001 0.799 0.000 0.836 0.000 0.854 0.000 0.857 0.000
as.factor(Region)10 0.979 0.372 0.956 0.068 0.972 0.243 0.953 0.052 1.044 0.094 0.863 0.000 0.907 0.000 1.032 0.202 0.958 0.092
Veteran 1.332 0.000 1.360 0.000 1.375 0.000 1.340 0.000 1.340 0.000 1.379 0.000 1.346 0.000 1.390 0.000 1.418 0.000
myprint2(do.call(cbind,t7), a1[7], mynames[17:34])
Arthritis
2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 0.058 0.000 0.061 0.000 0.059 0.000 0.054 0.000 0.051 0.000 0.054 0.000 0.056 0.000 0.058 0.000 0.058 0.000
as.factor(Age)2 2.799 0.000 2.806 0.000 2.799 0.000 3.300 0.000 2.826 0.000 2.785 0.000 2.962 0.000 2.723 0.000 2.729 0.000
as.factor(Age)3 6.060 0.000 6.110 0.000 6.145 0.000 7.285 0.000 6.265 0.000 6.614 0.000 6.335 0.000 6.155 0.000 5.759 0.000
as.factor(Age)4 13.090 0.000 12.759 0.000 12.951 0.000 14.828 0.000 13.411 0.000 13.460 0.000 13.405 0.000 12.542 0.000 11.402 0.000
as.factor(Age)5 22.209 0.000 22.393 0.000 21.656 0.000 24.639 0.000 21.903 0.000 22.906 0.000 22.025 0.000 21.196 0.000 19.098 0.000
as.factor(Age)6 27.638 0.000 27.175 0.000 27.651 0.000 31.984 0.000 27.740 0.000 28.508 0.000 27.297 0.000 26.642 0.000 24.318 0.000
Caucasian 1.168 0.000 1.168 0.000 1.161 0.000 1.169 0.000 1.262 0.000 1.268 0.000 1.219 0.000 1.220 0.000 1.259 0.000
Hispanic 0.703 0.000 0.712 0.000 0.650 0.000 0.677 0.000 0.712 0.000 0.697 0.000 0.651 0.000 0.740 0.000 0.695 0.000
Male 0.682 0.000 0.682 0.000 0.673 0.000 0.673 0.000 0.663 0.000 0.655 0.000 0.636 0.000 0.656 0.000 0.653 0.000
Married 0.860 0.000 0.852 0.000 0.862 0.000 0.838 0.000 0.856 0.000 0.853 0.000 0.839 0.000 0.844 0.000 0.844 0.000
Income.75K 0.789 0.000 0.806 0.000 0.771 0.000 0.797 0.000 0.794 0.000 0.773 0.000 0.790 0.000 0.822 0.000 0.848 0.000
College.Graduate 0.686 0.000 0.650 0.000 0.673 0.000 0.669 0.000 0.663 0.000 0.671 0.000 0.668 0.000 0.688 0.000 0.718 0.000
Employed.for.Wages 0.618 0.000 0.603 0.000 0.611 0.000 0.608 0.000 0.604 0.000 0.607 0.000 0.615 0.000 0.637 0.000 0.675 0.000
as.factor(Region)2 1.012 0.471 1.177 0.000 1.050 0.005 1.189 0.000 1.256 0.000 1.203 0.000 1.154 0.000 1.222 0.000 1.179 0.000
as.factor(Region)3 0.885 0.000 0.883 0.000 0.950 0.000 0.884 0.000 0.960 0.006 0.914 0.000 0.890 0.000 0.834 0.000 0.854 0.000
as.factor(Region)4 0.830 0.000 0.850 0.000 0.835 0.000 0.829 0.000 0.881 0.000 0.884 0.000 0.821 0.000 0.865 0.000 0.853 0.000
as.factor(Region)5 0.840 0.000 0.828 0.000 0.899 0.000 0.928 0.000 0.939 0.002 0.929 0.000 0.892 0.000 0.897 0.000 0.862 0.000
as.factor(Region)6 0.756 0.000 0.858 0.000 0.832 0.000 0.757 0.000 0.807 0.000 0.748 0.000 0.787 0.000 0.785 0.000 0.786 0.000
as.factor(Region)7 0.914 0.000 0.878 0.000 0.919 0.000 0.924 0.000 0.958 0.001 0.885 0.000 0.856 0.000 0.914 0.000 0.887 0.000
as.factor(Region)8 0.861 0.000 1.059 0.123 1.077 0.046 1.093 0.019 1.062 0.130 1.066 0.085 0.956 0.271 1.030 0.467 0.942 0.157
as.factor(Region)9 0.823 0.000 0.839 0.000 0.828 0.000 0.825 0.000 0.913 0.000 0.885 0.000 0.809 0.000 0.874 0.000 0.820 0.000
as.factor(Region)10 0.838 0.000 0.875 0.000 0.874 0.000 0.778 0.000 0.905 0.000 0.921 0.000 0.922 0.000 0.876 0.000 0.861 0.000
Veteran 1.160 0.000 1.232 0.000 1.201 0.000 1.192 0.000 1.246 0.000 1.300 0.000 1.312 0.000 1.327 0.000 1.292 0.000
myprint2(do.call(cbind,t8), a1[8],mynames)
Depression
2003 p 2004 p 2005 p 2006 p 2007 p 2008 p 2009 p 2010 p 2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 1.427 0.000 1.349 0.000 1.289 0.000 1.437 0 1.430 0.000 1.335 0.000 1.308 0.000 1.264 0.000 1.258 0.000 1.308 0.000 1.246 0.000 1.255 0.000 1.286 0.000 1.353 0.000 1.500 0.000 1.633 0.000 1.874 0.000
as.factor(Age)2 0.919 0.000 0.927 0.000 0.901 0.000 0.932 0 0.910 0.000 0.920 0.000 0.977 0.067 0.997 0.802 0.971 0.007 0.944 0.000 0.958 0.000 0.993 0.538 0.899 0.000 0.883 0.000 0.870 0.000 0.806 0.000 0.800 0.000
as.factor(Age)3 0.803 0.000 0.834 0.000 0.861 0.000 0.870 0 0.881 0.000 0.904 0.000 0.970 0.016 0.977 0.076 0.991 0.432 0.926 0.000 0.886 0.000 0.925 0.000 0.848 0.000 0.794 0.000 0.764 0.000 0.673 0.000 0.642 0.000
as.factor(Age)4 0.706 0.000 0.737 0.000 0.756 0.000 0.769 0 0.748 0.000 0.832 0.000 0.855 0.000 0.901 0.000 0.879 0.000 0.840 0.000 0.832 0.000 0.815 0.000 0.748 0.000 0.722 0.000 0.652 0.000 0.563 0.000 0.516 0.000
as.factor(Age)5 0.448 0.000 0.501 0.000 0.527 0.000 0.559 0 0.568 0.000 0.588 0.000 0.615 0.000 0.678 0.000 0.668 0.000 0.660 0.000 0.661 0.000 0.663 0.000 0.590 0.000 0.540 0.000 0.513 0.000 0.453 0.000 0.399 0.000
as.factor(Age)6 0.226 0.000 0.240 0.000 0.237 0.000 0.245 0 0.254 0.000 0.261 0.000 0.283 0.000 0.306 0.000 0.299 0.000 0.311 0.000 0.306 0.000 0.325 0.000 0.296 0.000 0.263 0.000 0.261 0.000 0.222 0.000 0.208 0.000
Caucasian 1.113 0.000 1.152 0.000 1.177 0.000 1.162 0 1.136 0.000 1.194 0.000 1.121 0.000 1.163 0.000 1.166 0.000 1.167 0.000 1.163 0.000 1.210 0.000 1.279 0.000 1.291 0.000 1.292 0.000 1.318 0.000 1.361 0.000
Hispanic 0.936 0.000 0.922 0.000 0.841 0.000 0.915 0 0.849 0.000 0.953 0.000 1.008 0.488 1.014 0.236 1.021 0.057 0.967 0.002 0.936 0.000 0.944 0.000 0.940 0.000 0.810 0.000 0.847 0.000 0.814 0.000 0.812 0.000
Male 0.635 0.000 0.616 0.000 0.623 0.000 0.636 0 0.653 0.000 0.633 0.000 0.648 0.000 0.656 0.000 0.669 0.000 0.656 0.000 0.658 0.000 0.657 0.000 0.654 0.000 0.633 0.000 0.635 0.000 0.622 0.000 0.614 0.000
Married 0.687 0.000 0.709 0.000 0.715 0.000 0.695 0 0.691 0.000 0.683 0.000 0.694 0.000 0.715 0.000 0.698 0.000 0.703 0.000 0.697 0.000 0.687 0.000 0.671 0.000 0.656 0.000 0.657 0.000 0.657 0.000 0.640 0.000
Income.75K 0.904 0.000 0.855 0.000 0.832 0.000 0.848 0 0.826 0.000 0.847 0.000 0.796 0.000 0.793 0.000 0.790 0.000 0.788 0.000 0.803 0.000 0.824 0.000 0.877 0.000 0.863 0.000 0.864 0.000 0.857 0.000 0.922 0.000
College.Graduate 0.901 0.000 0.907 0.000 0.914 0.000 0.909 0 0.940 0.000 0.921 0.000 0.934 0.000 0.880 0.000 0.924 0.000 0.899 0.000 0.936 0.000 0.923 0.000 0.930 0.000 0.991 0.260 0.980 0.013 1.011 0.176 1.041 0.000
Employed.for.Wages 0.795 0.000 0.818 0.000 0.796 0.000 0.790 0 0.784 0.000 0.777 0.000 0.778 0.000 0.740 0.000 0.720 0.000 0.741 0.000 0.727 0.000 0.758 0.000 0.756 0.000 0.751 0.000 0.792 0.000 0.797 0.000 0.812 0.000
as.factor(Region)2 0.794 0.000 0.764 0.000 0.791 0.000 0.741 0 0.708 0.000 0.872 0.000 0.802 0.000 0.744 0.000 0.839 0.000 0.864 0.000 0.834 0.000 0.845 0.000 0.858 0.000 0.949 0.001 0.907 0.000 0.964 0.023 1.010 0.549
as.factor(Region)3 0.996 0.814 1.017 0.255 0.978 0.103 0.841 0 0.889 0.000 0.885 0.000 0.896 0.000 0.919 0.000 0.943 0.000 0.930 0.000 0.944 0.000 0.889 0.000 0.897 0.000 0.968 0.008 0.953 0.000 0.907 0.000 0.848 0.000
as.factor(Region)4 1.041 0.043 0.994 0.763 1.011 0.510 0.854 0 0.979 0.161 0.942 0.000 0.954 0.002 0.960 0.006 0.951 0.000 0.963 0.009 0.976 0.087 0.875 0.000 0.966 0.021 1.022 0.126 0.949 0.000 0.990 0.487 1.036 0.018
as.factor(Region)5 0.983 0.425 0.945 0.006 0.913 0.000 0.816 0 0.931 0.000 0.941 0.001 0.929 0.000 0.874 0.000 0.966 0.032 0.976 0.145 0.993 0.664 0.928 0.000 0.958 0.013 1.023 0.178 0.974 0.132 1.002 0.920 0.881 0.000
as.factor(Region)6 1.236 0.000 1.081 0.000 1.219 0.000 1.144 0 1.058 0.000 1.066 0.000 1.041 0.001 1.057 0.000 1.054 0.000 1.087 0.000 1.079 0.000 0.932 0.000 1.029 0.015 1.075 0.000 0.957 0.000 1.084 0.000 0.994 0.645
as.factor(Region)7 0.844 0.000 0.849 0.000 0.914 0.000 0.806 0 0.820 0.000 0.834 0.000 0.782 0.000 0.817 0.000 0.871 0.000 0.866 0.000 0.865 0.000 0.808 0.000 0.856 0.000 0.839 0.000 0.827 0.000 0.896 0.000 0.874 0.000
as.factor(Region)8 0.441 0.000 0.684 0.000 0.943 0.104 0.696 0 0.723 0.000 0.682 0.000 0.435 0.000 0.470 0.000 0.487 0.000 0.473 0.000 0.595 0.000 0.485 0.000 0.517 0.000 0.630 0.000 0.576 0.000 0.664 0.000 0.552 0.000
as.factor(Region)9 0.854 0.000 0.844 0.000 0.893 0.000 0.757 0 0.782 0.000 0.820 0.000 0.850 0.000 0.802 0.000 0.808 0.000 0.803 0.000 0.804 0.000 0.753 0.000 0.783 0.000 0.834 0.000 0.794 0.000 0.889 0.000 0.884 0.000
as.factor(Region)10 0.881 0.000 0.885 0.000 0.730 0.000 0.721 0 0.790 0.000 0.778 0.000 0.889 0.000 0.897 0.000 0.910 0.000 0.873 0.000 0.764 0.000 0.731 0.000 0.788 0.000 0.890 0.000 0.897 0.000 0.919 0.000 0.921 0.000
Veteran 0.951 0.001 0.981 0.161 0.988 0.340 0.947 0 0.912 0.000 0.962 0.002 0.969 0.010 0.946 0.000 0.938 0.000 0.944 0.000 0.957 0.000 0.960 0.000 0.982 0.126 1.002 0.893 0.990 0.410 0.980 0.092 1.010 0.420
myprint2(do.call(cbind,t9), a1[9], mynames[17:34])
Kidney.Disease
2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 0.015 0.000 0.013 0.000 0.013 0.000 0.013 0.000 0.010 0.000 0.011 0.000 0.017 0.000 0.013 0.000 0.011 0.000
as.factor(Age)2 1.224 0.000 1.914 0.000 1.574 0.000 1.878 0.000 2.161 0.000 2.160 0.000 1.666 0.000 1.556 0.000 2.181 0.000
as.factor(Age)3 2.234 0.000 2.502 0.000 2.281 0.000 3.198 0.000 3.702 0.000 3.790 0.000 2.599 0.000 2.450 0.000 2.869 0.000
as.factor(Age)4 3.179 0.000 3.906 0.000 3.723 0.000 4.463 0.000 5.057 0.000 5.810 0.000 3.732 0.000 4.226 0.000 4.951 0.000
as.factor(Age)5 4.493 0.000 5.673 0.000 4.976 0.000 5.959 0.000 7.197 0.000 7.410 0.000 5.290 0.000 5.993 0.000 6.777 0.000
as.factor(Age)6 5.176 0.000 6.186 0.000 5.971 0.000 7.864 0.000 8.833 0.000 9.587 0.000 6.593 0.000 8.067 0.000 9.153 0.000
Caucasian 0.858 0.000 0.872 0.000 0.874 0.000 0.818 0.000 0.949 0.035 0.919 0.000 0.843 0.000 0.791 0.000 0.836 0.000
Hispanic 1.069 0.042 1.158 0.000 1.064 0.047 0.980 0.538 1.080 0.024 1.071 0.024 0.809 0.000 1.090 0.005 1.110 0.001
Male 0.944 0.005 0.974 0.202 0.899 0.000 0.884 0.000 0.981 0.357 0.950 0.007 0.900 0.000 1.042 0.036 1.077 0.000
Married 0.870 0.000 0.832 0.000 0.857 0.000 0.817 0.000 0.788 0.000 0.817 0.000 0.848 0.000 0.843 0.000 0.834 0.000
Income.75K 0.776 0.000 0.719 0.000 0.658 0.000 0.741 0.000 0.775 0.000 0.736 0.000 0.849 0.000 0.784 0.000 0.762 0.000
College.Graduate 0.846 0.000 0.719 0.000 0.814 0.000 0.829 0.000 0.791 0.000 0.772 0.000 0.735 0.000 0.823 0.000 0.821 0.000
Employed.for.Wages 0.504 0.000 0.460 0.000 0.480 0.000 0.480 0.000 0.485 0.000 0.449 0.000 0.446 0.000 0.489 0.000 0.463 0.000
as.factor(Region)2 0.875 0.003 0.897 0.016 0.974 0.563 1.061 0.165 1.000 0.995 1.022 0.590 1.069 0.101 1.181 0.000 1.090 0.043
as.factor(Region)3 0.959 0.224 0.817 0.000 0.891 0.001 0.807 0.000 0.778 0.000 0.743 0.000 0.812 0.000 0.939 0.076 0.840 0.000
as.factor(Region)4 1.093 0.027 1.085 0.047 1.148 0.001 1.082 0.055 1.036 0.394 0.899 0.007 1.042 0.303 0.981 0.638 0.977 0.579
as.factor(Region)5 0.801 0.000 0.810 0.000 0.903 0.049 1.012 0.813 0.817 0.000 0.791 0.000 0.959 0.381 0.924 0.116 0.789 0.000
as.factor(Region)6 0.916 0.009 0.967 0.325 1.158 0.000 0.926 0.024 0.897 0.002 0.914 0.004 1.054 0.099 0.934 0.044 0.940 0.066
as.factor(Region)7 0.955 0.138 1.022 0.478 1.081 0.013 1.070 0.028 0.935 0.036 0.895 0.000 0.921 0.006 1.011 0.728 1.094 0.003
as.factor(Region)8 0.587 0.000 0.620 0.000 0.791 0.008 0.765 0.003 0.696 0.000 0.786 0.003 0.954 0.584 0.744 0.001 0.752 0.001
as.factor(Region)9 0.793 0.000 0.887 0.008 0.839 0.000 0.886 0.008 0.842 0.000 0.868 0.001 0.895 0.010 0.858 0.001 0.888 0.009
as.factor(Region)10 1.042 0.247 1.007 0.838 1.190 0.000 1.038 0.294 0.966 0.344 0.909 0.005 1.028 0.422 1.192 0.000 1.123 0.001
Veteran 1.103 0.001 1.181 0.000 1.201 0.000 1.156 0.000 1.112 0.000 1.195 0.000 1.160 0.000 1.123 0.000 1.117 0.000
myprint2(do.call(cbind,t10), a1[10], mynames)
Diabetes
2003 p 2004 p 2005 p 2006 p 2007 p 2008 p 2009 p 2010 p 2011 p 2012 p 2013 p 2014 p 2015 p 2016 p 2017 p 2018 p 2019 p
(Intercept) 0.017 0.000 0.013 0.000 0.014 0.000 0.011 0.000 0.016 0.000 0.017 0.000 0.019 0.000 0.014 0.000 0.015 0.000 0.016 0.000 0.018 0.000 0.014 0.000 0.012 0.000 0.014 0.000 0.017 0.000 0.020 0.000 0.017 0.000
as.factor(Age)2 2.639 0.000 2.301 0.000 3.011 0.000 2.673 0.000 2.375 0.000 2.713 0.000 2.336 0.000 3.393 0.000 2.724 0.000 2.454 0.000 2.102 0.000 2.826 0.000 2.690 0.000 2.501 0.000 2.146 0.000 2.073 0.000 2.232 0.000
as.factor(Age)3 6.349 0.000 6.928 0.000 7.224 0.000 7.366 0.000 6.374 0.000 6.910 0.000 5.797 0.000 7.549 0.000 7.493 0.000 7.141 0.000 5.688 0.000 7.700 0.000 8.649 0.000 7.786 0.000 6.456 0.000 5.658 0.000 6.286 0.000
as.factor(Age)4 12.843 0.000 14.253 0.000 14.701 0.000 16.868 0.000 13.464 0.000 13.947 0.000 11.736 0.000 15.982 0.000 16.041 0.000 14.253 0.000 13.060 0.000 17.687 0.000 18.601 0.000 17.603 0.000 15.131 0.000 13.372 0.000 15.597 0.000
as.factor(Age)5 22.863 0.000 25.695 0.000 27.807 0.000 31.082 0.000 24.837 0.000 25.488 0.000 21.794 0.000 30.038 0.000 26.824 0.000 24.036 0.000 21.591 0.000 29.721 0.000 32.508 0.000 28.778 0.000 24.965 0.000 20.964 0.000 24.443 0.000
as.factor(Age)6 21.482 0.000 26.426 0.000 28.153 0.000 33.012 0.000 26.984 0.000 27.232 0.000 24.232 0.000 32.682 0.000 31.111 0.000 28.025 0.000 26.149 0.000 36.072 0.000 39.365 0.000 35.190 0.000 29.930 0.000 26.515 0.000 30.916 0.000
Caucasian 0.556 0.000 0.565 0.000 0.566 0.000 0.571 0.000 0.585 0.000 0.564 0.000 0.566 0.000 0.571 0.000 0.595 0.000 0.607 0.000 0.576 0.000 0.598 0.000 0.584 0.000 0.627 0.000 0.590 0.000 0.594 0.000 0.596 0.000
Hispanic 0.879 0.000 0.876 0.000 0.948 0.027 0.962 0.105 1.003 0.902 0.986 0.507 0.993 0.707 0.979 0.272 0.985 0.403 1.072 0.000 1.032 0.061 1.043 0.016 0.981 0.283 0.999 0.950 1.031 0.080 0.970 0.074 1.088 0.000
Male 1.177 0.000 1.289 0.000 1.253 0.000 1.353 0.000 1.233 0.000 1.286 0.000 1.290 0.000 1.289 0.000 1.173 0.000 1.227 0.000 1.147 0.000 1.243 0.000 1.244 0.000 1.214 0.000 1.228 0.000 1.212 0.000 1.219 0.000
Married 0.910 0.000 0.910 0.000 0.913 0.000 0.894 0.000 0.900 0.000 0.879 0.000 0.887 0.000 0.881 0.000 0.935 0.000 0.915 0.000 0.951 0.000 0.922 0.000 0.945 0.000 0.946 0.000 0.957 0.000 0.966 0.001 0.919 0.000
Income.75K 0.692 0.000 0.633 0.000 0.668 0.000 0.648 0.000 0.674 0.000 0.678 0.000 0.694 0.000 0.693 0.000 0.676 0.000 0.691 0.000 0.682 0.000 0.697 0.000 0.743 0.000 0.684 0.000 0.685 0.000 0.689 0.000 0.724 0.000
College.Graduate 0.692 0.000 0.716 0.000 0.662 0.000 0.674 0.000 0.715 0.000 0.689 0.000 0.664 0.000 0.693 0.000 0.669 0.000 0.657 0.000 0.651 0.000 0.675 0.000 0.665 0.000 0.673 0.000 0.633 0.000 0.668 0.000 0.658 0.000
Employed.for.Wages 0.641 0.000 0.650 0.000 0.664 0.000 0.704 0.000 0.686 0.000 0.666 0.000 0.715 0.000 0.663 0.000 0.675 0.000 0.679 0.000 0.668 0.000 0.673 0.000 0.647 0.000 0.647 0.000 0.661 0.000 0.661 0.000 0.687 0.000
as.factor(Region)2 1.135 0.000 1.123 0.001 1.133 0.000 1.302 0.000 1.186 0.000 1.103 0.000 1.146 0.000 1.150 0.000 1.129 0.000 1.111 0.000 1.199 0.000 1.212 0.000 1.279 0.000 1.232 0.000 1.204 0.000 1.249 0.000 1.224 0.000
as.factor(Region)3 0.902 0.000 1.000 0.989 0.951 0.040 0.963 0.125 0.913 0.000 0.881 0.000 0.917 0.000 0.960 0.042 0.958 0.023 0.901 0.000 0.966 0.067 0.936 0.001 0.897 0.000 0.927 0.000 0.928 0.000 0.965 0.071 0.891 0.000
as.factor(Region)4 0.734 0.000 0.846 0.000 0.833 0.000 0.936 0.039 0.828 0.000 0.798 0.000 0.825 0.000 0.815 0.000 0.864 0.000 0.839 0.000 0.855 0.000 0.820 0.000 0.816 0.000 0.828 0.000 0.802 0.000 0.845 0.000 0.814 0.000
as.factor(Region)5 0.827 0.000 0.925 0.055 0.877 0.000 0.916 0.016 0.881 0.000 0.838 0.000 0.868 0.000 0.850 0.000 0.908 0.000 0.882 0.000 0.895 0.000 0.922 0.004 0.879 0.000 0.895 0.000 0.892 0.000 0.865 0.000 0.807 0.000
as.factor(Region)6 0.872 0.000 0.962 0.166 0.825 0.000 0.985 0.543 0.801 0.000 0.844 0.000 0.905 0.000 0.830 0.000 0.801 0.000 0.883 0.000 0.923 0.000 0.887 0.000 0.903 0.000 0.888 0.000 0.877 0.000 0.892 0.000 0.805 0.000
as.factor(Region)7 0.950 0.049 1.039 0.129 1.012 0.594 1.079 0.001 0.979 0.270 0.968 0.101 1.005 0.778 0.992 0.656 1.011 0.539 0.971 0.090 1.015 0.373 0.976 0.160 0.988 0.496 1.001 0.951 0.927 0.000 1.040 0.027 0.963 0.035
as.factor(Region)8 0.995 0.933 1.051 0.425 1.052 0.353 1.026 0.646 0.932 0.150 0.868 0.005 0.932 0.139 0.922 0.086 0.968 0.461 1.023 0.594 0.935 0.122 0.972 0.519 1.057 0.229 0.945 0.195 0.981 0.661 0.925 0.097 0.916 0.059
as.factor(Region)9 0.826 0.000 0.921 0.021 0.896 0.001 0.925 0.015 0.819 0.000 0.882 0.000 0.849 0.000 0.890 0.000 0.940 0.011 0.919 0.001 0.898 0.000 0.968 0.187 0.943 0.024 0.926 0.001 0.909 0.000 0.971 0.243 0.922 0.001
as.factor(Region)10 0.994 0.850 1.084 0.006 1.021 0.418 1.098 0.000 1.111 0.000 1.046 0.048 1.030 0.177 1.029 0.174 1.067 0.001 1.034 0.095 1.065 0.001 1.043 0.036 1.120 0.000 1.051 0.010 1.095 0.000 1.154 0.000 1.062 0.003
Veteran 1.169 0.000 1.103 0.000 1.050 0.012 1.095 0.000 1.142 0.000 1.068 0.000 1.091 0.000 1.065 0.000 1.156 0.000 1.124 0.000 1.184 0.000 1.057 0.000 1.099 0.000 1.122 0.000 1.061 0.000 1.078 0.000 1.095 0.000
colnames(or1)=c(a1[1],"p", a1[2],"p", a1[3], "p", a1[4], "p", a1[5], "p",
               a1[6],"p", a1[7],"p", a1[8], "p", a1[9], "p", a1[10], "p")

myprint(or1)
Overweight.Obese p Heart.Disease p Stroke p Skin.Cancer p Cancer p COPD p Arthritis p Depression p Kidney.Disease p Diabetes p
(Intercept) 0.478 0.000 0.004 0.000 0.007 0.000 0.001 0.000 0.009 0.000 0.035 0.000 0.055 0.000 1.356 0.000 0.012 0.000 0.015 0.000
as.factor(Age)2 2.010 0.000 2.074 0.000 2.753 0.000 1.504 0.000 2.663 0.000 1.834 0.000 2.870 0.000 0.915 0.000 1.728 0.000 2.402 0.000
as.factor(Age)3 2.827 0.000 4.645 0.000 5.724 0.000 3.518 0.000 4.340 0.000 2.863 0.000 6.354 0.000 0.844 0.000 2.765 0.000 6.931 0.000
as.factor(Age)4 3.305 0.000 11.432 0.000 10.822 0.000 8.587 0.000 7.742 0.000 4.831 0.000 13.269 0.000 0.748 0.000 4.158 0.000 15.439 0.000
as.factor(Age)5 3.526 0.000 21.610 0.000 15.455 0.000 15.812 0.000 12.604 0.000 6.285 0.000 22.293 0.000 0.587 0.000 5.744 0.000 25.632 0.000
as.factor(Age)6 3.111 0.000 32.130 0.000 19.386 0.000 33.477 0.000 21.720 0.000 5.484 0.000 27.973 0.000 0.283 0.000 7.123 0.000 30.954 0.000
Caucasian 0.960 0.000 1.128 0.000 0.768 0.000 6.149 0.000 1.342 0.000 1.291 0.000 1.205 0.000 1.235 0.000 0.864 0.000 0.596 0.000
Hispanic 1.155 0.000 0.906 0.000 0.655 0.000 1.171 0.000 0.812 0.000 0.653 0.000 0.693 0.000 0.903 0.000 1.033 0.275 1.013 0.398
Male 1.889 0.000 1.676 0.000 1.093 0.000 1.058 0.000 0.618 0.000 0.816 0.000 0.664 0.000 0.648 0.000 0.948 0.000 1.212 0.000
Married 1.072 0.000 0.877 0.000 0.718 0.000 1.086 0.000 0.953 0.000 0.643 0.000 0.850 0.000 0.677 0.000 0.833 0.000 0.943 0.000
Income.75K 1.044 0.000 0.749 0.000 0.579 0.000 1.180 0.000 0.953 0.000 0.515 0.000 0.794 0.000 0.836 0.000 0.759 0.000 0.694 0.000
College.Graduate 0.689 0.000 0.739 0.000 0.639 0.000 1.249 0.000 1.016 0.060 0.473 0.000 0.671 0.000 0.949 0.000 0.789 0.000 0.661 0.000
Employed.for.Wages 1.186 0.000 0.472 0.000 0.351 0.000 0.780 0.000 0.651 0.000 0.460 0.000 0.613 0.000 0.756 0.000 0.473 0.000 0.663 0.000
as.factor(Region)2 1.076 0.000 1.217 0.000 1.283 0.000 1.480 0.000 1.037 0.011 1.297 0.000 1.157 0.000 0.882 0.000 1.016 0.483 1.203 0.000
as.factor(Region)3 0.817 0.000 0.999 0.976 0.869 0.000 0.977 0.132 0.986 0.321 0.887 0.000 0.899 0.000 0.928 0.000 0.839 0.000 0.934 0.000
as.factor(Region)4 0.787 0.000 0.802 0.000 0.885 0.000 1.476 0.000 0.994 0.660 0.840 0.000 0.849 0.000 0.961 0.000 1.039 0.072 0.832 0.000
as.factor(Region)5 0.790 0.000 0.890 0.000 0.840 0.000 1.062 0.000 1.039 0.004 0.870 0.000 0.894 0.000 0.976 0.005 0.877 0.000 0.892 0.000
as.factor(Region)6 0.743 0.000 0.843 0.000 0.845 0.000 1.417 0.000 1.027 0.102 0.766 0.000 0.790 0.000 1.035 0.000 0.968 0.209 0.883 0.000
as.factor(Region)7 0.883 0.000 1.022 0.131 1.027 0.134 1.562 0.000 0.996 0.713 1.022 0.092 0.905 0.000 0.853 0.000 0.982 0.344 0.990 0.370
as.factor(Region)8 0.986 0.358 1.975 0.000 0.597 0.000 0.890 0.048 0.902 0.002 0.785 0.000 1.025 0.163 0.553 0.000 0.745 0.000 0.976 0.264
as.factor(Region)9 0.961 0.000 0.913 0.000 0.976 0.179 1.067 0.000 0.998 0.887 0.849 0.000 0.849 0.000 0.808 0.000 0.860 0.000 0.935 0.000
as.factor(Region)10 0.980 0.061 1.052 0.017 1.107 0.000 1.306 0.000 0.979 0.251 0.962 0.039 0.873 0.000 0.845 0.000 1.044 0.138 1.080 0.000
Veteran 1.239 0.000 1.340 0.000 1.276 0.000 1.302 0.000 1.509 0.000 1.355 0.000 1.247 0.000 0.972 0.001 1.152 0.000 1.108 0.000
as.factor(Year)2012 1.009 0.324 1.048 0.008 0.988 0.584 0.961 0.007 0.976 0.110 1.029 0.078 1.069 0.000 1.011 0.254 1.061 0.022 1.040 0.003
as.factor(Year)2013 1.006 0.532 0.956 0.010 0.998 0.938 0.999 0.926 1.006 0.705 1.049 0.002 1.022 0.022 0.939 0.000 1.047 0.083 1.040 0.004
as.factor(Year)2014 0.974 0.003 0.992 0.629 1.045 0.041 0.979 0.139 0.974 0.086 1.069 0.000 1.054 0.000 0.941 0.000 1.093 0.000 1.066 0.000
as.factor(Year)2015 0.956 0.000 0.904 0.000 1.029 0.185 1.000 0.999 1.003 0.831 1.019 0.233 0.988 0.223 0.970 0.001 1.054 0.034 1.055 0.000
as.factor(Year)2016 0.960 0.000 0.965 0.046 1.076 0.001 0.990 0.490 0.988 0.418 1.067 0.000 1.034 0.000 0.982 0.048 1.163 0.000 1.088 0.000
as.factor(Year)2017 0.982 0.048 0.890 0.000 1.087 0.000 0.989 0.493 1.018 0.261 1.076 0.000 0.972 0.006 1.043 0.000 1.243 0.000 1.102 0.000
as.factor(Year)2018 1.000 0.992 0.958 0.028 1.156 0.000 1.013 0.413 1.020 0.195 1.142 0.000 1.056 0.000 1.076 0.000 1.224 0.000 1.165 0.000