#load data
library(car)
## Loading required package: carData
library(stargazer)
##
## Please cite as:
## Hlavac, Marek (2018). stargazer: Well-Formatted Regression and Summary Statistics Tables.
## R package version 5.2.2. https://CRAN.R-project.org/package=stargazer
library(survey)
## Loading required package: grid
## Loading required package: Matrix
## Loading required package: survival
##
## Attaching package: 'survey'
## The following object is masked from 'package:graphics':
##
## dotchart
library(questionr)
library(dplyr)
##
## Attaching package: 'dplyr'
## The following object is masked from 'package:car':
##
## recode
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
library(foreign)
brfss2017=read.xport("~/Desktop/DataFolder/brfss.XPT ")
#recode variables of interest
brfss2017$marital<-Recode(brfss2017$MARITAL,recodes="1=1;5=0;else=NA",as.factor=T)
brfss2017$sxorient<-Recode(brfss2017$SXORIENT,recodes="1='a) Straight';2='b) Gay';else=NA")
brfss2017$genh<-Recode(brfss2017$GENHLTH,recodes="1:3=0;4:5=1;else=NA")
brfss2017$int<-interaction(brfss2017$marital,brfss2017$sxorient)
brfss2017$age<-Recode(brfss2017$X_AGE_G, recodes="1='18-24';2='25-34';3='35-44';4='45-54';5='55-64';6='65+';else=NA",as.factor=T)
brfss2017$educ<-Recode(brfss2017$EDUCA, recodes="1:2='0Prim'; 3='1somehs'; 4='2hsgrad'; 5='3somecol'; 6='4colgrad';9=NA", as.factor=T)
brfss2017$educ<-relevel(brfss2017$educ, ref='0Prim')
#create a survey design
sub<-brfss2017%>%
select(educ,age,genh,marital,int,sxorient,X_LLCPWT,X_STRWT)%>%
filter(complete.cases(.))
options(survey.lonely.psu = "adjust")
des<-svydesign(ids = ~1,strata=~X_STRWT,weights = ~X_LLCPWT,data = sub)
library(ggplot2)
cat<-svyby(formula = ~genh, by = ~marital, design = des, FUN = svymean, na.rm=T)
svychisq(~genh+marital, design = des)
##
## Pearson's X^2: Rao & Scott adjustment
##
## data: svychisq(~genh + marital, design = des)
## F = 0.2595, ndf = 1, ddf = 129120, p-value = 0.6105
cat1<-svyby(formula = ~genh, by = ~sxorient, design = des, FUN = svymean, na.rm=T)
svychisq(~genh+sxorient, design = des)
##
## Pearson's X^2: Rao & Scott adjustment
##
## data: svychisq(~genh + sxorient, design = des)
## F = 1.1194, ndf = 1, ddf = 129120, p-value = 0.29
cat2<-svyby(formula = ~genh, by = ~age, design = des, FUN = svymean, na.rm=T)
svychisq(~genh+age, design = des)
##
## Pearson's X^2: Rao & Scott adjustment
##
## data: svychisq(~genh + age, design = des)
## F = 75.031, ndf = 4.8154e+00, ddf = 6.2177e+05, p-value < 2.2e-16
cat3<-svyby(formula = ~genh, by = ~educ, design = des, FUN = svymean, na.rm=T)
svychisq(~genh+educ, design = des)
##
## Pearson's X^2: Rao & Scott adjustment
##
## data: svychisq(~genh + educ, design = des)
## F = 395.2, ndf = 3.6917e+00, ddf = 4.7667e+05, p-value < 2.2e-16
cat4<-svyby(formula = ~marital, by = ~sxorient, design = des, FUN = svymean, na.rm=T)
svychisq(~marital+sxorient, design = des)
##
## Pearson's X^2: Rao & Scott adjustment
##
## data: svychisq(~marital + sxorient, design = des)
## F = 441.94, ndf = 1, ddf = 129120, p-value < 2.2e-16
cat5<-svyby(formula = ~genh, by = ~int+educ+age, design = des, FUN = svymean, na.rm=T)
cat5
## int educ age genh
## 0.a) Straight.0Prim.18-24 0.a) Straight 0Prim 18-24 0.21512696
## 1.a) Straight.0Prim.18-24 1.a) Straight 0Prim 18-24 0.09555236
## 0.a) Straight.1somehs.18-24 0.a) Straight 1somehs 18-24 0.15376637
## 1.a) Straight.1somehs.18-24 1.a) Straight 1somehs 18-24 0.03556219
## 0.b) Gay.1somehs.18-24 0.b) Gay 1somehs 18-24 0.21650135
## 1.b) Gay.1somehs.18-24 1.b) Gay 1somehs 18-24 0.00000000
## 0.a) Straight.2hsgrad.18-24 0.a) Straight 2hsgrad 18-24 0.09347302
## 1.a) Straight.2hsgrad.18-24 1.a) Straight 2hsgrad 18-24 0.11432728
## 0.b) Gay.2hsgrad.18-24 0.b) Gay 2hsgrad 18-24 0.06166635
## 1.b) Gay.2hsgrad.18-24 1.b) Gay 2hsgrad 18-24 0.79037196
## 0.a) Straight.3somecol.18-24 0.a) Straight 3somecol 18-24 0.07588447
## 1.a) Straight.3somecol.18-24 1.a) Straight 3somecol 18-24 0.07134385
## 0.b) Gay.3somecol.18-24 0.b) Gay 3somecol 18-24 0.05712565
## 1.b) Gay.3somecol.18-24 1.b) Gay 3somecol 18-24 0.27460452
## 0.a) Straight.4colgrad.18-24 0.a) Straight 4colgrad 18-24 0.05434934
## 1.a) Straight.4colgrad.18-24 1.a) Straight 4colgrad 18-24 0.04706812
## 0.b) Gay.4colgrad.18-24 0.b) Gay 4colgrad 18-24 0.02790748
## 1.b) Gay.4colgrad.18-24 1.b) Gay 4colgrad 18-24 0.00000000
## 0.a) Straight.0Prim.25-34 0.a) Straight 0Prim 25-34 0.39772903
## 1.a) Straight.0Prim.25-34 1.a) Straight 0Prim 25-34 0.27521409
## 0.b) Gay.0Prim.25-34 0.b) Gay 0Prim 25-34 0.00000000
## 1.b) Gay.0Prim.25-34 1.b) Gay 0Prim 25-34 0.85483226
## 0.a) Straight.1somehs.25-34 0.a) Straight 1somehs 25-34 0.33527984
## 1.a) Straight.1somehs.25-34 1.a) Straight 1somehs 25-34 0.23972564
## 0.b) Gay.1somehs.25-34 0.b) Gay 1somehs 25-34 0.27751577
## 1.b) Gay.1somehs.25-34 1.b) Gay 1somehs 25-34 0.00000000
## 0.a) Straight.2hsgrad.25-34 0.a) Straight 2hsgrad 25-34 0.16643305
## 1.a) Straight.2hsgrad.25-34 1.a) Straight 2hsgrad 25-34 0.11548115
## 0.b) Gay.2hsgrad.25-34 0.b) Gay 2hsgrad 25-34 0.23597933
## 1.b) Gay.2hsgrad.25-34 1.b) Gay 2hsgrad 25-34 0.46823948
## 0.a) Straight.3somecol.25-34 0.a) Straight 3somecol 25-34 0.11469667
## 1.a) Straight.3somecol.25-34 1.a) Straight 3somecol 25-34 0.08272997
## 0.b) Gay.3somecol.25-34 0.b) Gay 3somecol 25-34 0.06743122
## 1.b) Gay.3somecol.25-34 1.b) Gay 3somecol 25-34 0.01432897
## 0.a) Straight.4colgrad.25-34 0.a) Straight 4colgrad 25-34 0.05747497
## 1.a) Straight.4colgrad.25-34 1.a) Straight 4colgrad 25-34 0.02671382
## 0.b) Gay.4colgrad.25-34 0.b) Gay 4colgrad 25-34 0.09506727
## 1.b) Gay.4colgrad.25-34 1.b) Gay 4colgrad 25-34 0.01554524
## 0.a) Straight.0Prim.35-44 0.a) Straight 0Prim 35-44 0.35665769
## 1.a) Straight.0Prim.35-44 1.a) Straight 0Prim 35-44 0.37645416
## 0.b) Gay.0Prim.35-44 0.b) Gay 0Prim 35-44 0.10737216
## 0.a) Straight.1somehs.35-44 0.a) Straight 1somehs 35-44 0.36798594
## 1.a) Straight.1somehs.35-44 1.a) Straight 1somehs 35-44 0.27089875
## 0.b) Gay.1somehs.35-44 0.b) Gay 1somehs 35-44 0.00000000
## 1.b) Gay.1somehs.35-44 1.b) Gay 1somehs 35-44 0.42004439
## 0.a) Straight.2hsgrad.35-44 0.a) Straight 2hsgrad 35-44 0.21257575
## 1.a) Straight.2hsgrad.35-44 1.a) Straight 2hsgrad 35-44 0.13198914
## 0.b) Gay.2hsgrad.35-44 0.b) Gay 2hsgrad 35-44 0.08487684
## 1.b) Gay.2hsgrad.35-44 1.b) Gay 2hsgrad 35-44 0.07524950
## 0.a) Straight.3somecol.35-44 0.a) Straight 3somecol 35-44 0.20603275
## 1.a) Straight.3somecol.35-44 1.a) Straight 3somecol 35-44 0.12440275
## 0.b) Gay.3somecol.35-44 0.b) Gay 3somecol 35-44 0.16274581
## 1.b) Gay.3somecol.35-44 1.b) Gay 3somecol 35-44 0.40492372
## 0.a) Straight.4colgrad.35-44 0.a) Straight 4colgrad 35-44 0.10571111
## 1.a) Straight.4colgrad.35-44 1.a) Straight 4colgrad 35-44 0.05095177
## 0.b) Gay.4colgrad.35-44 0.b) Gay 4colgrad 35-44 0.13116690
## 1.b) Gay.4colgrad.35-44 1.b) Gay 4colgrad 35-44 0.05740687
## 0.a) Straight.0Prim.45-54 0.a) Straight 0Prim 45-54 0.49811856
## 1.a) Straight.0Prim.45-54 1.a) Straight 0Prim 45-54 0.49795271
## 0.b) Gay.0Prim.45-54 0.b) Gay 0Prim 45-54 0.38024134
## 0.a) Straight.1somehs.45-54 0.a) Straight 1somehs 45-54 0.35200492
## 1.a) Straight.1somehs.45-54 1.a) Straight 1somehs 45-54 0.38599364
## 0.b) Gay.1somehs.45-54 0.b) Gay 1somehs 45-54 0.30267806
## 1.b) Gay.1somehs.45-54 1.b) Gay 1somehs 45-54 0.39005177
## 0.a) Straight.2hsgrad.45-54 0.a) Straight 2hsgrad 45-54 0.30367155
## 1.a) Straight.2hsgrad.45-54 1.a) Straight 2hsgrad 45-54 0.20526277
## 0.b) Gay.2hsgrad.45-54 0.b) Gay 2hsgrad 45-54 0.62169472
## 1.b) Gay.2hsgrad.45-54 1.b) Gay 2hsgrad 45-54 0.11554980
## 0.a) Straight.3somecol.45-54 0.a) Straight 3somecol 45-54 0.24096674
## 1.a) Straight.3somecol.45-54 1.a) Straight 3somecol 45-54 0.13441386
## 0.b) Gay.3somecol.45-54 0.b) Gay 3somecol 45-54 0.23756717
## 1.b) Gay.3somecol.45-54 1.b) Gay 3somecol 45-54 0.06231890
## 0.a) Straight.4colgrad.45-54 0.a) Straight 4colgrad 45-54 0.14108931
## 1.a) Straight.4colgrad.45-54 1.a) Straight 4colgrad 45-54 0.05167209
## 0.b) Gay.4colgrad.45-54 0.b) Gay 4colgrad 45-54 0.25960130
## 1.b) Gay.4colgrad.45-54 1.b) Gay 4colgrad 45-54 0.03620881
## 0.a) Straight.0Prim.55-64 0.a) Straight 0Prim 55-64 0.64215166
## 1.a) Straight.0Prim.55-64 1.a) Straight 0Prim 55-64 0.64181598
## 0.b) Gay.0Prim.55-64 0.b) Gay 0Prim 55-64 0.09936658
## 0.a) Straight.1somehs.55-64 0.a) Straight 1somehs 55-64 0.48549688
## 1.a) Straight.1somehs.55-64 1.a) Straight 1somehs 55-64 0.37995245
## 0.b) Gay.1somehs.55-64 0.b) Gay 1somehs 55-64 0.28722740
## 1.b) Gay.1somehs.55-64 1.b) Gay 1somehs 55-64 1.00000000
## 0.a) Straight.2hsgrad.55-64 0.a) Straight 2hsgrad 55-64 0.36958103
## 1.a) Straight.2hsgrad.55-64 1.a) Straight 2hsgrad 55-64 0.23132547
## 0.b) Gay.2hsgrad.55-64 0.b) Gay 2hsgrad 55-64 0.33459130
## 1.b) Gay.2hsgrad.55-64 1.b) Gay 2hsgrad 55-64 0.15316561
## 0.a) Straight.3somecol.55-64 0.a) Straight 3somecol 55-64 0.27079094
## 1.a) Straight.3somecol.55-64 1.a) Straight 3somecol 55-64 0.15619267
## 0.b) Gay.3somecol.55-64 0.b) Gay 3somecol 55-64 0.11381011
## 1.b) Gay.3somecol.55-64 1.b) Gay 3somecol 55-64 0.08472192
## 0.a) Straight.4colgrad.55-64 0.a) Straight 4colgrad 55-64 0.12807383
## 1.a) Straight.4colgrad.55-64 1.a) Straight 4colgrad 55-64 0.07930944
## 0.b) Gay.4colgrad.55-64 0.b) Gay 4colgrad 55-64 0.22098379
## 1.b) Gay.4colgrad.55-64 1.b) Gay 4colgrad 55-64 0.12742092
## 0.a) Straight.0Prim.65+ 0.a) Straight 0Prim 65+ 0.50284074
## 1.a) Straight.0Prim.65+ 1.a) Straight 0Prim 65+ 0.56193429
## 0.b) Gay.0Prim.65+ 0.b) Gay 0Prim 65+ 0.56358122
## 1.b) Gay.0Prim.65+ 1.b) Gay 0Prim 65+ 0.00000000
## 0.a) Straight.1somehs.65+ 0.a) Straight 1somehs 65+ 0.31767947
## 1.a) Straight.1somehs.65+ 1.a) Straight 1somehs 65+ 0.36062406
## 0.b) Gay.1somehs.65+ 0.b) Gay 1somehs 65+ 0.64345732
## 1.b) Gay.1somehs.65+ 1.b) Gay 1somehs 65+ 0.57875537
## 0.a) Straight.2hsgrad.65+ 0.a) Straight 2hsgrad 65+ 0.30123949
## 1.a) Straight.2hsgrad.65+ 1.a) Straight 2hsgrad 65+ 0.24881887
## 0.b) Gay.2hsgrad.65+ 0.b) Gay 2hsgrad 65+ 0.19967970
## 1.b) Gay.2hsgrad.65+ 1.b) Gay 2hsgrad 65+ 0.48066807
## 0.a) Straight.3somecol.65+ 0.a) Straight 3somecol 65+ 0.33210326
## 1.a) Straight.3somecol.65+ 1.a) Straight 3somecol 65+ 0.17517737
## 0.b) Gay.3somecol.65+ 0.b) Gay 3somecol 65+ 0.26927958
## 1.b) Gay.3somecol.65+ 1.b) Gay 3somecol 65+ 0.36337454
## 0.a) Straight.4colgrad.65+ 0.a) Straight 4colgrad 65+ 0.15489159
## 1.a) Straight.4colgrad.65+ 1.a) Straight 4colgrad 65+ 0.11480154
## 0.b) Gay.4colgrad.65+ 0.b) Gay 4colgrad 65+ 0.20526149
## 1.b) Gay.4colgrad.65+ 1.b) Gay 4colgrad 65+ 0.07262925
## se
## 0.a) Straight.0Prim.18-24 0.087956642
## 1.a) Straight.0Prim.18-24 0.059456759
## 0.a) Straight.1somehs.18-24 0.027827682
## 1.a) Straight.1somehs.18-24 0.020681740
## 0.b) Gay.1somehs.18-24 0.127589027
## 1.b) Gay.1somehs.18-24 0.000000000
## 0.a) Straight.2hsgrad.18-24 0.009389558
## 1.a) Straight.2hsgrad.18-24 0.026646138
## 0.b) Gay.2hsgrad.18-24 0.021275247
## 1.b) Gay.2hsgrad.18-24 0.186975103
## 0.a) Straight.3somecol.18-24 0.008319055
## 1.a) Straight.3somecol.18-24 0.025751889
## 0.b) Gay.3somecol.18-24 0.023723630
## 1.b) Gay.3somecol.18-24 0.254054458
## 0.a) Straight.4colgrad.18-24 0.012001485
## 1.a) Straight.4colgrad.18-24 0.034097275
## 0.b) Gay.4colgrad.18-24 0.024787243
## 1.b) Gay.4colgrad.18-24 0.000000000
## 0.a) Straight.0Prim.25-34 0.101932150
## 1.a) Straight.0Prim.25-34 0.053916625
## 0.b) Gay.0Prim.25-34 0.000000000
## 1.b) Gay.0Prim.25-34 0.175495512
## 0.a) Straight.1somehs.25-34 0.042186596
## 1.a) Straight.1somehs.25-34 0.037760577
## 0.b) Gay.1somehs.25-34 0.208135589
## 1.b) Gay.1somehs.25-34 0.000000000
## 0.a) Straight.2hsgrad.25-34 0.014799731
## 1.a) Straight.2hsgrad.25-34 0.013320137
## 0.b) Gay.2hsgrad.25-34 0.087516213
## 1.b) Gay.2hsgrad.25-34 0.258122160
## 0.a) Straight.3somecol.25-34 0.012899018
## 1.a) Straight.3somecol.25-34 0.010209014
## 0.b) Gay.3somecol.25-34 0.024635408
## 1.b) Gay.3somecol.25-34 0.015045410
## 0.a) Straight.4colgrad.25-34 0.008369864
## 1.a) Straight.4colgrad.25-34 0.003861933
## 0.b) Gay.4colgrad.25-34 0.032607099
## 1.b) Gay.4colgrad.25-34 0.016203626
## 0.a) Straight.0Prim.35-44 0.070302242
## 1.a) Straight.0Prim.35-44 0.048799831
## 0.b) Gay.0Prim.35-44 0.135543010
## 0.a) Straight.1somehs.35-44 0.044592790
## 1.a) Straight.1somehs.35-44 0.041868474
## 0.b) Gay.1somehs.35-44 0.000000000
## 1.b) Gay.1somehs.35-44 0.283932661
## 0.a) Straight.2hsgrad.35-44 0.022319070
## 1.a) Straight.2hsgrad.35-44 0.012296821
## 0.b) Gay.2hsgrad.35-44 0.046933644
## 1.b) Gay.2hsgrad.35-44 0.073211947
## 0.a) Straight.3somecol.35-44 0.028544810
## 1.a) Straight.3somecol.35-44 0.010797559
## 0.b) Gay.3somecol.35-44 0.069591951
## 1.b) Gay.3somecol.35-44 0.183990072
## 0.a) Straight.4colgrad.35-44 0.013351305
## 1.a) Straight.4colgrad.35-44 0.006273928
## 0.b) Gay.4colgrad.35-44 0.092463061
## 1.b) Gay.4colgrad.35-44 0.027762286
## 0.a) Straight.0Prim.45-54 0.083532338
## 1.a) Straight.0Prim.45-54 0.044254741
## 0.b) Gay.0Prim.45-54 0.333270548
## 0.a) Straight.1somehs.45-54 0.045715544
## 1.a) Straight.1somehs.45-54 0.036465446
## 0.b) Gay.1somehs.45-54 0.202224841
## 1.b) Gay.1somehs.45-54 0.336457509
## 0.a) Straight.2hsgrad.45-54 0.025628325
## 1.a) Straight.2hsgrad.45-54 0.012908920
## 0.b) Gay.2hsgrad.45-54 0.120865785
## 1.b) Gay.2hsgrad.45-54 0.067636251
## 0.a) Straight.3somecol.45-54 0.028373985
## 1.a) Straight.3somecol.45-54 0.009461858
## 0.b) Gay.3somecol.45-54 0.078155803
## 1.b) Gay.3somecol.45-54 0.032933947
## 0.a) Straight.4colgrad.45-54 0.023167476
## 1.a) Straight.4colgrad.45-54 0.003926968
## 0.b) Gay.4colgrad.45-54 0.075277479
## 1.b) Gay.4colgrad.45-54 0.017291696
## 0.a) Straight.0Prim.55-64 0.083588898
## 1.a) Straight.0Prim.55-64 0.042185444
## 0.b) Gay.0Prim.55-64 0.114823857
## 0.a) Straight.1somehs.55-64 0.055474935
## 1.a) Straight.1somehs.55-64 0.034419100
## 0.b) Gay.1somehs.55-64 0.181233422
## 1.b) Gay.1somehs.55-64 0.000000000
## 0.a) Straight.2hsgrad.55-64 0.030008238
## 1.a) Straight.2hsgrad.55-64 0.012291214
## 0.b) Gay.2hsgrad.55-64 0.117562234
## 1.b) Gay.2hsgrad.55-64 0.093152193
## 0.a) Straight.3somecol.55-64 0.032278627
## 1.a) Straight.3somecol.55-64 0.007989706
## 0.b) Gay.3somecol.55-64 0.043036307
## 1.b) Gay.3somecol.55-64 0.049575337
## 0.a) Straight.4colgrad.55-64 0.014648743
## 1.a) Straight.4colgrad.55-64 0.004991557
## 0.b) Gay.4colgrad.55-64 0.061464397
## 1.b) Gay.4colgrad.55-64 0.051552418
## 0.a) Straight.0Prim.65+ 0.125442033
## 1.a) Straight.0Prim.65+ 0.046607227
## 0.b) Gay.0Prim.65+ 0.342854143
## 1.b) Gay.0Prim.65+ 0.000000000
## 0.a) Straight.1somehs.65+ 0.064238814
## 1.a) Straight.1somehs.65+ 0.025892000
## 0.b) Gay.1somehs.65+ 0.304505757
## 1.b) Gay.1somehs.65+ 0.339816760
## 0.a) Straight.2hsgrad.65+ 0.041384575
## 1.a) Straight.2hsgrad.65+ 0.009593694
## 0.b) Gay.2hsgrad.65+ 0.096701777
## 1.b) Gay.2hsgrad.65+ 0.200237599
## 0.a) Straight.3somecol.65+ 0.045486752
## 1.a) Straight.3somecol.65+ 0.008096467
## 0.b) Gay.3somecol.65+ 0.100182732
## 1.b) Gay.3somecol.65+ 0.137241362
## 0.a) Straight.4colgrad.65+ 0.021423877
## 1.a) Straight.4colgrad.65+ 0.005924703
## 0.b) Gay.4colgrad.65+ 0.066478915
## 1.b) Gay.4colgrad.65+ 0.030044430