Homework 5 Stats II

Author

Jacob Souch

library(viridis)

Loading required package: viridisLite

library(stringr)

library(purrr)
library(haven)
library(janitor)


Attaching package: 'janitor'

The following objects are masked from 'package:stats':

    chisq.test, fisher.test

library(car)

Loading required package: carData


Attaching package: 'car'

The following object is masked from 'package:purrr':

    some

library(usmap)
library(ggthemes)
library(tigris)

To enable caching of data, set `options(tigris_use_cache = TRUE)`
in your R script or .Rprofile.

library(tidyverse)

── Attaching packages ─────────────────────────────────────── tidyverse 1.3.2
──

✔ ggplot2 3.4.0      ✔ dplyr   1.0.10
✔ tibble  3.1.8      ✔ stringr 1.5.0 
✔ tidyr   1.3.0      ✔ forcats 0.5.2 
✔ readr   2.1.3      
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag()    masks stats::lag()
✖ dplyr::recode() masks car::recode()
✖ car::some()     masks purrr::some()

library(tidycensus)
library(lme4)

Loading required package: Matrix

Attaching package: 'Matrix'

The following objects are masked from 'package:tidyr':

    expand, pack, unpack

library(haven)
library(foreign)
library(FactoMineR)
library(ggrepel)
library(margins)
library(emmeans)
library(stargazer)


Please cite as: 

 Hlavac, Marek (2022). stargazer: Well-Formatted Regression and Summary Statistics Tables.
 R package version 5.2.3. https://CRAN.R-project.org/package=stargazer

library(stringr)

Ingest

setwd("C:/Users/jacob/OneDrive - University of Texas at San Antonio/Courses/Stats for Demographic Data 2/Homework 4")

brfss21 <- haven::read_xpt("LLCP2021.XPT_")

str(brfss21, list.len = ncol(brfss21))

