Final analytic model of psychiatric morbidity using Poisson regression (n = 705)

Variables within hierarchical modelsIRRSE95% CIPOverall model statistics
Pseudo R2AICLR test
Base covariates< 0.00010.0612,560Base
Age0.980.0020.98–0.99< 0.001
Gender (ref. female)
    Male0.900.080.76–1.060.20
    Transgender, non-binary, gender non-conforming, or other0.970.150.71–1.320.86
Education (ref. ≤ HS)
    Associate degree or some college0.930.110.75–1.160.54
    Bachelor’s degree0.960.110.77–1.190.70
    Master’s degree or above0.930.110.75–1.170.54
MS subtype (ref. RRMS)
    PPMS1.270.160.99–1.620.06
    SPMS1.101.000.92–1.300.30
    PRMS1.340.270.91–1.980.14
    Unsure1.240.180.94–1.640.14
Pain med count1.100.021.06–1.15< 0.001
Childhood stressors< 0.00010.0902,485< 0.0001
Child stressor count0.970.050.88–1.070.57
Child stressor severity1.020.011.00–1.040.11
Adult stressors< 0.00010.1162,420< 0.0001
Adult stressor count1.000.0040.99–1.0040.27
Adult stressor severity1.010.0021.006–1.014< 0.0001

≤ HS: high school equivalency or below; IRR: incident rate ratio; PPMS: primary progressive MS; PRMS: progressive-relapsing MS; ref.: reference; SPMS: secondary progressive MS. Blank cells indicate not applicable to individual variables