How I made the “Blows” variable
physical_abuse == “Yes” playing_sports == “Yes” other_injury == “Yes”
military_exposures == “Yes” If any of these conditions is true, the
overall condition evaluates to true, and the value of Blows is set to
“Yes”.
PTH <- PTH %>%
mutate(
TBIUTAH = case_when(TBI == "No TBI"&Blows=="No" ~ "No TBI", TBI == "No TBI"&Blows=="Yes" ~ "Sub-concussive",unconsciousness == "No"&TBI=="TBI"&dazed=="No"~"Mild TBI without LOC",longest_unconsciousness == "Less than 5 minutes"|longest_unconsciousness == "5-30 minutes"&longest_dazed_new=="1 to 7 days" ~ "Mild",
TRUE ~ "Moderate/Severe"
)
)
PTH <- PTH %>%
mutate(
TBIOSU = case_when(TBI == "No TBI"&Blows=="No" ~ "No TBI", TBI == "No TBI"&Blows=="Yes" ~ "Sub-concussive",unconsciousness == "No"&TBI=="TBI"~"Mild TBI without LOC",longest_unconsciousness == "Less than 5 minutes"|longest_unconsciousness == "5-30 minutes" ~ "Mild",
TRUE ~ "Moderate/Severe"
)
)
table(PTH$TBIUTAH)
##
## Mild Mild TBI without LOC Moderate/Severe
## 187 217 265
## No TBI Sub-concussive
## 233 99
table(PTH$TBIOSU)
##
## Mild Mild TBI without LOC Moderate/Severe
## 296 295 78
## No TBI Sub-concussive
## 233 99
PTH %>%
select(TBIUTAH, TBIOSU,longest_unconsciousness, longest_dazed_new, TBI ) %>% head()
## TBIUTAH TBIOSU longest_unconsciousness
## 1 No TBI No TBI No
## 2 Moderate/Severe Mild 5-30 minutes
## 3 Moderate/Severe Mild TBI without LOC No
## 4 No TBI No TBI No
## 5 No TBI No TBI No
## 6 Sub-concussive Sub-concussive No
## longest_dazed_new TBI
## 1 No No TBI
## 2 More TBI
## 3 More TBI
## 4 No No TBI
## 5 No No TBI
## 6 No No TBI
Descriptive
Stratified by TBI status
## Warning in compareGroups.fit(X = X, y = y, include.label = include.label, :
## variables waist specified in 'method' not found
## Warning: glm.fit: algorithm did not converge
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
##
## --------Summary descriptives table by 'TBI'---------
##
## ___________________________________________________________
## No TBI TBI p.overall
## N=332 N=669
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯
## age 49.4 (17.9) 47.7 (21.0) 0.167
## Race: 0.037
## Black 39 (11.7%) 77 (11.5%)
## Hispanic 16 (4.82%) 64 (9.57%)
## Mixed 24 (7.23%) 33 (4.93%)
## White 253 (76.2%) 495 (74.0%)
## education: 0.051
## High School 90 (27.1%) 146 (21.8%)
## Less than High School 24 (7.23%) 49 (7.32%)
## More than College 111 (33.4%) 201 (30.0%)
## Some College 107 (32.2%) 273 (40.8%)
## HIT 57.2 (9.60) 61.1 (8.78) <0.001
## pain_level 5.04 (2.54) 5.83 (2.40) <0.001
## sex: <0.001
## Female 213 (64.2%) 279 (41.7%)
## Male 119 (35.8%) 390 (58.3%)
## PHQ 1.79 (1.95) 2.64 (1.87) <0.001
## Cognetive 0.93 (1.03) 1.52 (1.13) <0.001
## Vestibular 0.90 (1.00) 1.50 (1.11) <0.001
## Alcohol: <0.001
## Heavy 57 (17.2%) 181 (27.1%)
## Moderate 135 (40.7%) 308 (46.0%)
## Non-drinker 140 (42.2%) 180 (26.9%)
## TBInumbers 0.00 (0.00) 4.72 (4.42) <0.001
## GAD 1.89 (1.96) 2.62 (1.94) <0.001
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯
Stratified by TBI severity
##
## --------Summary descriptives table by 'TBIOSU'---------
##
## _______________________________________________________________________________________________________________
## Mild Mild TBI without LOC Moderate/Severe No TBI Sub-concussive p.overall
## N=296 N=295 N=78 N=233 N=99
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯
## age 48.2 (22.9) 46.7 (18.1) 49.0 (23.2) 50.2 (17.0) 47.7 (19.7) 0.412
## Race: .
## Black 36 (12.2%) 31 (10.5%) 10 (12.8%) 26 (11.2%) 13 (13.1%)
## Hispanic 32 (10.8%) 24 (8.14%) 8 (10.3%) 11 (4.72%) 5 (5.05%)
## Mixed 13 (4.39%) 18 (6.10%) 2 (2.56%) 14 (6.01%) 10 (10.1%)
## White 215 (72.6%) 222 (75.3%) 58 (74.4%) 182 (78.1%) 71 (71.7%)
## education: 0.061
## High School 58 (19.6%) 73 (24.7%) 15 (19.2%) 71 (30.5%) 19 (19.2%)
## Less than High School 22 (7.43%) 22 (7.46%) 5 (6.41%) 12 (5.15%) 12 (12.1%)
## More than College 91 (30.7%) 88 (29.8%) 22 (28.2%) 79 (33.9%) 32 (32.3%)
## Some College 125 (42.2%) 112 (38.0%) 36 (46.2%) 71 (30.5%) 36 (36.4%)
## HIT 61.9 (8.75) 59.4 (8.55) 64.7 (8.35) 56.3 (9.20) 59.4 (10.2) <0.001
## sex: <0.001
## Female 111 (37.5%) 150 (50.8%) 18 (23.1%) 158 (67.8%) 55 (55.6%)
## Male 185 (62.5%) 145 (49.2%) 60 (76.9%) 75 (32.2%) 44 (44.4%)
## PHQ 2.83 (1.77) 2.30 (1.91) 3.23 (1.87) 1.50 (1.88) 2.47 (1.96) <0.001
## Cognetive 1.76 (1.11) 1.12 (0.98) 2.12 (1.19) 0.74 (0.86) 1.37 (1.24) <0.001
## Vestibular 1.67 (1.07) 1.15 (0.98) 2.16 (1.29) 0.69 (0.84) 1.39 (1.14) <0.001
## Alcohol: <0.001
## Heavy 87 (29.4%) 59 (20.0%) 35 (44.9%) 31 (13.3%) 26 (26.3%)
## Moderate 148 (50.0%) 133 (45.1%) 27 (34.6%) 92 (39.5%) 43 (43.4%)
## Non-drinker 61 (20.6%) 103 (34.9%) 16 (20.5%) 110 (47.2%) 30 (30.3%)
## TBInumbers 5.56 (4.22) 3.07 (3.25) 7.79 (6.16) 0.00 (0.00) 0.00 (0.00) <0.001
## GAD 2.89 (1.84) 2.23 (1.98) 3.08 (1.90) 1.68 (1.93) 2.36 (1.97) <0.001
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯
fit1=lm(HIT~TBIOSU+TBInumbers+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
fit3=lm(Cognetive~TBIOSU+TBInumbers+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
fit4=lm(Vestibular~TBIOSU+TBInumbers+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
fit5=lm(PHQ~TBIOSU+TBInumbers+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
tab_model(fit1, fit3, fit4, fit5, show.reflvl = T)
|
|
HIT
|
Cognetive
|
Vestibular
|
PHQ
|
|
Predictors
|
Estimates
|
CI
|
p
|
Estimates
|
CI
|
p
|
Estimates
|
CI
|
p
|
Estimates
|
CI
|
p
|
|
(Intercept)
|
53.36
|
51.20 – 55.53
|
<0.001
|
0.19
|
-0.04 – 0.42
|
0.113
|
0.15
|
-0.09 – 0.38
|
0.229
|
0.45
|
0.14 – 0.77
|
0.005
|
|
age
|
-0.03
|
-0.05 – -0.00
|
0.021
|
-0.00
|
-0.00 – 0.00
|
0.342
|
-0.00
|
-0.00 – 0.00
|
0.893
|
-0.00
|
-0.01 – 0.00
|
0.054
|
|
No TBI
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Mild
|
1.11
|
-0.60 – 2.82
|
0.205
|
0.18
|
-0.00 – 0.37
|
0.051
|
0.10
|
-0.08 – 0.29
|
0.279
|
0.05
|
-0.20 – 0.30
|
0.706
|
|
Mild TBI without LOC
|
0.54
|
-0.99 – 2.07
|
0.489
|
-0.07
|
-0.23 – 0.10
|
0.416
|
-0.04
|
-0.21 – 0.13
|
0.664
|
0.15
|
-0.08 – 0.37
|
0.200
|
|
BlowsYes
|
1.11
|
-0.27 – 2.50
|
0.115
|
0.15
|
0.00 – 0.30
|
0.047
|
0.25
|
0.10 – 0.40
|
0.001
|
0.10
|
-0.10 – 0.30
|
0.343
|
|
Moderate/Severe
|
2.86
|
0.44 – 5.28
|
0.020
|
0.32
|
0.06 – 0.58
|
0.017
|
0.34
|
0.08 – 0.61
|
0.012
|
0.15
|
-0.21 – 0.50
|
0.413
|
|
Sub-concussive
|
0.91
|
-1.36 – 3.18
|
0.431
|
0.23
|
-0.02 – 0.47
|
0.069
|
0.22
|
-0.03 – 0.47
|
0.081
|
0.31
|
-0.02 – 0.64
|
0.068
|
|
Male
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Female
|
2.29
|
1.20 – 3.37
|
<0.001
|
0.01
|
-0.11 – 0.13
|
0.887
|
-0.01
|
-0.13 – 0.11
|
0.921
|
-0.08
|
-0.24 – 0.08
|
0.321
|
|
GAD
|
1.88
|
1.62 – 2.15
|
<0.001
|
0.27
|
0.24 – 0.30
|
<0.001
|
0.21
|
0.18 – 0.24
|
<0.001
|
0.74
|
0.71 – 0.78
|
<0.001
|
|
White
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Black
|
-2.38
|
-3.94 – -0.83
|
0.003
|
0.01
|
-0.16 – 0.17
|
0.949
|
-0.03
|
-0.20 – 0.14
|
0.750
|
0.01
|
-0.22 – 0.24
|
0.929
|
|
Hispanic
|
-1.12
|
-2.96 – 0.71
|
0.231
|
-0.21
|
-0.41 – -0.01
|
0.036
|
0.01
|
-0.19 – 0.21
|
0.947
|
-0.05
|
-0.32 – 0.22
|
0.717
|
|
Mixed
|
-1.99
|
-4.09 – 0.11
|
0.063
|
-0.26
|
-0.49 – -0.03
|
0.024
|
-0.16
|
-0.39 – 0.07
|
0.178
|
-0.07
|
-0.38 – 0.23
|
0.642
|
|
High School
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Less than High School
|
-0.57
|
-2.65 – 1.51
|
0.590
|
0.18
|
-0.05 – 0.40
|
0.123
|
0.31
|
0.08 – 0.54
|
0.009
|
0.31
|
0.01 – 0.61
|
0.044
|
|
TBInumbers
|
0.35
|
0.19 – 0.50
|
<0.001
|
0.06
|
0.04 – 0.07
|
<0.001
|
0.06
|
0.04 – 0.08
|
<0.001
|
0.05
|
0.02 – 0.07
|
<0.001
|
|
More than College
|
-0.24
|
-1.56 – 1.07
|
0.718
|
0.06
|
-0.08 – 0.21
|
0.375
|
0.05
|
-0.10 – 0.19
|
0.502
|
-0.03
|
-0.22 – 0.16
|
0.746
|
|
Some College
|
-0.02
|
-1.35 – 1.30
|
0.973
|
0.14
|
0.00 – 0.29
|
0.049
|
0.22
|
0.08 – 0.37
|
0.003
|
0.08
|
-0.12 – 0.27
|
0.441
|
|
Non-drinker
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Heavy
|
0.89
|
-0.52 – 2.30
|
0.217
|
0.36
|
0.21 – 0.52
|
<0.001
|
0.43
|
0.27 – 0.58
|
<0.001
|
0.16
|
-0.05 – 0.36
|
0.130
|
|
Moderate
|
0.36
|
-0.83 – 1.54
|
0.551
|
0.16
|
0.04 – 0.29
|
0.011
|
0.16
|
0.03 – 0.29
|
0.016
|
-0.04
|
-0.21 – 0.13
|
0.654
|
|
Observations
|
1001
|
1001
|
1001
|
1001
|
|
R2 / R2 adjusted
|
0.313 / 0.301
|
0.471 / 0.462
|
0.426 / 0.416
|
0.670 / 0.664
|
fit1=lm(HIT~TBIUTAH+TBInumbers+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
fit3=lm(Cognetive~TBIUTAH+TBInumbers+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
fit4=lm(Vestibular~TBIUTAH+TBInumbers+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
fit5=lm(PHQ~TBIUTAH+TBInumbers+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
tab_model(fit1, fit3, fit4, fit5, show.reflvl = T)
|
|
HIT
|
Cognetive
|
Vestibular
|
PHQ
|
|
Predictors
|
Estimates
|
CI
|
p
|
Estimates
|
CI
|
p
|
Estimates
|
CI
|
p
|
Estimates
|
CI
|
p
|
|
(Intercept)
|
53.45
|
51.29 – 55.61
|
<0.001
|
0.20
|
-0.03 – 0.43
|
0.094
|
0.16
|
-0.08 – 0.39
|
0.197
|
0.45
|
0.14 – 0.77
|
0.005
|
|
age
|
-0.03
|
-0.06 – -0.01
|
0.015
|
-0.00
|
-0.00 – 0.00
|
0.229
|
-0.00
|
-0.00 – 0.00
|
0.761
|
-0.00
|
-0.01 – 0.00
|
0.059
|
|
No TBI
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Mild
|
0.59
|
-1.25 – 2.44
|
0.528
|
0.09
|
-0.11 – 0.28
|
0.394
|
-0.00
|
-0.20 – 0.20
|
0.987
|
0.01
|
-0.26 – 0.28
|
0.930
|
|
Mild TBI without LOC
|
0.07
|
-1.53 – 1.67
|
0.930
|
-0.16
|
-0.33 – 0.01
|
0.063
|
-0.10
|
-0.27 – 0.08
|
0.284
|
0.18
|
-0.05 – 0.41
|
0.130
|
|
BlowsYes
|
1.10
|
-0.29 – 2.48
|
0.121
|
0.14
|
-0.00 – 0.29
|
0.055
|
0.25
|
0.10 – 0.40
|
0.001
|
0.10
|
-0.10 – 0.30
|
0.321
|
|
Moderate/Severe
|
2.01
|
0.29 – 3.72
|
0.022
|
0.28
|
0.10 – 0.47
|
0.003
|
0.21
|
0.02 – 0.40
|
0.031
|
0.06
|
-0.19 – 0.31
|
0.623
|
|
Sub-concussive
|
0.90
|
-1.36 – 3.17
|
0.433
|
0.23
|
-0.01 – 0.47
|
0.064
|
0.22
|
-0.03 – 0.47
|
0.085
|
0.30
|
-0.03 – 0.63
|
0.074
|
|
Male
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Female
|
2.28
|
1.19 – 3.36
|
<0.001
|
0.01
|
-0.10 – 0.13
|
0.841
|
-0.01
|
-0.13 – 0.11
|
0.872
|
-0.09
|
-0.24 – 0.07
|
0.282
|
|
GAD
|
1.89
|
1.63 – 2.16
|
<0.001
|
0.27
|
0.24 – 0.30
|
<0.001
|
0.21
|
0.18 – 0.24
|
<0.001
|
0.75
|
0.71 – 0.78
|
<0.001
|
|
White
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Black
|
-2.34
|
-3.90 – -0.79
|
0.003
|
0.01
|
-0.15 – 0.18
|
0.866
|
-0.02
|
-0.19 – 0.15
|
0.805
|
0.01
|
-0.22 – 0.23
|
0.941
|
|
Hispanic
|
-1.17
|
-3.00 – 0.67
|
0.212
|
-0.21
|
-0.41 – -0.01
|
0.035
|
0.00
|
-0.20 – 0.20
|
0.999
|
-0.06
|
-0.32 – 0.21
|
0.674
|
|
Mixed
|
-2.15
|
-4.24 – -0.05
|
0.045
|
-0.29
|
-0.52 – -0.07
|
0.011
|
-0.18
|
-0.41 – 0.05
|
0.117
|
-0.07
|
-0.37 – 0.24
|
0.661
|
|
High School
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Less than High School
|
-0.59
|
-2.66 – 1.49
|
0.579
|
0.17
|
-0.05 – 0.39
|
0.130
|
0.30
|
0.07 – 0.53
|
0.009
|
0.31
|
0.01 – 0.62
|
0.042
|
|
TBInumbers
|
0.36
|
0.21 – 0.51
|
<0.001
|
0.06
|
0.04 – 0.07
|
<0.001
|
0.06
|
0.04 – 0.08
|
<0.001
|
0.05
|
0.03 – 0.07
|
<0.001
|
|
More than College
|
-0.28
|
-1.60 – 1.03
|
0.671
|
0.06
|
-0.08 – 0.20
|
0.421
|
0.04
|
-0.10 – 0.19
|
0.544
|
-0.03
|
-0.22 – 0.16
|
0.755
|
|
Some College
|
-0.05
|
-1.38 – 1.28
|
0.940
|
0.14
|
0.00 – 0.28
|
0.050
|
0.22
|
0.07 – 0.37
|
0.003
|
0.07
|
-0.12 – 0.27
|
0.458
|
|
Non-drinker
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
Reference
|
|
|
|
Heavy
|
1.02
|
-0.38 – 2.43
|
0.153
|
0.39
|
0.24 – 0.54
|
<0.001
|
0.45
|
0.30 – 0.60
|
<0.001
|
0.16
|
-0.05 – 0.36
|
0.127
|
|
Moderate
|
0.37
|
-0.81 – 1.56
|
0.535
|
0.17
|
0.05 – 0.30
|
0.007
|
0.17
|
0.04 – 0.29
|
0.013
|
-0.04
|
-0.21 – 0.13
|
0.640
|
|
Observations
|
1001
|
1001
|
1001
|
1001
|
|
R2 / R2 adjusted
|
0.315 / 0.303
|
0.479 / 0.470
|
0.429 / 0.419
|
0.670 / 0.664
|
Chisq test
result <- chisq.test(table(PTH$TBIUTAH, PTH$TBIOSU))
result
##
## Pearson's Chi-squared test
##
## data: table(PTH$TBIUTAH, PTH$TBIOSU)
## X-squared = 2893.9, df = 16, p-value < 2.2e-16