res2 <- compareGroups(sex ~ kind_headache+ education, data = PTH)
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
createTable(res2, show.ratio = TRUE)
## 
## --------Summary descriptives table by 'sex'---------
## 
## ___________________________________________________________________________________________________________________________ 
##                                                                    Female       Male            OR        p.ratio p.overall 
##                                                                     N=492       N=512                                       
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯ 
## kind_headache:                                                                                                      0.021   
##     0Not Answer                                                  29 (5.89%)  44 (8.59%)        Ref.        Ref.             
##     Band like                                                    74 (15.0%)  101 (19.7%) 0.90 [0.51;1.57]  0.714            
##     One-sided                                                    113 (23.0%) 121 (23.6%) 0.71 [0.41;1.20]  0.204            
##     Other (specify):                                             13 (2.64%)  11 (2.15%)  0.56 [0.22;1.44]  0.229            
##     Pressing or tightening (non-pulsating) quality               115 (23.4%) 81 (15.8%)  0.47 [0.27;0.80]  0.006            
##     Pulsating                                                    148 (30.1%) 154 (30.1%) 0.69 [0.40;1.15]  0.157            
## education:                                                                                                            .     
##     College 1 year to 3 years (some college or technical school) 184 (37.4%) 129 (25.2%)       Ref.        Ref.             
##     College 4 years or more (college graduate)                   121 (24.6%) 261 (51.0%) 3.07 [2.25;4.20] <0.001            
##     Grade 12 or GED (high school graduate)                       156 (31.7%) 80 (15.6%)  0.73 [0.51;1.04]  0.081            
##     Grades 1-8 (elementary – middle school)                     5 (1.02%)  11 (2.15%)  3.08 [1.08;10.2]  0.036            
##     Grades 9-11 (some high school)                               23 (4.67%)  25 (4.88%)  1.55 [0.84;2.87]  0.162            
##     Never attended school or only attended kindergarten           3 (0.61%)   6 (1.17%)  2.77 [0.70;14.1]  0.150            
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯
compareGroups(sex ~ kind_headache+ education, data = PTH)
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## 
## 
## -------- Summary of results by groups of 'sex'---------
## 
## 
##   var           N    p.value method      selection
## 1 kind_headache 1004 0.021** categorical ALL      
## 2 education     1004 .       categorical ALL      
## -----
## Signif. codes:  0 '**' 0.05 '*' 0.1 ' ' 1

TOTAL INJURY