tibble [438,693 × 303] (S3: tbl_df/tbl/data.frame)
 $ _STATE  : num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "STATE FIPS CODE"
 $ FMONTH  : num [1:438693] 1 1 1 1 1 1 1 1 2 2 ...
  ..- attr(*, "label")= chr "FILE MONTH"
 $ IDATE   : chr [1:438693] "01192021" "01212021" "01212021" "01172021" ...
  ..- attr(*, "label")= chr "INTERVIEW DATE"
 $ IMONTH  : chr [1:438693] "01" "01" "01" "01" ...
  ..- attr(*, "label")= chr "INTERVIEW MONTH"
 $ IDAY    : chr [1:438693] "19" "21" "21" "17" ...
  ..- attr(*, "label")= chr "INTERVIEW DAY"
 $ IYEAR   : chr [1:438693] "2021" "2021" "2021" "2021" ...
  ..- attr(*, "label")= chr "INTERVIEW YEAR"
 $ DISPCODE: num [1:438693] 1100 1100 1100 1100 1100 1100 1100 1100 1100 1100 ...
  ..- attr(*, "label")= chr "FINAL DISPOSITION"
 $ SEQNO   : chr [1:438693] "2021000001" "2021000002" "2021000003" "2021000004" ...
  ..- attr(*, "label")= chr "ANNUAL SEQUENCE NUMBER"
 $ _PSU    : num [1:438693] 2.02e+09 2.02e+09 2.02e+09 2.02e+09 2.02e+09 ...
  ..- attr(*, "label")= chr "PRIMARY SAMPLING UNIT"
 $ CTELENM1: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "CORRECT TELEPHONE NUMBER?"
 $ PVTRESD1: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "PRIVATE RESIDENCE?"
 $ COLGHOUS: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DO YOU LIVE IN COLLEGE HOUSING?"
 $ STATERE1: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "RESIDENT OF STATE"
 $ CELPHON1: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "CELLULAR TELEPHONE"
 $ LADULT1 : num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "ARE YOU 18 YEARS OF AGE OR OLDER?"
 $ COLGSEX : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "ARE YOU MALE OR FEMALE?"
 $ NUMADULT: num [1:438693] 2 2 2 2 2 1 2 1 1 2 ...
  ..- attr(*, "label")= chr "NUMBER OF ADULTS IN HOUSEHOLD"
 $ LANDSEX : num [1:438693] NA NA NA NA NA 1 NA 1 2 NA ...
  ..- attr(*, "label")= chr "ARE YOU MALE OR FEMALE?"
 $ NUMMEN  : num [1:438693] 1 1 1 1 1 NA 1 NA NA 1 ...
  ..- attr(*, "label")= chr "NUMBER OF ADULT MEN IN HOUSEHOLD"
 $ NUMWOMEN: num [1:438693] 1 1 1 1 1 NA 1 NA NA 1 ...
  ..- attr(*, "label")= chr "NUMBER OF ADULT WOMEN IN HOUSEHOLD"
 $ RESPSLCT: num [1:438693] 2 2 2 2 1 NA 1 NA NA 2 ...
  ..- attr(*, "label")= chr "RESPONDENT SELECTION"
 $ SAFETIME: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "SAFE TIME TO TALK?"
 $ CTELNUM1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "CORRECT PHONE NUMBER?"
 $ CELLFON5: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "IS THIS A CELL PHONE?"
 $ CADULT1 : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "ARE YOU 18 YEARS OF AGE OR OLDER?"
 $ CELLSEX : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "ARE YOU MALE OR FEMALE?"
 $ PVTRESD3: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DO YOU LIVE IN A PRIVATE RESIDENCE?"
 $ CCLGHOUS: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DO YOU LIVE IN COLLEGE HOUSING?"
 $ CSTATE1 : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DO YOU CURRENTLY LIVE IN  ____(STATE)___"
 $ LANDLINE: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DO YOU ALSO HAVE A LANDLINE TELEPHONE?"
 $ HHADULT : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "NUMBER OF ADULTS IN HOUSEHOLD"
 $ SEXVAR  : num [1:438693] 2 2 2 2 1 1 1 1 2 2 ...
  ..- attr(*, "label")= chr "SEX OF RESPONDENT"
 $ GENHLTH : num [1:438693] 5 3 2 2 5 3 3 4 2 3 ...
  ..- attr(*, "label")= chr "GENERAL HEALTH"
 $ PHYSHLTH: num [1:438693] 20 88 88 88 30 88 30 88 88 25 ...
  ..- attr(*, "label")= chr "NUMBER OF DAYS PHYSICAL HEALTH NOT GOOD"
 $ MENTHLTH: num [1:438693] 10 88 88 10 88 88 88 88 88 5 ...
  ..- attr(*, "label")= chr "NUMBER OF DAYS MENTAL HEALTH NOT GOOD"
 $ POORHLTH: num [1:438693] 88 NA NA 88 30 NA 2 NA NA 5 ...
  ..- attr(*, "label")= chr "POOR PHYSICAL OR MENTAL HEALTH"
 $ PRIMINSR: num [1:438693] 3 1 2 2 3 3 1 2 3 3 ...
  ..- attr(*, "label")= chr "WHAT IS PRIMARY SOURCE OF HEALTH INSURAN"
 $ PERSDOC3: num [1:438693] 1 2 2 1 1 1 1 1 1 2 ...
  ..- attr(*, "label")= chr "HAVE PERSONAL HEALTH CARE PROVIDER?"
 $ MEDCOST1: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "COULD NOT AFFORD TO SEE DOCTOR"
 $ CHECKUP1: num [1:438693] 2 1 1 1 1 1 1 1 2 1 ...
  ..- attr(*, "label")= chr "LENGTH OF TIME SINCE LAST ROUTINE CHECKU"
 $ EXERANY2: num [1:438693] 2 1 2 1 1 2 2 1 2 1 ...
  ..- attr(*, "label")= chr "EXERCISE IN PAST 30 DAYS"
 $ BPHIGH6 : num [1:438693] 3 1 1 1 4 3 1 1 3 1 ...
  ..- attr(*, "label")= chr "EVER TOLD BLOOD PRESSURE HIGH"
 $ BPMEDS  : num [1:438693] NA 1 1 1 NA NA 1 1 NA 1 ...
  ..- attr(*, "label")= chr "CURRENTLY TAKING BLOOD PRESSURE MEDICATI"
 $ CHOLCHK3: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "HOW LONG SINCE CHOLESTEROL CHECKED"
 $ TOLDHI3 : num [1:438693] 1 1 2 1 1 2 2 1 2 2 ...
  ..- attr(*, "label")= chr "EVER TOLD CHOLESTEROL IS HIGH"
 $ CHOLMED3: num [1:438693] 1 1 2 2 1 2 2 1 2 1 ...
  ..- attr(*, "label")= chr "CURRENTLY TAKING MEDICINE FOR HIGH CHOLE"
 $ CVDINFR4: num [1:438693] 2 2 2 2 1 2 2 1 2 2 ...
  ..- attr(*, "label")= chr "EVER DIAGNOSED WITH HEART ATTACK"
 $ CVDCRHD4: num [1:438693] 2 1 1 2 7 2 2 1 2 2 ...
  ..- attr(*, "label")= chr "EVER DIAGNOSED WITH ANGINA OR CORONARY H"
 $ CVDSTRK3: num [1:438693] 2 2 2 2 1 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "EVER DIAGNOSED WITH A STROKE"
 $ ASTHMA3 : num [1:438693] 1 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "EVER TOLD HAD ASTHMA"
 $ ASTHNOW : num [1:438693] 1 NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "STILL HAVE ASTHMA"
 $ CHCSCNCR: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "(EVER TOLD) YOU HAD SKIN CANCER?"
 $ CHCOCNCR: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "(EVER TOLD) YOU HAD ANY OTHER TYPES OF C"
 $ CHCCOPD3: num [1:438693] 1 2 2 2 2 1 2 1 2 2 ...
  ..- attr(*, "label")= chr "EVER TOLD YOU HAD C.O.P.D. EMPHYSEMA OR"
 $ ADDEPEV3: num [1:438693] 2 2 2 2 2 2 1 2 2 2 ...
  ..- attr(*, "label")= chr "(EVER TOLD) YOU HAD A DEPRESSIVE DISORDE"
 $ CHCKDNY2: num [1:438693] 2 1 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "EVER TOLD YOU HAVE KIDNEY DISEASE?"
 $ DIABETE4: num [1:438693] 3 1 1 1 1 3 3 3 3 3 ...
  ..- attr(*, "label")= chr "(EVER TOLD) YOU HAD DIABETES"
 $ DIABAGE3: num [1:438693] NA 98 98 56 65 NA NA NA NA NA ...
  ..- attr(*, "label")= chr "AGE WHEN TOLD DIABETES"
 $ HAVARTH5: num [1:438693] 1 1 2 2 2 2 1 1 2 1 ...
  ..- attr(*, "label")= chr "TOLD HAVE ARTHRITIS"
 $ ARTHEXER: num [1:438693] 2 1 NA NA NA NA 1 2 NA 1 ...
  ..- attr(*, "label")= chr "DR. SUGGEST USE OF PHYSICAL ACTIVITY OR"
 $ ARTHEDU : num [1:438693] 2 2 NA NA NA NA 2 2 NA 2 ...
  ..- attr(*, "label")= chr "EVER TAKEN CLASS IN MANAGING ARTHRITIS O"
 $ LMTJOIN3: num [1:438693] 2 1 NA NA NA NA 1 2 NA 2 ...
  ..- attr(*, "label")= chr "LIMITED BECAUSE OF JOINT SYMPTOMS"
 $ ARTHDIS2: num [1:438693] 1 1 NA NA NA NA 2 2 NA 2 ...
  ..- attr(*, "label")= chr "DOES ARTHRITIS AFFECT WHETHER YOU WORK"
 $ JOINPAI2: num [1:438693] 8 10 NA NA NA NA 8 8 NA 4 ...
  ..- attr(*, "label")= chr "HOW BAD WAS JOINT PAIN"
 $ MARITAL : num [1:438693] 1 9 3 1 1 1 1 2 2 1 ...
  ..- attr(*, "label")= chr "MARITAL STATUS"
 $ EDUCA   : num [1:438693] 4 6 4 4 3 5 6 3 3 4 ...
  ..- attr(*, "label")= chr "EDUCATION LEVEL"
 $ RENTHOM1: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "OWN OR RENT HOME"
 $ NUMHHOL3: num [1:438693] 1 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "HOUSEHOLD TELEPHONES"
 $ NUMPHON3: num [1:438693] 1 NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "RESIDENTIAL PHONES"
 $ CPDEMO1B: num [1:438693] 1 1 1 1 8 1 2 8 1 2 ...
  ..- attr(*, "label")= chr "DO YOU HAVE A CELL PHONE FOR PERSONAL US"
 $ VETERAN3: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "ARE YOU A VETERAN"
 $ EMPLOY1 : num [1:438693] 7 8 7 7 8 7 8 2 7 7 ...
  ..- attr(*, "label")= chr "EMPLOYMENT STATUS"
 $ CHILDREN: num [1:438693] 88 88 88 88 88 88 88 88 88 88 ...
  ..- attr(*, "label")= chr "NUMBER OF CHILDREN IN HOUSEHOLD"
 $ INCOME3 : num [1:438693] 5 77 3 7 4 6 77 99 77 8 ...
  ..- attr(*, "label")= chr "INCOME LEVEL"
 $ PREGNANT: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "PREGNANCY STATUS"
 $ WEIGHT2 : num [1:438693] 72 7777 170 195 206 ...
  ..- attr(*, "label")= chr "REPORTED WEIGHT IN POUNDS"
 $ HEIGHT3 : num [1:438693] 411 506 505 504 511 603 600 509 503 505 ...
  ..- attr(*, "label")= chr "REPORTED HEIGHT IN FEET AND INCHES"
 $ DEAF    : num [1:438693] 2 2 2 2 1 2 1 2 2 2 ...
  ..- attr(*, "label")= chr "ARE YOU DEAF OR DO YOU HAVE SERIOUS DIFF"
 $ BLIND   : num [1:438693] 2 1 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "BLIND OR DIFFICULTY SEEING"
 $ DECIDE  : num [1:438693] 2 1 2 2 2 2 2 1 2 2 ...
  ..- attr(*, "label")= chr "DIFFICULTY CONCENTRATING OR REMEMBERING"
 $ DIFFWALK: num [1:438693] 2 1 2 2 1 1 1 2 2 1 ...
  ..- attr(*, "label")= chr "DIFFICULTY WALKING OR CLIMBING STAIRS"
 $ DIFFDRES: num [1:438693] 2 2 2 2 2 2 1 2 2 2 ...
  ..- attr(*, "label")= chr "DIFFICULTY DRESSING OR BATHING"
 $ DIFFALON: num [1:438693] 1 1 2 2 2 2 1 2 1 2 ...
  ..- attr(*, "label")= chr "DIFFICULTY DOING ERRANDS ALONE"
 $ SMOKE100: num [1:438693] 1 2 2 2 2 1 2 1 2 1 ...
  ..- attr(*, "label")= chr "SMOKED AT LEAST 100 CIGARETTES"
 $ SMOKDAY2: num [1:438693] 3 NA NA NA NA 3 NA 3 NA 3 ...
  ..- attr(*, "label")= chr "FREQUENCY OF DAYS NOW SMOKING"
 $ USENOW3 : num [1:438693] 3 3 3 3 2 2 3 3 3 3 ...
  ..- attr(*, "label")= chr "USE OF SMOKELESS TOBACCO PRODUCTS"
 $ ECIGNOW1: num [1:438693] 3 3 3 3 3 3 3 3 3 3 ...
  ..- attr(*, "label")= chr "DO YOU NOW USE E-CIGARETTES, EVERY DAY,"
 $ ALCDAY5 : num [1:438693] 888 888 888 101 888 888 888 888 888 203 ...
  ..- attr(*, "label")= chr "DAYS IN PAST 30 HAD ALCOHOLIC BEVERAGE"
 $ AVEDRNK3: num [1:438693] NA NA NA 3 NA NA NA NA NA 2 ...
  ..- attr(*, "label")= chr "AVG ALCOHOLIC DRINKS PER DAY IN PAST 30"
 $ DRNK3GE5: num [1:438693] NA NA NA 1 NA NA NA NA NA 88 ...
  ..- attr(*, "label")= chr "BINGE DRINKING"
 $ MAXDRNKS: num [1:438693] NA NA NA 6 NA NA NA NA NA 2 ...
  ..- attr(*, "label")= chr "MOST DRINKS ON SINGLE OCCASION PAST 30 D"
 $ FLUSHOT7: num [1:438693] 1 2 2 1 1 1 1 1 1 2 ...
  ..- attr(*, "label")= chr "ADULT FLU SHOT/SPRAY PAST 12 MOS"
 $ FLSHTMY3: num [1:438693] 92020 NA NA 102020 92020 ...
  ..- attr(*, "label")= chr "WHEN RECEIVED MOST RECENT SEASONAL FLU S"
 $ IMFVPLA2: num [1:438693] 1 NA NA 1 1 5 1 1 1 NA ...
  ..- attr(*, "label")= chr "WHERE DID YOU GET YOUR LAST FLU SHOT/VAC"
 $ PNEUVAC4: num [1:438693] 1 2 2 2 1 1 2 1 2 1 ...
  ..- attr(*, "label")= chr "PNEUMONIA SHOT EVER"
 $ HIVTST7 : num [1:438693] 2 2 2 2 1 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "EVER TESTED H.I.V."
 $ HIVTSTD3: num [1:438693] NA NA NA NA 777777 ...
  ..- attr(*, "label")= chr "MONTH AND YEAR OF LAST HIV TEST"
 $ FRUIT2  : num [1:438693] 101 101 101 203 101 202 312 204 202 101 ...
  ..- attr(*, "label")= chr "HOW MANY TIMES DID YOU EAT FRUIT?"
 $ FRUITJU2: num [1:438693] 555 555 555 205 555 555 301 555 101 555 ...
  ..- attr(*, "label")= chr "HOW MANY TIMES DID YOU DRINK 100 PERCENT"
 $ FVGREEN1: num [1:438693] 204 201 555 303 101 201 203 202 555 202 ...
  ..- attr(*, "label")= chr "HOW MANY TIMES DID YOU EAT DARK GREEN VE"
 $ FRENCHF1: num [1:438693] 203 555 201 204 202 555 203 555 777 202 ...
  ..- attr(*, "label")= chr "HOW OFTEN DO YOU EAT FRENCH FRIES OR FRI"
 $ POTATOE1: num [1:438693] 201 201 201 308 202 201 302 204 555 203 ...
  ..- attr(*, "label")= chr "HOW OFTEN DO YOU EAT FRENCH FRIES OR FRI"
 $ VEGETAB2: num [1:438693] 101 207 203 205 101 201 202 204 202 204 ...
  ..- attr(*, "label")= chr "HOW OFTEN DO YOU EAT FRENCH FRIES OR FRI"
 $ PDIABTST: num [1:438693] 2 NA NA NA NA 2 1 1 1 1 ...
  ..- attr(*, "label")= chr "HAD A TEST FOR HIGH BLOOD SUGAR IN PAST"
 $ PREDIAB1: num [1:438693] 3 NA NA NA NA 3 3 3 3 3 ...
  ..- attr(*, "label")= chr "EVER BEEN TOLD YOU HAVE PRE-DIABETES OR"
 $ INSULIN1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "NOW TAKING INSULIN"
 $ BLDSUGAR: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW OFTEN CHECK BLOOD FOR GLUCOSE"
 $ FEETCHK3: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW OFTEN CHECK FEET FOR SORES OR IRRITA"
 $ DOCTDIAB: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TIMES SEEN HEALTH PROFESSIONAL FOR DIABE"
 $ CHKHEMO3: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TIMES CHECKED FOR GLYCOSYLATED HEMOGLOBI"
 $ FEETCHK : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TIMES FEET CHECK FOR SORES/IRRITATIONS"
 $ EYEEXAM1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "LAST EYE EXAM WHERE PUPILS WERE DILATED"
 $ DIABEYE : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER TOLD DIABETES HAS AFFECTED EYES"
 $ DIABEDU : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER TAKEN CLASS IN MANAGING DIABETES"
 $ TOLDCFS : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TOLD HAD CHRONIC FATIGUE SYNDROME OR MYA"
 $ HAVECFS : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "STILL HAVE CHRONIC FATIGUE SYNDROME OR M"
 $ WORKCFS : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW MANY HOURS A WEEK ARE YOU BEEN ABLE"
 $ TOLDHEPC: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TOLD HAD HEPATITIS C"
 $ TRETHEPC: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TREATED FOR HEPATITIS C"
 $ PRIRHEPC: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "WERE YOU TREATED FOR HEPATITIS C PRIOR T"
 $ HAVEHEPC: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "STILL HAVE HEPATITIS C"
 $ HAVEHEPB: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TOLD  HAD HEPATITIS B"
 $ MEDSHEPB: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "CURRENTLY TAKING MEDICINE FOR HEPATITIS"
 $ HPVADVC4: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER HAD AN H.P.V. VACCINATION?"
 $ HPVADSHT: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "NUMBER OF HPV SHOTS RECEIVED"
 $ TETANUS1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "RECEIVED TETANUS SHOT SINCE 2005?"
 $ SHINGLE2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HAVE YOU EVER HAD THE SHINGLES OR ZOSTER"
 $ LCSFIRST: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW OLD WHEN YOU FIRST STARTED SMOKING?"
 $ LCSLAST : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW OLD WHEN YOU LAST SMOKED?"
 $ LCSNUMCG: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "ON AVERAGE, HOW MANY CIGARETTES DO YOU S"
 $ LCSCTSCN: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DID YOU HAVE A CT OR CAT SCAN?"
 $ HADMAM  : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HAVE YOU EVER HAD A MAMMOGRAM"
 $ HOWLONG : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW LONG SINCE LAST MAMMOGRAM"
 $ CERVSCRN: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HAVE YOU EVER HAD A CERVICAL CANCER SCRE"
 $ CRVCLCNC: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TIME SINCE LAST CERVICAL CANCER SCREENIN"
 $ CRVCLPAP: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HAVE A PAP TEST AND RECENT CERVICAL CANC"
 $ CRVCLHPV: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HAVE AN H.P.V. TEST AND RECENT CERVICAL"
 $ HADHYST2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HAD HYSTERECTOMY"
 $ PSATEST1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER HAD PSA TEST"
 $ PSATIME1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TIME SINCE MOST RECENT PSA TEST"
 $ PCPSARS2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "WHAT WAS THE MAIN REASON YOU HAD THIS PS"
 $ PCSTALK : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DID YOU TALK ABOUT THE ADVANTAGES OR DIS"
 $ HADSIGM4: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER HAD SIGMOIDOSCOPY/COLONOSCOPY"
 $ COLNSIGM: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER HAD A COLONOSCOPY, SIGMOIDOSCOPY, O"
 $ COLNTES1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW LONG HAS IT BEEN SINCE YOU HAD COLON"
 $ SIGMTES1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW LONG HAS IT BEEN SINCE YOU HAD SIGMO"
 $ LASTSIG4: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TIME SINCE LAST SIGMOIDOSCOPY/COLONOSCOP"
 $ COLNCNCR: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER HAD ANY OTHER KIND OF TEST FOR COLO"
 $ VIRCOLO1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HAVE YOU EVER HAD A VIRTUAL COLONOSCOPY?"
 $ VCLNTES1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW LONG HAS IT BEEN SINCE YOU HAD VIRTU"
 $ SMALSTOL: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER HAD STOOL TEST?"
 $ STOLTEST: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW LONG SINCE YOU HAD STOOL TEST?"
 $ STOOLDN1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER HAD STOOL DNA TEST?"
 $ BLDSTFIT: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "WAS TEST PART OF COLOGUARD TEST?"
 $ SDNATES1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW LONG SINCE YOU HAD STOOL DNA?"
 $ CNCRDIFF: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW MANY TYPES OF CANCER?"
 $ CNCRAGE : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "AGE TOLD HAD CANCER"
 $ CNCRTYP1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TYPE OF CANCER"
 $ CSRVTRT3: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "CURRENTLY RECEIVING TREATMENT FOR CANCER"
 $ CSRVDOC1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "WHAT TYPE OF DOCTOR PROVIDES MAJORITY OF"
 $ CSRVSUM : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DID YOU RECEIVE A SUMMARY OF CANCER TREA"
 $ CSRVRTRN: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER RECEIVE INSTRUCTIONS FROM A DOCTOR"
 $ CSRVINST: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "INSTRUCTIONS WRITTEN OR PRINTED"
 $ CSRVINSR: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DID HEALTH INSURANCE PAY FOR ALL OF YOUR"
 $ CSRVDEIN: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "EVER DENIED INSURANCE COVERAGE BECAUSE O"
 $ CSRVCLIN: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "PARTICIPATE IN CLINICAL TRIAL AS PART OF"
 $ CSRVPAIN: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "CURRENTLY HAVE PHYSICAL PAIN FROM CANCER"
 $ CSRVCTL2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "IS PAIN UNDER CONTROL?"
 $ HOMBPCHK: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "TOLD CHECK BLOOD PRESSURE AT HOME"
 $ HOMRGCHK: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "REGULARLY CHECK BLOOD PRESSURE AT HOME"
 $ WHEREBP : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "WHERE DO YOU GET BLOOD PRESSURE TAKEN"
 $ SHAREBP : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HOW DO YOU SHARE YOUR BLOOD PRESSURE NUM"
 $ WTCHSALT: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "WATCHING SODIUM OR SALT INTAKE"
 $ DRADVISE: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DOCTOR ADVISED REDUCED SODIUM/SALT INTAK"
 $ CIMEMLOS: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HAVE YOU EXPERIENCED CONFUSION OR MEMORY"
 $ CDHOUSE : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "GIVEN UP DAY-TO-DAY CHORES DUE TO CONFUS"
 $ CDASSIST: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "NEED ASSISTANCE WITH DAY-TO_DAY ACTIVITI"
 $ CDHELP  : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "WHEN YOU NEED HELP WITH DAY-TO-DAY ACTIV"
 $ CDSOCIAL: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DOES CONFUSION OR MEMORY LOSS INTERFERE"
 $ CDDISCUS: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HAVE YOU DISCUSSED YOUR CONFUSION OR MEM"
 $ CAREGIV1: num [1:438693] 2 2 2 2 2 2 2 2 2 1 ...
  ..- attr(*, "label")= chr "PROVIDED REGULAR CARE FOR FAMILY OR FRIE"
 $ CRGVREL4: num [1:438693] NA NA NA NA NA NA NA NA NA 15 ...
  ..- attr(*, "label")= chr "RELATIONSHIP OF PERSON TO WHOM YOU ARE G"
 $ CRGVLNG1: num [1:438693] NA NA NA NA NA NA NA NA NA 2 ...
  ..- attr(*, "label")= chr "HOW LONG PROVIDED CARE FOR PERSON."
 $ CRGVHRS1: num [1:438693] NA NA NA NA NA NA NA NA NA 1 ...
  ..- attr(*, "label")= chr "HOW MANY HOURS DO YOU PROVIDE CARE FOR P"
 $ CRGVPRB3: num [1:438693] NA NA NA NA NA NA NA NA NA 15 ...
  ..- attr(*, "label")= chr "WHAT IS THE MAJOR HEALTH PROBLEM, ILLNES"
 $ CRGVALZD: num [1:438693] NA NA NA NA NA NA NA NA NA 2 ...
  ..- attr(*, "label")= chr "DOES PERSON BEING CARED FOR HAVE ALZHEIM"
 $ CRGVPER1: num [1:438693] NA NA NA NA NA NA NA NA NA 2 ...
  ..- attr(*, "label")= chr "MANAGED PERSONAL CARE"
 $ CRGVHOU1: num [1:438693] NA NA NA NA NA NA NA NA NA 2 ...
  ..- attr(*, "label")= chr "MANAGED HOUSEHOLD TASKS"
 $ CRGVEXPT: num [1:438693] 2 2 2 7 2 2 2 2 7 NA ...
  ..- attr(*, "label")= chr "DO YOU EXPECT TO HAVE A RELATIVE YOU WIL"
 $ ACEDEPRS: num [1:438693] 2 2 2 1 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "LIVE WITH ANYONE DEPRESSED, MENTALLY ILL"
 $ ACEDRINK: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "LIVE WITH A PROBLEM DRINKER/ALCOHOLIC?"
 $ ACEDRUGS: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "LIVE WITH ANYONE WHO USED ILLEGAL DRUGS"
 $ ACEPRISN: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "LIVE WITH ANYONE WHO SERVED TIME IN PRIS"
 $ ACEDIVRC: num [1:438693] 1 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "WERE YOUR PARENTS DIVORCED/SEPERATED?"
 $ ACEPUNCH: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "HOW OFTEN DID YOUR PARENTS BEAT EACH OTH"
 $ ACEHURT1: num [1:438693] 1 1 1 3 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "HOW OFTEN DID A PARENT PHYSICALLY HURT Y"
 $ ACESWEAR: num [1:438693] 1 1 1 1 1 1 1 1 1 3 ...
  ..- attr(*, "label")= chr "HOW OFTEN DID A PARENT SWEAR AT YOU?"
 $ ACETOUCH: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "HOW OFTEN DID ANYONE EVER TOUCH YOU SEXU"
 $ ACETTHEM: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "HOW OFTEN DID ANYONE MAKE YOU TOUCH THEM"
 $ ACEHVSEX: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "HOW OFTEN DID ANYONE EVER FORCE YOU TO H"
 $ ACEADSAF: num [1:438693] 5 5 5 5 5 5 5 5 5 4 ...
  ..- attr(*, "label")= chr "DID AN ADULT MAKE YOU FEEL SAFE AND PROT"
 $ ACEADNED: num [1:438693] 3 5 5 5 5 5 5 5 5 4 ...
  ..- attr(*, "label")= chr "DID AN ADULT MAKE SURE BASIC NEEDS WERE"
 $ MARIJAN1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DURING THE PAST 30 DAYS, ON HOW MANY DAY"
 $ USEMRJN3: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DURING THE PAST 30 DAYS, HOW DID YOU PRI"
 $ RSNMRJN2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "WHAT WAS THE REASON YOU USED MARIJUANA?"
 $ LASTSMK2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "INTERVAL SINCE LAST SMOKED"
 $ STOPSMK2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "STOPPED SMOKING IN PAST 12 MONTHS"
 $ FIREARM5: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "ANY FIREARMS IN HOME"
 $ GUNLOAD : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "ANY FIREARMS LOADED"
 $ LOADULK2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "ANY LOADED FIREARMS ALSO UNLOCKED"
 $ RCSGENDR: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "GENDER OF CHILD"
 $ RCSRLTN2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "RELATIONSHIP TO CHILD"
 $ CASTHDX2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "HLTH PRO EVER SAID CHILD HAS ASTHMA"
 $ CASTHNO2: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "CHILD STILL HAVE ASTHMA?"
 $ BIRTHSEX: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "ARE YOU MALE OR FEMALE?"
 $ SOMALE  : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "SEXUAL ORIENTATION"
 $ SOFEMALE: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "SEXUAL ORIENTATION"
 $ TRNSGNDR: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "DO YOU CONSIDER YOURSELF TO BE TRANSGEND"
 $ QSTVER  : num [1:438693] 10 10 10 10 10 10 10 10 10 10 ...
  ..- attr(*, "label")= chr "QUESTIONNAIRE VERSION IDENTIFIER"
 $ QSTLANG : num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "LANGUAGE IDENTIFIER"
 $ _METSTAT: num [1:438693] 1 1 1 1 2 2 2 1 1 1 ...
  ..- attr(*, "label")= chr "METROPOLITAN STATUS"
 $ _URBSTAT: num [1:438693] 1 1 1 1 1 2 1 1 1 1 ...
  ..- attr(*, "label")= chr "URBAN/RURAL STATUS"
 $ MSCODE  : num [1:438693] 1 2 1 3 2 5 5 3 1 1 ...
  ..- attr(*, "label")= chr "METROPOLITAN STATUS CODE"
 $ _STSTR  : num [1:438693] 11011 11011 11011 11011 11011 ...
  ..- attr(*, "label")= chr "SAMPLE DESIGN STRATIFICATION VARIABLE"
 $ _STRWT  : num [1:438693] 39.8 39.8 39.8 39.8 39.8 ...
  ..- attr(*, "label")= chr "STRATUM WEIGHT"
 $ _RAWRAKE: num [1:438693] 2 2 2 2 2 1 2 1 1 2 ...
  ..- attr(*, "label")= chr "RAW WEIGHTING FACTOR USED IN RAKING"
 $ _WT2RAKE: num [1:438693] 79.5 79.5 79.5 79.5 79.5 ...
  ..- attr(*, "label")= chr "DESIGN WEIGHT USED IN RAKING"
 $ _IMPRACE: num [1:438693] 1 2 2 1 6 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "IMPUTED RACE/ETHNICITY VALUE"
 $ _CHISPNC: num [1:438693] 9 9 9 9 9 9 9 9 9 9 ...
  ..- attr(*, "label")= chr "CHILD HISPANIC, LATINO/A, OR SPANISH ORI"
 $ _CRACE1 : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "CHILD NON-HISPANIC RACE INCLUDING MULTIR"
 $ _CPRACE1: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "PREFERRED CHILD RACE CATEGORIES"
 $ CAGEG   : num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "FOUR LEVEL CHILD AGE"
 $ _CLLCPWT: num [1:438693] NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "label")= chr "FINAL CHILD WEIGHT: LAND-LINE AND CELL-P"
 $ _DUALUSE: num [1:438693] 1 1 1 1 9 1 1 9 1 1 ...
  ..- attr(*, "label")= chr "DUAL PHONE USE CATEGORIES"
 $ _DUALCOR: num [1:438693] 0.519 0.519 0.519 0.519 NA ...
  ..- attr(*, "label")= chr "DUAL PHONE USE CORRECTION FACTOR"
 $ _LLCPWT2: num [1:438693] 874 874 874 874 2145 ...
  ..- attr(*, "label")= chr "TRUNCATED DESIGN WEIGHT USED IN ADULT CO"
 $ _LLCPWT : num [1:438693] 745 299 588 1100 1712 ...
  ..- attr(*, "label")= chr "FINAL WEIGHT: LAND-LINE AND CELL-PHONE D"
 $ _RFHLTH : num [1:438693] 2 1 1 1 2 1 1 2 1 1 ...
  ..- attr(*, "label")= chr "ADULTS WITH GOOD OR BETTER HEALTH"
 $ _PHYS14D: num [1:438693] 3 1 1 1 3 1 3 1 1 3 ...
  ..- attr(*, "label")= chr "COMPUTED PHYSICAL HEALTH STATUS"
 $ _MENT14D: num [1:438693] 2 1 1 2 1 1 1 1 1 2 ...
  ..- attr(*, "label")= chr "COMPUTED MENTAL HEALTH STATUS"
 $ _HLTHPLN: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "HAVE ANY HEALTH INSURANCE"
 $ _HCVU652: num [1:438693] 9 9 9 1 9 9 1 1 9 9 ...
  ..- attr(*, "label")= chr "RESPONDENTS AGED 18-64 WITH HEALTH INSUR"
 $ _TOTINDA: num [1:438693] 2 1 2 1 1 2 2 1 2 1 ...
  ..- attr(*, "label")= chr "LEISURE TIME PHYSICAL ACTIVITY CALCULATE"
 $ _RFHYPE6: num [1:438693] 1 2 2 2 1 1 2 2 1 2 ...
  ..- attr(*, "label")= chr "HIGH BLOOD PRESSURE CALCULATED VARIABLE"
 $ _CHOLCH3: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "CHOLESTEROL CHECKED CALCULATED VARIABLE"
 $ _RFCHOL3: num [1:438693] 2 2 1 2 2 1 1 2 1 1 ...
  ..- attr(*, "label")= chr "HIGH CHOLESTEROL CALCULATED VARIABLE"
 $ _MICHD  : num [1:438693] 2 1 1 2 1 2 2 1 2 2 ...
  ..- attr(*, "label")= chr "RESPONDENTS THAT HAVE EVER REPORTED HAVI"
 $ _LTASTH1: num [1:438693] 2 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "LIFETIME ASTHMA CALCULATED VARIABLE"
 $ _CASTHM1: num [1:438693] 2 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "CURRENT ASTHMA CALCULATED VARIABLE"
 $ _ASTHMS1: num [1:438693] 1 3 3 3 3 3 3 3 3 3 ...
  ..- attr(*, "label")= chr "COMPUTED ASTHMA STATUS"
 $ _DRDXAR3: num [1:438693] 1 1 2 2 2 2 1 1 2 1 ...
  ..- attr(*, "label")= chr "RESPONDENTS DIAGNOSED WITH ARTHRITIS"
 $ _LMTACT3: num [1:438693] 2 1 3 3 3 3 1 2 3 2 ...
  ..- attr(*, "label")= chr "LIMITED USUAL ACTIVITIES"
 $ _LMTWRK3: num [1:438693] 1 1 3 3 3 3 2 2 3 2 ...
  ..- attr(*, "label")= chr "LIMITED WORK ACTIVITIES"
 $ _PRACE1 : num [1:438693] 1 2 2 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "COMPUTED PREFERRED RACE"
 $ _MRACE1 : num [1:438693] 1 2 2 1 7 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "CALCULATED NON-HISPANIC RACE INCLUDING M"
 $ _HISPANC: num [1:438693] 2 2 2 2 2 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "HISPANIC, LATINO/A, OR SPANISH ORIGIN CA"
 $ _RACE   : num [1:438693] 1 2 2 1 7 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "COMPUTED RACE-ETHNICITY GROUPING"
 $ _RACEG21: num [1:438693] 1 2 2 1 2 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "COMPUTED NON-HISPANIC WHITES/ALL OTHERS"
 $ _RACEGR3: num [1:438693] 1 2 2 1 4 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "COMPUTED FIVE LEVEL RACE/ETHNICITY CATEG"
 $ _RACEPRV: num [1:438693] 1 2 2 1 7 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "COMPUTED RACE GROUPS USED FOR INTERNET P"
 $ _SEX    : num [1:438693] 2 2 2 2 1 1 1 1 2 2 ...
  ..- attr(*, "label")= chr "CALCULATED SEX VARIABLE"
 $ _AGEG5YR: num [1:438693] 11 10 11 9 12 13 9 9 12 10 ...
  ..- attr(*, "label")= chr "REPORTED AGE IN FIVE-YEAR AGE CATEGORIES"
 $ _AGE65YR: num [1:438693] 2 2 2 1 2 2 1 1 2 2 ...
  ..- attr(*, "label")= chr "REPORTED AGE IN TWO AGE GROUPS CALCULATE"
 $ _AGE80  : num [1:438693] 70 67 72 62 76 80 63 62 78 65 ...
  ..- attr(*, "label")= chr "IMPUTED AGE VALUE COLLAPSED ABOVE 80"
 $ _AGE_G  : num [1:438693] 6 6 6 5 6 6 5 5 6 6 ...
  ..- attr(*, "label")= chr "IMPUTED AGE IN SIX GROUPS"
 $ HTIN4   : num [1:438693] 59 66 65 64 71 75 72 69 63 65 ...
  ..- attr(*, "label")= chr "COMPUTED HEIGHT IN INCHES"
 $ HTM4    : num [1:438693] 150 168 165 163 180 191 183 175 160 165 ...
  ..- attr(*, "label")= chr "COMPUTED HEIGHT IN METERS"
 $ WTKG3   : num [1:438693] 3266 NA 7711 8845 9344 ...
  ..- attr(*, "label")= chr "COMPUTED WEIGHT IN KILOGRAMS"
 $ _BMI5   : num [1:438693] 1454 NA 2829 3347 2873 ...
  ..- attr(*, "label")= chr "COMPUTED BODY MASS INDEX"
 $ _BMI5CAT: num [1:438693] 1 NA 3 4 3 2 4 2 NA 4 ...
  ..- attr(*, "label")= chr "COMPUTED BODY MASS INDEX CATEGORIES"
 $ _RFBMI5 : num [1:438693] 1 9 2 2 2 1 2 1 9 2 ...
  ..- attr(*, "label")= chr "OVERWEIGHT OR OBESE CALCULATED VARIABLE"
 $ _CHLDCNT: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "COMPUTED NUMBER OF CHILDREN IN HOUSEHOLD"
 $ _EDUCAG : num [1:438693] 2 4 2 2 1 3 4 1 1 2 ...
  ..- attr(*, "label")= chr "COMPUTED LEVEL OF EDUCATION COMPLETED CA"
 $ _INCOMG1: num [1:438693] 3 9 2 5 2 4 9 9 9 5 ...
  ..- attr(*, "label")= chr "COMPUTED INCOME CATEGORIES"
 $ _SMOKER3: num [1:438693] 3 4 4 4 4 3 4 3 4 3 ...
  ..- attr(*, "label")= chr "COMPUTED SMOKING STATUS"
 $ _RFSMOK3: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "CURRENT SMOKING CALCULATED VARIABLE"
 $ _CURECI1: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "CURRENT E-CIGARETTE USER CALCULATED VARI"
 $ DRNKANY5: num [1:438693] 2 2 2 1 2 2 2 2 2 1 ...
  ..- attr(*, "label")= chr "DRINK ANY ALCOHOLIC BEVERAGES IN PAST 30"
 $ DROCDY3_: num [1:438693] 0 0 0 14 0 0 0 0 0 10 ...
  ..- attr(*, "label")= chr "COMPUTED DRINK-OCCASIONS-PER-DAY"
 $ _RFBING5: num [1:438693] 1 1 1 2 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "BINGE DRINKING CALCULATED VARIABLE"
 $ _DRNKWK1: num [1:438693] 0 0 0 300 0 0 0 0 0 140 ...
  ..- attr(*, "label")= chr "COMPUTED NUMBER OF DRINKS OF ALCOHOL BEV"
 $ _RFDRHV7: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "HEAVY ALCOHOL CONSUMPTION  CALCULATED VA"
 $ _FLSHOT7: num [1:438693] 1 2 2 NA 1 1 NA NA 1 2 ...
  ..- attr(*, "label")= chr "FLU SHOT CALCULATED VARIABLE"
 $ _PNEUMO3: num [1:438693] 1 2 2 NA 1 1 NA NA 2 1 ...
  ..- attr(*, "label")= chr "PNEUMONIA VACCINATION CALCULATED VARIABL"
 $ _AIDTST4: num [1:438693] 2 2 2 2 1 2 2 2 2 2 ...
  ..- attr(*, "label")= chr "EVER BEEN TESTED FOR HIV CALCULATED VARI"
 $ FTJUDA2_: num [1:438693] 0 0 0 71 0 0 3 0 100 0 ...
  ..- attr(*, "label")= chr "COMPUTED FRUIT JUICE INTAKE IN TIMES PER"
 $ FRUTDA2_: num [1:438693] 100 100 100 43 100 29 40 57 29 100 ...
  ..- attr(*, "label")= chr "COMPUTED FRUIT INTAKE IN TIMES PER DAY"
 $ GRENDA1_: num [1:438693] 57 14 0 10 100 14 43 29 0 29 ...
  ..- attr(*, "label")= chr "COMPUTED DARK GREEN VEGETABLE INTAKE IN"
 $ FRNCHDA_: num [1:438693] 43 0 14 57 29 0 43 0 NA 29 ...
  ..- attr(*, "label")= chr "FRENCH FRY INTAKE IN TIMES PER DAY"
 $ POTADA1_: num [1:438693] 14 14 14 27 29 14 7 57 0 43 ...
  ..- attr(*, "label")= chr "COMPUTED POTATO SERVINGS PER DAY"
 $ VEGEDA2_: num [1:438693] 100 100 43 71 100 14 29 57 29 57 ...
  ..- attr(*, "label")= chr "COMPUTED OTHER VEGETABLE INTAKE IN TIMES"
 $ _MISFRT1: num [1:438693] 0 0 0 0 0 0 0 0 0 0 ...
  ..- attr(*, "label")= chr "THE NUMBER OF MISSING FRUIT RESPONSES"
 $ _MISVEG1: num [1:438693] 0 0 0 0 0 0 0 0 1 0 ...
  ..- attr(*, "label")= chr "THE NUMBER OF MISSING VEGETABLE RESPONSE"
 $ _FRTRES1: num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "MISSING ANY FRUIT RESPONSES"
 $ _VEGRES1: num [1:438693] 1 1 1 1 1 1 1 1 0 1 ...
  ..- attr(*, "label")= chr "MISSING ANY VEGETABLE RESPONSES"
 $ _FRUTSU1: num [1:438693] 100 100 100 114 100 29 43 57 129 100 ...
  ..- attr(*, "label")= chr "TOTAL FRUITS CONSUMED PER DAY"
 $ _VEGESU1: num [1:438693] 214 128 71 165 258 42 122 143 NA 158 ...
  ..- attr(*, "label")= chr "TOTAL VEGETABLES CONSUMED PER DAY"
 $ _FRTLT1A: num [1:438693] 1 1 1 1 1 2 2 2 1 1 ...
  ..- attr(*, "label")= chr "CONSUME FRUIT 1 OR MORE TIMES PER DAY"
 $ _VEGLT1A: num [1:438693] 1 1 2 1 1 2 1 1 9 1 ...
  ..- attr(*, "label")= chr "CONSUME VEGETABLES 1 OR MORE TIMES PER D"
 $ _FRT16A : num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "REPORTED CONSUMING FRUIT >16/DAY"
 $ _VEG23A : num [1:438693] 1 1 1 1 1 1 1 1 1 1 ...
  ..- attr(*, "label")= chr "REPORTED CONSUMING VEGETABLES >23/DAY"
 $ _FRUITE1: num [1:438693] 0 0 0 0 0 0 0 0 0 0 ...
  ..- attr(*, "label")= chr "FRUIT EXCLUSION FROM ANALYSES"
 $ _VEGETE1: num [1:438693] 0 0 0 0 0 0 0 0 1 0 ...
  ..- attr(*, "label")= chr "VEGETABLE EXCLUSION FROM ANALYSES"

