Here are some descriptive statistics for the variables you used in your analyses. These statistics should match the values in SPSS (within rounding).
These are the dependent variables.
faren[, varsDep]
7 Variables 3721 Observations
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
LabResWBCimp Format:%12.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2802 919 237 6.539 3.700 4.100 5.100 6.199 7.600 9.299 10.500
lowest : 1.8 1.9 2.2 2.3 2.4, highest: 16.5 16.9 17.0 18.5 20.1
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
LabResCRP Format:%9.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2557 1164 333 4.95 0.2000 0.3000 0.7999 2.0996 5.5996 11.6992 18.6180
lowest : 4.999e-03 2.000e-02 3.000e-02 3.999e-02 5.000e-02, highest: 1.160e+02 1.410e+02 1.520e+02 1.530e+02 1.630e+02
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
NLratio Format:%12.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2745 976 2629 2.124 0.8066 0.9964 1.3300 1.8451 2.5639 3.3805 4.1123
lowest : 0.1642 0.1744 0.2353 0.2691 0.2781, highest: 10.7143 11.2857 16.8537 26.4438 90.0000
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
IL1b : Aushon hIL1b (pg/ml) Format:%12.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
244 3477 59 2.562 0.200 0.230 0.400 0.800 1.900 4.570 8.535
lowest : 0.03 0.10 0.20 0.30 0.40, highest: 13.30 18.30 40.40 46.50 111.30
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
IL6 : Aushon hIL6 (pg/ml) Format:%12.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
259 3462 158 13.75 2.70 3.70 5.40 8.50 14.15 30.56 43.96
lowest : 1.4 1.7 1.8 1.9 2.0, highest: 79.6 86.1 89.2 92.3 136.1
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
IL10 : Aushon hIL10 (pg/ml) Format:%12.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
259 3462 55 2.707 0.90 1.00 1.20 1.60 2.30 3.64 6.00
lowest : 0.2 0.3 0.5 0.6 0.7, highest: 13.4 16.0 24.6 39.3 100.0
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
TNFa : Aushon hTNFa (pg/ml) Format:%12.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
244 3477 153 31.13 0.900 1.400 2.975 5.850 16.575 45.470 74.500
lowest : 0.4 0.5 0.7 0.8 0.9, highest: 301.8 304.6 666.0 699.7 1736.1
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
These are the independent variables (CES-D subscales).
faren[, varsInd]
5 Variables 3721 Observations
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
CESimp Format:%12.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2802 919 123 15.13 1.0 3.0 6.0 13.0 21.0 31.9 38.0
lowest : 0 1 2 3 4, highest: 53 54 55 57 59
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
pa Format:%8.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2731 990 13 9.511 4 5 8 10 12 12 12
0 1 2 3 4 5 6 7 8 9 10 11 12
Frequency 11 12 14 50 80 111 137 168 243 286 331 372 916
% 0 0 1 2 3 4 5 6 9 10 12 14 34
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
sc Format:%8.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2736 985 22 6.88 0 1 3 6 10 13 15
lowest : 0 1 2 3 4, highest: 17 18 19 20 21
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
ip Format:%8.0g
n missing unique Mean
2735 986 7 1.035
0 1 2 3 4 5 6
Frequency 1433 479 422 207 116 45 33
% 52 18 15 8 4 2 1
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
na Format:%8.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2733 988 22 4.749 0 0 0 3 8 12 15
lowest : 0 1 2 3 4, highest: 17 18 19 20 21
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
These are the covariates.
faren[, varsCov]
12 Variables 3721 Observations
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Age0 Format:%8.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2802 919 35 48.06 32 35 41 48 56 60 62
lowest : 30 31 32 33 34, highest: 60 61 62 63 64
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Sex Format:%8.0g
n missing unique Sum Mean
2802 919 2 1222 0.4361
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Race Format:%8.0g
n missing unique Sum Mean
2802 919 2 1617 0.5771
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
PovStat Format:%8.0g
n missing unique Sum Mean
2802 919 2 1142 0.4076
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
CVDclusterdich Format:%8.0g
n missing unique Sum Mean
2802 919 2 154 0.05496
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
DMMclusterdich Format:%8.0g
n missing unique Sum Mean
2802 919 2 793 0.283
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
PhysBMI Format:%9.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2802 919 1515 30.07 19.71 21.30 24.44 28.63 34.54 41.03 45.48
lowest : 14.36 14.70 14.77 14.87 15.21, highest: 59.94 60.55 60.89 61.89 62.46
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
MedHxAlcStatusNFC Format:%8.0g
n missing unique Mean
2802 919 3 1.294
0 (660, 24%), 1 (658, 23%), 2 (1484, 53%)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
MedHxCigaretteStatusNFC Format:%8.0g
n missing unique Mean
2802 919 3 1.214
0 (702, 25%), 1 (799, 29%), 2 (1301, 46%)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
MedHxHardDrugNFC Format:%8.0g
n missing unique Mean
2802 919 3 0.3026
0 (2109, 75%), 1 (538, 19%), 2 (155, 6%)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
PhysMeanSBPsitimp Format:%12.0g
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2802 919 199 120.7 96.0 100.0 108.2 119.0 131.0 144.5 152.0
lowest : 79.0 80.0 81.0 82.5 83.0, highest: 187.0 189.0 190.0 194.0 194.5
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
ThesisDrug Format:%8.0g
n missing unique Sum Mean
2652 1069 2 1116 0.4208
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
First examine the complete models separately for each dependent variable for each depression subscales plus covariates.
CES-D (imputed):
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.137478 *** 0.065976 *** 0.0297856 *** 0.036081 0.11172 * 0.085363 0.04828
CESimp 0.002313 0.009406 -0.0009309 -0.030579 0.11327 0.004405 0.24362
Age0 -0.020486 *** 0.016092 -0.0038427 0.023249 0.09051 -0.004067 0.32591
Sex -0.175148 * -0.820300 * 0.0073667 -0.555032 2.26798 -0.023833 0.37363
Race -0.963470 *** 0.852616 * -0.5058570 *** -0.068909 0.37034 -0.876341 -11.33052
PovStat 0.184905 * 0.627629 0.1700536 -0.195168 0.41443 0.134319 -2.98122
CVDclusterdich 0.749587 *** 3.127713 *** 1.0062671 *** -0.465275 1.38296 0.010935 -6.42834
DMMclusterdich 0.159134 -0.114836 0.1897409 1.026302 0.20707 -0.288212 -5.14544
PhysBMI 0.038749 *** 0.271544 *** 0.0005589 -0.015476 0.07791 -0.021701 0.17765
MedHxAlcStatusNFC -0.027671 -0.124476 0.0750011 0.451984 -0.98435 0.134008 1.09644
MedHxCigaretteStatusNFC 0.670772 *** 0.702424 ** 0.0429583 0.555476 3.08119 * -0.048145 1.31487
MedHxHardDrugNFC -0.103235 -0.241081 -0.1608947 * -0.057115 4.17025 * 1.633428 ** 21.22386 *
PhysMeanSBPsitimp 0.005916 * -0.001091 0.0005254 -0.007591 -0.06646 0.006743 0.11956
ThesisDrug 0.134964 0.331546 -0.0879350 -0.427625 2.44931 0.102784 -5.30282
PA:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.1381689 *** 0.066798 *** 0.0298064 *** 0.04339 0.11651 * 0.092293 0.05636
pa -0.0005724 -0.034276 -0.0022900 0.15809 -0.13569 0.080902 -1.04170
Age0 -0.0206837 *** 0.012027 -0.0032337 0.02787 0.09699 -0.001367 0.38755
Sex -0.1399624 -0.733253 0.0060224 -0.45917 1.96996 -0.047942 1.41498
Race -0.9857406 *** 0.750673 -0.5027818 *** 0.05087 0.08135 -0.864921 -6.39963
PovStat 0.1911382 * 0.747056 0.1842661 * -0.53406 0.83625 0.089924 -7.69643
CVDclusterdich 0.7925026 *** 3.372134 *** 1.0188734 *** -0.57489 1.11585 -0.097942 -5.78136
DMMclusterdich 0.1578554 -0.190817 0.2017780 * 1.13394 0.09344 -0.237127 -4.59205
PhysBMI 0.0385471 *** 0.276067 *** -0.0006598 -0.02093 0.06162 -0.026035 0.23786
MedHxAlcStatusNFC -0.0422856 -0.181727 0.0740855 0.48151 -1.27657 0.126932 -0.77794
MedHxCigaretteStatusNFC 0.6775566 *** 0.734256 ** 0.0407247 0.63527 3.39081 * 0.016802 3.76884
MedHxHardDrugNFC -0.0821097 -0.217621 -0.1641722 * -0.01075 5.04601 * 1.786988 *** 22.94611 *
PhysMeanSBPsitimp 0.0049923 * -0.002095 0.0002964 -0.01106 -0.06703 0.005457 -0.03404
ThesisDrug 0.1477656 0.367675 -0.0954957 -0.28893 3.30786 0.247849 -2.20998
SC:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.140095 *** 0.06655 *** 0.0298773 *** 0.043592 0.12157 * 0.090782 0.05563
sc 0.015096 0.05197 -0.0046764 -0.086924 0.20036 -0.006117 -0.07748
Age0 -0.020115 *** 0.01315 -0.0035031 0.026241 0.10158 -0.001592 0.39632
Sex -0.132929 -0.71064 0.0010659 -0.572307 2.31816 -0.043191 1.10868
Race -0.985077 *** 0.74866 -0.5046415 *** 0.117818 0.42520 -0.734672 -8.28380
PovStat 0.166046 * 0.69333 0.1906610 * -0.496544 0.66555 0.059215 -7.15151
CVDclusterdich 0.781522 *** 3.34484 *** 1.0237875 *** -0.476627 0.96353 -0.053675 -6.55158
DMMclusterdich 0.159474 -0.20926 0.2063970 * 1.136207 -0.14084 -0.288738 -3.92006
PhysBMI 0.039490 *** 0.27472 *** -0.0005381 -0.015689 0.02976 -0.028872 0.27748
MedHxAlcStatusNFC -0.040287 -0.17116 0.0727848 0.430807 -1.09614 0.136883 -0.94097
MedHxCigaretteStatusNFC 0.670026 *** 0.69867 ** 0.0446665 0.688091 3.07548 * -0.023367 4.47038
MedHxHardDrugNFC -0.090722 -0.22853 -0.1609817 0.009295 4.78220 * 1.732965 ** 23.90052 *
PhysMeanSBPsitimp 0.005043 * -0.00199 0.0003414 -0.010093 -0.07357 0.004464 -0.01633
ThesisDrug 0.125833 0.34629 -0.0935062 -0.303125 3.56839 0.281860 -2.56314
IP:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.137154 *** 0.066748 *** 0.0303086 *** 0.04623 0.13411 ** 0.133637 ** 0.05647
ip -0.008523 -0.056198 0.0295217 0.39108 1.66833 * 0.628764 ** 2.17142
Age0 -0.020896 *** 0.010050 -0.0028413 0.02633 0.09682 -0.001622 0.38783
Sex -0.141219 -0.745588 0.0033020 -0.31121 2.39082 0.130106 1.85051
Race -0.983007 *** 0.730515 -0.5035991 *** 0.17231 -0.23996 -0.827757 -7.75641
PovStat 0.191827 * 0.770064 0.1758104 -0.81149 -0.05656 -0.327288 -8.42760
CVDclusterdich 0.775558 *** 3.345017 *** 1.0306013 *** -0.20044 2.25022 0.414743 -5.03271
DMMclusterdich 0.157647 -0.184606 0.2031684 * 1.16023 0.45651 -0.146345 -3.87786
PhysBMI 0.038627 *** 0.276396 *** -0.0005942 -0.01726 0.06152 -0.023032 0.21954
MedHxAlcStatusNFC -0.040787 -0.183932 0.0767796 0.53527 -1.17308 0.181769 -0.64822
MedHxCigaretteStatusNFC 0.674141 *** 0.744812 ** 0.0408090 0.63560 3.62770 ** 0.069399 4.28996
MedHxHardDrugNFC -0.080335 -0.197550 -0.1695413 * -0.41851 4.20073 * 1.362894 ** 22.49179 *
PhysMeanSBPsitimp 0.005191 * -0.001405 0.0001606 -0.01531 -0.07268 0.001988 -0.03238
ThesisDrug 0.146612 0.377555 -0.0968781 -0.42984 2.91230 0.064584 -2.65127
NA:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.1390581 *** 6.706e-02 *** 0.0300196 *** 0.04227 0.11886 * 0.089411 0.05724
na 0.0007375 4.957e-05 -0.0042782 -0.06787 0.15879 -0.003758 0.63573
Age0 -0.0208993 *** 1.080e-02 -0.0035756 0.02706 0.09834 -0.001510 0.39449
Sex -0.1447305 -7.502e-01 -0.0006358 -0.50117 2.14289 -0.040601 1.98267
Race -0.9815047 *** 7.372e-01 -0.5071111 *** 0.12008 0.18447 -0.774247 -6.35112
PovStat 0.1830126 * 7.431e-01 0.1876524 * -0.48116 0.59019 0.048249 -8.49644
CVDclusterdich 0.8203489 *** 3.447e+00 *** 1.0287200 *** -0.54737 1.12986 -0.048167 -6.18520
DMMclusterdich 0.1671707 -1.690e-01 0.2072532 * 1.10199 0.11983 -0.264849 -4.41645
PhysBMI 0.0384835 *** 2.765e-01 *** -0.0006417 -0.01737 0.05850 -0.024172 0.22713
MedHxAlcStatusNFC -0.0411079 -1.809e-01 0.0729537 0.45346 -1.18385 0.134053 -0.31465
MedHxCigaretteStatusNFC 0.6793986 *** 7.513e-01 ** 0.0450834 0.65479 3.27719 * -0.002091 3.34789
MedHxHardDrugNFC -0.0840928 -2.065e-01 -0.1602083 -0.04024 4.87411 * 1.728970 *** 22.62728 *
PhysMeanSBPsitimp 0.0050809 * -1.785e-03 0.0003519 -0.01256 -0.06646 0.004657 -0.02808
ThesisDrug 0.1430323 3.719e-01 -0.0915830 -0.18543 3.02468 0.234436 -3.43697
As a first cut, eliminate covariates that are nonsignificant in any of the analyses: sex, MedHxAlcStatusNFC, ThesisDrug.
CES-D (imputed):
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.127264 *** 0.064324 *** 0.0290107 *** 0.026513 0.07322 * 0.03770 0.03142
CESimp 0.004829 0.016767 -0.0009213 -0.019318 0.07093 -0.02988 -0.42634
Age0 -0.017226 *** 0.021638 -0.0054686 -0.084055 0.09253 -0.03764 -0.98668
Race -0.937622 *** 0.894698 * -0.5196067 *** 0.939949 1.45610 -0.65369 -3.23248
PovStat 0.238888 ** 0.662280 0.1523641 0.713380 2.09515 0.97166 16.65367
CVDclusterdich 0.763569 *** 3.169588 *** 0.9728182 *** -0.929840 0.57702 -0.69081 -20.49824
DMMclusterdich 0.187433 * -0.010532 0.1579286 1.022803 -0.10909 -0.77713 -17.41271
PhysBMI 0.040055 *** 0.284232 *** -0.0018820 -0.098325 0.03659 -0.05469 -1.07263
MedHxCigaretteStatusNFC 0.607507 *** 0.627486 ** 0.0498288 -0.617946 3.12756 * 0.56913 6.53363
MedHxHardDrugNFC -0.137863 -0.407089 -0.1391544 -0.727530 3.29662 1.09140 7.51102
PhysMeanSBPsitimp 0.005699 * -0.002619 0.0010483 -0.004503 -0.01313 0.05216 1.02330
PA:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.127734 *** 0.065931 *** 0.0289502 *** 0.04097 0.07823 * 0.04235 0.02817
pa -0.005809 -0.048915 -0.0024789 -0.44307 0.05961 0.23628 -1.20707
Age0 -0.017519 *** 0.017188 -0.0049966 -0.08606 0.09126 -0.03663 -0.96255
Race -0.958626 *** 0.788499 * -0.5166438 *** 1.71391 0.98869 -0.70733 4.66688
PovStat 0.252155 ** 0.816236 * 0.1656577 0.18458 2.32954 0.90755 9.84843
CVDclusterdich 0.821104 *** 3.461995 *** 0.9839249 *** -0.78612 0.59941 -0.89391 -19.84341
DMMclusterdich 0.190759 * -0.071360 0.1686822 1.00782 0.13187 -0.68921 -17.04881
PhysBMI 0.040167 *** 0.290893 *** -0.0031228 -0.10006 0.04791 -0.06620 -1.10288
MedHxCigaretteStatusNFC 0.613179 *** 0.664660 ** 0.0460010 -0.73721 3.50338 * 0.56653 6.19034
MedHxHardDrugNFC -0.114196 -0.381543 -0.1410681 -1.21292 3.96174 1.18870 5.19312
PhysMeanSBPsitimp 0.004875 * -0.003675 0.0008181 -0.01234 -0.01166 0.05534 0.90487
SC:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.130362 *** 0.065949 *** 0.0289901 *** 0.030779 0.08005 * 0.03702 0.02852
sc 0.020713 * 0.068106 -0.0040891 -0.142271 0.16455 -0.02703 -1.08501
Age0 -0.017261 *** 0.018048 -0.0052032 -0.078396 0.09362 -0.03569 -0.92427
Race -0.961024 *** 0.776682 * -0.5179066 *** 1.198618 1.44699 -0.44706 2.66002
PovStat 0.219941 ** 0.742555 0.1722905 * 0.701047 2.02726 0.78825 12.50223
CVDclusterdich 0.800834 *** 3.417430 *** 0.9886404 *** -0.955671 0.61669 -0.74913 -20.73491
DMMclusterdich 0.185783 * -0.101214 0.1734547 1.261888 -0.11680 -0.79731 -16.13782
PhysBMI 0.040662 *** 0.289141 *** -0.0029364 -0.108987 0.02132 -0.06629 -1.14829
MedHxCigaretteStatusNFC 0.605335 *** 0.628538 * 0.0488621 -0.550516 3.21166 * 0.53604 7.45049
MedHxHardDrugNFC -0.120994 -0.386977 -0.1390344 -0.646397 3.72783 1.04301 7.26096
PhysMeanSBPsitimp 0.004905 * -0.003541 0.0008524 -0.005766 -0.01624 0.05310 0.94024
IP:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.126751 *** 0.065773 *** 0.0293485 *** 0.02683 0.08517 * 0.03798 0.03179
ip 0.004128 -0.040223 0.0230748 -0.09479 1.15461 0.23736 -7.46029
Age0 -0.017713 *** 0.015093 -0.0047446 -0.08256 0.08773 -0.03721 -0.89910
Race -0.958968 *** 0.760381 * -0.5183260 *** 1.26343 0.95600 -0.48284 3.56162
PovStat 0.251924 ** 0.845156 * 0.1586378 0.50336 1.63228 0.59927 14.79419
CVDclusterdich 0.804306 *** 3.438319 *** 0.9954000 *** -1.05973 1.38307 -0.60281 -25.30530
DMMclusterdich 0.190514 * -0.060717 0.1704192 1.18092 0.23399 -0.77080 -17.48727
PhysBMI 0.040345 *** 0.291653 *** -0.0030794 -0.10476 0.04287 -0.06375 -1.07974
MedHxCigaretteStatusNFC 0.611108 *** 0.675491 ** 0.0470495 -0.66223 3.63725 * 0.56553 5.59793
MedHxHardDrugNFC -0.113248 -0.365156 -0.1450448 -0.76462 3.37810 0.90365 9.83934
PhysMeanSBPsitimp 0.005034 * -0.003007 0.0007096 -0.00669 -0.01561 0.05220 0.95452
NA:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.128682 *** 0.066184 *** 0.0291543 *** 0.029278 0.07893 * 0.03912 0.02904
na 0.006401 0.018617 -0.0041373 -0.091012 0.09319 -0.07639 -1.01833
Age0 -0.017747 *** 0.016273 -0.0052545 -0.079192 0.08874 -0.03438 -0.91940
Race -0.954217 *** 0.777365 * -0.5202044 *** 1.112806 1.18646 -0.57515 1.68193
PovStat 0.239236 ** 0.797906 * 0.1701500 0.643423 2.08906 0.89473 12.90349
CVDclusterdich 0.843587 *** 3.525558 *** 0.9945156 *** -1.024987 0.64463 -0.76551 -20.76334
DMMclusterdich 0.196298 * -0.056615 0.1749402 1.272368 0.03539 -0.74069 -15.79515
PhysBMI 0.040055 *** 0.291280 *** -0.0029766 -0.100931 0.04467 -0.05911 -1.09326
MedHxCigaretteStatusNFC 0.614040 *** 0.675819 ** 0.0494564 -0.570089 3.41952 * 0.59897 7.27153
MedHxHardDrugNFC -0.117581 -0.378401 -0.1384946 -0.704502 3.80754 1.10500 7.30684
PhysMeanSBPsitimp 0.004935 * -0.003454 0.0008595 -0.007567 -0.01225 0.05287 0.92385
Now look for covariates that are were significant before but now are all nonsignificant and covariates that have few significant but small effects on the outcomes: MedHxHardDrugNFC (all nonsignificant); Age0, DMMclusterdich, PhysMeanSBPsitimp (significant only for WBC).
CES-D (imputed):
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.120421 *** 0.06341 *** 0.0267421 *** 0.01694 0.061784 * 0.01915 0.01464
CESimp 0.005655 0.01438 -0.0006556 -0.02669 0.098614 -0.02303 -0.40949
Race -0.932531 *** 0.84645 * -0.5292395 *** 0.64722 2.045286 -0.19441 3.58284
PovStat 0.252368 ** 0.60472 0.1498711 0.84727 1.940389 1.13888 20.42244
CVDclusterdich 0.709587 *** 3.33203 *** 0.9797630 *** -1.21163 0.713984 -0.98090 -27.26742
PhysBMI 0.045205 *** 0.28687 *** 0.0010261 -0.08947 -0.004028 -0.04382 -0.71756
MedHxCigaretteStatusNFC 0.589134 *** 0.56438 * 0.0319658 -0.58115 3.648315 ** 0.68670 6.09131
PA:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.121355 *** 0.06523 *** 0.0266549 *** 0.02836 0.062502 * 0.02269 0.01438
pa -0.009891 -0.04064 -0.0036816 -0.39845 -0.070401 0.19356 -1.35603
Race -0.951663 *** 0.73749 -0.5268260 *** 1.28345 1.663492 -0.16039 10.95103
PovStat 0.268344 *** 0.75819 0.1633762 0.25605 2.288373 1.08466 12.72428
CVDclusterdich 0.768315 *** 3.57606 *** 1.0002064 *** -1.07485 0.743601 -1.12119 -26.33028
PhysBMI 0.044468 *** 0.29208 *** -0.0002799 -0.09179 0.004418 -0.05252 -0.78222
MedHxCigaretteStatusNFC 0.597930 *** 0.60411 * 0.0280764 -0.78042 4.179416 ** 0.69729 5.59464
SC:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.12405 *** 0.06520 *** 0.0266490 *** 0.02135 0.06584 * 0.018970 0.01447
sc 0.02247 * 0.06386 -0.0032968 -0.15474 0.22511 -0.009996 -0.99548
Race -0.95671 *** 0.73008 -0.5283948 *** 0.90578 1.99336 0.046595 9.39369
PovStat 0.23440 ** 0.68439 0.1699939 0.84471 1.93226 0.945694 15.55234
CVDclusterdich 0.74799 *** 3.53447 *** 1.0041489 *** -1.17599 0.64957 -1.010221 -27.02894
PhysBMI 0.04502 *** 0.29012 *** -0.0000699 -0.10056 -0.02577 -0.053929 -0.81871
MedHxCigaretteStatusNFC 0.58904 *** 0.56622 * 0.0312126 -0.50006 3.81270 ** 0.629777 7.15958
IP:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.12023 *** 0.06517 *** 0.0269963 *** 0.01687 0.073135 ** 0.02134 0.01697
ip 0.01149 -0.05881 0.0223855 -0.18774 1.402769 0.31650 -6.43804
Race -0.95384 *** 0.71779 -0.5299951 *** 0.97595 1.353776 -0.04108 10.68610
PovStat 0.26806 *** 0.79374 * 0.1565248 0.65569 1.522073 0.73704 17.50063
CVDclusterdich 0.75230 *** 3.53673 *** 1.0129123 *** -1.37800 1.679851 -0.79554 -31.25236
PhysBMI 0.04480 *** 0.29293 *** -0.0002115 -0.09560 0.002443 -0.05018 -0.78730
MedHxCigaretteStatusNFC 0.59692 *** 0.61632 ** 0.0289149 -0.64406 4.264046 ** 0.65734 5.80364
NA:
LabResWBCimp LabResCRP NLratio IL1b IL6 IL10 TNFa
Rsq 0.122072 *** 0.06551 *** 2.680e-02 *** 0.01934 0.064412 * 0.02095 0.01544
na 0.007921 0.01468 -3.650e-03 -0.09822 0.147940 -0.07135 -1.12047
Race -0.948777 *** 0.72982 -5.304e-01 *** 0.78679 1.796755 -0.08049 8.12234
PovStat 0.255229 ** 0.74328 1.685e-01 0.77140 2.024394 1.07883 16.31908
CVDclusterdich 0.790778 *** 3.64101 *** 1.011e+00 *** -1.25508 0.697325 -1.01169 -27.00807
PhysBMI 0.044502 *** 0.29264 *** -6.765e-05 -0.09173 -0.002421 -0.04710 -0.77107
MedHxCigaretteStatusNFC 0.598819 *** 0.61539 ** 3.214e-02 -0.52657 4.041076 ** 0.71901 7.25185
Based on the preceding analyses it doesn't look like overall symptoms of depression or any of the subscales are associated with inflammatory markers. However, there are widely different numbers of participants with complete data for WBC, CRP, and NLratio (1st set: N ~ 2400) compared with the IL markers and TNFa (2nd set: N ~ 190). Rather than restricting the entire sample to the fewest numbers of participants, we can divide the outcomes into two sets and examine them separately. For comparability, I will use only those with complete data on WBC, CRP, and NLratio in the first set and only those with complete data on the IL markers and TNFa in the second set.
But before we look at these associations, we should check the distributions of the inflammatory marks.
Looks like these distributions are pretty skewed. Let's see what a log transform does.
Much better! These distributions pass the eyeball test. Let's do the analyses on complete data for the 1st and 2nd set with log transformed inflammatory markers.
f1
23 Variables 2407 Observations
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
LabResWBCimp
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 131 6.568 3.700 4.100 5.000 6.199 7.800 9.398 10.699
lowest : 1.9 2.2 2.3 2.4 2.5, highest: 15.7 16.5 16.9 18.5 20.1
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
LabResCRP
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 324 4.957 0.2000 0.3000 0.7999 2.1997 5.5996 10.9391 18.5977
lowest : 4.999e-03 2.000e-02 3.000e-02 3.999e-02 5.000e-02, highest: 1.160e+02 1.410e+02 1.520e+02 1.530e+02 1.630e+02
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
NLratio
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 2318 2.121 0.8045 0.9955 1.3308 1.8362 2.5564 3.3397 4.0944
lowest : 0.1642 0.1744 0.2353 0.2691 0.2781, highest: 10.6738 10.7143 11.2857 26.4438 90.0000
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
CESimp
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 57 15.12 1 3 6 13 22 32 38
lowest : 0 1 2 3 4, highest: 52 53 54 55 59
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
pa
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 13 9.501 4 5 8 10 12 12 12
0 1 2 3 4 5 6 7 8 9 10 11 12
Frequency 10 11 12 43 72 98 122 153 214 245 294 329 804
% 0 0 0 2 3 4 5 6 9 10 12 14 33
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
sc
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 22 6.857 0 1 4 6 10 13 15
lowest : 0 1 2 3 4, highest: 17 18 19 20 21
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
ip
n missing unique Mean
2407 0 7 1.034
0 1 2 3 4 5 6
Frequency 1261 416 375 190 101 36 28
% 52 17 16 8 4 1 1
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
na
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 22 4.731 0 0 0 3 8 12 15
lowest : 0 1 2 3 4, highest: 17 18 19 20 21
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Age0
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 35 48.07 32 35 41 48 56 61 62
lowest : 30 31 32 33 34, highest: 60 61 62 63 64
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Sex
n missing unique Sum Mean
2407 0 2 1056 0.4387
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Race
n missing unique Sum Mean
2407 0 2 1366 0.5675
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
PovStat
n missing unique Sum Mean
2407 0 2 945 0.3926
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
CVDclusterdich
n missing unique Sum Mean
2407 0 2 139 0.05775
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
DMMclusterdich
n missing unique Sum Mean
2407 0 2 703 0.2921
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
PhysBMI
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 1375 30.13 19.82 21.31 24.54 28.70 34.61 41.02 45.37
lowest : 14.36 14.77 14.87 15.62 15.79, highest: 56.88 57.94 58.45 59.94 60.89
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
MedHxAlcStatusNFC
n missing unique Mean
2407 0 3 1.375
0 (453, 19%), 1 (599, 25%), 2 (1355, 56%)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
MedHxCigaretteStatusNFC
n missing unique Mean
2407 0 3 1.22
0 (583, 24%), 1 (712, 30%), 2 (1112, 46%)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
MedHxHardDrugNFC
n missing unique Mean
2407 0 3 0.3074
0 (1803, 75%), 1 (468, 19%), 2 (136, 6%)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
PhysMeanSBPsitimp
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 177 120.6 97.0 100.0 108.0 119.0 131.0 144.0 151.7
lowest : 80.0 81.0 82.5 83.0 84.0, highest: 181.5 183.0 185.0 189.0 194.0
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
ThesisDrug
n missing unique Sum Mean
2407 0 2 1012 0.4204
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
logLabResWBCimp : LabResWBCimp
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 131 1.829 1.308 1.411 1.609 1.824 2.054 2.241 2.370
lowest : 0.6418 0.7883 0.8328 0.8754 0.9163, highest: 2.7536 2.8034 2.8272 2.9178 3.0006
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
logLabResCRP : LabResCRP
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 324 0.7235 -1.6095 -1.2040 -0.2232 0.7883 1.7227 2.3923 2.9230
lowest : -5.298 -3.912 -3.507 -3.219 -2.996, highest: 4.754 4.949 5.024 5.030 5.094
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
logNLratio
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
2407 0 2318 0.6075 -0.217529 -0.004496 0.285783 0.607684 0.938598 1.205873 1.409619
lowest : -1.807 -1.746 -1.447 -1.312 -1.280, highest: 2.368 2.372 2.424 3.275 4.500
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
f2
25 Variables 193 Observations
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
IL1b : Aushon hIL1b (pg/ml)
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 53 2.236 0.20 0.30 0.40 0.90 1.90 4.06 8.88
lowest : 0.03 0.10 0.20 0.30 0.40, highest: 11.60 13.30 18.30 40.40 46.50
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
IL6 : Aushon hIL6 (pg/ml)
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 125 13.1 2.70 3.70 5.40 8.40 14.10 23.54 40.84
lowest : 1.4 1.7 1.9 2.0 2.1, highest: 62.6 72.0 79.6 86.1 92.3
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
IL10 : Aushon hIL10 (pg/ml)
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 49 2.473 0.90 0.92 1.20 1.60 2.40 3.88 5.82
lowest : 0.2 0.5 0.7 0.8 0.9, highest: 10.9 13.4 16.0 24.6 39.3
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
TNFa : Aushon hTNFa (pg/ml)
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 135 22.28 0.90 1.52 2.90 5.80 16.30 42.34 62.74
lowest : 0.4 0.5 0.7 0.8 0.9, highest: 200.3 288.6 301.8 304.6 666.0
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
CESimp
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 45 15.61 2.0 3.0 7.0 13.0 23.0 32.0 37.4
lowest : 0 1 2 3 4, highest: 40 44 46 47 49
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
pa
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 11 9.782 5 6 8 11 12 12 12
2 3 4 5 6 7 8 9 10 11 12
Frequency 1 1 7 8 9 10 19 16 22 25 75
% 1 1 4 4 5 5 10 8 11 13 39
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
sc
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 20 7.244 1.0 1.2 4.0 7.0 10.0 13.0 15.0
lowest : 0 1 2 3 4, highest: 15 16 17 18 19
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
ip
n missing unique Mean
193 0 7 0.9482
0 1 2 3 4 5 6
Frequency 101 39 30 13 7 1 2
% 52 20 16 7 4 1 1
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
na
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 20 5.197 0.0 0.0 1.0 3.0 8.0 13.8 16.0
lowest : 0 1 2 3 4, highest: 15 16 17 18 19
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Age0
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 35 49.77 33 36 43 50 57 61 63
lowest : 30 31 32 33 34, highest: 60 61 62 63 64
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Sex
n missing unique Sum Mean
193 0 2 67 0.3472
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Race
n missing unique Sum Mean
193 0 2 138 0.715
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
PovStat
n missing unique Sum Mean
193 0 2 99 0.513
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
CVDclusterdich
n missing unique Sum Mean
193 0 2 31 0.1606
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
DMMclusterdich
n missing unique Sum Mean
193 0 2 69 0.3575
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
PhysBMI
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 183 31.2 20.28 21.49 24.54 29.38 36.33 45.13 47.67
lowest : 16.61 16.98 17.21 17.36 17.97, highest: 54.33 54.54 56.08 56.57 56.75
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
MedHxAlcStatusNFC
n missing unique Mean
193 0 3 1.244
0 (48, 25%), 1 (50, 26%), 2 (95, 49%)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
MedHxCigaretteStatusNFC
n missing unique Mean
193 0 3 1.285
0 (39, 20%), 1 (60, 31%), 2 (94, 49%)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
MedHxHardDrugNFC
n missing unique Mean
193 0 3 0.3057
0 (144, 75%), 1 (39, 20%), 2 (10, 5%)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
PhysMeanSBPsitimp
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 72 121.3 95.0 99.2 108.0 120.0 133.0 147.6 152.0
lowest : 90 91 92 94 95, highest: 159 160 161 165 169
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
ThesisDrug
n missing unique Sum Mean
193 0 2 101 0.5233
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
logIL1b : Aushon hIL1b (pg/ml)
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 53 -0.03387 -1.6094 -1.2040 -0.9163 -0.1054 0.6419 1.4010 2.1687
lowest : -3.5066 -2.3026 -1.6094 -1.2040 -0.9163, highest: 2.4510 2.5878 2.9069 3.6988 3.8395
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
logIL6 : Aushon hIL6 (pg/ml)
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 125 2.204 0.9933 1.3083 1.6864 2.1282 2.6462 3.1587 3.7096
lowest : 0.3365 0.5306 0.6419 0.6931 0.7419, highest: 4.1368 4.2767 4.3770 4.4555 4.5250
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
logIL10 : Aushon hIL10 (pg/ml)
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 49 0.5968 -0.10536 -0.08429 0.18232 0.47000 0.87547 1.35578 1.76098
lowest : -1.6094 -0.6931 -0.3567 -0.2231 -0.1054, highest: 2.3888 2.5953 2.7726 3.2027 3.6712
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
logTNFa : Aushon hTNFa (pg/ml)
n missing unique Mean .05 .10 .25 .50 .75 .90 .95
193 0 135 1.958 -0.1054 0.4184 1.0647 1.7579 2.7912 3.7457 4.1390
lowest : -0.9163 -0.6931 -0.3567 -0.2231 -0.1054, highest: 5.2998 5.6650 5.7098 5.7190 6.5013
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Examine the first set (I ran out of time before I got to the 2nd set).
CES-D (imputed):
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1475814 *** 0.2255966 *** 0.110158 ***
CESimp 0.0004431 0.0004098 -0.001230
Age0 -0.0029259 *** 0.0083282 ** -0.002925 *
Sex -0.0158984 -0.2249405 *** 0.007394
Race -0.1516819 *** -0.0921001 -0.322132 ***
PovStat 0.0168071 0.0290064 0.041113
CVDclusterdich 0.1218163 *** 0.3352016 ** 0.146470 ***
DMMclusterdich 0.0230160 -0.1518231 ** 0.027500
PhysBMI 0.0067381 *** 0.0745787 *** 0.000864
MedHxAlcStatusNFC -0.0044868 -0.0464175 0.010081
MedHxCigaretteStatusNFC 0.1065518 *** 0.1907321 *** 0.007440
MedHxHardDrugNFC -0.0150288 -0.1163799 * -0.068964 ***
PhysMeanSBPsitimp 0.0008249 * 0.0020143 0.001368 *
ThesisDrug 0.0180042 0.0930947 0.006384
PA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1473723 *** 0.225623 *** 0.1103481 ***
pa -0.0003717 -0.003218 0.0058235
Age0 -0.0029785 *** 0.008369 ** -0.0029283 *
Sex -0.0169225 -0.224938 *** 0.0086537
Race -0.1521761 *** -0.090553 -0.3241019 ***
PovStat 0.0179686 0.029042 0.0396208
CVDclusterdich 0.1225226 *** 0.335251 ** 0.1455162 ***
DMMclusterdich 0.0233764 -0.152086 ** 0.0274937
PhysBMI 0.0067330 *** 0.074556 *** 0.0009085
MedHxAlcStatusNFC -0.0047248 -0.046493 0.0105015
MedHxCigaretteStatusNFC 0.1071781 *** 0.190347 *** 0.0073092
MedHxHardDrugNFC -0.0144382 -0.116519 * -0.0694594 ***
PhysMeanSBPsitimp 0.0008317 * 0.002009 0.0013690 *
ThesisDrug 0.0188296 0.093025 0.0054815
SC:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1483507 *** 0.226555 *** 0.1103390 ***
sc 0.0024369 0.010098 -0.0035720
Age0 -0.0028800 *** 0.008725 ** -0.0029083 *
Sex -0.0150674 -0.217806 *** 0.0076775
Race -0.1518981 *** -0.090571 -0.3208290 ***
PovStat 0.0150636 0.017611 0.0419723 *
CVDclusterdich 0.1205187 *** 0.327300 ** 0.1473445 ***
DMMclusterdich 0.0225116 -0.155321 ** 0.0276663
PhysBMI 0.0067442 *** 0.074613 *** 0.0008587
MedHxAlcStatusNFC -0.0043398 -0.044982 0.0102019
MedHxCigaretteStatusNFC 0.1055941 *** 0.184347 *** 0.0078607
MedHxHardDrugNFC -0.0153188 -0.119768 * -0.0694279 ***
PhysMeanSBPsitimp 0.0008267 * 0.001995 0.0013542 *
ThesisDrug 0.0169610 0.085769 0.0066908
IP:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1473807 *** 0.225972 *** 0.1094742 ***
ip -0.0010024 -0.019895 -0.0006335
Age0 -0.0030061 *** 0.007950 ** -0.0027571 *
Sex -0.0170012 -0.225128 *** 0.0106043
Race -0.1525014 *** -0.094110 -0.3200835 ***
PovStat 0.0183515 0.035140 0.0376745
CVDclusterdich 0.1225004 *** 0.333937 ** 0.1442293 ***
DMMclusterdich 0.0234899 -0.150697 ** 0.0263087
PhysBMI 0.0067318 *** 0.074505 *** 0.0008694
MedHxAlcStatusNFC -0.0048873 -0.049509 0.0107027
MedHxCigaretteStatusNFC 0.1073823 *** 0.193006 *** 0.0054065
MedHxHardDrugNFC -0.0141598 -0.111986 * -0.0707288 ***
PhysMeanSBPsitimp 0.0008374 * 0.002103 0.0013471 *
ThesisDrug 0.0189949 0.095099 0.0038302
NA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1473935 *** 0.225900 *** 0.1097634 ***
na 0.0003730 -0.005053 -0.0018215
Age0 -0.0029705 *** 0.008004 ** -0.0028424 *
Sex -0.0164894 -0.233531 *** 0.0078623
Race -0.1521680 *** -0.096418 -0.3213468 ***
PovStat 0.0177046 0.035600 0.0394621
CVDclusterdich 0.1223002 *** 0.339999 ** 0.1457604 ***
DMMclusterdich 0.0233203 -0.149614 * 0.0269363
PhysBMI 0.0067374 *** 0.074549 *** 0.0008619
MedHxAlcStatusNFC -0.0046614 -0.047766 0.0103930
MedHxCigaretteStatusNFC 0.1070868 *** 0.194352 *** 0.0064106
MedHxHardDrugNFC -0.0146134 -0.112177 * -0.0695602 ***
PhysMeanSBPsitimp 0.0008304 * 0.002062 0.0013589 *
ThesisDrug 0.0185839 0.098747 0.0055202
In the first set we can eliminate MedHxAlcStatusNFC and ThesisDrug
CES-D (imputed):
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1469105 *** 0.2240967 *** 0.1099234 ***
CESimp 0.0005150 0.0008441 -0.0012358
Age0 -0.0026633 *** 0.0098531 *** -0.0029157 *
Sex -0.0171179 -0.2349653 *** 0.0088314
Race -0.1508327 *** -0.0850872 -0.3231496 ***
PovStat 0.0170301 0.0329367 0.0397959
CVDclusterdich 0.1270074 *** 0.3647321 *** 0.1469583 ***
DMMclusterdich 0.0277119 -0.1249375 * 0.0278555
PhysBMI 0.0069817 *** 0.0759797 *** 0.0008791
MedHxCigaretteStatusNFC 0.1054233 *** 0.1806469 *** 0.0091768
MedHxHardDrugNFC -0.0167074 -0.1292355 ** -0.0674593 ***
PhysMeanSBPsitimp 0.0008258 * 0.0019797 0.0013875 *
PA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1466323 *** 0.224114 *** 0.110103 ***
pa -0.0005431 -0.004163 0.005824
Age0 -0.0027076 *** 0.009870 *** -0.002932 *
Sex -0.0183511 -0.235843 *** 0.010212
Race -0.1512875 *** -0.083582 -0.325168 ***
PovStat 0.0183674 0.033933 0.038239
CVDclusterdich 0.1280860 *** 0.365539 *** 0.145698 ***
DMMclusterdich 0.0283599 * -0.124799 * 0.027576
PhysBMI 0.0069879 *** 0.075957 *** 0.000910
MedHxCigaretteStatusNFC 0.1060544 *** 0.180625 *** 0.009122
MedHxHardDrugNFC -0.0161361 -0.129009 ** -0.067849 ***
PhysMeanSBPsitimp 0.0008332 * 0.001978 0.001389 *
SC:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1477498 *** 0.225224 *** 0.1100968 ***
sc 0.0026041 0.011063 -0.0035685
Age0 -0.0026339 *** 0.010132 *** -0.0028953 *
Sex -0.0162769 -0.227679 *** 0.0091436
Race -0.1511628 *** -0.084328 -0.3218447 ***
PovStat 0.0152789 0.021490 0.0406148
CVDclusterdich 0.1254031 *** 0.354686 ** 0.1478822 ***
DMMclusterdich 0.0269538 -0.130249 * 0.0280694
PhysBMI 0.0069747 *** 0.075919 *** 0.0008766
MedHxCigaretteStatusNFC 0.1045062 *** 0.174646 *** 0.0096032
MedHxHardDrugNFC -0.0168851 -0.131789 ** -0.0679252 ***
PhysMeanSBPsitimp 0.0008281 * 0.001964 0.0013743 *
IP:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1466215 *** 0.224342 *** 0.1092321 ***
ip -0.0006770 -0.017214 -0.0010591
Age0 -0.0027320 *** 0.009510 ** -0.0027892 *
Sex -0.0185288 -0.236993 *** 0.0122646
Race -0.1516946 *** -0.087299 -0.3212144 ***
PovStat 0.0187459 0.040032 0.0363923
CVDclusterdich 0.1282001 *** 0.365606 *** 0.1439161 ***
DMMclusterdich 0.0285583 * -0.122477 * 0.0260031
PhysBMI 0.0069921 *** 0.075964 *** 0.0008488
MedHxCigaretteStatusNFC 0.1062657 *** 0.182882 *** 0.0072975
MedHxHardDrugNFC -0.0159089 -0.125326 ** -0.0689422 ***
PhysMeanSBPsitimp 0.0008381 * 0.002062 0.0013681 *
NA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1466771 *** 0.224254 *** 0.1095229 ***
na 0.0005407 -0.004039 -0.0018345
Age0 -0.0027012 *** 0.009598 *** -0.0028459 *
Sex -0.0176961 -0.243607 *** 0.0093804
Race -0.1512935 *** -0.089190 -0.3224087 ***
PovStat 0.0179772 0.039851 0.0381076
CVDclusterdich 0.1276582 *** 0.371149 *** 0.1459797 ***
DMMclusterdich 0.0281843 * -0.121167 * 0.0270485
PhysBMI 0.0069894 *** 0.076029 *** 0.0008646
MedHxCigaretteStatusNFC 0.1059441 *** 0.184103 *** 0.0082104
MedHxHardDrugNFC -0.0163572 -0.125573 ** -0.0679631 ***
PhysMeanSBPsitimp 0.0008312 * 0.002028 0.0013788 *
It's a little disappointing that there's only one significant association between a depression measure and an inflammatory marker (SC with WBC).
CES-D distributions are notoriously skewed. Let's check the distributions for the 1st set.
Maybe a log transform will make the distributions a little more symmetric?
Except for PA, the distributions look a little better. Let's re-run the analyses.
Re-examine the first set.
CES-D (imputed):
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1475980 *** 0.225856 *** 0.1104751 ***
CESimp 0.0059970 0.026940 -0.0193991
Age0 -0.0029225 *** 0.008573 ** -0.0029656 *
Sex -0.0162498 -0.222428 *** 0.0080037
Race -0.1518141 *** -0.090007 -0.3220509 ***
PovStat 0.0166433 0.023648 0.0422386 *
CVDclusterdich 0.1217929 *** 0.332298 ** 0.1469057 ***
DMMclusterdich 0.0229661 -0.153608 ** 0.0278623
PhysBMI 0.0067310 *** 0.074556 *** 0.0008859
MedHxAlcStatusNFC -0.0046095 -0.046053 0.0103603
MedHxCigaretteStatusNFC 0.1064003 *** 0.187373 *** 0.0082754
MedHxHardDrugNFC -0.0150074 -0.118707 * -0.0687206 ***
PhysMeanSBPsitimp 0.0008251 * 0.001985 0.0013712 *
ThesisDrug 0.0181270 0.090317 0.0064153
PA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1474281 *** 0.225635 *** 0.1107216 ***
pa -0.0077737 0.028334 0.0534976
Age0 -0.0029607 *** 0.008162 ** -0.0029488 *
Sex -0.0166713 -0.227364 *** 0.0080035
Race -0.1517465 *** -0.095307 -0.3247816 ***
PovStat 0.0177667 0.031430 0.0398304
CVDclusterdich 0.1224711 *** 0.336399 ** 0.1451843 ***
DMMclusterdich 0.0232257 -0.150588 * 0.0278537
PhysBMI 0.0067318 *** 0.074589 *** 0.0008959
MedHxAlcStatusNFC -0.0046582 -0.046965 0.0102071
MedHxCigaretteStatusNFC 0.1069657 *** 0.192655 *** 0.0076755
MedHxHardDrugNFC -0.0145740 -0.114933 * -0.0693037 ***
PhysMeanSBPsitimp 0.0008285 * 0.002039 0.0013773 *
ThesisDrug 0.0186406 0.095039 0.0058361
SC:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1486964 *** 0.227233 *** 0.1101342 ***
sc 0.0176433 0.082018 * -0.0194535
Age0 -0.0028591 *** 0.008878 ** -0.0028913 *
Sex -0.0155687 -0.219137 *** 0.0089485
Race -0.1521455 *** -0.091450 -0.3203612 ***
PovStat 0.0146229 0.014028 0.0413544
CVDclusterdich 0.1203195 *** 0.325322 ** 0.1468079 ***
DMMclusterdich 0.0222498 -0.157014 ** 0.0276129
PhysBMI 0.0067123 *** 0.074469 *** 0.0008971
MedHxAlcStatusNFC -0.0046054 -0.046013 0.0106413
MedHxCigaretteStatusNFC 0.1050558 *** 0.180982 *** 0.0078337
MedHxHardDrugNFC -0.0156344 -0.121725 ** -0.0694320 ***
PhysMeanSBPsitimp 0.0008234 * 0.001976 0.0013554 *
ThesisDrug 0.0168677 0.084336 0.0060759
IP:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1474464 *** 0.226127 *** 0.109551 ***
ip -0.0050946 -0.054453 -0.007808
Age0 -0.0030301 *** 0.007841 ** -0.002808 *
Sex -0.0169630 -0.225122 *** 0.010703
Race -0.1526178 *** -0.094748 -0.320321 ***
PovStat 0.0186338 0.035880 0.038331
CVDclusterdich 0.1223276 *** 0.333008 ** 0.143874 **
DMMclusterdich 0.0234839 -0.151094 ** 0.026332
PhysBMI 0.0067271 *** 0.074488 *** 0.000859
MedHxAlcStatusNFC -0.0050317 -0.049736 0.010352
MedHxCigaretteStatusNFC 0.1074867 *** 0.193393 *** 0.005638
MedHxHardDrugNFC -0.0139553 -0.111538 * -0.070245 ***
PhysMeanSBPsitimp 0.0008422 * 0.002117 0.001358 *
ThesisDrug 0.0190712 0.095388 0.003999
NA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1474389 *** 0.225805 *** 0.1100251 ***
na 0.0030268 -0.021586 -0.0128289
Age0 -0.0029595 *** 0.008051 ** -0.0028764 *
Sex -0.0162774 -0.231477 *** 0.0073224
Race -0.1520225 *** -0.095740 -0.3217910 ***
PovStat 0.0174501 0.034860 0.0402838
CVDclusterdich 0.1222554 *** 0.338389 ** 0.1457561 ***
DMMclusterdich 0.0231880 -0.149525 * 0.0274112
PhysBMI 0.0067393 *** 0.074548 *** 0.0008551
MedHxAlcStatusNFC -0.0046701 -0.047178 0.0104824
MedHxCigaretteStatusNFC 0.1069653 *** 0.193833 *** 0.0067863
MedHxHardDrugNFC -0.0148058 -0.112498 * -0.0689148 ***
PhysMeanSBPsitimp 0.0008282 * 0.002058 0.0013662 *
ThesisDrug 0.0183955 0.097822 0.0060906
Eliminate MedHxAlcStatusNFC and ThesisDrug
CES-D (imputed):
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1469112 *** 0.224413 *** 0.1102288 ***
CESimp 0.0067378 0.031052 -0.0193559
Age0 -0.0026587 *** 0.010050 *** -0.0029569 *
Sex -0.0175883 -0.232908 *** 0.0095095
Race -0.1509920 *** -0.083361 -0.3230839 ***
PovStat 0.0169137 0.027840 0.0408569
CVDclusterdich 0.1270760 *** 0.361372 *** 0.1473468 ***
DMMclusterdich 0.0277314 -0.127227 * 0.0281785
PhysBMI 0.0069765 *** 0.075921 *** 0.0008995
MedHxCigaretteStatusNFC 0.1052536 *** 0.177470 *** 0.0100457
MedHxHardDrugNFC -0.0166908 -0.131206 ** -0.0671810 ***
PhysMeanSBPsitimp 0.0008261 * 0.001952 0.0013913 *
PA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1467042 *** 0.224075 *** 0.1104865 ***
pa -0.0091942 0.019853 0.0536751
Age0 -0.0026927 *** 0.009692 *** -0.0029471 *
Sex -0.0180687 -0.238317 *** 0.0094898
Race -0.1508622 *** -0.088170 -0.3258256 ***
PovStat 0.0181581 0.036334 0.0384911
CVDclusterdich 0.1279861 *** 0.367322 *** 0.1455018 ***
DMMclusterdich 0.0281485 * -0.122836 * 0.0280619
PhysBMI 0.0069844 *** 0.076018 *** 0.0009037
MedHxCigaretteStatusNFC 0.1058534 *** 0.182782 *** 0.0094387
MedHxHardDrugNFC -0.0162506 -0.127604 ** -0.0677568 ***
PhysMeanSBPsitimp 0.0008299 * 0.002008 0.0013970 *
SC:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1480870 *** 0.225895 *** 0.1098788 ***
sc 0.0184973 * 0.086533 * -0.0192773
Age0 -0.0026130 *** 0.010268 *** -0.0028880 *
Sex -0.0168883 -0.229567 *** 0.0105236
Race -0.1514026 *** -0.085236 -0.3214312 ***
PovStat 0.0149153 0.018376 0.0399151
CVDclusterdich 0.1252558 *** 0.352775 ** 0.1471029 ***
DMMclusterdich 0.0267178 -0.132023 * 0.0278036
PhysBMI 0.0069421 *** 0.075759 *** 0.0009044
MedHxCigaretteStatusNFC 0.1039258 *** 0.171190 *** 0.0096425
MedHxHardDrugNFC -0.0172444 -0.133842 ** -0.0678353 ***
PhysMeanSBPsitimp 0.0008243 * 0.001943 0.0013758 *
IP:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1466745 *** 0.224484 *** 0.1093228 ***
ip -0.0043555 -0.048574 -0.0086664
Age0 -0.0027534 *** 0.009410 ** -0.0028359 *
Sex -0.0185160 -0.237021 *** 0.0122951
Race -0.1517888 *** -0.087870 -0.3214161 ***
PovStat 0.0190504 0.040858 0.0370748
CVDclusterdich 0.1280728 *** 0.364835 *** 0.1436435 ***
DMMclusterdich 0.0285951 * -0.122711 * 0.0260954
PhysBMI 0.0069893 *** 0.075952 *** 0.0008427
MedHxCigaretteStatusNFC 0.1063417 *** 0.183217 *** 0.0074639
MedHxHardDrugNFC -0.0157297 -0.124888 ** -0.0685390 ***
PhysMeanSBPsitimp 0.0008428 * 0.002077 0.0013785 *
NA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1467332 *** 0.224190 *** 0.109776 ***
na 0.0037990 -0.017124 -0.012765
Age0 -0.0026931 *** 0.009624 *** -0.002872 *
Sex -0.0175365 -0.241824 *** 0.008872
Race -0.1511623 *** -0.088693 -0.322838 ***
PovStat 0.0177400 0.039177 0.038886
CVDclusterdich 0.1276190 *** 0.369586 *** 0.146092 ***
DMMclusterdich 0.0280073 * -0.121344 * 0.027632
PhysBMI 0.0069895 *** 0.076017 *** 0.000864
MedHxCigaretteStatusNFC 0.1058309 *** 0.183767 *** 0.008582
MedHxHardDrugNFC -0.0165448 -0.125732 ** -0.067344 ***
PhysMeanSBPsitimp 0.0008289 * 0.002026 0.001386 *
Try some interactions
CES-D (imputed):
logLabResWBCimp logLabResCRP logNLratio
Rsq 1.502e-01 *** 0.227740 *** 0.1194693 ***
CESimp 9.024e-02 0.404660 0.0716465
Age0 8.556e-04 0.036642 * -0.0008463
Sex -2.227e-02 0.939778 -0.1501388
Race 3.466e-01 1.901486 0.2757633
PovStat 1.619e-02 0.023769 0.0399259
CVDclusterdich 1.251e-01 *** 0.370886 *** 0.1451501 ***
DMMclusterdich 2.729e-02 -0.122880 * 0.0298757
PhysBMI 6.821e-03 *** 0.075817 *** 0.0004693
MedHxCigaretteStatusNFC 1.053e-01 *** 0.180492 *** 0.0116700
MedHxHardDrugNFC -1.574e-02 -0.126777 ** -0.0650982 ***
PhysMeanSBPsitimp 8.873e-04 * 0.001992 0.0014779 *
CESimp:Age0 -1.457e-03 -0.007919 -0.0012646
CESimp:Sex -3.485e-03 -0.430277 -0.0658616
Age0:Sex -2.417e-05 -0.027769 0.0060528
CESimp:Race -1.852e-01 -0.586552 -0.1864508
Age0:Race -8.695e-03 -0.037435 -0.0077382
Sex:Race -1.332e-01 -1.433896 -0.0202791
CESimp:Age0:Sex 2.403e-04 0.010693 0.0007136
CESimp:Age0:Race 3.290e-03 0.010924 0.0023460
CESimp:Sex:Race 3.163e-02 0.582705 0.0731036
Age0:Sex:Race 2.683e-03 0.026803 -0.0070606
CESimp:Age0:Sex:Race -8.383e-04 -0.011545 0.0007139
PA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 1.484e-01 *** 0.227409 *** 0.1192491 ***
pa -9.203e-02 -0.528010 0.0708242
Age0 -7.804e-03 -0.014578 -0.0014679
Sex -4.501e-01 -1.652103 0.5785158
Race -7.299e-02 -3.098452 -0.5294154
PovStat 1.830e-02 0.030699 0.0375842
CVDclusterdich 1.258e-01 *** 0.368132 *** 0.1430427 ***
DMMclusterdich 2.823e-02 * -0.121126 * 0.0304003
PhysBMI 6.933e-03 *** 0.075654 *** 0.0004492
MedHxCigaretteStatusNFC 1.067e-01 *** 0.186353 *** 0.0105113
MedHxHardDrugNFC -1.593e-02 -0.123394 ** -0.0643398 ***
PhysMeanSBPsitimp 8.658e-04 * 0.002184 0.0014377 *
pa:Age0 2.163e-03 0.013597 -0.0011702
pa:Sex 1.756e-01 0.648701 -0.3827145
Age0:Sex 1.113e-02 0.046887 -0.0164846
pa:Race -2.116e-02 1.514368 0.1337237
Age0:Race -6.027e-05 0.067011 0.0010694
Sex:Race -7.380e-01 4.057002 -1.9829422
pa:Age0:Sex -4.488e-03 -0.020824 0.0104291
pa:Age0:Race -1.741e-04 -0.033160 -0.0010251
pa:Sex:Race 2.982e-01 -1.716159 0.9271693
Age0:Sex:Race 1.312e-02 -0.098042 0.0518079
pa:Age0:Sex:Race -5.399e-03 0.041053 -0.0246565
SC:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1532995 *** 0.229977 *** 0.1190072 ***
sc 0.1633999 0.963610 * 0.1429646
Age0 0.0017571 0.048815 ** 0.0006962
Sex 0.0356197 1.561464 -0.0697280
Race 0.4038721 2.175095 * 0.2785156
PovStat 0.0135346 0.012315 0.0377590
CVDclusterdich 0.1232695 *** 0.364493 *** 0.1433043 ***
DMMclusterdich 0.0263134 -0.129429 * 0.0302823
PhysBMI 0.0067176 *** 0.075257 *** 0.0004179
MedHxCigaretteStatusNFC 0.1036906 *** 0.171697 *** 0.0105223
MedHxHardDrugNFC -0.0163889 -0.127455 ** -0.0658733 ***
PhysMeanSBPsitimp 0.0008854 * 0.001947 0.0014746 *
sc:Age0 -0.0024957 -0.017208 * -0.0025534
sc:Sex -0.0383001 -0.907570 -0.1345874
Age0:Sex -0.0011856 -0.037006 0.0043004
sc:Race -0.2805168 * -0.946340 -0.2551668
Age0:Race -0.0093574 -0.041275 * -0.0080311
Sex:Race -0.3156403 -1.226822 -0.2078460
sc:Age0:Sex 0.0010296 0.019290 0.0019595
sc:Age0:Race 0.0047997 0.016989 0.0033958
sc:Sex:Race 0.1436777 0.679468 0.2076132
Age0:Sex:Race 0.0064006 0.021703 -0.0029657
sc:Age0:Sex:Race -0.0032192 -0.013014 -0.0014604
IP:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1480506 *** 0.2287911 *** 0.1173474 ***
ip -0.0061793 -0.7359699 0.1048243
Age0 -0.0028776 0.0105314 -0.0032142
Sex -0.0366131 -0.2688381 -0.3219735
Race -0.1273014 -0.0052491 -0.1194859
PovStat 0.0194708 0.0395930 0.0367887
CVDclusterdich 0.1252780 *** 0.3614245 *** 0.1385806 **
DMMclusterdich 0.0289485 * -0.1144574 * 0.0287061
PhysBMI 0.0068828 *** 0.0760302 *** 0.0003392
MedHxCigaretteStatusNFC 0.1068252 *** 0.1836820 *** 0.0088879
MedHxHardDrugNFC -0.0151712 -0.1205471 ** -0.0644595 ***
PhysMeanSBPsitimp 0.0008563 * 0.0019222 0.0014760 *
ip:Age0 -0.0001741 0.0113487 -0.0017484
ip:Sex -0.0328664 0.2783817 -0.0353774
Age0:Sex 0.0002155 -0.0008585 0.0080201
ip:Race 0.0029293 0.8354787 -0.1733088
Age0:Race -0.0004172 -0.0022813 -0.0024514
Sex:Race 0.0123706 -0.0816276 0.2462886
ip:Age0:Sex 0.0014565 -0.0011489 0.0006688
ip:Age0:Race 0.0001241 -0.0152072 0.0018605
ip:Sex:Race -0.0742666 0.1817882 -0.0812657
Age0:Sex:Race -0.0001812 0.0010942 -0.0078940
ip:Age0:Sex:Race 0.0006151 -0.0057069 0.0033719
NA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1495919 *** 0.228242 *** 0.1174961 ***
na 0.0657020 0.350506 0.0834083
Age0 -0.0013389 0.027134 ** -0.0016806
Sex 0.0129231 0.590029 -0.1484815
Race 0.0850836 1.533850 * 0.0373854
PovStat 0.0182487 0.036268 0.0389041
CVDclusterdich 0.1260640 *** 0.373492 *** 0.1426299 **
DMMclusterdich 0.0275992 -0.116763 * 0.0294736
PhysBMI 0.0068604 *** 0.075995 *** 0.0005111
MedHxCigaretteStatusNFC 0.1066076 *** 0.187688 *** 0.0112270
MedHxHardDrugNFC -0.0154930 -0.121397 ** -0.0649188 ***
PhysMeanSBPsitimp 0.0008892 * 0.002091 0.0014833 *
na:Age0 -0.0010309 -0.007423 -0.0016323
na:Sex -0.0300970 -0.538845 -0.1308436
Age0:Sex -0.0003382 -0.018181 0.0047982
na:Race -0.1453087 -0.795353 * -0.1741760
Age0:Race -0.0040260 -0.031279 * -0.0052022
Sex:Race -0.0900291 -1.361907 -0.0745551
na:Age0:Sex 0.0006315 0.012644 0.0024350
na:Age0:Race 0.0025518 0.015320 * 0.0026158
na:Sex:Race 0.0083477 1.053303 0.1981335
Age0:Sex:Race 0.0012415 0.026968 -0.0022719
na:Age0:Sex:Race -0.0001113 -0.022262 -0.0027636
Looks like we can eliminate sex (and some other interactions):
CES-D (imputed):
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1489515 *** 0.219610 *** 0.1124197 ***
CESimp 0.0908333 0.223525 0.0698341
Age0 0.0009340 0.023841 0.0032390
Race 0.2842699 1.166749 0.3063738
PovStat 0.0176008 0.042286 0.0384663
CVDclusterdich 0.1259838 *** 0.354446 ** 0.1504577 ***
DMMclusterdich 0.0270617 -0.125396 * 0.0286054
PhysBMI 0.0071136 *** 0.078757 *** 0.0006835
MedHxCigaretteStatusNFC 0.1042124 *** 0.166848 *** 0.0102249
MedHxHardDrugNFC -0.0188056 -0.161430 *** -0.0641900 ***
PhysMeanSBPsitimp 0.0008229 * 0.001390 0.0014911 *
CESimp:Age0 -0.0013950 -0.003152 -0.0015898
CESimp:Race -0.1703860 * -0.307754 -0.1739990
Age0:Race -0.0074727 -0.022839 -0.0117363
CESimp:Age0:Race 0.0029229 0.005177 0.0030845
PA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1462405 *** 0.218926 *** 0.1124953 ***
pa -0.0563926 -0.436088 -0.0634603
Age0 -0.0052474 -0.004224 -0.0066805
Race -0.3090583 -1.781773 -1.1593879
PovStat 0.0199518 0.052376 0.0358166
CVDclusterdich 0.1269577 *** 0.359699 *** 0.1467992 ***
DMMclusterdich 0.0283141 * -0.119860 * 0.0289682
PhysBMI 0.0072131 *** 0.078925 *** 0.0006907
MedHxCigaretteStatusNFC 0.1051566 *** 0.173257 *** 0.0100202
MedHxHardDrugNFC -0.0189738 -0.159452 *** -0.0654768 ***
PhysMeanSBPsitimp 0.0007813 * 0.001432 0.0014897 *
pa:Age0 0.0011649 0.008753 0.0025795
pa:Race 0.0713481 0.945167 0.4391051
Age0:Race 0.0043982 0.034443 0.0182541
pa:Age0:Race -0.0019586 -0.019239 -0.0095129
SC:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1518351 *** 0.222216 *** 0.1120160 ***
sc 0.1482977 * 0.545732 0.1170880
Age0 0.0012449 0.030537 ** 0.0039276
Race 0.2792327 1.438747 0.2399125
PovStat 0.0148876 0.031187 0.0371349
CVDclusterdich 0.1242597 *** 0.347349 ** 0.1506481 ***
DMMclusterdich 0.0261340 -0.129738 * 0.0286106
PhysBMI 0.0069890 *** 0.078318 *** 0.0006466
MedHxCigaretteStatusNFC 0.1026832 *** 0.160261 *** 0.0097619
MedHxHardDrugNFC -0.0192670 -0.163844 *** -0.0649566 ***
PhysMeanSBPsitimp 0.0008269 * 0.001430 0.0014644 *
sc:Age0 -0.0020800 -0.008003 -0.0025101
sc:Race -0.2264545 * -0.570014 -0.1969145
Age0:Race -0.0068432 -0.026871 -0.0104252
sc:Age0:Race 0.0035860 0.009302 0.0034289
IP:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1461974 *** 0.2204517 *** 0.1114442 ***
ip -0.0208843 -0.6279948 0.0804300
Age0 -0.0028168 * 0.0098321 0.0001179
Race -0.1253030 -0.0961420 -0.0220770
PovStat 0.0205679 0.0623641 0.0345808
CVDclusterdich 0.1270833 *** 0.3541665 ** 0.1460921 ***
DMMclusterdich 0.0288150 * -0.1175777 * 0.0269729
PhysBMI 0.0072081 *** 0.0791099 *** 0.0005990
MedHxCigaretteStatusNFC 0.1055719 *** 0.1720273 *** 0.0078350
MedHxHardDrugNFC -0.0182296 -0.1573504 *** -0.0650695 ***
PhysMeanSBPsitimp 0.0007943 * 0.0013506 0.0014539 *
ip:Age0 0.0004761 0.0111683 -0.0012798
ip:Race -0.0271119 0.9601798 * -0.1930020
Age0:Race -0.0004232 -0.0006817 -0.0056785 *
ip:Age0:Race 0.0003386 -0.0186052 * 0.0030155
NA:
logLabResWBCimp logLabResCRP logNLratio
Rsq 0.1482148 *** 0.219377 *** 0.1114181 ***
na 0.0556425 0.144045 0.0577675
Age0 -0.0014610 0.018854 ** 0.0010692
Race 0.0348268 0.854071 0.0145404
PovStat 0.0194486 0.055288 0.0374789
CVDclusterdich 0.1264021 *** 0.360463 *** 0.1495607 ***
DMMclusterdich 0.0273337 -0.120977 * 0.0280699
PhysBMI 0.0071639 *** 0.078989 *** 0.0007041
MedHxCigaretteStatusNFC 0.1052971 *** 0.173689 *** 0.0090948
MedHxHardDrugNFC -0.0187606 -0.157800 *** -0.0647943 ***
PhysMeanSBPsitimp 0.0008302 * 0.001490 0.0014825 *
na:Age0 -0.0008136 -0.002270 -0.0013078
na:Race -0.1346968 -0.340956 -0.1109883
Age0:Race -0.0032927 -0.017696 -0.0066036
na:Age0:Race 0.0023894 0.005761 0.0020039
It appears as though there's a weak relationship between depressive symptoms and inflammatory markers, but the restricted range and funky distributions seem to interfere with straightforward analyses. It seems like there's something here, but it may take some incremental model building to see what's going on. It might be valuable to start with NA (for example), age, sex, race, and poverty status with all their interactions to see whether some combination of interactions reveals a pattern of relationships with inflammatory markers. I would add the other covariates afterwards to see which moderators kill the association.
These scatterplots illustrate part of the problem.
Here's untransformed CRP with NA:
log(CRP) with NA:
log(CRP) with log(NA):