##       TBI crash fall fight blast
## 1  No TBI    No   No    No    No
## 2     TBI   Yes  Yes   Yes   Yes
## 3     TBI    No  Yes    No    No
## 4  No TBI    No   No    No    No
## 5  No TBI    No   No    No    No
## 6  No TBI    No   No    No    No
## 7  No TBI    No   No    No    No
## 8     TBI   Yes  Yes   Yes   Yes
## 9  No TBI    No   No    No    No
## 10    TBI   Yes  Yes    No    No
## 11 No TBI    No   No    No    No
## 12    TBI   Yes   No    No    No
## 13    TBI   Yes  Yes    No    No
## 14 No TBI    No   No    No   Yes
## 15    TBI    No  Yes   Yes    No
## 16    TBI   Yes  Yes   Yes   Yes
## 17 No TBI    No   No    No    No
## 18    TBI    No  Yes   Yes    No
## 19 No TBI    No   No    No    No
## 20    TBI   Yes  Yes   Yes   Yes
## 21    TBI   Yes  Yes   Yes   Yes
## 22    TBI   Yes  Yes   Yes   Yes
## 23    TBI   Yes  Yes   Yes   Yes
## 24 No TBI    No   No    No    No
## 25    TBI    No  Yes   Yes   Yes
## 26    TBI    No  Yes    No    No
## 27 No TBI    No   No    No    No
## 28    TBI   Yes   No    No    No
## 29    TBI   Yes  Yes   Yes   Yes
## 30    TBI    No  Yes    No    No
##     total_injury_cat total_injury crash_num fall_num fight_num blast_num
## 1          2 or more            4         1        1         1         1
## 2          2 or more            4         1        1         1         1
## 3          2 or more            2         1        1         0         0
## 4          2 or more            2         1        1         0         0
## 5          2 or more            2         0        1         1         0
## 6          2 or more            4         1        1         1         1
## 7          2 or more            2         0        1         1         0
## 8          2 or more            4         1        1         1         1
## 9          2 or more            4         1        1         1         1
## 10         2 or more            4         1        1         1         1
## 11         2 or more            4         1        1         1         1
## 12         2 or more            3         0        1         1         1
## 13         2 or more            4         1        1         1         1
## 14         2 or more            2         1        1         0         0
## 15         2 or more            2         1        1         0         0
## 16         2 or more            2         1        0         1         0
## 17         2 or more            3         1        1         1         0
## 18         2 or more            4         1        1         1         1
## 19         2 or more            2         1        1         0         0
## 20         2 or more            2         0        1         0         1
## 21         2 or more            3         1        1         1         0
## 22         2 or more            2         1        1         0         0
## 23         2 or more            3         1        1         0         1
## 24         2 or more            2         0        1         1         0
## 25         2 or more            2         0        1         1         0
## 26         2 or more            4         1        1         1         1
## 27         2 or more            2         0        0         1         1
## 28         2 or more            3         1        1         0         1
## 29         2 or more            2         1        1         0         0
## 30         2 or more            2         1        1         0         0
## 31         2 or more            2         0        1         1         0
## 32         2 or more            2         1        1         0         0
## 33         2 or more            2         0        1         1         0
## 34         2 or more            4         1        1         1         1
## 35         2 or more            2         0        1         0         1
## 36         2 or more            4         1        1         1         1
## 37         2 or more            3         1        1         1         0
## 38         2 or more            4         1        1         1         1
## 39         2 or more            4         1        1         1         1
## 40         2 or more            3         1        1         1         0
## 41         2 or more            4         1        1         1         1
## 42         2 or more            2         1        1         0         0
## 43         2 or more            3         1        1         1         0
## 44         2 or more            2         1        0         0         1
## 45         2 or more            3         1        1         1         0
## 46         2 or more            2         1        1         0         0
## 47         2 or more            3         1        1         1         0
## 48         2 or more            2         1        1         0         0
## 49         2 or more            2         1        1         0         0
## 50         2 or more            2         1        1         0         0
## 51         2 or more            2         1        0         1         0
## 52         2 or more            2         1        1         0         0
## 53         2 or more            2         1        1         0         0
## 54         2 or more            2         1        1         0         0
## 55         2 or more            3         1        1         1         0
## 56         2 or more            3         1        1         1         0
## 57         2 or more            2         0        1         1         0
## 58         2 or more            2         1        1         0         0
## 59         2 or more            4         1        1         1         1
## 60         2 or more            2         1        1         0         0
## 61         2 or more            2         0        1         1         0
## 62         2 or more            3         1        1         1         0
## 63         2 or more            3         1        1         0         1
## 64         2 or more            3         1        1         0         1
## 65         2 or more            3         1        1         0         1
## 66         2 or more            3         1        1         1         0
## 67         2 or more            4         1        1         1         1
## 68         2 or more            4         1        1         1         1
## 69         2 or more            2         1        1         0         0
## 70         2 or more            2         1        1         0         0
## 71         2 or more            3         1        1         1         0
## 72         2 or more            3         1        1         1         0
## 73         2 or more            4         1        1         1         1
## 74         2 or more            4         1        1         1         1
## 75         2 or more            4         1        1         1         1
## 76         2 or more            4         1        1         1         1
## 77         2 or more            2         1        0         0         1
## 78         2 or more            4         1        1         1         1
## 79         2 or more            3         1        1         0         1
## 80         2 or more            3         1        1         1         0
## 81         2 or more            3         1        1         0         1
## 82         2 or more            4         1        1         1         1
## 83         2 or more            3         0        1         1         1
## 84         2 or more            4         1        1         1         1
## 85         2 or more            3         1        1         0         1
## 86         2 or more            2         0        1         1         0
## 87         2 or more            4         1        1         1         1
## 88         2 or more            4         1        1         1         1
## 89         2 or more            2         0        1         0         1
## 90         2 or more            4         1        1         1         1
## 91         2 or more            4         1        1         1         1
## 92         2 or more            4         1        1         1         1
## 93         2 or more            3         1        1         1         0
## 94         2 or more            3         1        1         1         0
## 95         2 or more            3         1        1         0         1
## 96         2 or more            3         1        1         1         0
## 97         2 or more            4         1        1         1         1
## 98         2 or more            3         1        1         0         1
## 99         2 or more            3         1        1         1         0
## 100        2 or more            2         0        1         0         1
## 101        2 or more            3         1        1         1         0
## 102        2 or more            2         1        1         0         0
## 103        2 or more            2         1        1         0         0
## 104        2 or more            3         0        1         1         1
## 105        2 or more            3         1        0         1         1
## 106        2 or more            4         1        1         1         1
## 107        2 or more            2         1        1         0         0
## 108        2 or more            3         1        1         0         1
## 109        2 or more            3         0        1         1         1
## 110        2 or more            4         1        1         1         1
## 111        2 or more            4         1        1         1         1
## 112        2 or more            4         1        1         1         1
## 113        2 or more            4         1        1         1         1
## 114        2 or more            4         1        1         1         1
## 115        2 or more            4         1        1         1         1
## 116        2 or more            4         1        1         1         1
## 117        2 or more            2         1        1         0         0
## 118        2 or more            4         1        1         1         1
## 119        2 or more            2         1        1         0         0
## 120        2 or more            2         0        1         1         0
## 121        2 or more            4         1        1         1         1
## 122        2 or more            2         0        1         1         0
## 123        2 or more            3         1        1         1         0
## 124        2 or more            3         1        1         1         0
## 125        2 or more            2         1        1         0         0
## 126        2 or more            3         1        1         1         0
## 127        2 or more            4         1        1         1         1
## 128        2 or more            4         1        1         1         1
## 129        2 or more            4         1        1         1         1
## 130        2 or more            4         1        1         1         1
## 131        2 or more            2         0        1         1         0
## 132        2 or more            3         1        1         1         0
## 133        2 or more            3         1        1         1         0
## 134        2 or more            2         1        1         0         0
## 135        2 or more            4         1        1         1         1
## 136        2 or more            3         1        1         1         0
## 137        2 or more            4         1        1         1         1
## 138        2 or more            2         1        1         0         0
## 139        2 or more            2         0        1         1         0
## 140        2 or more            3         1        1         1         0
## 141        2 or more            3         1        1         0         1
## 142        2 or more            3         1        1         1         0
## 143        2 or more            2         0        1         1         0
## 144        2 or more            2         1        1         0         0
## 145        2 or more            2         0        1         1         0
## 146        2 or more            2         1        1         0         0
## 147        2 or more            2         1        0         0         1
## 148        2 or more            2         1        1         0         0
## 149        2 or more            2         1        1         0         0
## 150        2 or more            2         0        1         1         0
## 151        2 or more            2         0        1         1         0
## 152        2 or more            3         0        1         1         1
## 153        2 or more            2         0        1         1         0
## 154        2 or more            4         1        1         1         1
## 155        2 or more            2         0        1         1         0
## 156        2 or more            4         1        1         1         1
## 157        2 or more            2         1        1         0         0
## 158        2 or more            4         1        1         1         1
## 159        2 or more            2         0        1         1         0
## 160        2 or more            3         1        1         1         0
## 161        2 or more            2         1        1         0         0
## 162        2 or more            2         1        1         0         0
## 163        2 or more            4         1        1         1         1
## 164        2 or more            4         1        1         1         1
## 165        2 or more            4         1        1         1         1
## 166        2 or more            3         1        1         1         0
## 167        2 or more            4         1        1         1         1
## 168        2 or more            4         1        1         1         1
## 169        2 or more            4         1        1         1         1
## 170        2 or more            4         1        1         1         1
## 171        2 or more            2         1        1         0         0
## 172        2 or more            2         1        1         0         0
## 173        2 or more            3         1        1         1         0
## 174        2 or more            3         1        1         0         1
## 175        2 or more            3         1        1         1         0
## 176        2 or more            4         1        1         1         1
## 177        2 or more            2         1        1         0         0
## 178        2 or more            3         1        1         1         0
## 179        2 or more            2         1        1         0         0
## 180        2 or more            4         1        1         1         1
## 181        2 or more            2         0        0         1         1
## 182        2 or more            4         1        1         1         1
## 183        2 or more            2         0        1         1         0
## 184        2 or more            4         1        1         1         1
## 185        2 or more            3         1        1         0         1
## 186        2 or more            2         0        1         1         0
## 187        2 or more            3         1        1         1         0
## 188        2 or more            2         1        0         1         0
## 189        2 or more            2         1        1         0         0
## 190        2 or more            3         1        1         1         0
## 191        2 or more            3         1        1         0         1
## 192        2 or more            2         1        1         0         0
## 193        2 or more            2         1        1         0         0
## 194        2 or more            4         1        1         1         1
## 195        2 or more            2         1        0         1         0
## 196        2 or more            2         1        1         0         0
## 197        2 or more            3         1        1         1         0
## 198        2 or more            2         1        0         1         0
## 199        2 or more            2         0        1         1         0
## 200        2 or more            3         1        1         0         1
## 201        2 or more            2         0        1         1         0
## 202        2 or more            4         1        1         1         1
## 203        2 or more            2         0        1         1         0
## 204        2 or more            4         1        1         1         1
## 205        2 or more            2         1        1         0         0
## 206        2 or more            4         1        1         1         1
## 207        2 or more            2         1        1         0         0
## 208        2 or more            2         1        0         0         1
## 209        2 or more            4         1        1         1         1
## 210        2 or more            4         1        1         1         1
## 211        2 or more            3         1        1         1         0
## 212        2 or more            3         1        1         0         1
## 213        2 or more            2         1        1         0         0
## 214        2 or more            4         1        1         1         1
## 215        2 or more            2         0        1         1         0
## 216        2 or more            3         1        1         1         0
## 217        2 or more            2         1        1         0         0
## 218        2 or more            4         1        1         1         1
## 219        2 or more            2         1        1         0         0
## 220        2 or more            2         1        1         0         0
## 221        2 or more            4         1        1         1         1
## 222        2 or more            3         1        1         1         0
## 223        2 or more            2         1        1         0         0
## 224        2 or more            2         1        1         0         0
## 225        2 or more            3         1        1         1         0
## 226        2 or more            2         0        1         0         1
## 227        2 or more            4         1        1         1         1
## 228        2 or more            3         1        1         1         0
## 229        2 or more            3         1        1         1         0
## 230        2 or more            2         1        1         0         0
## 231        2 or more            3         1        0         1         1
## 232        2 or more            2         1        1         0         0
## 233        2 or more            2         0        1         1         0
## 234        2 or more            3         1        1         1         0
## 235        2 or more            2         1        1         0         0
## 236        2 or more            2         1        1         0         0
## 237        2 or more            2         1        1         0         0
## 238        2 or more            2         1        1         0         0
## 239        2 or more            4         1        1         1         1
## 240        2 or more            2         1        1         0         0
## 241        2 or more            2         1        1         0         0
## 242        2 or more            2         1        1         0         0
## 243        2 or more            2         0        1         1         0
## 244        2 or more            2         1        1         0         0
## 245        2 or more            4         1        1         1         1
## 246        2 or more            3         1        1         1         0
## 247        2 or more            4         1        1         1         1
## 248        2 or more            2         1        0         1         0
## 249        2 or more            2         1        1         0         0
## 250        2 or more            2         0        1         0         1
## 251        2 or more            3         1        1         1         0
## 252        2 or more            3         1        1         0         1
## 253        2 or more            2         1        1         0         0
## 254        2 or more            3         1        0         1         1
## 255        2 or more            3         1        1         1         0
## 256        2 or more            2         1        0         1         0
## 257        2 or more            3         1        1         1         0
## 258        2 or more            3         1        1         1         0
## 259        2 or more            3         1        1         1         0
## 260        2 or more            2         1        1         0         0
## 261        2 or more            3         1        1         1         0
## 262        2 or more            3         1        1         1         0
## 263        2 or more            2         0        1         1         0
## 264        2 or more            3         1        1         1         0
## 265        2 or more            3         1        1         0         1
## 266        2 or more            2         0        1         1         0
## 267        2 or more            2         1        1         0         0
## 268        2 or more            3         1        1         1         0
## 269        2 or more            3         1        1         1         0
## 270        2 or more            3         1        1         1         0
## 271        2 or more            4         1        1         1         1
## 272        2 or more            3         1        1         1         0
## 273        2 or more            4         1        1         1         1
## 274        2 or more            3         1        1         1         0
## 275        2 or more            3         1        1         1         0
## 276        2 or more            3         1        1         1         0
## 277        2 or more            2         0        1         0         1
## 278        2 or more            2         1        1         0         0
## 279        2 or more            2         0        1         1         0
## 280        2 or more            2         1        0         1         0
## 281        2 or more            3         1        1         1         0
## 282        2 or more            3         1        1         1         0
## 283        2 or more            4         1        1         1         1
## 284        2 or more            2         1        0         1         0
## 285        2 or more            2         1        1         0         0
## 286        2 or more            2         0        1         1         0
## 287        2 or more            3         1        1         0         1
## 288        2 or more            2         1        1         0         0
## 289        2 or more            4         1        1         1         1
## 290        2 or more            2         1        0         0         1
## 291        2 or more            2         1        1         0         0
## 292        2 or more            3         1        1         1         0
## 293        2 or more            4         1        1         1         1
## 294        2 or more            2         0        1         1         0
## 295        2 or more            2         1        1         0         0
## 296        2 or more            2         1        1         0         0
## 297        2 or more            2         1        1         0         0
## 298        2 or more            3         1        1         1         0
## 299        2 or more            2         1        1         0         0
## 300        2 or more            2         0        1         1         0
## 301        2 or more            2         1        1         0         0
## 302        2 or more            3         1        1         1         0
## 303        2 or more            2         0        1         1         0
## 304        2 or more            2         1        0         0         1
## 305        2 or more            2         1        1         0         0
## 306        2 or more            2         1        1         0         0
## 307        2 or more            3         1        0         1         1
## 308        2 or more            3         1        1         1         0
## 309        2 or more            2         1        1         0         0
## 310        2 or more            4         1        1         1         1
## 311        2 or more            3         1        1         1         0
## 312        2 or more            3         1        1         1         0
## 313        2 or more            4         1        1         1         1
## 314        2 or more            2         0        1         1         0
## 315        2 or more            2         1        1         0         0
## 316        2 or more            3         1        1         1         0
## 317        2 or more            2         1        1         0         0
## 318        2 or more            3         1        1         0         1
## 319        2 or more            2         1        1         0         0
## 320        2 or more            3         1        0         1         1
## 321        2 or more            2         1        1         0         0
## 322        2 or more            2         0        1         0         1
## 323        2 or more            3         1        1         1         0
## 324        2 or more            4         1        1         1         1
## 325        2 or more            2         1        1         0         0
## 326        2 or more            3         1        1         0         1
## 327        2 or more            3         1        1         0         1
## 328        2 or more            2         1        0         1         0
## 329        2 or more            4         1        1         1         1
## 330        2 or more            4         1        1         1         1
## 331        2 or more            2         0        1         0         1
## 332        2 or more            3         1        1         1         0
## 333        2 or more            3         1        1         1         0
## 334        2 or more            3         1        1         0         1
## 335        2 or more            4         1        1         1         1
## 336        2 or more            2         1        0         1         0
## 337        2 or more            4         1        1         1         1
## 338        2 or more            4         1        1         1         1
## 339        2 or more            4         1        1         1         1
## 340        2 or more            4         1        1         1         1
## 341        2 or more            4         1        1         1         1
## 342        2 or more            4         1        1         1         1
## 343        2 or more            4         1        1         1         1
## 344        2 or more            2         1        1         0         0
## 345        2 or more            4         1        1         1         1
## 346        2 or more            2         0        1         1         0
## 347        2 or more            2         1        1         0         0
## 348        2 or more            4         1        1         1         1
## 349        2 or more            4         1        1         1         1
## 350        2 or more            4         1        1         1         1
## 351        2 or more            4         1        1         1         1
## 352        2 or more            2         0        1         1         0
## 353        2 or more            2         1        0         1         0
## 354        2 or more            2         1        1         0         0
## 355        2 or more            4         1        1         1         1
## 356        2 or more            4         1        1         1         1
## 357        2 or more            2         1        1         0         0
## 358        2 or more            4         1        1         1         1
## 359        2 or more            2         0        1         1         0
## 360        2 or more            3         1        1         0         1
## 361        2 or more            4         1        1         1         1
## 362        2 or more            3         1        1         1         0
## 363        2 or more            4         1        1         1         1
## 364        2 or more            2         1        1         0         0
## 365        2 or more            4         1        1         1         1
## 366        2 or more            2         1        1         0         0
## 367        2 or more            3         1        1         1         0
## 368        2 or more            4         1        1         1         1
## 369        2 or more            4         1        1         1         1
## 370        2 or more            4         1        1         1         1
## 371        2 or more            4         1        1         1         1
## 372        2 or more            2         1        1         0         0
## 373        2 or more            2         1        1         0         0
## 374        2 or more            3         1        1         1         0
## 375        2 or more            3         1        1         1         0
## 376        2 or more            4         1        1         1         1
## 377        2 or more            2         1        1         0         0
## 378        2 or more            3         1        1         0         1
## 379        2 or more            4         1        1         1         1
## 380        2 or more            2         1        1         0         0
## 381        2 or more            2         1        1         0         0
## 382        2 or more            2         0        1         0         1
## 383        2 or more            3         1        1         1         0
## 384        2 or more            2         1        1         0         0
## 385        2 or more            4         1        1         1         1
## 386        2 or more            2         1        0         1         0
## 387        2 or more            4         1        1         1         1
## 388        2 or more            3         1        1         0         1
## 389        2 or more            3         0        1         1         1
## 390        2 or more            4         1        1         1         1
## 391        2 or more            3         1        1         1         0
## 392        2 or more            4         1        1         1         1
## 393        2 or more            2         0        1         0         1
## 394        2 or more            2         1        1         0         0
## 395        2 or more            2         1        1         0         0
## 396        2 or more            2         1        1         0         0
## 397        2 or more            4         1        1         1         1
## 398        2 or more            3         1        1         1         0
## 399        2 or more            4         1        1         1         1
## 400        2 or more            2         1        1         0         0
## 401        2 or more            3         0        1         1         1
## 402        2 or more            4         1        1         1         1
## 403        2 or more            2         1        1         0         0
## 404        2 or more            2         1        1         0         0
## 405        2 or more            4         1        1         1         1
## 406        2 or more            2         1        0         0         1
## 407        2 or more            4         1        1         1         1
## 408        2 or more            3         1        1         1         0
## 409        2 or more            4         1        1         1         1
## 410        2 or more            4         1        1         1         1
## 411        2 or more            4         1        1         1         1
## 412        2 or more            4         1        1         1         1
## 413        2 or more            4         1        1         1         1
## 414        2 or more            2         1        1         0         0
## 415        2 or more            2         1        0         0         1
## 416        2 or more            2         1        1         0         0
## 417        2 or more            3         1        1         1         0
## 418        2 or more            4         1        1         1         1
## 419        2 or more            3         1        1         1         0
## 420        2 or more            3         0        1         1         1
## 421        2 or more            2         1        1         0         0
## 422        2 or more            3         1        1         1         0
## 423        2 or more            2         0        1         1         0
## 424        2 or more            2         0        1         1         0
## 425        2 or more            3         1        1         1         0
## 426        2 or more            4         1        1         1         1
## 427        2 or more            4         1        1         1         1
## 428        2 or more            3         1        1         1         0
## 429        2 or more            2         1        1         0         0
## 430        2 or more            3         1        1         1         0
## 431        2 or more            2         0        1         1         0
## 432        2 or more            3         1        1         1         0
## 433        2 or more            2         0        1         1         0
## 434        2 or more            2         1        1         0         0