attr(brfss21, which = "labels")

NULL

names(brfss21) <- str_replace(names(brfss21),"_","")

tabyl(brfss21$SEXVAR)

 brfss21$SEXVAR      n   percent
              1 203810 0.4645846
              2 234883 0.5354154

tabyl(brfss21$EDUCAG)

 brfss21$EDUCAG      n     percent
              1  25991 0.059246443
              2 111545 0.254266651
              3 120102 0.273772319
              4 178577 0.407065989
              9   2478 0.005648597

tabyl(brfss21$RFSMOK3 )

 brfss21$RFSMOK3      n    percent
               1 359891 0.82037097
               2  53832 0.12270996
               9  24970 0.05691908

tabyl(brfss21$CHCSCNCR)

 brfss21$CHCSCNCR      n      percent valid_percent
                1  41112 9.371474e-02  0.0937151663
                2 396140 9.030005e-01  0.9030046206
                7   1151 2.623703e-03  0.0026237146
                9    288 6.564955e-04  0.0006564985
               NA      2 4.558997e-06            NA

brfss21.c <- brfss21 %>% filter(EDUCAG %in% c(1:4))

brfss21.c <-filter(brfss21.c, AGEG5YR != 14)

brfss21.c <- filter(brfss21.c, !CHCSCNCR %in% c(7:9))

