I revised the analyses slightly. This probably is not a final analysis, but perhaps it is close.
I revised cigarette use to ever v. never. I combined alcohol and illicit drug use into ever (either one) v never (both).
If we can find a justification for age stratified analyses, then these results look pretty interesting. Note that I have yet to detect an age interaction so there is no was to justify stratification based on these results. Nevertheless, perhaps we can find a separate justification.
In the following, I analyzed Disc02b, TotDisc5, numsources. For each analysis, the complete group analysis is first, followed by age (median) stratified young, then old.
Call:
lm(formula = Formula, data = discrimIMT, na.action = na.omit)
Residuals:
Min 1Q Median 3Q Max
-0.36813321 -0.08489957 -0.01213427 0.07452405 0.68051165
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.3613062563 0.0451353283 8.00495 3.3521e-15
Disc02b -0.0079563177 0.0070156252 -1.13409 0.257035
CESimp -0.0014982633 0.0008128932 -1.84312 0.065611
Age0 0.0037514366 0.0004836137 7.75709 2.1664e-14
ses01 -0.0088666937 0.0083983990 -1.05576 0.291337
Sex 0.0199143008 0.0085093553 2.34028 0.019468
WRATtotal 0.0005961526 0.0005287673 1.12744 0.259832
everCigEver 0.0201038072 0.0110852986 1.81356 0.070051
everImbEver -0.0075841340 0.0156798270 -0.48369 0.628715
HTN 0.0157819383 0.0090252815 1.74864 0.080666
DIABETES 0.0211364043 0.0119728234 1.76537 0.077813
Cholimp 0.0001414332 0.0001008911 1.40184 0.161278
MedsChol -0.0122627525 0.0131951412 -0.92934 0.352942
RevCVDclusterdich 0.0158345379 0.0134250962 1.17947 0.238495
PhysBMI 0.0025797343 0.0005916417 4.36030 1.4355e-05
Disc02b:CESimp 0.0006996169 0.0003575792 1.95654 0.050684
Residual standard error: 0.1221447 on 984 degrees of freedom
(8 observations deleted due to missingness)
Multiple R-squared: 0.1419073, Adjusted R-squared: 0.1288266
F-statistic: 10.84862 on 15 and 984 DF, p-value: < 2.2204e-16
Call:
lm(formula = Formula, data = yIMT, na.action = na.omit)
Residuals:
Min 1Q Median 3Q Max
-0.23312431 -0.07349269 -0.01712227 0.06942937 0.43628462
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.4110307311 0.0613945164 6.69491 0.000000000055964
Disc02b -0.0132670122 0.0091715416 -1.44654 0.148627
CESimp -0.0016460792 0.0009426662 -1.74620 0.081366
Age0 0.0038646175 0.0009251636 4.17723 0.000034596225721
ses01 -0.0211216724 0.0102665373 -2.05733 0.040151
Sex 0.0083280115 0.0103028920 0.80832 0.419277
WRATtotal -0.0001046910 0.0006658808 -0.15722 0.875131
everCigEver 0.0175236063 0.0122580246 1.42956 0.153442
everImbEver -0.0101538201 0.0190060002 -0.53424 0.593401
HTN 0.0055601898 0.0115559410 0.48115 0.630609
DIABETES 0.0166179773 0.0178626693 0.93032 0.352637
Cholimp 0.0001201900 0.0001281804 0.93766 0.348852
MedsChol 0.0079916958 0.0262078731 0.30493 0.760537
RevCVDclusterdich 0.0441902126 0.0195936145 2.25534 0.024526
PhysBMI 0.0031580717 0.0006968141 4.53216 0.000007248055257
Disc02b:CESimp 0.0008888475 0.0004282654 2.07546 0.038434
Residual standard error: 0.1078847 on 521 degrees of freedom
(5 observations deleted due to missingness)
Multiple R-squared: 0.120014, Adjusted R-squared: 0.09467851
F-statistic: 4.736992 on 15 and 521 DF, p-value: 0.00000001276742
Call:
lm(formula = Formula, data = oIMT, na.action = na.omit)
Residuals:
Min 1Q Median 3Q Max
-0.35818271 -0.09874219 -0.00367312 0.07519110 0.67097446
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.3093120672 0.1022434760 3.02525 0.0026275
Disc02b -0.0057117132 0.0109385445 -0.52216 0.6018149
CESimp -0.0018065973 0.0014463130 -1.24911 0.2122803
Age0 0.0036805227 0.0013909079 2.64613 0.0084291
ses01 0.0052074480 0.0138712054 0.37541 0.7075305
Sex 0.0381909992 0.0143894485 2.65410 0.0082356
WRATtotal 0.0013630677 0.0008430980 1.61674 0.1066405
everCigEver 0.0260635127 0.0212713856 1.22529 0.2211130
everImbEver -0.0143884399 0.0272059136 -0.52887 0.5971569
HTN 0.0275973449 0.0142994618 1.92996 0.0542446
DIABETES 0.0286976188 0.0171762853 1.67077 0.0954671
Cholimp 0.0002153153 0.0001605725 1.34092 0.1806264
MedsChol -0.0153420721 0.0168542763 -0.91028 0.3631666
RevCVDclusterdich -0.0059921159 0.0190867173 -0.31394 0.7537118
PhysBMI 0.0016661802 0.0010201766 1.63323 0.1031253
Disc02b:CESimp 0.0006333105 0.0006180186 1.02474 0.3060385
Residual standard error: 0.1368211 on 447 degrees of freedom
(3 observations deleted due to missingness)
Multiple R-squared: 0.06838304, Adjusted R-squared: 0.03712073
F-statistic: 2.187395 on 15 and 447 DF, p-value: 0.006202495
Call:
lm(formula = Formula, data = discrimIMT, na.action = na.omit)
Residuals:
Min 1Q Median 3Q Max
-0.37898522 -0.08356048 -0.01131626 0.07206452 0.65464272
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.32971758193 0.04535362370 7.26993 7.3451e-13
TotDisc5 0.00075729577 0.00074427899 1.01749 0.309172
CESimp -0.00032608462 0.00065594425 -0.49712 0.619214
Age0 0.00381647545 0.00047919894 7.96428 4.5868e-15
ses01 -0.00774449607 0.00841692103 -0.92011 0.357742
Sex 0.01989214448 0.00829433202 2.39828 0.016658
WRATtotal 0.00061538694 0.00052465411 1.17294 0.241106
everCigEver 0.01841759468 0.01106691426 1.66420 0.096391
everImbEver -0.00554448254 0.01571603837 -0.35279 0.724321
HTN 0.01780640787 0.00904561082 1.96851 0.049290
DIABETES 0.01986049071 0.01197445379 1.65857 0.097522
Cholimp 0.00013729705 0.00010107051 1.35843 0.174640
MedsChol -0.01230935574 0.01324562552 -0.92931 0.352955
RevCVDclusterdich 0.01551250731 0.01341205131 1.15661 0.247714
PhysBMI 0.00265543450 0.00059417965 4.46908 8.7719e-06
TotDisc5:CESimp 0.00001118564 0.00003854120 0.29023 0.771705
Residual standard error: 0.1220363 on 980 degrees of freedom
(12 observations deleted due to missingness)
Multiple R-squared: 0.1419613, Adjusted R-squared: 0.128828
F-statistic: 10.80931 on 15 and 980 DF, p-value: < 2.2204e-16
Call:
lm(formula = Formula, data = yIMT, na.action = na.omit)
Residuals:
Min 1Q Median 3Q Max
-0.23344491 -0.07608254 -0.01668302 0.06548131 0.41632587
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.39172415009 0.06151525453 6.36792 0.00000000042232
TotDisc5 -0.00071295937 0.00094978727 -0.75065 0.453203
CESimp -0.00078091499 0.00082216061 -0.94983 0.342639
Age0 0.00382303811 0.00092361573 4.13921 0.00004066407614
ses01 -0.02003833063 0.01027923013 -1.94940 0.051786
Sex 0.00806802309 0.01011942885 0.79728 0.425653
WRATtotal -0.00011376503 0.00066677951 -0.17062 0.864590
everCigEver 0.01530368744 0.01225618352 1.24865 0.212356
everImbEver -0.00701898991 0.01917461861 -0.36606 0.714472
HTN 0.00897685034 0.01162482005 0.77221 0.440339
DIABETES 0.01538647718 0.01790539509 0.85932 0.390560
Cholimp 0.00012934848 0.00012866623 1.00530 0.315220
MedsChol 0.00855329246 0.02630581409 0.32515 0.745200
RevCVDclusterdich 0.04052068772 0.01963201751 2.06401 0.039513
PhysBMI 0.00315266896 0.00070164145 4.49328 0.00000865367789
TotDisc5:CESimp 0.00005973564 0.00004652457 1.28396 0.199729
Residual standard error: 0.108072 on 519 degrees of freedom
(7 observations deleted due to missingness)
Multiple R-squared: 0.11293, Adjusted R-squared: 0.08729212
F-statistic: 4.404813 on 15 and 519 DF, p-value: 0.00000007754818
Call:
lm(formula = Formula, data = oIMT, na.action = na.omit)
Residuals:
Min 1Q Median 3Q Max
-0.37265862 -0.09797528 -0.00069189 0.07862809 0.62800996
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.26621043189 0.10256154168 2.59562 0.0097543
TotDisc5 0.00206556412 0.00118489879 1.74324 0.0819822
CESimp -0.00012527533 0.00106021953 -0.11816 0.9059944
Age0 0.00383751491 0.00138605396 2.76866 0.0058638
ses01 0.00503072559 0.01386945272 0.36272 0.7169864
Sex 0.03755638718 0.01390122809 2.70166 0.0071630
WRATtotal 0.00131928621 0.00082762215 1.59407 0.1116307
everCigEver 0.02641646107 0.02121867654 1.24496 0.2138007
everImbEver -0.01334626365 0.02710426317 -0.49240 0.6226762
HTN 0.02891709196 0.01429378036 2.02305 0.0436650
DIABETES 0.02480350843 0.01717347270 1.44429 0.1493606
Cholimp 0.00020005504 0.00015997015 1.25058 0.2117459
MedsChol -0.01360272057 0.01690667842 -0.80458 0.4214936
RevCVDclusterdich -0.00437311264 0.01899965927 -0.23017 0.8180670
PhysBMI 0.00191063599 0.00102176342 1.86994 0.0621486
TotDisc5:CESimp -0.00004084492 0.00006637913 -0.61533 0.5386529
Residual standard error: 0.1363049 on 445 degrees of freedom
(5 observations deleted due to missingness)
Multiple R-squared: 0.0748858, Adjusted R-squared: 0.04370218
F-statistic: 2.401446 on 15 and 445 DF, p-value: 0.002349365
Call:
lm(formula = Formula, data = discrimIMT, na.action = na.omit)
Residuals:
Min 1Q Median 3Q Max
-0.37900456 -0.08461656 -0.01110389 0.07265546 0.67112272
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.36471994071 0.04706539703 7.74922 2.3009e-14
numsources -0.00108652108 0.00112131026 -0.96897 0.332796
CESimp -0.00235590933 0.00105006056 -2.24359 0.025081
Age0 0.00371030067 0.00048147029 7.70619 3.1640e-14
ses01 -0.00721345048 0.00840070361 -0.85867 0.390731
Sex 0.01928605593 0.00832109451 2.31773 0.020669
WRATtotal 0.00056467947 0.00052608374 1.07336 0.283371
everCigEver 0.02013392549 0.01104305334 1.82322 0.068574
everImbEver -0.00712519841 0.01564616214 -0.45540 0.648925
HTN 0.01623263091 0.00901555755 1.80051 0.072086
DIABETES 0.01949462164 0.01196137035 1.62980 0.103465
Cholimp 0.00015496693 0.00010103969 1.53372 0.125420
MedsChol -0.01154987179 0.01318324492 -0.87610 0.381189
RevCVDclusterdich 0.01266676170 0.01346215499 0.94092 0.346979
PhysBMI 0.00263201447 0.00059108831 4.45283 9.4466e-06
numsources:CESimp 0.00012666014 0.00005548987 2.28258 0.022668
Residual standard error: 0.1219813 on 982 degrees of freedom
(10 observations deleted due to missingness)
Multiple R-squared: 0.1449411, Adjusted R-squared: 0.1318801
F-statistic: 11.09726 on 15 and 982 DF, p-value: < 2.2204e-16
Call:
lm(formula = Formula, data = yIMT, na.action = na.omit)
Residuals:
Min 1Q Median 3Q Max
-0.23204372 -0.07510456 -0.01618633 0.06783925 0.43846928
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.43038222447 0.06298126079 6.83350 0.000000000023218
numsources -0.00238373127 0.00137121189 -1.73841 0.0827290
CESimp -0.00327464141 0.00118163223 -2.77129 0.0057831
Age0 0.00378548293 0.00091598117 4.13271 0.000041769740857
ses01 -0.01994769954 0.01017914935 -1.95966 0.0505680
Sex 0.00682877476 0.01005483544 0.67915 0.4973423
WRATtotal -0.00016274787 0.00066041682 -0.24643 0.8054449
everCigEver 0.01827433209 0.01213535653 1.50588 0.1327052
everImbEver -0.00960191500 0.01887175300 -0.50880 0.6111090
HTN 0.00710235852 0.01147579594 0.61890 0.5362534
DIABETES 0.01423234937 0.01773342340 0.80257 0.4225881
Cholimp 0.00012435442 0.00012745429 0.97568 0.3296766
MedsChol 0.01245090981 0.02609754580 0.47709 0.6334972
RevCVDclusterdich 0.04135635624 0.01944436355 2.12691 0.0338970
PhysBMI 0.00321472389 0.00069139108 4.64965 0.000004219527490
numsources:CESimp 0.00019071886 0.00006302277 3.02619 0.0025990
Residual standard error: 0.107246 on 521 degrees of freedom
(5 observations deleted due to missingness)
Multiple R-squared: 0.1304027, Adjusted R-squared: 0.1053664
F-statistic: 5.208528 on 15 and 521 DF, p-value: 0.0000000009746393
Call:
lm(formula = Formula, data = oIMT, na.action = na.omit)
Residuals:
Min 1Q Median 3Q Max
-0.36921284 -0.09749887 -0.00143991 0.07754219 0.66012505
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.30215258698 0.10480706894 2.88294 0.0041307
numsources -0.00017713739 0.00186680618 -0.09489 0.9244466
CESimp -0.00188162420 0.00195276677 -0.96357 0.3357857
Age0 0.00360088208 0.00139665575 2.57822 0.0102516
ses01 0.00646012277 0.01394647937 0.46321 0.6434416
Sex 0.03796054844 0.01406539334 2.69886 0.0072225
WRATtotal 0.00133244246 0.00084066619 1.58498 0.1136803
everCigEver 0.02619455026 0.02129874648 1.22986 0.2193980
everImbEver -0.01417144800 0.02720459649 -0.52092 0.6026809
HTN 0.02749561789 0.01436510641 1.91406 0.0562541
DIABETES 0.02720184435 0.01723069443 1.57869 0.1151188
Cholimp 0.00023831899 0.00016114162 1.47894 0.1398635
MedsChol -0.01515439281 0.01687478131 -0.89805 0.3696447
RevCVDclusterdich -0.00849654992 0.01927418982 -0.44083 0.6595536
PhysBMI 0.00172997101 0.00102554573 1.68688 0.0923275
numsources:CESimp 0.00007574543 0.00010183056 0.74384 0.4573671
Residual standard error: 0.1369742 on 445 degrees of freedom
(5 observations deleted due to missingness)
Multiple R-squared: 0.06804541, Adjusted R-squared: 0.03663121
F-statistic: 2.166072 on 15 and 445 DF, p-value: 0.006824846
/prj/hnd/git/dataMgr/Projects/Beatty/2016-06-06-discrimIMT/03-R/2017-03-01-discIMT.Rmd