Overall
(N=1004)
total_injury_cat
2 or more 434 (43.2%)
Just one 251 (25.0%)
no injury 319 (31.8%)
total_injury
Mean (SD) 1.49 (1.37)
Median [Min, Max] 1.00 [0, 4.00]
## 
##  No Yes 
## 553 451
## 
##   0   1 
## 553 451

How I made TBI

##       TBI crash fall fight
## 1  No TBI    No   No    No
## 2     TBI   Yes  Yes   Yes
## 3     TBI    No  Yes    No
## 4  No TBI    No   No    No
## 5  No TBI    No   No    No
## 6  No TBI    No   No    No
## 7  No TBI    No   No    No
## 8     TBI   Yes  Yes   Yes
## 9  No TBI    No   No    No
## 10    TBI   Yes  Yes    No
## 11 No TBI    No   No    No
## 12    TBI   Yes   No    No
## 13    TBI   Yes  Yes    No
## 14 No TBI    No   No    No
## 15    TBI    No  Yes   Yes
## 16    TBI   Yes  Yes   Yes
## 17 No TBI    No   No    No
## 18    TBI    No  Yes   Yes
## 19 No TBI    No   No    No
## 20    TBI   Yes  Yes   Yes
## 21    TBI   Yes  Yes   Yes
## 22    TBI   Yes  Yes   Yes
## 23    TBI   Yes  Yes   Yes
## 24 No TBI    No   No    No
## 25    TBI    No  Yes   Yes
## 26    TBI    No  Yes    No
## 27 No TBI    No   No    No
## 28    TBI   Yes   No    No
## 29    TBI   Yes  Yes   Yes
## 30    TBI    No  Yes    No

#Plot TBI severity

## Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
## `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.

## Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
## `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.

TBI numbers

## [1] TBI            TBInumbers     TBInumbers_cat crash          fall          
## [6] fight          blast         
## <0 rows> (or 0-length row.names)

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

##       TBI fall_times_num crash_times_num fight_times_num TBInumbers
## 1  No TBI             NA              NA              NA          0
## 2     TBI              2               3               6         11
## 3     TBI              1              NA              NA          1
## 4  No TBI             NA              NA              NA          0
## 5  No TBI             NA              NA              NA          0
## 6  No TBI             NA              NA              NA          0
## 7  No TBI             NA              NA              NA          0
## 8     TBI              4               4               4         12
## 9  No TBI             NA              NA              NA          0
## 10    TBI              5               1              NA          6
## 11 No TBI             NA              NA              NA          0
## 12    TBI             NA               1              NA          1
## 13    TBI              1               1              NA          2
## 14 No TBI             NA              NA              NA          0
## 15    TBI              1              NA               1          2
## 16    TBI              2               2               2          6
## 17 No TBI             NA              NA              NA          0
## 18    TBI              2              NA               9         11
## 19 No TBI             NA              NA              NA          0
## 20    TBI              1               2               4          7
## 21    TBI              1               3               2          6
## 22    TBI              2               2               2          6
## 23    TBI              3               2               1          6
## 24 No TBI             NA              NA              NA          0
## 25    TBI              2              NA               2          4
## 26    TBI              1              NA              NA          1
## 27 No TBI             NA              NA              NA          0
## 28    TBI             NA               3              NA          3
## 29    TBI              1               3               1          5
## 30    TBI              1              NA              NA          1
## 31 No TBI             NA              NA              NA          0
## 32    TBI              2               2              NA          4
## 33    TBI              1               1              NA          2
## 34 No TBI             NA              NA              NA          0
## 35    TBI              3              NA              NA          3
## 36    TBI             NA               1               1          2
## 37 No TBI             NA              NA              NA          0
## 38 No TBI             NA              NA              NA          0
## 39    TBI              1              NA              NA          1
## 40 No TBI             NA              NA              NA          0
## 41 No TBI             NA              NA              NA          0
## 42    TBI              1               1               1          3
## 43 No TBI             NA              NA              NA          0
## 44    TBI              1               2               3          6
## 45    TBI              1               1              NA          2
## 46    TBI              1              NA              NA          1
## 47    TBI              3               1               1          5
## 48    TBI              1              NA              NA          1
## 49    TBI              1               1              NA          2
## 50 No TBI             NA              NA              NA          0

experienced multiple repeated blows or impacts to your head