brfss21.c <- filter(brfss21.c, RFSMOK3 !=9)

nrow(brfss21.c)

[1] 403589

Recode

brfss21.c$SEXVAR.female  <- brfss21.c$SEXVAR-1 # 0 = male; 1 =  female

brfss21.c$CHCSCNCR.skincancer<- car::recode(brfss21.c$CHCSCNCR, recodes = "1=1;2=0; else=NA")

brfss21.c$RFSMOK3.smoke <- car::recode(brfss21.c$RFSMOK3, recodes=  "1=0;2=1; else=NA")

brfss21.c$EDUCAG.graduate <- car::recode(brfss21.c$EDUCAG, recodes=  "1:3 = 0; 4=1")

Means

mean(brfss21.c$AGE80)

[1] 54.73399

mean(brfss21.c$SEXVAR.female)

[1] 0.5360602

mean(brfss21.c$CHCSCNCR.skincancer, na.rm=T)

[1] 0.09562946

mean(brfss21.c$EDUCAG.graduate)

[1] 0.4101475

Age

brfss21.c$AGE80SQ<- brfss21.c$AGE80*brfss21.c$AGE80


brfss21.c$AGE80SQ<- brfss21.c$AGE80*brfss21.c$AGE80

ageprobs<- brfss21.c %>%
    group_by(AGE80) %>%
    summarize(p=mean(CHCSCNCR.skincancer),n=n())
              
