library(car)## Loading required package: carDatalibrary(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=stargazerlibrary(survey)## Loading required package: grid## Loading required package: Matrix## Loading required package: survival##
## Attaching package: 'survey'## The following object is masked from 'package:graphics':
##
## dotchartlibrary(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, unionlibrary(tidyverse)## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --## v ggplot2 3.3.5 v purrr 0.3.4
## v tibble 3.1.6 v stringr 1.4.0
## v tidyr 1.1.4 v forcats 0.5.1
## v readr 2.1.1## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x tidyr::expand() masks Matrix::expand()
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()
## x tidyr::pack() masks Matrix::pack()
## x dplyr::recode() masks car::recode()
## x purrr::some() masks car::some()
## x tidyr::unpack() masks Matrix::unpack()library(broom)
library(emmeans)
library(table1)##
## Attaching package: 'table1'## The following objects are masked from 'package:base':
##
## units, units<-library(ggplot2)
library(mice)## Warning: package 'mice' was built under R version 4.1.3##
## Attaching package: 'mice'## The following object is masked from 'package:stats':
##
## filter## The following objects are masked from 'package:base':
##
## cbind, rbindlibrary(boot)## Warning: package 'boot' was built under R version 4.1.3##
## Attaching package: 'boot'## The following object is masked from 'package:survival':
##
## aml## The following object is masked from 'package:car':
##
## logitlibrary(tableone)library(haven)
NSDUH_2019 <- read_sav("NSDUH_2019.SAV")
View(NSDUH_2019)nams<-names(NSDUH_2019)
head(nams, n=10)## [1] "QUESTID2" "FILEDATE" "CIGEVER" "CIGOFRSM" "CIGWILYR" "CIGTRY"
## [7] "CIGYFU" "CIGMFU" "CIGREC" "CIG30USE"newnames<-tolower(gsub(pattern = "_",replacement = "",x = nams))
names(NSDUH_2019)<-newnames## attempted suicide
NSDUH_2019$attempt_suicide<-Recode(NSDUH_2019$adwrsatp, recodes="1=1; 2=0;else=NA", as.factor=T)
summary(NSDUH_2019$attempt_suicide, na.rm = TRUE)## 0 1 NA's
## 2650 887 52599## education recodes
NSDUH_2019$educ<-Recode(NSDUH_2019$ireduhighst2, recodes="1:7='LssThnHgh'; 8='highschool'; 9:10='someCollege';11='colgrad';else=NA", as.factor=T)
NSDUH_2019$educ<-relevel(NSDUH_2019$educ, ref='colgrad')
## gender recodes
NSDUH_2019$male<-as.factor(ifelse(NSDUH_2019$irsex==1, "Male", "Female"))
## Race recoded items
NSDUH_2019$black<-Recode(NSDUH_2019$newrace2, recodes="2=1; 9=NA; else=0")
NSDUH_2019$white<-Recode(NSDUH_2019$newrace2, recodes="1=1; 9=NA; else=0")
NSDUH_2019$other<-Recode(NSDUH_2019$newrace2, recodes="3:4=1; 5:6=1; 9=NA; else=0")
NSDUH_2019$hispanic<-Recode(NSDUH_2019$newrace2, recodes="7=1; 9=NA; else=0")
NSDUH_2019$race_eth<-Recode(NSDUH_2019$newrace2,
recodes="1='white'; 2='black'; 3:6='other'; 7='hispanic'; else=NA",
as.factor = T)
NSDUH_2019$served2001<-Recode(NSDUH_2019$actd2001, recodes="1=1; 2=0;else=NA", as.factor=T)
summary(NSDUH_2019$served2001, na.rm = TRUE)## 0 1 NA's
## 1122 677 54337NSDUH_2019$service2<-Recode(NSDUH_2019$service, recodes="1=1; 2=0;else=NA", as.factor=T)
summary(NSDUH_2019$service2, na.rm = TRUE)## 0 1 NA's
## 42608 2336 11192NSDUH_2019$prior2001<-Recode(NSDUH_2019$actd9001, recodes="1=1; 2=0;else=NA", as.factor=T)
summary(NSDUH_2019$prior2001, na.rm = TRUE)## 0 1 NA's
## 1364 435 54337NSDUH_2019$impsocm2<-Recode(NSDUH_2019$impsocm, recodes="1=1; 2=0;else=NA", as.factor=T)
summary(NSDUH_2019$impsocm2, na.rm = TRUE)## 0 1 NA's
## 644 1079 54413NSDUH_2019$imfam<-Recode(NSDUH_2019$imprespm, recodes="1=1; 2=0;else=NA", as.factor=T)
summary(NSDUH_2019$imfam, na.rm = TRUE)## 0 1 NA's
## 1120 411 54605NSDUH_2019$alcserpb2<-Recode(NSDUH_2019$alcserpb, recodes="1=1; 2=0;else=NA", as.factor=T)
summary(NSDUH_2019$alcserpb2, na.rm = TRUE)## 0 1 NA's
## 27180 789 28167NSDUH_2019$alcfmfpb2<-Recode(NSDUH_2019$alcfmfpb, recodes="1=1; 2=0;else=NA", as.factor=T)
summary(NSDUH_2019$alcfmfpb2, na.rm = TRUE)## 0 1 NA's
## 26627 1350 28159## sexuality recodes
NSDUH_2019$sexuality<-Recode(NSDUH_2019$sexident, recodes="1='Heterosexual'; 2='Les/Gay'; 3='Bisexual';else=NA", as.factor=T)
NSDUH_2019$sexuality<-relevel(NSDUH_2019$sexuality, ref='Heterosexual')#First we tell R our survey design
options(survey.lonely.psu = "adjust")
des<-svydesign(ids= ~1,
strata= ~vestr,
weights= ~analwtc
, data = NSDUH_2019 )
imp<-mice(data= NSDUH_2019[,c( "male", "sexuality", "race_eth", "educ", "served2001", "prior2001", "impsocm2", "imfam", "alcserpb2", "alcfmfpb2", "attempt_suicide")], m=5)##
## iter imp variable
## 1 1 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 1 2 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 1 3 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 1 4 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 1 5 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 2 1 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 2 2 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 2 3 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 2 4 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 2 5 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 3 1 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 3 2 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 3 3 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 3 4 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 3 5 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 4 1 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 4 2 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 4 3 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 4 4 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 4 5 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 5 1 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 5 2 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 5 3 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 5 4 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicide
## 5 5 sexuality served2001 prior2001 impsocm2 imfam alcserpb2 alcfmfpb2 attempt_suicidedat.imp<-complete(imp, action = 1)
head(dat.imp, n=10)## male sexuality race_eth educ served2001 prior2001 impsocm2 imfam
## 1 Male Heterosexual white colgrad 1 0 1 0
## 2 Female Heterosexual other colgrad 0 1 1 0
## 3 Male Heterosexual white colgrad 1 0 1 0
## 4 Female Heterosexual white LssThnHgh 0 0 1 0
## 5 Male Heterosexual white LssThnHgh 0 0 0 0
## 6 Male Heterosexual black someCollege 1 0 0 0
## 7 Female Heterosexual white someCollege 0 1 1 0
## 8 Male Heterosexual white someCollege 0 0 0 0
## 9 Female Heterosexual other colgrad 0 1 1 0
## 10 Female Heterosexual black colgrad 0 1 1 0
## alcserpb2 alcfmfpb2 attempt_suicide
## 1 0 0 0
## 2 0 0 0
## 3 0 0 0
## 4 0 0 0
## 5 0 0 1
## 6 0 0 0
## 7 0 0 1
## 8 0 0 1
## 9 0 0 0
## 10 0 0 1library(table1)
t1<-CreateTableOne(vars = c("male", "sexuality", "race_eth", "educ", "served2001", "prior2001", "impsocm2", "imfam", "alcserpb2", "alcfmfpb2"), strata = "attempt_suicide", test = T, data = dat.imp)
#st1<-print(st1, format="p")
print(t1, format="p")## Stratified by attempt_suicide
## 0 1 p test
## n 41066 15070
## male = Male (%) 48.9 44.6 <0.001
## sexuality (%) <0.001
## Heterosexual 92.1 88.8
## Bisexual 5.7 8.8
## Les/Gay 2.2 2.4
## race_eth (%) <0.001
## black 13.0 12.8
## hispanic 19.0 20.2
## other 9.1 14.8
## white 59.0 52.2
## educ (%) <0.001
## colgrad 23.8 12.5
## highschool 20.9 19.9
## LssThnHgh 28.7 44.5
## someCollege 26.7 23.1
## served2001 = 1 (%) 39.7 28.9 <0.001
## prior2001 = 1 (%) 24.3 23.4 0.017
## impsocm2 = 1 (%) 58.7 52.4 <0.001
## imfam = 1 (%) 26.0 31.9 <0.001
## alcserpb2 = 1 (%) 3.2 4.3 <0.001
## alcfmfpb2 = 1 (%) 5.3 6.4 <0.001table1(~ male + sexuality + educ + race_eth + served2001 + prior2001 + alcserpb2 + alcfmfpb2 + impsocm2 + imfam| attempt_suicide, data=dat.imp, overall="Total")| 0 (N=41066) |
1 (N=15070) |
Total (N=56136) |
|
|---|---|---|---|
| male | |||
| Female | 20998 (51.1%) | 8350 (55.4%) | 29348 (52.3%) |
| Male | 20068 (48.9%) | 6720 (44.6%) | 26788 (47.7%) |
| sexuality | |||
| Heterosexual | 37809 (92.1%) | 13383 (88.8%) | 51192 (91.2%) |
| Bisexual | 2335 (5.7%) | 1324 (8.8%) | 3659 (6.5%) |
| Les/Gay | 922 (2.2%) | 363 (2.4%) | 1285 (2.3%) |
| educ | |||
| colgrad | 9756 (23.8%) | 1880 (12.5%) | 11636 (20.7%) |
| highschool | 8580 (20.9%) | 2997 (19.9%) | 11577 (20.6%) |
| LssThnHgh | 11768 (28.7%) | 6707 (44.5%) | 18475 (32.9%) |
| someCollege | 10962 (26.7%) | 3486 (23.1%) | 14448 (25.7%) |
| race_eth | |||
| black | 5325 (13.0%) | 1931 (12.8%) | 7256 (12.9%) |
| hispanic | 7797 (19.0%) | 3051 (20.2%) | 10848 (19.3%) |
| other | 3718 (9.1%) | 2225 (14.8%) | 5943 (10.6%) |
| white | 24226 (59.0%) | 7863 (52.2%) | 32089 (57.2%) |
| served2001 | |||
| 0 | 24770 (60.3%) | 10716 (71.1%) | 35486 (63.2%) |
| 1 | 16296 (39.7%) | 4354 (28.9%) | 20650 (36.8%) |
| prior2001 | |||
| 0 | 31076 (75.7%) | 11551 (76.6%) | 42627 (75.9%) |
| 1 | 9990 (24.3%) | 3519 (23.4%) | 13509 (24.1%) |
| alcserpb2 | |||
| 0 | 39758 (96.8%) | 14423 (95.7%) | 54181 (96.5%) |
| 1 | 1308 (3.2%) | 647 (4.3%) | 1955 (3.5%) |
| alcfmfpb2 | |||
| 0 | 38907 (94.7%) | 14105 (93.6%) | 53012 (94.4%) |
| 1 | 2159 (5.3%) | 965 (6.4%) | 3124 (5.6%) |
| impsocm2 | |||
| 0 | 16972 (41.3%) | 7168 (47.6%) | 24140 (43.0%) |
| 1 | 24094 (58.7%) | 7902 (52.4%) | 31996 (57.0%) |
| imfam | |||
| 0 | 30408 (74.0%) | 10257 (68.1%) | 40665 (72.4%) |
| 1 | 10658 (26.0%) | 4813 (31.9%) | 15471 (27.6%) |
glm1<-glm(attempt_suicide~factor(served2001),
data=dat.imp,
family = binomial)
summary(glm1)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001), family = binomial,
## data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.8479 -0.8479 -0.6882 1.5475 1.7644
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.83790 0.01156 -72.47 <2e-16 ***
## factor(served2001)1 -0.48193 0.02061 -23.38 <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: 65310 on 56135 degrees of freedom
## Residual deviance: 64745 on 56134 degrees of freedom
## AIC: 64749
##
## Number of Fisher Scoring iterations: 4glm2<-glm(attempt_suicide~factor(served2001)+factor(race_eth),
data=dat.imp,
family = binomial)
summary(glm2)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.0593 -0.8026 -0.7089 1.3003 1.8431
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.82987 0.02754 -30.134 < 2e-16 ***
## factor(served2001)1 -0.52879 0.02087 -25.340 < 2e-16 ***
## factor(race_eth)hispanic 0.10573 0.03432 3.081 0.00206 **
## factor(race_eth)other 0.54548 0.03805 14.334 < 2e-16 ***
## factor(race_eth)white -0.13783 0.02976 -4.632 3.62e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 64227 on 56131 degrees of freedom
## AIC: 64237
##
## Number of Fisher Scoring iterations: 4glm3<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male),
data=dat.imp,
family = binomial)
summary(glm3)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male), family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.1060 -0.8399 -0.7378 1.3460 1.9000
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.72533 0.02904 -24.979 < 2e-16 ***
## factor(served2001)1 -0.54707 0.02096 -26.100 < 2e-16 ***
## factor(race_eth)hispanic 0.11036 0.03436 3.212 0.00132 **
## factor(race_eth)other 0.55493 0.03811 14.560 < 2e-16 ***
## factor(race_eth)white -0.13532 0.02979 -4.543 5.56e-06 ***
## factor(male)Male -0.21764 0.01940 -11.220 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 64100 on 56130 degrees of freedom
## AIC: 64112
##
## Number of Fisher Scoring iterations: 4glm4<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male)+factor(educ),
data=dat.imp,
family = binomial)
summary(glm4)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male) + factor(educ), family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.2674 -0.8408 -0.6842 1.2912 2.1264
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.46777 0.03890 -37.731 <2e-16 ***
## factor(served2001)1 -0.40670 0.02192 -18.554 <2e-16 ***
## factor(race_eth)hispanic 0.04468 0.03476 1.285 0.199
## factor(race_eth)other 0.63926 0.03877 16.488 <2e-16 ***
## factor(race_eth)white -0.02331 0.03024 -0.771 0.441
## factor(male)Male -0.25282 0.01968 -12.848 <2e-16 ***
## factor(educ)highschool 0.63606 0.03348 19.000 <2e-16 ***
## factor(educ)LssThnHgh 1.03759 0.03051 34.008 <2e-16 ***
## factor(educ)someCollege 0.56504 0.03227 17.510 <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: 65310 on 56135 degrees of freedom
## Residual deviance: 62821 on 56127 degrees of freedom
## AIC: 62839
##
## Number of Fisher Scoring iterations: 4glm5<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality),
data=dat.imp,
family = binomial)
summary(glm5)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality), family = binomial,
## data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.4855 -0.8275 -0.6546 1.2039 2.1508
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.48242 0.03900 -38.010 < 2e-16 ***
## factor(served2001)1 -0.47301 0.02256 -20.969 < 2e-16 ***
## factor(race_eth)hispanic 0.04911 0.03483 1.410 0.158539
## factor(race_eth)other 0.63418 0.03886 16.319 < 2e-16 ***
## factor(race_eth)white -0.03148 0.03031 -1.039 0.298965
## factor(male)Male -0.22177 0.01986 -11.167 < 2e-16 ***
## factor(educ)highschool 0.61808 0.03355 18.423 < 2e-16 ***
## factor(educ)LssThnHgh 1.00785 0.03060 32.931 < 2e-16 ***
## factor(educ)someCollege 0.55536 0.03233 17.177 < 2e-16 ***
## factor(sexuality)Bisexual 0.54075 0.03804 14.216 < 2e-16 ***
## factor(sexuality)Les/Gay 0.22615 0.06473 3.494 0.000476 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 62619 on 56125 degrees of freedom
## AIC: 62641
##
## Number of Fisher Scoring iterations: 4glm6<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(alcserpb2),
data=dat.imp,
family = binomial)
summary(glm6)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(alcserpb2),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.4799 -0.8267 -0.6538 1.2082 2.1520
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.48284 0.03900 -38.022 < 2e-16 ***
## factor(served2001)1 -0.47417 0.02256 -21.016 < 2e-16 ***
## factor(race_eth)hispanic 0.04733 0.03484 1.359 0.174286
## factor(race_eth)other 0.63171 0.03888 16.250 < 2e-16 ***
## factor(race_eth)white -0.03088 0.03031 -1.019 0.308290
## factor(male)Male -0.22376 0.01988 -11.258 < 2e-16 ***
## factor(educ)highschool 0.61574 0.03356 18.347 < 2e-16 ***
## factor(educ)LssThnHgh 1.00278 0.03066 32.708 < 2e-16 ***
## factor(educ)someCollege 0.55350 0.03234 17.116 < 2e-16 ***
## factor(sexuality)Bisexual 0.53615 0.03808 14.079 < 2e-16 ***
## factor(sexuality)Les/Gay 0.22226 0.06476 3.432 0.000599 ***
## factor(alcserpb2)1 0.14209 0.05059 2.809 0.004971 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 62611 on 56124 degrees of freedom
## AIC: 62635
##
## Number of Fisher Scoring iterations: 4glm7<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(alcfmfpb2),
data=dat.imp,
family = binomial)
summary(glm7)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(alcfmfpb2),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.4781 -0.8264 -0.6488 1.2104 2.1605
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.47797 0.03901 -37.891 < 2e-16 ***
## factor(served2001)1 -0.49189 0.02286 -21.522 < 2e-16 ***
## factor(race_eth)hispanic 0.04633 0.03484 1.330 0.183585
## factor(race_eth)other 0.63119 0.03888 16.235 < 2e-16 ***
## factor(race_eth)white -0.03652 0.03033 -1.204 0.228610
## factor(male)Male -0.22560 0.01988 -11.349 < 2e-16 ***
## factor(educ)highschool 0.61542 0.03356 18.339 < 2e-16 ***
## factor(educ)LssThnHgh 0.99513 0.03069 32.421 < 2e-16 ***
## factor(educ)someCollege 0.55515 0.03234 17.168 < 2e-16 ***
## factor(sexuality)Bisexual 0.53539 0.03807 14.062 < 2e-16 ***
## factor(sexuality)Les/Gay 0.22098 0.06479 3.411 0.000648 ***
## factor(alcfmfpb2)1 0.22731 0.04187 5.428 5.69e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 62590 on 56124 degrees of freedom
## AIC: 62614
##
## Number of Fisher Scoring iterations: 4glm8<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(alcserpb2)+factor(alcfmfpb2),
data=dat.imp,
family = binomial)
summary(glm8)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(alcserpb2) +
## factor(alcfmfpb2), family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.4770 -0.8262 -0.6489 1.2111 2.1603
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.47830 0.03901 -37.897 < 2e-16 ***
## factor(served2001)1 -0.49115 0.02288 -21.465 < 2e-16 ***
## factor(race_eth)hispanic 0.04601 0.03484 1.321 0.186643
## factor(race_eth)other 0.63071 0.03889 16.220 < 2e-16 ***
## factor(race_eth)white -0.03608 0.03034 -1.189 0.234328
## factor(male)Male -0.22592 0.01989 -11.361 < 2e-16 ***
## factor(educ)highschool 0.61494 0.03357 18.321 < 2e-16 ***
## factor(educ)LssThnHgh 0.99448 0.03071 32.384 < 2e-16 ***
## factor(educ)someCollege 0.55466 0.03235 17.148 < 2e-16 ***
## factor(sexuality)Bisexual 0.53448 0.03810 14.029 < 2e-16 ***
## factor(sexuality)Les/Gay 0.22024 0.06480 3.399 0.000677 ***
## factor(alcserpb2)1 0.03711 0.05549 0.669 0.503658
## factor(alcfmfpb2)1 0.21477 0.04591 4.678 2.9e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 62589 on 56123 degrees of freedom
## AIC: 62615
##
## Number of Fisher Scoring iterations: 4glm9<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(impsocm2),
data=dat.imp,
family = binomial)
summary(glm9)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(impsocm2),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.5614 -0.8280 -0.6525 1.2302 2.1997
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.360816 0.040377 -33.703 < 2e-16 ***
## factor(served2001)1 -0.460318 0.022610 -20.359 < 2e-16 ***
## factor(race_eth)hispanic 0.084249 0.035007 2.407 0.0161 *
## factor(race_eth)other 0.664843 0.039016 17.040 < 2e-16 ***
## factor(race_eth)white -0.004807 0.030439 -0.158 0.8745
## factor(male)Male -0.264794 0.020253 -13.074 < 2e-16 ***
## factor(educ)highschool 0.617658 0.033581 18.393 < 2e-16 ***
## factor(educ)LssThnHgh 0.984214 0.030691 32.069 < 2e-16 ***
## factor(educ)someCollege 0.572793 0.032402 17.678 < 2e-16 ***
## factor(sexuality)Bisexual 0.580335 0.038236 15.178 < 2e-16 ***
## factor(sexuality)Les/Gay 0.271024 0.064889 4.177 2.96e-05 ***
## factor(impsocm2)1 -0.235376 0.020555 -11.451 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 62488 on 56124 degrees of freedom
## AIC: 62512
##
## Number of Fisher Scoring iterations: 4glm10<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(imfam),
data=dat.imp,
family = binomial)
summary(glm10)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(imfam),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.6032 -0.8336 -0.6553 1.2233 2.2081
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.55110 0.03933 -39.436 < 2e-16 ***
## factor(served2001)1 -0.55062 0.02313 -23.805 < 2e-16 ***
## factor(race_eth)hispanic 0.04008 0.03493 1.147 0.251
## factor(race_eth)other 0.69408 0.03915 17.730 < 2e-16 ***
## factor(race_eth)white -0.01614 0.03041 -0.531 0.596
## factor(male)Male -0.22849 0.01992 -11.471 < 2e-16 ***
## factor(educ)highschool 0.59201 0.03366 17.590 < 2e-16 ***
## factor(educ)LssThnHgh 0.97817 0.03071 31.854 < 2e-16 ***
## factor(educ)someCollege 0.53286 0.03243 16.431 < 2e-16 ***
## factor(sexuality)Bisexual 0.45586 0.03848 11.847 < 2e-16 ***
## factor(sexuality)Les/Gay 0.28140 0.06503 4.327 1.51e-05 ***
## factor(imfam)1 0.38438 0.02238 17.174 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 62327 on 56124 degrees of freedom
## AIC: 62351
##
## Number of Fisher Scoring iterations: 4glm11<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(impsocm2)+factor(imfam),
data=dat.imp,
family = binomial)
summary(glm11)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(impsocm2) +
## factor(imfam), family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.6070 -0.8283 -0.6374 1.1835 2.3319
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.35368 0.04051 -33.418 < 2e-16 ***
## factor(served2001)1 -0.56568 0.02324 -24.344 < 2e-16 ***
## factor(race_eth)hispanic 0.10235 0.03519 2.909 0.00363 **
## factor(race_eth)other 0.78378 0.03959 19.799 < 2e-16 ***
## factor(race_eth)white 0.04238 0.03065 1.383 0.16667
## factor(male)Male -0.31370 0.02046 -15.334 < 2e-16 ***
## factor(educ)highschool 0.57940 0.03377 17.157 < 2e-16 ***
## factor(educ)LssThnHgh 0.91899 0.03091 29.728 < 2e-16 ***
## factor(educ)someCollege 0.55565 0.03257 17.059 < 2e-16 ***
## factor(sexuality)Bisexual 0.48639 0.03865 12.585 < 2e-16 ***
## factor(sexuality)Les/Gay 0.39656 0.06555 6.049 1.45e-09 ***
## factor(impsocm2)1 -0.45998 0.02294 -20.050 < 2e-16 ***
## factor(imfam)1 0.59430 0.02495 23.818 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 61921 on 56123 degrees of freedom
## AIC: 61947
##
## Number of Fisher Scoring iterations: 4glm12<-glm(attempt_suicide~factor(served2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+ factor(alcserpb2)+factor(alcfmfpb2)+factor(impsocm2)+factor(imfam),
data=dat.imp,
family = binomial)
summary(glm12)##
## Call:
## glm(formula = attempt_suicide ~ factor(served2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(alcserpb2) +
## factor(alcfmfpb2) + factor(impsocm2) + factor(imfam), family = binomial,
## data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.5994 -0.8307 -0.6345 1.1730 2.3379
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.34846 0.04052 -33.277 < 2e-16 ***
## factor(served2001)1 -0.57836 0.02351 -24.600 < 2e-16 ***
## factor(race_eth)hispanic 0.10060 0.03520 2.858 0.004264 **
## factor(race_eth)other 0.78062 0.03962 19.705 < 2e-16 ***
## factor(race_eth)white 0.03834 0.03067 1.250 0.211349
## factor(male)Male -0.31651 0.02048 -15.453 < 2e-16 ***
## factor(educ)highschool 0.57817 0.03378 17.116 < 2e-16 ***
## factor(educ)LssThnHgh 0.91030 0.03101 29.358 < 2e-16 ***
## factor(educ)someCollege 0.55618 0.03258 17.070 < 2e-16 ***
## factor(sexuality)Bisexual 0.48458 0.03868 12.529 < 2e-16 ***
## factor(sexuality)Les/Gay 0.39211 0.06562 5.976 2.29e-09 ***
## factor(alcserpb2)1 -0.01171 0.05592 -0.209 0.834164
## factor(alcfmfpb2)1 0.17154 0.04626 3.708 0.000209 ***
## factor(impsocm2)1 -0.46069 0.02295 -20.077 < 2e-16 ***
## factor(imfam)1 0.58659 0.02508 23.388 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 61906 on 56121 degrees of freedom
## AIC: 61936
##
## Number of Fisher Scoring iterations: 4exp(coefficients(glm1))## (Intercept) factor(served2001)1
## 0.4326201 0.6175906exp(coefficients(glm2))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.4361045 0.5893168 1.1115268
## factor(race_eth)other factor(race_eth)white
## 1.7254404 0.8712425exp(coefficients(glm3))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.4841634 0.5786401 1.1166807
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.7418153 0.8734321 0.8044187exp(coefficients(glm4))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.2304387 0.6658465 1.0456925
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8950740 0.9769559 0.7766101
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8890273 2.8224172 1.7595267exp(coefficients(glm5))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.2270866 0.6231209 1.0503319
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8854801 0.9690088 0.8011003
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8553549 2.7397045 1.7425635
## factor(sexuality)Bisexual factor(sexuality)Les/Gay
## 1.7173026 1.2537680exp(coefficients(glm6))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.2269914 0.6223995 1.0484662
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8808264 0.9695889 0.7995068
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8510339 2.7258358 1.7393245
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(alcserpb2)1
## 1.7094205 1.2488907 1.1526817exp(coefficients(glm7))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.2281012 0.6114709 1.0474198
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8798430 0.9641420 0.7980360
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8504317 2.7050677 1.7421960
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(alcfmfpb2)1
## 1.7081150 1.2473036 1.2552129exp(coefficients(glm8))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.2280255 0.6119225 1.0470888
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8789389 0.9645626 0.7977820
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8495521 2.7033277 1.7413453
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(alcserpb2)1
## 1.7065548 1.2463742 1.0378052
## factor(alcfmfpb2)1
## 1.2395806exp(coefficients(glm9))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.2564515 0.6310827 1.0878997
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.9441854 0.9952042 0.7673644
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8545794 2.6757090 1.7732128
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(impsocm2)1
## 1.7866365 1.3113070 0.7902735exp(coefficients(glm10))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.2120147 0.5765936 1.0408938
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 2.0018728 0.9839939 0.7957354
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8076197 2.6595813 1.7037995
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(imfam)1
## 1.5775221 1.3249885 1.4687070exp(coefficients(glm11))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.2582877 0.5679750 1.1077748
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 2.1897390 1.0432960 0.7307389
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.7849608 2.5067517 1.7430773
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(impsocm2)1
## 1.6264313 1.4867027 0.6312987
## factor(imfam)1
## 1.8117555exp(coefficients(glm12))## (Intercept) factor(served2001)1 factor(race_eth)hispanic
## 0.2596388 0.5608169 1.1058330
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 2.1828217 1.0390808 0.7286894
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.7827783 2.4850627 1.7439920
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(alcserpb2)1
## 1.6234853 1.4800971 0.9883608
## factor(alcfmfpb2)1 factor(impsocm2)1 factor(imfam)1
## 1.1871370 0.6308452 1.7978562glm13<-glm(attempt_suicide~factor(prior2001),
data=dat.imp,
family = binomial)
summary(glm13)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001), family = binomial,
## data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.7951 -0.7951 -0.7951 1.6160 1.6402
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.98966 0.01090 -90.817 <2e-16 ***
## factor(prior2001)1 -0.05374 0.02243 -2.396 0.0166 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 65305 on 56134 degrees of freedom
## AIC: 65309
##
## Number of Fisher Scoring iterations: 4glm14<-glm(attempt_suicide~factor(prior2001)+factor(race_eth),
data=dat.imp,
family = binomial)
summary(glm14)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.9779 -0.7916 -0.7557 1.3909 1.7043
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.99478 0.02715 -36.647 < 2e-16 ***
## factor(prior2001)1 -0.07831 0.02258 -3.468 0.000525 ***
## factor(race_eth)hispanic 0.07327 0.03410 2.149 0.031653 *
## factor(race_eth)other 0.50561 0.03777 13.388 < 2e-16 ***
## factor(race_eth)white -0.11255 0.02957 -3.806 0.000141 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 64881 on 56131 degrees of freedom
## AIC: 64891
##
## Number of Fisher Scoring iterations: 4glm15<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male),
data=dat.imp,
family = binomial)
summary(glm15)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male), family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.0178 -0.7874 -0.7514 1.3939 1.7597
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.90217 0.02871 -31.421 < 2e-16 ***
## factor(prior2001)1 -0.10796 0.02280 -4.734 2.2e-06 ***
## factor(race_eth)hispanic 0.07528 0.03413 2.206 0.027399 *
## factor(race_eth)other 0.51433 0.03781 13.601 < 2e-16 ***
## factor(race_eth)white -0.11010 0.02960 -3.720 0.000199 ***
## factor(male)Male -0.18905 0.01939 -9.750 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 64786 on 56130 degrees of freedom
## AIC: 64798
##
## Number of Fisher Scoring iterations: 4glm16<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male)+factor(educ),
data=dat.imp,
family = binomial)
summary(glm16)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male) + factor(educ), family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.2508 -0.8201 -0.7051 1.3593 2.0248
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.694287 0.039232 -43.186 < 2e-16 ***
## factor(prior2001)1 0.109576 0.023913 4.582 4.6e-06 ***
## factor(race_eth)hispanic 0.013439 0.034598 0.388 0.698
## factor(race_eth)other 0.602316 0.038585 15.610 < 2e-16 ***
## factor(race_eth)white 0.008201 0.030130 0.272 0.785
## factor(male)Male -0.217905 0.019705 -11.059 < 2e-16 ***
## factor(educ)highschool 0.659565 0.033596 19.632 < 2e-16 ***
## factor(educ)LssThnHgh 1.153328 0.030776 37.475 < 2e-16 ***
## factor(educ)someCollege 0.529282 0.032113 16.482 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 63151 on 56127 degrees of freedom
## AIC: 63169
##
## Number of Fisher Scoring iterations: 4glm17<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality),
data=dat.imp,
family = binomial)
summary(glm17)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality), family = binomial,
## data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.4028 -0.8534 -0.7029 1.2555 2.0287
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.730067 0.039479 -43.822 < 2e-16 ***
## factor(prior2001)1 0.122870 0.023969 5.126 2.96e-07 ***
## factor(race_eth)hispanic 0.013102 0.034631 0.378 0.7052
## factor(race_eth)other 0.594366 0.038640 15.382 < 2e-16 ***
## factor(race_eth)white 0.005747 0.030171 0.190 0.8489
## factor(male)Male -0.191130 0.019911 -9.599 < 2e-16 ***
## factor(educ)highschool 0.648975 0.033640 19.292 < 2e-16 ***
## factor(educ)LssThnHgh 1.144974 0.030809 37.164 < 2e-16 ***
## factor(educ)someCollege 0.518089 0.032159 16.110 < 2e-16 ***
## factor(sexuality)Bisexual 0.383537 0.036995 10.367 < 2e-16 ***
## factor(sexuality)Les/Gay 0.138413 0.064245 2.154 0.0312 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 63044 on 56125 degrees of freedom
## AIC: 63066
##
## Number of Fisher Scoring iterations: 4glm18<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(alcserpb2),
data=dat.imp,
family = binomial)
summary(glm18)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(alcserpb2),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.3976 -0.8513 -0.7016 1.2593 2.0296
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.730495 0.039479 -43.834 < 2e-16 ***
## factor(prior2001)1 0.122365 0.023971 5.105 3.31e-07 ***
## factor(race_eth)hispanic 0.011490 0.034640 0.332 0.7401
## factor(race_eth)other 0.592192 0.038654 15.320 < 2e-16 ***
## factor(race_eth)white 0.006252 0.030173 0.207 0.8358
## factor(male)Male -0.192853 0.019926 -9.678 < 2e-16 ***
## factor(educ)highschool 0.646935 0.033652 19.224 < 2e-16 ***
## factor(educ)LssThnHgh 1.140746 0.030861 36.964 < 2e-16 ***
## factor(educ)someCollege 0.516427 0.032166 16.055 < 2e-16 ***
## factor(sexuality)Bisexual 0.379301 0.037045 10.239 < 2e-16 ***
## factor(sexuality)Les/Gay 0.134954 0.064275 2.100 0.0358 *
## factor(alcserpb2)1 0.119084 0.050365 2.364 0.0181 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 63039 on 56124 degrees of freedom
## AIC: 63063
##
## Number of Fisher Scoring iterations: 4glm19<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(alcfmfpb2),
data=dat.imp,
family = binomial)
summary(glm19)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(alcfmfpb2),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.3979 -0.8514 -0.7022 1.2580 2.0291
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.729790 0.039477 -43.817 < 2e-16 ***
## factor(prior2001)1 0.120218 0.024021 5.005 5.59e-07 ***
## factor(race_eth)hispanic 0.011727 0.034641 0.339 0.7350
## factor(race_eth)other 0.593130 0.038648 15.347 < 2e-16 ***
## factor(race_eth)white 0.004423 0.030181 0.147 0.8835
## factor(male)Male -0.192391 0.019926 -9.655 < 2e-16 ***
## factor(educ)highschool 0.647834 0.033646 19.254 < 2e-16 ***
## factor(educ)LssThnHgh 1.141576 0.030871 36.979 < 2e-16 ***
## factor(educ)someCollege 0.517455 0.032160 16.090 < 2e-16 ***
## factor(sexuality)Bisexual 0.379609 0.037070 10.240 < 2e-16 ***
## factor(sexuality)Les/Gay 0.135609 0.064279 2.110 0.0349 *
## factor(alcfmfpb2)1 0.070396 0.041148 1.711 0.0871 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 63041 on 56124 degrees of freedom
## AIC: 63065
##
## Number of Fisher Scoring iterations: 4glm20<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(alcserpb2)+factor(alcfmfpb2),
data=dat.imp,
family = binomial)
summary(glm20)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(alcserpb2) +
## factor(alcfmfpb2), family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.4033 -0.8506 -0.7014 1.2600 2.0297
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.730269 0.039479 -43.828 < 2e-16 ***
## factor(prior2001)1 0.121039 0.024025 5.038 4.7e-07 ***
## factor(race_eth)hispanic 0.011010 0.034645 0.318 0.7506
## factor(race_eth)other 0.591875 0.038656 15.311 < 2e-16 ***
## factor(race_eth)white 0.005474 0.030188 0.181 0.8561
## factor(male)Male -0.193257 0.019933 -9.696 < 2e-16 ***
## factor(educ)highschool 0.646642 0.033653 19.215 < 2e-16 ***
## factor(educ)LssThnHgh 1.139596 0.030892 36.890 < 2e-16 ***
## factor(educ)someCollege 0.516343 0.032166 16.052 < 2e-16 ***
## factor(sexuality)Bisexual 0.377877 0.037086 10.189 < 2e-16 ***
## factor(sexuality)Les/Gay 0.134002 0.064290 2.084 0.0371 *
## factor(alcserpb2)1 0.100564 0.055160 1.823 0.0683 .
## factor(alcfmfpb2)1 0.037246 0.045090 0.826 0.4088
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 65310 on 56135 degrees of freedom
## Residual deviance: 63038 on 56123 degrees of freedom
## AIC: 63064
##
## Number of Fisher Scoring iterations: 4glm21<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(impsocm2),
data=dat.imp,
family = binomial)
summary(glm21)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(impsocm2),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.4354 -0.8164 -0.6869 1.2414 2.1070
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.59842 0.04069 -39.283 < 2e-16 ***
## factor(prior2001)1 0.15190 0.02412 6.299 3.01e-10 ***
## factor(race_eth)hispanic 0.05502 0.03483 1.580 0.1142
## factor(race_eth)other 0.62877 0.03881 16.202 < 2e-16 ***
## factor(race_eth)white 0.03639 0.03031 1.200 0.2300
## factor(male)Male -0.23757 0.02027 -11.718 < 2e-16 ***
## factor(educ)highschool 0.65286 0.03369 19.380 < 2e-16 ***
## factor(educ)LssThnHgh 1.12214 0.03089 36.331 < 2e-16 ***
## factor(educ)someCollege 0.54009 0.03225 16.747 < 2e-16 ***
## factor(sexuality)Bisexual 0.43618 0.03727 11.702 < 2e-16 ***
## factor(sexuality)Les/Gay 0.19334 0.06444 3.000 0.0027 **
## factor(impsocm2)1 -0.26865 0.02056 -13.069 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 62874 on 56124 degrees of freedom
## AIC: 62898
##
## Number of Fisher Scoring iterations: 4glm22<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(imfam),
data=dat.imp,
family = binomial)
summary(glm22)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(imfam),
## family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.4581 -0.8422 -0.6829 1.2723 2.0552
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.786243 0.039827 -44.850 < 2e-16 ***
## factor(prior2001)1 0.093401 0.024123 3.872 0.000108 ***
## factor(race_eth)hispanic 0.001667 0.034695 0.048 0.961681
## factor(race_eth)other 0.633633 0.038837 16.315 < 2e-16 ***
## factor(race_eth)white 0.018690 0.030231 0.618 0.536420
## factor(male)Male -0.196777 0.019943 -9.867 < 2e-16 ***
## factor(educ)highschool 0.627351 0.033718 18.606 < 2e-16 ***
## factor(educ)LssThnHgh 1.129041 0.030855 36.592 < 2e-16 ***
## factor(educ)someCollege 0.496940 0.032236 15.416 < 2e-16 ***
## factor(sexuality)Bisexual 0.303313 0.037639 8.059 7.72e-16 ***
## factor(sexuality)Les/Gay 0.165775 0.064405 2.574 0.010055 *
## factor(imfam)1 0.266143 0.021857 12.177 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 62898 on 56124 degrees of freedom
## AIC: 62922
##
## Number of Fisher Scoring iterations: 4glm23<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+
factor(impsocm2)+factor(imfam),
data=dat.imp,
family = binomial)
summary(glm23)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(impsocm2) +
## factor(imfam), family = binomial, data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.5187 -0.8263 -0.6736 1.2213 2.2031
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.61020 0.04078 -39.487 < 2e-16 ***
## factor(prior2001)1 0.11982 0.02426 4.939 7.83e-07 ***
## factor(race_eth)hispanic 0.06209 0.03494 1.777 0.0756 .
## factor(race_eth)other 0.71809 0.03925 18.297 < 2e-16 ***
## factor(race_eth)white 0.07777 0.03048 2.552 0.0107 *
## factor(male)Male -0.27638 0.02044 -13.524 < 2e-16 ***
## factor(educ)highschool 0.61946 0.03382 18.315 < 2e-16 ***
## factor(educ)LssThnHgh 1.08113 0.03101 34.861 < 2e-16 ***
## factor(educ)someCollege 0.51868 0.03237 16.025 < 2e-16 ***
## factor(sexuality)Bisexual 0.33143 0.03780 8.767 < 2e-16 ***
## factor(sexuality)Les/Gay 0.27613 0.06488 4.256 2.08e-05 ***
## factor(impsocm2)1 -0.44789 0.02285 -19.605 < 2e-16 ***
## factor(imfam)1 0.46505 0.02432 19.123 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 62510 on 56123 degrees of freedom
## AIC: 62536
##
## Number of Fisher Scoring iterations: 4glm24<-glm(attempt_suicide~factor(prior2001)+factor(race_eth)+factor(male)+factor(educ)+factor(sexuality)+ factor(alcserpb2)+factor(alcfmfpb2)+factor(impsocm2)+factor(imfam),
data=dat.imp,
family = binomial)
summary(glm24)##
## Call:
## glm(formula = attempt_suicide ~ factor(prior2001) + factor(race_eth) +
## factor(male) + factor(educ) + factor(sexuality) + factor(alcserpb2) +
## factor(alcfmfpb2) + factor(impsocm2) + factor(imfam), family = binomial,
## data = dat.imp)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.5459 -0.8257 -0.6732 1.2236 2.2036
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.60975 0.04078 -39.471 < 2e-16 ***
## factor(prior2001)1 0.12031 0.02430 4.951 7.38e-07 ***
## factor(race_eth)hispanic 0.06144 0.03495 1.758 0.0788 .
## factor(race_eth)other 0.71656 0.03927 18.246 < 2e-16 ***
## factor(race_eth)white 0.07828 0.03050 2.567 0.0103 *
## factor(male)Male -0.27733 0.02046 -13.557 < 2e-16 ***
## factor(educ)highschool 0.61853 0.03383 18.283 < 2e-16 ***
## factor(educ)LssThnHgh 1.07908 0.03109 34.711 < 2e-16 ***
## factor(educ)someCollege 0.51794 0.03237 15.999 < 2e-16 ***
## factor(sexuality)Bisexual 0.33018 0.03784 8.725 < 2e-16 ***
## factor(sexuality)Les/Gay 0.27415 0.06492 4.223 2.41e-05 ***
## factor(alcserpb2)1 0.07839 0.05549 1.413 0.1578
## factor(alcfmfpb2)1 -0.01288 0.04545 -0.283 0.7770
## factor(impsocm2)1 -0.44841 0.02285 -19.625 < 2e-16 ***
## factor(imfam)1 0.46262 0.02449 18.892 < 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: 65310 on 56135 degrees of freedom
## Residual deviance: 62508 on 56121 degrees of freedom
## AIC: 62538
##
## Number of Fisher Scoring iterations: 4exp(coefficients(glm13))## (Intercept) factor(prior2001)1
## 0.3717016 0.9476747exp(coefficients(glm14))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.3698053 0.9246779 1.0760162
## factor(race_eth)other factor(race_eth)white
## 1.6579992 0.8935513exp(coefficients(glm15))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.4056895 0.8976669 1.0781823
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.6725231 0.8957402 0.8277443exp(coefficients(glm16))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.1837301 1.1158045 1.0135293
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8263432 1.0082350 0.8042021
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.9339516 3.1687216 1.6977130exp(coefficients(glm17))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.1772725 1.1307378 1.0131882
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8118814 1.0057631 0.8260252
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.9135785 3.1423585 1.6788164
## factor(sexuality)Bisexual factor(sexuality)Les/Gay
## 1.4674664 1.1484503exp(coefficients(glm18))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.1771966 1.1301666 1.0115559
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8079463 1.0062720 0.8246032
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.9096784 3.1291021 1.6760284
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(alcserpb2)1
## 1.4612623 1.1444846 1.1264640exp(coefficients(glm19))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.1773217 1.1277422 1.0117957
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8096442 1.0044332 0.8249840
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.9113955 3.1316993 1.6777519
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(alcfmfpb2)1
## 1.4617133 1.1452335 1.0729330exp(coefficients(glm20))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.1772367 1.1286689 1.0110705
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8073750 1.0054888 0.8242700
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.9091196 3.1255053 1.6758874
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(alcserpb2)1
## 1.4591831 1.1433953 1.1057943
## factor(alcfmfpb2)1
## 1.0379482exp(coefficients(glm21))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.2022160 1.1640467 1.0565629
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8752959 1.0370554 0.7885402
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.9210252 3.0714135 1.7161684
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(impsocm2)1
## 1.5467826 1.2132930 0.7644116exp(coefficients(glm22))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.1675887 1.0979019 1.0016683
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 1.8844439 1.0188658 0.8213735
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8726429 3.0926886 1.6436846
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(imfam)1
## 1.3543383 1.1803070 1.3049217exp(coefficients(glm23))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.1998476 1.1272950 1.0640557
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 2.0505051 1.0808741 0.7585246
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8579199 2.9480219 1.6798163
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(impsocm2)1
## 1.3929568 1.3180257 0.6389744
## factor(imfam)1
## 1.5920972exp(coefficients(glm24))## (Intercept) factor(prior2001)1 factor(race_eth)hispanic
## 0.1999373 1.1278516 1.0633693
## factor(race_eth)other factor(race_eth)white factor(male)Male
## 2.0473836 1.0814214 0.7578042
## factor(educ)highschool factor(educ)LssThnHgh factor(educ)someCollege
## 1.8562033 2.9419738 1.6785594
## factor(sexuality)Bisexual factor(sexuality)Les/Gay factor(alcserpb2)1
## 1.3912116 1.3154091 1.0815438
## factor(alcfmfpb2)1 factor(impsocm2)1 factor(imfam)1
## 0.9872071 0.6386446 1.5882290