##    physical_abuse playing_sports other_injury military_exposures Blows
## 1              No             No           No                 No    No
## 2             Yes            Yes          Yes                 No   Yes
## 3              No             No           No                 No    No
## 4              No             No           No                 No    No
## 5              No             No           No                 No    No
## 6              No            Yes          Yes                 No   Yes
## 7              No             No           No                 No    No
## 8              No             No           No                 No    No
## 9              No             No           No                 No    No
## 10            Yes             No           No                 No   Yes
## 11             No             No           No                 No    No
## 12             No             No           No                 No    No
## 13             No             No           No                 No    No
## 14            Yes             No          Yes                Yes   Yes
## 15            Yes            Yes           No                 No   Yes
## 16            Yes             No           No                Yes   Yes
## 17             No            Yes           No                 No   Yes
## 18             No            Yes          Yes                 No   Yes
## 19             No             No           No                 No    No
## 20            Yes            Yes          Yes                 No   Yes
## 21            Yes            Yes          Yes                Yes   Yes
## 22            Yes            Yes           No                Yes   Yes
## 23             No             No           No                 No    No
## 24             No             No           No                 No    No
## 25            Yes            Yes          Yes                 No   Yes
## 26             No            Yes           No                 No   Yes
## 27             No             No           No                 No    No
## 28             No            Yes           No                 No   Yes
## 29            Yes            Yes           No                 No   Yes
## 30             No             No           No                 No    No

Interpreting Findings

##                          TBI_severity Worst first_unconsciousness_age First
## 1                              No TBI    No               0Not Answer    No
## 2  Mild with LOC from 5 to 30 minutes    No                        63    No
## 3                               Dazed    No               0Not Answer    No
## 4                              No TBI    No               0Not Answer    No
## 5                              No TBI    No               0Not Answer    No
## 6                              No TBI    No               0Not Answer    No
## 7                              No TBI    No               0Not Answer    No
## 8  Mild with LOC from 5 to 30 minutes    No                        18    No
## 9                              No TBI    No               0Not Answer    No
## 10               Mild TBI without LOC    No               0Not Answer    No
## 11                             No TBI    No               0Not Answer    No
## 12               Mild TBI without LOC    No               0Not Answer    No
## 13               Mild TBI without LOC    No               0Not Answer    No
## 14                             No TBI    No                        21    No
## 15               Mild TBI without LOC    No               0Not Answer    No
## 16            Mild with LOC<5 minutes    No                        25    No
## 17                             No TBI    No               0Not Answer    No
## 18 Mild with LOC from 5 to 30 minutes    No                        55    No
## 19                             No TBI    No               0Not Answer    No
## 20 Mild with LOC from 5 to 30 minutes    No                        64    No
## 21 Mild with LOC from 5 to 30 minutes    No                       100    No
## 22 Mild with LOC from 5 to 30 minutes    No                        86    No
## 23            Mild with LOC<5 minutes    No                        46    No
## 24                             No TBI    No               0Not Answer    No
## 25  Moderate [LOC from\n30min to 24h]   Yes                         9   Yes
## 26               Mild TBI without LOC    No               0Not Answer    No
## 27                             No TBI    No               0Not Answer    No
## 28               Mild TBI without LOC    No               0Not Answer    No
## 29  Moderate [LOC from\n30min to 24h]   Yes                        21    No
## 30               Mild TBI without LOC    No               0Not Answer    No
##    unconsciousness_drug unconsciousness_choked Other_TBI
## 1                    No                     No        No
## 2                   Yes                    Yes       Yes
## 3                    No                     No        No
## 4                    No                     No        No
## 5                    No                     No        No
## 6                    No                     No        No
## 7                    No                     No        No
## 8                    No                     No        No
## 9                    No                     No        No
## 10                   No                     No        No
## 11                   No                     No        No
## 12                   No                     No        No
## 13                   No                     No        No
## 14                   No                     No        No
## 15                   No                     No        No
## 16                  Yes                     No       Yes
## 17                   No                     No        No
## 18                   No                     No        No
## 19                   No                     No        No
## 20                   No                     No        No
## 21                  Yes                    Yes       Yes
## 22                  Yes                    Yes       Yes
## 23                  Yes                    Yes       Yes
## 24                   No                     No        No
## 25                   No                     No        No
## 26                   No                     No        No
## 27                   No                     No        No
## 28                   No                     No        No
## 29                  Yes                     No       Yes
## 30                   No                     No        No

Hospitalization

##    carsh_treatment fall_treatment fight_treatment blast_treatment
## 1      0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 2              Yes            Yes             Yes             Yes
## 3      0Not Answer            Yes     0Not Answer     0Not Answer
## 4      0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 5      0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 6      0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 7      0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 8              Yes            Yes             Yes             Yes
## 9      0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 10             Yes             No     0Not Answer     0Not Answer
## 11     0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 12              No    0Not Answer     0Not Answer     0Not Answer
## 13             Yes             No     0Not Answer     0Not Answer
## 14     0Not Answer    0Not Answer     0Not Answer             Yes
## 15     0Not Answer            Yes             Yes     0Not Answer
## 16             Yes            Yes             Yes             Yes
## 17     0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 18     0Not Answer            Yes              No     0Not Answer
## 19     0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 20             Yes            Yes             Yes             Yes
## 21             Yes            Yes             Yes             Yes
## 22              No             No              No             Yes
## 23             Yes            Yes             Yes             Yes
## 24     0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 25     0Not Answer            Yes             Yes             Yes
## 26     0Not Answer             No     0Not Answer     0Not Answer
## 27     0Not Answer    0Not Answer     0Not Answer     0Not Answer
## 28             Yes    0Not Answer     0Not Answer     0Not Answer
## 29             Yes            Yes             Yes             Yes
## 30     0Not Answer             No     0Not Answer     0Not Answer
##    Hospitalizatied
## 1               No
## 2              Yes
## 3              Yes
## 4               No
## 5               No
## 6               No
## 7               No
## 8              Yes
## 9               No
## 10             Yes
## 11              No
## 12              No
## 13             Yes
## 14             Yes
## 15             Yes
## 16             Yes
## 17              No
## 18             Yes
## 19              No
## 20             Yes
## 21             Yes
## 22             Yes
## 23             Yes
## 24              No
## 25             Yes
## 26              No
## 27              No
## 28             Yes
## 29             Yes
## 30              No

Plots

TBI Distribution

table1(~fall+crash+fight+physical_abuse+Hospitalizatied, data=PTH)
Overall
(N=1004)
fall
No 458 (45.6%)
Yes 546 (54.4%)
crash
No 553 (55.1%)
Yes 451 (44.9%)
fight
No 708 (70.5%)
Yes 296 (29.5%)
physical_abuse
No 702 (69.9%)
Yes 302 (30.1%)
Hospitalizatied
No 505 (50.3%)
Yes 499 (49.7%)

Stratified by TBI status