ageprobs$num <- ageprobs$p*(1-ageprobs$p)

ageprobs$sep <- sqrt(ageprobs$num/ageprobs$n)

ageprobs$me <- 2*ageprobs$sep
  
ggplot(data =ageprobs, aes(x = AGE80, y = p, ymin=p-me, ymax=p+me)) +
     geom_line() + ylim(0,1) + geom_ribbon(alpha=0.3,aes(color=NULL)) +
labs(title="Skin Cancer by Age", 
       subtitle="Respondents in Sample", 
       caption="Source: 2021 BRFSS Data") +
       xlab(label="Age at Survey") +
  ylab(label="Proportion of Respondents with Skin Cancer")

model0 <- glm(CHCSCNCR.skincancer ~ AGE80, family = "binomial", data = brfss21.c)
summary(model0)


Call:
glm(formula = CHCSCNCR.skincancer ~ AGE80, family = "binomial", 
    data = brfss21.c)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-0.8070  -0.5076  -0.3096  -0.1554   3.3220  

Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept) -6.8371771  0.0333722  -204.9   <2e-16 ***
AGE80        0.0735282  0.0004825   152.4   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 254562  on 403588  degrees of freedom
Residual deviance: 219754  on 403587  degrees of freedom
AIC: 219758

Number of Fisher Scoring iterations: 6

model0$coefficients = exp(model0$coefficients)


summary(model0)


Call:
glm(formula = CHCSCNCR.skincancer ~ AGE80, family = "binomial", 
    data = brfss21.c)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-0.8070  -0.5076  -0.3096  -0.1554   3.3220  

Coefficients:
             Estimate Std. Error  z value Pr(>|z|)    
(Intercept) 0.0010731  0.0333722    0.032    0.974    
AGE80       1.0762989  0.0004825 2230.726   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 254562  on 403588  degrees of freedom
Residual deviance: 219754  on 403587  degrees of freedom
AIC: 219758

Number of Fisher Scoring iterations: 6

#MERs<-margins(model0,variables="AGE80", at=list(AGE80=c(18,20,30,40,50,60,70,80)))

#summary(MERs)

ageprobs<- brfss21.c %>%
    group_by(EDUCAG.graduate, SEXVAR.female, AGE80) %>%
    summarize(pskincancer=mean(CHCSCNCR.skincancer),n=n())

`summarise()` has grouped output by 'EDUCAG.graduate', 'SEXVAR.female'. You can
override using the `.groups` argument.

ggplot(data =ageprobs, aes(x = AGE80, y = pskincancer, color = factor(EDUCAG.graduate), group = EDUCAG.graduate)) + geom_line() + facet_wrap(~ SEXVAR.female) + ylim(0,.5)