No TBI
(N=335)
TBI
(N=669)
Overall
(N=1004)
age
Mean (SD) 49.4 (17.9) 47.7 (21.0) 48.2 (20.0)
Median [Min, Max] 48.0 [18.0, 120] 41.0 [18.0, 120] 43.0 [18.0, 120]
Race
White 255 (76.1%) 495 (74.0%) 750 (74.7%)
Black 39 (11.6%) 77 (11.5%) 116 (11.6%)
Hispanic 16 (4.8%) 64 (9.6%) 80 (8.0%)
Mixed 25 (7.5%) 33 (4.9%) 58 (5.8%)
Education
High School 90 (26.9%) 146 (21.8%) 236 (23.5%)
Less than High School 24 (7.2%) 49 (7.3%) 73 (7.3%)
More than College 112 (33.4%) 201 (30.0%) 313 (31.2%)
Some College 109 (32.5%) 273 (40.8%) 382 (38.0%)
HIT
Mean (SD) 57.2 (9.56) 61.1 (8.78) 59.8 (9.23)
Median [Min, Max] 58.0 [36.0, 78.0] 61.0 [36.0, 78.0] 60.0 [36.0, 78.0]
Headache Severity
Mean (SD) 5.03 (2.54) 5.83 (2.40) 5.56 (2.47)
Median [Min, Max] 5.00 [0, 10.0] 6.00 [0, 10.0] 6.00 [0, 10.0]
Sex
Female 213 (63.6%) 279 (41.7%) 492 (49.0%)
Male 122 (36.4%) 390 (58.3%) 512 (51.0%)
PHQ_2
No Depression 228 (68.1%) 333 (49.8%) 561 (55.9%)
Depression 107 (31.9%) 336 (50.2%) 443 (44.1%)
PHQ
Mean (SD) 1.81 (1.95) 2.64 (1.87) 2.36 (1.94)
Median [Min, Max] 1.00 [0, 6.00] 3.00 [0, 6.00] 2.00 [0, 6.00]
Cognetive
Mean (SD) 0.938 (1.03) 1.52 (1.13) 1.33 (1.13)
Median [Min, Max] 0.500 [0, 4.00] 1.50 [0, 4.00] 1.25 [0, 4.00]
Vestibular
Mean (SD) 0.909 (1.00) 1.50 (1.11) 1.30 (1.11)
Median [Min, Max] 0.667 [0, 4.00] 1.33 [0, 4.00] 1.00 [0, 4.00]
Alcohol
Non-drinker 143 (42.7%) 180 (26.9%) 323 (32.2%)
Heavy 57 (17.0%) 181 (27.1%) 238 (23.7%)
Moderate 135 (40.3%) 308 (46.0%) 443 (44.1%)
TBInumbers_cat
0 335 (100%) 0 (0%) 335 (33.4%)
1 0 (0%) 136 (20.3%) 136 (13.5%)
2 0 (0%) 152 (22.7%) 152 (15.1%)
3 0 (0%) 91 (13.6%) 91 (9.1%)
4 0 (0%) 44 (6.6%) 44 (4.4%)
4-8 0 (0%) 139 (20.8%) 139 (13.8%)
9+ 0 (0%) 107 (16.0%) 107 (10.7%)
GAD_2
No Depression 225 (67.2%) 340 (50.8%) 565 (56.3%)
Depression 110 (32.8%) 329 (49.2%) 439 (43.7%)
## Warning in compareGroups.fit(X = X, y = y, include.label = include.label, :
## variables waist specified in 'method' not found
## 
## --------Summary descriptives table by 'TBI'---------
## 
## ___________________________________________________________ 
##                             No TBI        TBI     p.overall 
##                              N=335       N=669              
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯ 
## age                       49.4 (17.9) 47.7 (21.0)   0.176   
## Race:                                               0.029   
##     White                 255 (76.1%) 495 (74.0%)           
##     Black                 39 (11.6%)  77 (11.5%)            
##     Hispanic              16 (4.78%)  64 (9.57%)            
##     Mixed                 25 (7.46%)  33 (4.93%)            
## Education:                                          0.064   
##     High School           90 (26.9%)  146 (21.8%)           
##     Less than High School 24 (7.16%)  49 (7.32%)            
##     More than College     112 (33.4%) 201 (30.0%)           
##     Some College          109 (32.5%) 273 (40.8%)           
## HIT                       57.2 (9.56) 61.1 (8.78)  <0.001   
## Headache Severity         5.03 (2.54) 5.83 (2.40)  <0.001   
## Sex:                                               <0.001   
##     Female                213 (63.6%) 279 (41.7%)           
##     Male                  122 (36.4%) 390 (58.3%)           
## PHQ                       1.81 (1.95) 2.64 (1.87)  <0.001   
## Cognetive                 0.94 (1.03) 1.52 (1.13)  <0.001   
## Vestibular                0.91 (1.00) 1.50 (1.11)  <0.001   
## Alcohol:                                           <0.001   
##     Non-drinker           143 (42.7%) 180 (26.9%)           
##     Heavy                 57 (17.0%)  181 (27.1%)           
##     Moderate              135 (40.3%) 308 (46.0%)           
## TBInumbers_cat:                                    <0.001   
##     0                     335 (100%)   0 (0.00%)            
##     1                      0 (0.00%)  136 (20.3%)           
##     2                      0 (0.00%)  152 (22.7%)           
##     3                      0 (0.00%)  91 (13.6%)            
##     4                      0 (0.00%)  44 (6.58%)            
##     4-8                    0 (0.00%)  139 (20.8%)           
##     9+                     0 (0.00%)  107 (16.0%)           
## GAD_2:                                             <0.001   
##     No Depression         225 (67.2%) 340 (50.8%)           
##     Depression            110 (32.8%) 329 (49.2%)           
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯

Stratified by TBI severity

Mild
(N=225)
Mild TBI without LOC
(N=295)
Moderate/Severe
(N=149)
No TBI
(N=335)
Overall
(N=1004)
age
Mean (SD) 47.9 (20.5) 46.7 (18.1) 49.0 (26.3) 49.4 (17.9) 48.2 (20.0)
Median [Min, Max] 42.0 [18.0, 120] 41.0 [18.0, 112] 40.0 [21.0, 120] 48.0 [18.0, 120] 43.0 [18.0, 120]
Race
White 172 (76.4%) 222 (75.3%) 101 (67.8%) 255 (76.1%) 750 (74.7%)
Black 26 (11.6%) 31 (10.5%) 20 (13.4%) 39 (11.6%) 116 (11.6%)
Hispanic 18 (8.0%) 24 (8.1%) 22 (14.8%) 16 (4.8%) 80 (8.0%)
Mixed 9 (4.0%) 18 (6.1%) 6 (4.0%) 25 (7.5%) 58 (5.8%)
Education
High School 49 (21.8%) 73 (24.7%) 24 (16.1%) 90 (26.9%) 236 (23.5%)
Less than High School 11 (4.9%) 22 (7.5%) 16 (10.7%) 24 (7.2%) 73 (7.3%)
More than College 77 (34.2%) 88 (29.8%) 36 (24.2%) 112 (33.4%) 313 (31.2%)
Some College 88 (39.1%) 112 (38.0%) 73 (49.0%) 109 (32.5%) 382 (38.0%)
HIT
Mean (SD) 60.8 (8.48) 59.4 (8.55) 65.0 (8.54) 57.2 (9.56) 59.8 (9.23)
Median [Min, Max] 62.0 [36.0, 78.0] 60.0 [36.0, 78.0] 65.0 [36.0, 78.0] 58.0 [36.0, 78.0] 60.0 [36.0, 78.0]
Headache Severity
Mean (SD) 5.74 (2.34) 5.36 (2.30) 6.89 (2.37) 5.03 (2.54) 5.56 (2.47)
Median [Min, Max] 6.00 [0, 10.0] 5.00 [0, 10.0] 7.00 [1.00, 10.0] 5.00 [0, 10.0] 6.00 [0, 10.0]
Sex
Female 92 (40.9%) 150 (50.8%) 37 (24.8%) 213 (63.6%) 492 (49.0%)
Male 133 (59.1%) 145 (49.2%) 112 (75.2%) 122 (36.4%) 512 (51.0%)
PHQ_2
No Depression 119 (52.9%) 168 (56.9%) 46 (30.9%) 228 (68.1%) 561 (55.9%)
Depression 106 (47.1%) 127 (43.1%) 103 (69.1%) 107 (31.9%) 443 (44.1%)
PHQ
Mean (SD) 2.55 (1.74) 2.30 (1.91) 3.47 (1.75) 1.81 (1.95) 2.36 (1.94)
Median [Min, Max] 2.00 [0, 6.00] 2.00 [0, 6.00] 3.00 [0, 6.00] 1.00 [0, 6.00] 2.00 [0, 6.00]
Cognetive
Mean (SD) 1.55 (1.07) 1.12 (0.985) 2.27 (1.10) 0.938 (1.03) 1.33 (1.13)
Median [Min, Max] 1.50 [0, 4.00] 1.00 [0, 4.00] 2.25 [0, 4.00] 0.500 [0, 4.00] 1.25 [0, 4.00]
Vestibular
Mean (SD) 1.41 (0.966) 1.15 (0.978) 2.31 (1.16) 0.909 (1.00) 1.30 (1.11)
Median [Min, Max] 1.33 [0, 4.00] 1.00 [0, 4.00] 2.33 [0, 4.00] 0.667 [0, 4.00] 1.00 [0, 4.00]
Alcohol
Non-drinker 59 (26.2%) 103 (34.9%) 18 (12.1%) 143 (42.7%) 323 (32.2%)
Heavy 49 (21.8%) 59 (20.0%) 73 (49.0%) 57 (17.0%) 238 (23.7%)
Moderate 117 (52.0%) 133 (45.1%) 58 (38.9%) 135 (40.3%) 443 (44.1%)
## Warning in compareGroups.fit(X = X, y = y, include.label = include.label, :
## variables waist specified in 'method' not found
## 
## --------Summary descriptives table by 'TBI_severity_finalcc'---------
## 
## ________________________________________________________________________________________________ 
##                              Mild     Mild TBI without LOC Moderate/Severe   No TBI    p.overall 
##                              N=225           N=295              N=149         N=335              
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯ 
## age                       47.9 (20.5)     46.7 (18.1)        49.0 (26.3)   49.4 (17.9)   0.389   
## Race:                                                                                    0.029   
##     White                 172 (76.4%)     222 (75.3%)        101 (67.8%)   255 (76.1%)           
##     Black                 26 (11.6%)       31 (10.5%)        20 (13.4%)    39 (11.6%)            
##     Hispanic              18 (8.00%)       24 (8.14%)        22 (14.8%)    16 (4.78%)            
##     Mixed                  9 (4.00%)       18 (6.10%)         6 (4.03%)    25 (7.46%)            
## Education:                                                                               0.013   
##     High School           49 (21.8%)       73 (24.7%)        24 (16.1%)    90 (26.9%)            
##     Less than High School 11 (4.89%)       22 (7.46%)        16 (10.7%)    24 (7.16%)            
##     More than College     77 (34.2%)       88 (29.8%)        36 (24.2%)    112 (33.4%)           
##     Some College          88 (39.1%)      112 (38.0%)        73 (49.0%)    109 (32.5%)           
## HIT                       60.8 (8.48)     59.4 (8.55)        65.0 (8.54)   57.2 (9.56)  <0.001   
## Headache Severity         5.74 (2.34)     5.36 (2.30)        6.89 (2.37)   5.03 (2.54)  <0.001   
## Sex:                                                                                    <0.001   
##     Female                92 (40.9%)      150 (50.8%)        37 (24.8%)    213 (63.6%)           
##     Male                  133 (59.1%)     145 (49.2%)        112 (75.2%)   122 (36.4%)           
## PHQ                       2.55 (1.74)     2.30 (1.91)        3.47 (1.75)   1.81 (1.95)  <0.001   
## Cognetive                 1.55 (1.07)     1.12 (0.98)        2.27 (1.10)   0.94 (1.03)  <0.001   
## Vestibular                1.41 (0.97)     1.15 (0.98)        2.31 (1.16)   0.91 (1.00)  <0.001   
## Alcohol:                                                                                <0.001   
##     Non-drinker           59 (26.2%)      103 (34.9%)        18 (12.1%)    143 (42.7%)           
##     Heavy                 49 (21.8%)       59 (20.0%)        73 (49.0%)    57 (17.0%)            
##     Moderate              117 (52.0%)     133 (45.1%)        58 (38.9%)    135 (40.3%)           
## TBInumbers_cat:                                                                         <0.001   
##     0                      0 (0.00%)       0 (0.00%)          0 (0.00%)    335 (100%)            
##     1                     30 (13.3%)       95 (32.2%)        11 (7.38%)     0 (0.00%)            
##     2                     50 (22.2%)       94 (31.9%)         8 (5.37%)     0 (0.00%)            
##     3                     35 (15.6%)       38 (12.9%)        18 (12.1%)     0 (0.00%)            
##     4                     17 (7.56%)       16 (5.42%)        11 (7.38%)     0 (0.00%)            
##     4-8                   61 (27.1%)       36 (12.2%)        42 (28.2%)     0 (0.00%)            
##     9+                    32 (14.2%)       16 (5.42%)        59 (39.6%)     0 (0.00%)            
## GAD_2:                                                                                  <0.001   
##     No Depression         114 (50.7%)     176 (59.7%)        50 (33.6%)    225 (67.2%)           
##     Depression            111 (49.3%)     119 (40.3%)        99 (66.4%)    110 (32.8%)           
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯

TBI variables

Overall
(N=1004)
TBInumbers
Mean (SD) 3.15 (4.24)
Median [Min, Max] 2.00 [0, 27.0]

Headache Charactersitics

No TBI
(N=335)
TBI
(N=669)
Overall
(N=1004)
family_history_two
No 121 (36.1%) 134 (20.0%) 255 (25.4%)
Yes 214 (63.9%) 535 (80.0%) 749 (74.6%)
headache_hypersensitivity
1-4 times 137 (40.9%) 315 (47.1%) 452 (45.0%)
10+ times 43 (12.8%) 89 (13.3%) 132 (13.1%)
5-9 times 36 (10.7%) 110 (16.4%) 146 (14.5%)
Never 119 (35.5%) 155 (23.2%) 274 (27.3%)
any_headache_time
at least 3 months 76 (22.7%) 233 (34.8%) 309 (30.8%)
more than 6 months 31 (9.3%) 110 (16.4%) 141 (14.0%)
Never 138 (41.2%) 217 (32.4%) 355 (35.4%)
No 90 (26.9%) 109 (16.3%) 199 (19.8%)
days_migraine_num
Mean (SD) 3.92 (5.14) 7.06 (7.75) 6.01 (7.14)
Median [Min, Max] 3.00 [0, 30.0] 5.00 [0, 30.0] 3.00 [0, 30.0]
days_other_headache_num
Mean (SD) 4.09 (5.33) 7.39 (7.80) 6.29 (7.24)
Median [Min, Max] 3.00 [0, 30.0] 5.00 [0, 30.0] 4.00 [0, 30.0]
Mild
Mild TBI without LOC
Moderate/Severe
No TBI
Overall
TBI
(N=225)
TBI
(N=295)
TBI
(N=149)
No TBI
(N=335)
No TBI
(N=335)
TBI
(N=669)
family_history_two
No 40 (17.8%) 79 (26.8%) 15 (10.1%) 121 (36.1%) 121 (36.1%) 134 (20.0%)
Yes 185 (82.2%) 216 (73.2%) 134 (89.9%) 214 (63.9%) 214 (63.9%) 535 (80.0%)
headache_hypersensitivity
1-4 times 112 (49.8%) 128 (43.4%) 75 (50.3%) 137 (40.9%) 137 (40.9%) 315 (47.1%)
10+ times 26 (11.6%) 53 (18.0%) 10 (6.7%) 43 (12.8%) 43 (12.8%) 89 (13.3%)
5-9 times 34 (15.1%) 31 (10.5%) 45 (30.2%) 36 (10.7%) 36 (10.7%) 110 (16.4%)
Never 53 (23.6%) 83 (28.1%) 19 (12.8%) 119 (35.5%) 119 (35.5%) 155 (23.2%)
any_headache_time
at least 3 months 85 (37.8%) 81 (27.5%) 67 (45.0%) 76 (22.7%) 76 (22.7%) 233 (34.8%)
more than 6 months 31 (13.8%) 35 (11.9%) 44 (29.5%) 31 (9.3%) 31 (9.3%) 110 (16.4%)
Never 75 (33.3%) 121 (41.0%) 21 (14.1%) 138 (41.2%) 138 (41.2%) 217 (32.4%)
No 34 (15.1%) 58 (19.7%) 17 (11.4%) 90 (26.9%) 90 (26.9%) 109 (16.3%)
days_migraine_num
Mean (SD) 6.33 (6.94) 5.14 (6.05) 12.0 (9.64) 3.92 (5.14) 3.92 (5.14) 7.06 (7.75)
Median [Min, Max] 5.00 [0, 30.0] 3.00 [0, 30.0] 10.0 [0, 30.0] 3.00 [0, 30.0] 3.00 [0, 30.0] 5.00 [0, 30.0]
days_other_headache_num
Mean (SD) 6.84 (7.07) 5.65 (6.22) 11.7 (9.85) 4.09 (5.33) 4.09 (5.33) 7.39 (7.80)
Median [Min, Max] 5.00 [0, 30.0] 3.00 [0, 30.0] 9.00 [0, 30.0] 3.00 [0, 30.0] 3.00 [0, 30.0] 5.00 [0, 30.0]

Injuries reported

table1(~ crash+fall+ fight+ blast+physical_abuse+ playing_sports+military_exposures, data=PTH)
Overall
(N=1004)
crash
No 553 (55.1%)
Yes 451 (44.9%)
fall
No 458 (45.6%)
Yes 546 (54.4%)
fight
No 708 (70.5%)
Yes 296 (29.5%)
blast
No 801 (79.8%)
Yes 203 (20.2%)
physical_abuse
No 702 (69.9%)
Yes 302 (30.1%)
playing_sports
No 669 (66.6%)
Yes 335 (33.4%)
military_exposures
No 869 (86.6%)
Yes 135 (13.4%)
hist(PTH$TBInumbers)

PTH$TBInumbers<-as.numeric(PTH$TBInumbers)

Analysis of Covariance(ANCOVA)