model1 <- glm(CHCSCNCR.skincancer ~ AGE80 + EDUCAG.graduate + SEXVAR.female, family = "binomial", data = brfss21.c)
summary(model1)


Call:
glm(formula = CHCSCNCR.skincancer ~ AGE80 + EDUCAG.graduate + 
    SEXVAR.female, family = "binomial", data = brfss21.c)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-0.9288  -0.5119  -0.3077  -0.1532   3.4174  

Coefficients:
                 Estimate Std. Error z value Pr(>|z|)    
(Intercept)     -6.997992   0.034749 -201.39   <2e-16 ***
AGE80            0.074927   0.000488  153.56   <2e-16 ***
EDUCAG.graduate  0.386418   0.011304   34.18   <2e-16 ***
SEXVAR.female   -0.187187   0.011282  -16.59   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 254562  on 403588  degrees of freedom
Residual deviance: 218253  on 403585  degrees of freedom
AIC: 218261

Number of Fisher Scoring iterations: 6

model1$coefficients = exp(model1$coefficients)
summary(model1)


Call:
glm(formula = CHCSCNCR.skincancer ~ AGE80 + EDUCAG.graduate + 
    SEXVAR.female, family = "binomial", data = brfss21.c)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-0.9288  -0.5119  -0.3077  -0.1532   3.4174  