## Loading required package: mvtnorm
## Loading required package: survival
## Loading required package: TH.data
## Loading required package: MASS
## 
## Attaching package: 'MASS'
## The following object is masked from 'package:gtsummary':
## 
##     select
## The following object is masked from 'package:dplyr':
## 
##     select
## 
## Attaching package: 'TH.data'
## The following object is masked from 'package:MASS':
## 
##     geyser
ANCOVA for HIT
Sum Sq Df F value Pr(>F)
(Intercept) 163994.0058 1 2324.6404928 0.0000000
TBI_severity_finalcc 725.0563 3 3.4259287 0.0167103
TBInumbers 1833.6878 1 25.9928092 0.0000004
Hospitalizatied 250.6943 1 3.5536300 0.0597089
Blows 1067.5855 1 15.1331901 0.0001069
Gender 2984.7931 1 42.3099049 0.0000000
Race 877.1239 3 4.1444560 0.0062407
education 108.9495 3 0.5147922 0.6721643
Alcohol 240.7144 2 1.7060820 0.1821108
age 1076.7602 1 15.2632425 0.0000999
Residuals 69628.8670 987 NA NA
ANCOVA for Headaches Severity
Sum Sq Df F value Pr(>F)
(Intercept) 787.44369 1 149.7473563 0.0000000
TBI_severity_finalcc 60.50727 3 3.8355350 0.0095439
TBInumbers 89.26027 1 16.9745328 0.0000411
Hospitalizatied 47.85875 1 9.1012493 0.0026198
Blows 56.60095 1 10.7637451 0.0010713
Gender 206.73548 1 39.3146731 0.0000000
Race 10.17805 3 0.6451834 0.5860839
education 16.37394 3 1.0379386 0.3749018
Alcohol 51.46551 2 4.8935715 0.0076774
age 28.07361 1 5.3387294 0.0210621
Residuals 5190.12118 987 NA NA
ANCOVA for Cognetive
Sum Sq Df F value Pr(>F)
(Intercept) 48.6394327 1 52.9137932 0.0000000
TBI_severity_finalcc 30.3723935 3 11.0138246 0.0000004
TBInumbers 43.6770570 1 47.5153313 0.0000000
Hospitalizatied 0.5060583 1 0.5505300 0.4582778
Blows 27.6405976 1 30.0696119 0.0000001
Gender 9.7085014 1 10.5616700 0.0011936
Race 3.6384125 3 1.3193836 0.2666683
education 1.0262826 3 0.3721569 0.7731155
Alcohol 22.5337081 2 12.2569683 0.0000055
age 8.3873113 1 9.1243757 0.0025872
Residuals 907.2704333 987 NA NA
ANCOVA for Vestibular
Sum Sq Df F value Pr(>F)
(Intercept) 17.667382 1 20.894594 0.0000055
TBI_severity_finalcc 28.686815 3 11.308964 0.0000003
TBInumbers 38.106379 1 45.067080 0.0000000
Hospitalizatied 4.719861 1 5.582015 0.0183390
Blows 35.523859 1 42.012824 0.0000000
Gender 5.651988 1 6.684408 0.0098684
Race 2.700152 3 1.064458 0.3632723
education 6.471348 3 2.551146 0.0543647
Alcohol 23.668349 2 13.995864 0.0000010
age 3.159342 1 3.736443 0.0535234
Residuals 834.555856 987 NA NA
## 
##   Simultaneous Tests for General Linear Hypotheses
## 
## Fit: aov(formula = Vestibular ~ TBI_severity_finalcc + TBInumbers + 
##     Hospitalizatied + Blows + Gender + Race + education + Alcohol + 
##     age, data = PTH)
## 
## Linear Hypotheses:
##                                                Estimate Std. Error t value
## (Intercept) == 0                               0.634060   0.138712   4.571
## TBI_severity_finalccMild TBI without LOC == 0 -0.066961   0.085294  -0.785
## TBI_severity_finalccModerate/Severe == 0       0.472924   0.103129   4.586
## TBI_severity_finalccNo TBI == 0                0.137377   0.110904   1.239
## TBInumbers == 0                                0.062711   0.009341   6.713
## HospitalizatiedYes == 0                        0.207438   0.087800   2.363
## BlowsYes == 0                                  0.440707   0.067992   6.482
## GenderMale == 0                               -0.167839   0.064917  -2.585
## RaceBlack == 0                                -0.058852   0.094337  -0.624
## RaceHispanic == 0                              0.118714   0.110874   1.071
## RaceMixed == 0                                -0.148798   0.126327  -1.178
## educationLess than High School == 0            0.218923   0.125948   1.738
## educationMore than College == 0                0.014684   0.079665   0.184
## educationSome College == 0                     0.170787   0.080384   2.125
## AlcoholHeavy == 0                              0.451202   0.085954   5.249
## AlcoholModerate == 0                           0.157824   0.071521   2.207
## age == 0                                      -0.002884   0.001492  -1.933
##                                               Pr(>|t|)    
## (Intercept) == 0                                <0.001 ***
## TBI_severity_finalccMild TBI without LOC == 0    1.000    
## TBI_severity_finalccModerate/Severe == 0        <0.001 ***
## TBI_severity_finalccNo TBI == 0                  0.970    
## TBInumbers == 0                                 <0.001 ***
## HospitalizatiedYes == 0                          0.250    
## BlowsYes == 0                                   <0.001 ***
## GenderMale == 0                                  0.144    
## RaceBlack == 0                                   1.000    
## RaceHispanic == 0                                0.992    
## RaceMixed == 0                                   0.981    
## educationLess than High School == 0              0.720    
## educationMore than College == 0                  1.000    
## educationSome College == 0                       0.408    
## AlcoholHeavy == 0                               <0.001 ***
## AlcoholModerate == 0                             0.349    
## age == 0                                         0.562    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
ANCOVA for PHQ
Sum Sq Df F value Pr(>F)
(Intercept) 278.4030319 1 90.2782749 0.0000000
TBI_severity_finalcc 22.3499198 3 2.4158169 0.0650726
TBInumbers 76.4328886 1 24.7850366 0.0000008
Hospitalizatied 0.1116799 1 0.0362147 0.8491120
Blows 78.6039454 1 25.4890493 0.0000005
Gender 46.9181222 1 15.2142278 0.0001025
Race 15.0622392 3 1.6280869 0.1812070
education 5.2997015 3 0.5728481 0.6329573
Alcohol 28.7033804 2 4.6538496 0.0097349
age 62.6256625 1 20.3077414 0.0000074
Residuals 3043.7421700 987 NA NA
fit1=lm(HIT~TBI_severity_finalcc+TBInumbers_cat+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
#fit2=lm(pain_level~TBI_severity_finalcc+TBInumbers+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
fit3=lm(Cognetive~TBI_severity_finalcc+TBInumbers_cat+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
fit4=lm(Vestibular~TBI_severity_finalcc+TBInumbers_cat+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)
fit5=lm(PHQ~TBI_severity_finalcc+TBInumbers_cat+Blows+Gender+Race+education+Alcohol+age+GAD,PTH)


tab_model(fit1,  fit3, fit4, fit5)
## Warning: Model matrix is rank deficient. Parameters TBInumbers_cat9+ were not
## estimable.

## Warning: Model matrix is rank deficient. Parameters TBInumbers_cat9+ were not
## estimable.

## Warning: Model matrix is rank deficient. Parameters TBInumbers_cat9+ were not
## estimable.

## Warning: Model matrix is rank deficient. Parameters TBInumbers_cat9+ were not
## estimable.
  HIT Cognetive Vestibular PHQ
Predictors Estimates CI p Estimates CI p Estimates CI p Estimates CI p
(Intercept) 60.53 57.64 – 63.43 <0.001 0.88 0.57 – 1.19 <0.001 0.76 0.44 – 1.08 <0.001 0.88 0.45 – 1.30 <0.001
TBI_severity_finalccMild TBI without LOC -0.15 -1.56 – 1.26 0.832 -0.16 -0.31 – -0.00 0.043 -0.01 -0.17 – 0.14 0.885 0.15 -0.06 – 0.35 0.160
TBI_severity_finalccModerate/Severe 2.14 0.43 – 3.85 0.014 0.29 0.10 – 0.47 0.002 0.44 0.26 – 0.63 <0.001 0.19 -0.06 – 0.44 0.135
TBI_severity_finalccNo TBI -4.65 -6.77 – -2.52 <0.001 -0.64 -0.87 – -0.41 <0.001 -0.56 -0.79 – -0.32 <0.001 -0.45 -0.76 – -0.14 0.005
TBInumbers_cat1 -4.44 -6.71 – -2.16 <0.001 -0.59 -0.83 – -0.34 <0.001 -0.70 -0.94 – -0.45 <0.001 -0.52 -0.85 – -0.19 0.002
TBInumbers_cat2 -2.88 -5.05 – -0.72 0.009 -0.55 -0.78 – -0.31 <0.001 -0.46 -0.70 – -0.23 <0.001 -0.44 -0.75 – -0.12 0.007
TBInumbers_cat3 -3.02 -5.31 – -0.74 0.010 -0.46 -0.70 – -0.21 <0.001 -0.53 -0.78 – -0.28 <0.001 -0.38 -0.72 – -0.05 0.024
TBInumbers_cat4 -2.10 -4.91 – 0.72 0.144 -0.25 -0.55 – 0.05 0.106 -0.24 -0.55 – 0.06 0.120 -0.55 -0.96 – -0.14 0.009
TBInumbers_cat4-8 -2.45 -4.46 – -0.43 0.017 -0.23 -0.45 – -0.02 0.034 -0.24 -0.46 – -0.02 0.033 -0.32 -0.61 – -0.02 0.035
BlowsYes 1.35 0.20 – 2.50 0.021 0.22 0.10 – 0.34 0.001 0.31 0.18 – 0.43 <0.001 0.21 0.04 – 0.37 0.017
Gender: Male -2.38 -3.46 – -1.29 <0.001 -0.01 -0.13 – 0.10 0.803 -0.00 -0.12 – 0.12 0.972 0.07 -0.09 – 0.23 0.390
Race: Black -2.45 -4.00 – -0.89 0.002 -0.01 -0.17 – 0.16 0.944 -0.04 -0.21 – 0.13 0.623 -0.00 -0.23 – 0.23 0.986
Race: Hispanic -1.14 -2.98 – 0.70 0.225 -0.22 -0.42 – -0.02 0.031 -0.00 -0.20 – 0.20 0.976 -0.06 -0.33 – 0.21 0.661
Race: Mixed -2.15 -4.23 – -0.06 0.044 -0.24 -0.46 – -0.01 0.037 -0.13 -0.35 – 0.10 0.275 -0.09 -0.39 – 0.22 0.568
educationLess than High School -0.52 -2.60 – 1.56 0.624 0.21 -0.02 – 0.43 0.071 0.32 0.09 – 0.55 0.006 0.33 0.03 – 0.64 0.031
educationMore than College -0.19 -1.51 – 1.13 0.777 0.09 -0.06 – 0.23 0.236 0.06 -0.08 – 0.21 0.383 -0.01 -0.20 – 0.18 0.924
educationSome College -0.07 -1.40 – 1.26 0.922 0.16 0.01 – 0.30 0.032 0.23 0.09 – 0.38 0.002 0.07 -0.13 – 0.26 0.499
Alcohol: Heavy 0.74 -0.69 – 2.16 0.312 0.33 0.18 – 0.49 <0.001 0.37 0.21 – 0.52 <0.001 0.13 -0.07 – 0.34 0.206
Alcohol: Moderate 0.24 -0.94 – 1.42 0.688 0.14 0.01 – 0.27 0.032 0.12 -0.01 – 0.25 0.067 -0.05 -0.23 – 0.12 0.541
age -0.03 -0.05 – -0.00 0.028 -0.00 -0.00 – 0.00 0.434 -0.00 -0.00 – 0.00 0.985 -0.00 -0.01 – 0.00 0.086
GAD 1.88 1.61 – 2.14 <0.001 0.27 0.24 – 0.30 <0.001 0.21 0.18 – 0.24 <0.001 0.74 0.70 – 0.78 <0.001
Observations 1004 1004 1004 1004
R2 / R2 adjusted 0.313 / 0.299 0.472 / 0.461 0.435 / 0.423 0.667 / 0.661
# install.packages("gtsummary")

#tbl_regression(fit1)

PTH$GAD
##    [1] 0 3 0 2 6 0 0 1 0 4 6 1 0 6 0 1 5 0 6 5 4 5 2 0 1 2 0 1 6 1 1 4 0 4 4 3 2
##   [38] 4 6 0 0 1 0 4 4 6 6 3 2 2 0 0 1 2 0 4 2 4 1 0 4 0 3 3 5 1 0 0 4 6 4 4 0 2
##   [75] 3 0 2 4 0 0 2 0 4 0 2 2 0 0 6 2 1 0 0 2 0 2 4 4 6 2 4 4 5 1 3 2 3 2 0 4 0
##  [112] 0 4 0 4 4 3 3 4 5 0 0 1 2 0 0 2 1 0 0 0 0 2 0 0 4 2 3 3 4 0 6 1 3 3 4 0 0
##  [149] 0 6 0 0 1 0 2 0 0 4 3 2 6 3 3 5 4 4 0 0 1 6 0 4 6 4 2 4 0 0 0 4 2 0 0 0 6
##  [186] 6 3 4 2 2 0 2 2 6 4 6 4 2 2 0 3 0 0 1 4 3 4 0 0 0 2 4 3 5 1 2 2 2 0 6 2 5
##  [223] 2 6 5 3 6 4 3 1 6 6 3 4 0 2 3 4 1 0 4 1 4 0 2 6 3 4 0 5 0 0 0 2 3 4 6 3 2
##  [260] 5 0 4 3 5 3 6 5 3 1 1 2 5 2 0 2 2 4 1 1 1 3 1 4 0 1 0 0 1 0 0 0 4 2 5 3 3
##  [297] 0 2 1 0 3 2 0 1 0 5 0 1 1 0 2 4 2 0 6 1 4 0 3 5 0 5 5 3 6 2 3 0 4 2 0 1 1
##  [334] 2 4 2 1 6 3 6 4 0 0 3 2 4 3 3 4 0 1 0 4 6 0 4 2 1 2 1 2 1 0 0 2 4 1 0 4 0
##  [371] 0 5 2 1 2 6 4 0 0 0 1 0 0 1 1 6 3 2 2 1 5 2 3 5 2 2 0 3 0 2 5 4 4 4 0 0 5
##  [408] 3 3 4 3 2 3 6 1 1 2 2 0 1 2 0 1 2 0 2 0 0 3 0 5 1 0 0 3 0 3 0 2 4 0 2 2 3
##  [445] 2 5 2 4 0 4 2 3 2 2 2 0 3 3 3 2 4 0 6 2 3 1 2 5 0 6 3 1 6 4 3 6 3 3 0 3 6
##  [482] 2 1 6 5 6 2 2 0 4 1 2 0 0 4 2 0 2 2 2 6 0 6 3 2 2 2 0 3 5 2 6 6 5 0 0 4 3
##  [519] 6 0 6 4 2 0 0 4 5 0 1 3 0 1 4 0 3 3 1 2 1 6 2 5 2 6 0 2 0 0 5 1 0 4 3 0 2
##  [556] 4 2 0 0 4 0 1 6 1 4 2 6 6 3 1 4 3 0 0 4 3 6 2 6 1 2 1 5 4 4 4 2 1 6 0 6 2
##  [593] 4 2 1 1 4 6 3 0 0 0 2 2 2 1 3 2 3 2 5 3 6 0 3 4 2 2 1 4 0 0 6 2 0 0 2 0 0
##  [630] 6 4 5 2 3 1 0 6 3 6 4 1 3 3 1 0 5 2 6 0 6 3 6 4 3 0 4 6 0 0 0 2 0 0 4 3 4
##  [667] 2 3 3 4 4 2 6 3 0 0 2 0 4 2 2 6 1 0 0 1 6 1 1 0 6 4 0 2 4 3 3 0 1 0 4 6 6
##  [704] 4 6 4 4 2 0 2 0 2 4 6 4 4 5 3 2 0 6 6 2 6 2 2 2 6 2 6 2 6 6 1 6 0 4 0 4 5
##  [741] 3 6 0 0 4 1 0 4 6 1 0 2 2 0 1 3 2 3 4 6 5 5 0 1 1 6 0 4 2 2 0 0 2 0 4 4 1
##  [778] 2 6 2 0 0 3 3 3 0 0 1 1 0 2 0 0 2 0 0 2 1 2 4 0 0 0 0 4 4 4 2 2 0 4 3 2 4
##  [815] 5 6 2 6 2 0 0 4 6 6 5 0 4 1 2 5 4 4 4 6 2 2 3 5 3 3 6 0 0 3 2 3 2 4 2 3 4
##  [852] 3 0 1 2 3 4 4 0 4 3 5 0 0 4 4 3 0 2 3 0 1 3 4 3 2 0 2 0 0 1 0 2 0 1 3 4 0
##  [889] 2 1 6 3 4 1 2 6 4 3 5 4 2 2 0 4 1 1 0 2 6 0 1 3 1 2 0 2 2 2 0 0 1 0 0 0 6
##  [926] 4 2 2 0 3 1 0 2 4 0 4 3 0 2 2 1 0 0 6 3 1 4 2 3 3 0 2 2 2 2 3 2 4 1 5 0 3
##  [963] 4 3 6 1 3 2 2 6 0 6 0 5 3 6 6 6 2 0 2 4 4 0 5 0 0 0 6 0 3 0 0 0 0 4 3 1 2
## [1000] 0 3 6 3 2
#fit1 %>%
  #tbl_regression(
    #pvalue_fun = ~style_pvalue(.x, digits = 2),
  #) %>% 
  #bold_p(t = 0.10) %>%
  #bold_labels() %>%
  #italicize_levels()


#tab_model(fit1, fit2, fit3, fit4, fit5)
#set_gtsummary_theme(theme_gtsummary_compact(set_theme = TRUE))
#t1 <- tbl_regression(fit1) %>% bold_p() #t1 is just an object containing the model output.
#t2 <- tbl_regression(fit2) %>% bold_p() 
#theme_gtsummary_compact()

#t1



fit6=lm(HIT~TBI_severity_finalcc+TBInumbers+Blows,PTH)
fit7=lm(pain_level~TBI_severity_finalcc+TBInumbers+Blows,PTH)
fit8=lm(Cognetive~TBI_severity_finalcc+TBInumbers+Blows,PTH)
fit9=lm(Vestibular~TBI_severity_finalcc+TBInumbers+Blows,PTH)
fit10=lm(PHQ~TBI_severity_finalcc+TBInumbers+Blows,PTH)
tab_model(fit6, fit7, fit8, fit9, fit10)
  HIT Headache Severity Cognetive Vestibular PHQ
Predictors Estimates CI p Estimates CI p Estimates CI p Estimates CI p Estimates CI p
(Intercept) 57.47 55.96 – 58.98 <0.001 4.91 4.50 – 5.32 <0.001 0.92 0.75 – 1.09 <0.001 0.73 0.56 – 0.89 <0.001 1.66 1.35 – 1.97 <0.001
TBI severity finalcc
[Mild TBI without LOC]
-0.58 -2.12 – 0.96 0.458 -0.16 -0.57 – 0.26 0.466 -0.26 -0.43 – -0.09 0.003 -0.09 -0.26 – 0.07 0.277 -0.03 -0.35 – 0.29 0.848
TBI severity finalcc
[Moderate/Severe]
2.44 0.57 – 4.32 0.011 0.74 0.23 – 1.25 0.005 0.42 0.21 – 0.64 <0.001 0.58 0.38 – 0.78 <0.001 0.50 0.11 – 0.89 0.011
TBI severity finalcc [No
TBI]
-0.92 -2.60 – 0.77 0.287 -0.03 -0.49 – 0.43 0.891 -0.12 -0.31 – 0.07 0.229 0.02 -0.16 – 0.20 0.830 -0.07 -0.42 – 0.28 0.705
TBInumbers 0.40 0.23 – 0.57 <0.001 0.10 0.06 – 0.15 <0.001 0.07 0.05 – 0.09 <0.001 0.07 0.05 – 0.09 <0.001 0.09 0.05 – 0.12 <0.001
Blows [Yes] 2.22 1.02 – 3.43 <0.001 0.50 0.17 – 0.82 0.003 0.46 0.32 – 0.59 <0.001 0.53 0.40 – 0.66 <0.001 0.71 0.46 – 0.96 <0.001
Observations 1004 1004 1004 1004 1004
R2 / R2 adjusted 0.116 / 0.111 0.095 / 0.090 0.253 / 0.249 0.285 / 0.281 0.141 / 0.137
# https://rpubs.com/lindangulopez/708163
##                          Blows TBI_severity_finalcc TBInumbers
## Blows                1.0000000            0.3786597  0.4099976
## TBI_severity_finalcc 0.3786597            1.0000000  0.6309605
## TBInumbers           0.4099976            0.6309605  1.0000000
## 
## Attaching package: 'reshape'
## The following objects are masked from 'package:tidyr':
## 
##     expand, smiths
## The following object is masked from 'package:dplyr':
## 
##     rename
## The following object is masked from 'package:data.table':
## 
##     melt
## Warning in type.convert.default(X[[i]], ...): 'as.is' should be specified by the
## caller; using TRUE

## Warning in type.convert.default(X[[i]], ...): 'as.is' should be specified by the
## caller; using TRUE
##                     X1                   X2     value
## 1                Blows                Blows 1.0000000
## 2 TBI_severity_finalcc                Blows 0.3786597
## 3           TBInumbers                Blows 0.4099976
## 4                Blows TBI_severity_finalcc 0.3786597
## 5 TBI_severity_finalcc TBI_severity_finalcc 1.0000000
## 6           TBInumbers TBI_severity_finalcc 0.6309605

## Warning in type.convert.default(X[[i]], ...): 'as.is' should be specified by the
## caller; using TRUE

## Warning in type.convert.default(X[[i]], ...): 'as.is' should be specified by the
## caller; using TRUE

```

https://biostats.w.uib.no/creating-a-multiple-pie-donut-chart/