Coefficients:
                 Estimate Std. Error  z value Pr(>|z|)    
(Intercept)     0.0009137  0.0347492    0.026    0.979    
AGE80           1.0778060  0.0004880 2208.833   <2e-16 ***
EDUCAG.graduate 1.4717004  0.0113043  130.189   <2e-16 ***
SEXVAR.female   0.8292884  0.0112820   73.505   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 254562  on 403588  degrees of freedom
Residual deviance: 218253  on 403585  degrees of freedom
AIC: 218261

Number of Fisher Scoring iterations: 6

emmeans(model1, specs = "EDUCAG.graduate",
regrid = "response")

 EDUCAG.graduate prob       SE  df asymp.LCL asymp.UCL
               0    1 2.30e-18 Inf         1         1
               1    1 2.28e-18 Inf         1         1

Results are averaged over the levels of: SEXVAR.female 
Confidence level used: 0.95

emmeans(model1, specs = "EDUCAG.graduate",
regrid = "response") |>
contrast(method = "revpairwise") |>
confint()

 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0        0 2.51e-18 Inf -4.92e-18  4.92e-18

Results are averaged over the levels of: SEXVAR.female 
Confidence level used: 0.95

#MERs: marginal effects at representative values

for (i in c(20,40,50,60,80))
  
  {
  print(paste("At age = ", i))
emmeans(model1, specs = "EDUCAG.graduate", by = "SEXVAR.female", at = list(AGE80=i),
regrid = "response") |>
contrast(method = "revpairwise") |>
confint() %>% print()
  
}

[1] "At age =  20"
SEXVAR.female = 0:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0 3.35e-10 8.93e-12 Inf  3.17e-10  3.52e-10

SEXVAR.female = 1:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0 1.46e-10 3.93e-12 Inf  1.38e-10  1.54e-10

Confidence level used: 0.95 
[1] "At age =  40"
SEXVAR.female = 0:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0        0 2.51e-18 Inf -4.92e-18  4.92e-18

SEXVAR.female = 1:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0        0 2.51e-18 Inf -4.92e-18  4.92e-18

Confidence level used: 0.95 
[1] "At age =  50"
SEXVAR.female = 0:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0        0 2.51e-18 Inf -4.92e-18  4.92e-18

SEXVAR.female = 1:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0        0 2.51e-18 Inf -4.92e-18  4.92e-18

Confidence level used: 0.95 
[1] "At age =  60"
SEXVAR.female = 0:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0        0 2.51e-18 Inf -4.92e-18  4.92e-18

SEXVAR.female = 1:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0        0 2.51e-18 Inf -4.92e-18  4.92e-18

Confidence level used: 0.95 
[1] "At age =  80"
SEXVAR.female = 0:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0        0 2.51e-18 Inf -4.92e-18  4.92e-18

SEXVAR.female = 1:
 contrast                            estimate       SE  df asymp.LCL asymp.UCL
 EDUCAG.graduate1 - EDUCAG.graduate0        0 2.51e-18 Inf -4.92e-18  4.92e-18

Confidence level used: 0.95

EDUCAG_eff.grad <- predict(model1, newdata = data.frame(EDUCAG.graduate=1, SEXVAR.female=brfss21.c$SEXVAR.female,AGE80=brfss21.c$AGE80), 
type = "response")

EDUCAG_eff.nograd <- predict(model1, newdata = data.frame(EDUCAG.graduate=0, SEXVAR.female=brfss21.c$SEXVAR.female,AGE80=brfss21.c$AGE80), 
type = "response")

mean(EDUCAG_eff.grad - EDUCAG_eff.nograd)

[1] 2.587651e-11