PA ~ RE Pre-COVID
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ roaming_entropy + time_of_day + precipi + mean_temp +
## distance + dow + (1 | subject)
## data: prelockdown_df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 4944 -9805.14 19652.28 981(5880) 2.92e+02 1.2e+08
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1 1
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.06305758 0.02880767 2.189 0.028603 *
## time_of_day 0.00364605 0.00142483 2.559 0.010499 *
## precipi 0.09090892 0.16537613 0.550 0.582518
## mean_temp -0.01247655 0.00277387 -4.498 6.86e-06 ***
## distance 0.00013884 0.00006843 2.029 0.042472 *
## dowMonday -0.32445115 0.09834804 -3.299 0.000970 ***
## dowSaturday 0.44504380 0.10092955 4.409 1.04e-05 ***
## dowSunday -0.08811485 0.09845691 -0.895 0.370809
## dowThursday -0.36338187 0.09377795 -3.875 0.000107 ***
## dowTuesday -0.71914671 0.08994121 -7.996 1.29e-15 ***
## dowWednesday -0.39540160 0.09488587 -4.167 3.08e-05 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.5000 0.2514 -17.900
## (10,20]|(20,30] -3.6872 0.2436 -15.138
## (20,30]|(30,40] -2.8736 0.2398 -11.985
## (30,40]|(40,50] -2.0979 0.2382 -8.809
## (40,50]|(50,60] -1.3257 0.2373 -5.585
## (50,60]|(60,70] -0.3181 0.2368 -1.343
## (60,70]|(70,80] 0.4550 0.2368 1.921
## (70,80]|(80,90] 1.3511 0.2381 5.675
## (80,90]|(90,100] 2.5289 0.2447 10.334
## (6738 observations deleted due to missingness)
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ roaming_entropy + dow + (1 | subject)
## data: prelockdown_df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 4944 -9724.91 19483.83 3042(18428) 3.69e-03 4.9e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.266 1.125
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.09371 0.02919 3.210 0.00133 **
## dowMonday -0.28565 0.09722 -2.938 0.00330 **
## dowSaturday 0.48452 0.10405 4.657 0.00000322 ***
## dowSunday -0.01228 0.09942 -0.124 0.90170
## dowThursday -0.45616 0.09561 -4.771 0.00000183 ***
## dowTuesday -0.79533 0.09171 -8.673 < 2e-16 ***
## dowWednesday -0.44376 0.09667 -4.590 0.00000442 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.2721 0.1559 -27.403
## (10,20]|(20,30] -3.3948 0.1408 -24.119
## (20,30]|(30,40] -2.4886 0.1328 -18.736
## (30,40]|(40,50] -1.5976 0.1292 -12.366
## (40,50]|(50,60] -0.6771 0.1276 -5.307
## (50,60]|(60,70] 0.5594 0.1275 4.389
## (60,70]|(70,80] 1.5218 0.1293 11.767
## (70,80]|(80,90] 2.6060 0.1345 19.382
## (80,90]|(90,100] 3.9622 0.1511 26.224
## (6738 observations deleted due to missingness)
## df BIC
## complex_PA_RE 21 19788.91
## simple_PA_RE 17 19594.43
PA ~ NL Pre_COVID
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ novel_locations + time_of_day + precipi + mean_temp +
## distance + dow + (1 | subject)
## data: prelockdown_df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 4722 -9335.41 18712.82 958(5742) 9.76e+02 1.3e+08
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1 1
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.00243244 0.00045949 5.294 1.2e-07 ***
## time_of_day 0.00454691 0.00145994 3.114 0.001843 **
## precipi 0.05836599 0.16729864 0.349 0.727185
## mean_temp -0.01063658 0.00285510 -3.725 0.000195 ***
## distance 0.00004159 0.00006925 0.601 0.548166
## dowMonday -0.26681706 0.10312354 -2.587 0.009672 **
## dowSaturday 0.39145406 0.10288377 3.805 0.000142 ***
## dowSunday -0.11200762 0.10002073 -1.120 0.262780
## dowThursday -0.36596759 0.09727633 -3.762 0.000168 ***
## dowTuesday -0.69028981 0.09384723 -7.355 1.9e-13 ***
## dowWednesday -0.37591435 0.09964589 -3.773 0.000162 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.4104 0.2552 -17.284
## (10,20]|(20,30] -3.5993 0.2468 -14.581
## (20,30]|(30,40] -2.7668 0.2428 -11.397
## (30,40]|(40,50] -1.9798 0.2411 -8.212
## (40,50]|(50,60] -1.2016 0.2403 -5.000
## (50,60]|(60,70] -0.1847 0.2398 -0.770
## (60,70]|(70,80] 0.6020 0.2400 2.509
## (70,80]|(80,90] 1.4898 0.2413 6.174
## (80,90]|(90,100] 2.6670 0.2481 10.748
## (6960 observations deleted due to missingness)
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ novel_locations + dow + (1 | subject)
## data: prelockdown_df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 4722 -9253.81 18541.62 1905(11725) 4.83e-01 7.0e+05
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.268 1.126
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.0035111 0.0004594 7.642 2.14e-14 ***
## dowMonday -0.2198405 0.1014147 -2.168 0.0302 *
## dowSaturday 0.4162784 0.1058343 3.933 8.38e-05 ***
## dowSunday -0.0598573 0.1006830 -0.595 0.5522
## dowThursday -0.4459820 0.0988669 -4.511 6.45e-06 ***
## dowTuesday -0.7509536 0.0952865 -7.881 3.25e-15 ***
## dowWednesday -0.3964624 0.1010946 -3.922 8.79e-05 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.3414 0.1450 -29.935
## (10,20]|(20,30] -3.4664 0.1272 -27.262
## (20,30]|(30,40] -2.5404 0.1175 -21.627
## (30,40]|(40,50] -1.6365 0.1129 -14.493
## (40,50]|(50,60] -0.7104 0.1108 -6.409
## (50,60]|(60,70] 0.5370 0.1107 4.850
## (60,70]|(70,80] 1.5186 0.1131 13.431
## (70,80]|(80,90] 2.5952 0.1192 21.780
## (80,90]|(90,100] 3.9564 0.1384 28.594
## (6960 observations deleted due to missingness)
## df BIC
## complex_PA_NL 21 18848.48
## simple_PA_NL 17 18651.44
RE ~ Lockdown
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: roaming_entropy ~ lockdown + dow + (1 | subject)
## Data: df
##
## REML criterion at convergence: 56046.6
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.7418 -0.6409 0.0226 0.6227 8.2953
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 0.1129 0.3360
## Residual 0.7187 0.8477
## Number of obs: 22077, groups: subject, 236
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.47168 0.02708 483.70444 91.258 < 2e-16 ***
## lockdown -1.26665 0.01167 21958.74460 -108.519 < 2e-16 ***
## dowMonday -0.16938 0.02139 21839.60711 -7.919 2.51e-15 ***
## dowSaturday -0.13628 0.02165 21840.00786 -6.296 3.12e-10 ***
## dowSunday -0.36768 0.02153 21840.60334 -17.074 < 2e-16 ***
## dowThursday -0.04689 0.02108 21839.48620 -2.224 0.0262 *
## dowTuesday -0.13551 0.02116 21839.99873 -6.404 1.54e-10 ***
## dowWednesday -0.10346 0.02118 21839.61492 -4.885 1.04e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) lckdwn dwMndy dwStrd dwSndy dwThrs dwTsdy
## lockdown -0.195
## dowMonday -0.389 -0.007
## dowSaturday -0.381 -0.022 0.489
## dowSunday -0.390 0.010 0.491 0.485
## dowThursday -0.394 -0.011 0.502 0.496 0.498
## dowTuesday -0.395 0.001 0.500 0.494 0.497 0.507
## dowWednesdy -0.395 0.003 0.499 0.493 0.496 0.507 0.505
NL ~ Lockdown
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: novel_locations ~ lockdown + dow + (1 | subject)
## Data: df
##
## REML criterion at convergence: 222894
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.0679 -0.5095 -0.2037 0.2143 17.5727
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 293.7 17.14
## Residual 3074.6 55.45
## Number of obs: 20462, groups: subject, 236
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 65.848 1.563 713.324 42.126 < 2e-16 ***
## lockdown -30.876 0.792 20400.351 -38.986 < 2e-16 ***
## dowMonday -14.443 1.454 20228.137 -9.935 < 2e-16 ***
## dowSaturday 4.634 1.472 20228.912 3.148 0.00165 **
## dowSunday -6.405 1.464 20230.085 -4.375 1.22e-05 ***
## dowThursday -5.808 1.431 20227.895 -4.059 4.95e-05 ***
## dowTuesday -12.516 1.437 20228.888 -8.710 < 2e-16 ***
## dowWednesday -11.184 1.438 20228.141 -7.778 7.75e-15 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) lckdwn dwMndy dwStrd dwSndy dwThrs dwTsdy
## lockdown -0.247
## dowMonday -0.457 -0.008
## dowSaturday -0.447 -0.025 0.488
## dowSunday -0.458 0.009 0.490 0.484
## dowThursday -0.464 -0.011 0.502 0.496 0.498
## dowTuesday -0.465 0.001 0.500 0.493 0.496 0.508
## dowWednesdy -0.465 0.004 0.499 0.493 0.496 0.507 0.505
PA ~ Lockdown
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ lockdown + dow + (1 | subject)
## data: df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 8257 -16125.26 32284.53 3845(23427) 4.02e-03 4.7e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.378 1.174
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## lockdown -0.66039 0.04122 -16.021 < 2e-16 ***
## dowMonday -0.26437 0.07353 -3.595 0.000324 ***
## dowSaturday 0.35310 0.07733 4.566 4.97e-06 ***
## dowSunday 0.06728 0.07647 0.880 0.378984
## dowThursday -0.29577 0.07341 -4.029 5.60e-05 ***
## dowTuesday -0.50770 0.07285 -6.970 3.18e-12 ***
## dowWednesday -0.25541 0.07520 -3.397 0.000682 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.34675 0.11457 -37.941
## (10,20]|(20,30] -3.46858 0.10523 -32.963
## (20,30]|(30,40] -2.62277 0.10059 -26.075
## (30,40]|(40,50] -1.74588 0.09808 -17.801
## (40,50]|(50,60] -0.79273 0.09672 -8.196
## (50,60]|(60,70] 0.45717 0.09655 4.735
## (60,70]|(70,80] 1.43757 0.09809 14.655
## (70,80]|(80,90] 2.52669 0.10272 24.597
## (80,90]|(90,100] 3.87137 0.11825 32.738
## (20153 observations deleted due to missingness)
NA ~ Lockdown
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: NA_avg ~ lockdown + dow + (1 | subject)
## data: df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 8257 -15958.66 31951.33 3418(29264) 8.15e-03 6.7e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.384 1.544
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## lockdown 0.41073 0.04096 10.028 < 2e-16 ***
## dowMonday 0.15017 0.07370 2.038 0.04158 *
## dowSaturday -0.27849 0.07773 -3.583 0.00034 ***
## dowSunday -0.09953 0.07656 -1.300 0.19360
## dowThursday 0.21614 0.07395 2.923 0.00347 **
## dowTuesday 0.42769 0.07289 5.868 0.00000000441 ***
## dowWednesday 0.14854 0.07494 1.982 0.04746 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -2.9591 0.1231 -24.039
## (10,20]|(20,30] -1.7781 0.1188 -14.969
## (20,30]|(30,40] -0.8360 0.1173 -7.129
## (30,40]|(40,50] 0.0188 0.1168 0.161
## (40,50]|(50,60] 0.9595 0.1172 8.186
## (50,60]|(60,70] 2.1605 0.1191 18.140
## (60,70]|(70,80] 2.9993 0.1218 24.621
## (70,80]|(80,90] 4.0785 0.1288 31.666
## (80,90]|(90,100] 5.0634 0.1420 35.661
## (20153 observations deleted due to missingness)
PA ~ RE * Lockdown
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ roaming_entropy * lockdown + dow + (1 | subject)
## data: df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 8257 -16114.47 32266.93 4273(25852) 9.58e-03 6.9e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.379 1.174
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.07658 0.02821 2.715 0.006638 **
## lockdown -0.64080 0.08611 -7.441 9.97e-14 ***
## dowMonday -0.24756 0.07367 -3.361 0.000778 ***
## dowSaturday 0.34597 0.07772 4.452 8.52e-06 ***
## dowSunday 0.10057 0.07742 1.299 0.193943
## dowThursday -0.29347 0.07344 -3.996 6.45e-05 ***
## dowTuesday -0.50362 0.07293 -6.905 5.01e-12 ***
## dowWednesday -0.24348 0.07530 -3.234 0.001223 **
## roaming_entropy:lockdown 0.09906 0.05218 1.898 0.057635 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.1632 0.1337 -31.129
## (10,20]|(20,30] -3.2851 0.1258 -26.110
## (20,30]|(30,40] -2.4386 0.1220 -19.994
## (30,40]|(40,50] -1.5604 0.1200 -13.003
## (40,50]|(50,60] -0.6057 0.1190 -5.089
## (50,60]|(60,70] 0.6459 0.1191 5.425
## (60,70]|(70,80] 1.6280 0.1205 13.506
## (70,80]|(80,90] 2.7190 0.1246 21.830
## (80,90]|(90,100] 4.0652 0.1379 29.482
## (20153 observations deleted due to missingness)
PA ~ NL * Lockdown
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ novel_locations * lockdown + dow + (1 | subject)
## data: df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 8035 -15640.43 31318.86 2865(17164) 5.13e-01 7.0e+05
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.387 1.178
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.0036404 0.0004548 8.004 1.21e-15 ***
## lockdown -0.5325563 0.0521671 -10.209 < 2e-16 ***
## dowMonday -0.2078100 0.0754538 -2.754 0.005885 **
## dowSaturday 0.3237868 0.0786209 4.118 3.82e-05 ***
## dowSunday 0.0896366 0.0779748 1.150 0.250326
## dowThursday -0.2789739 0.0749592 -3.722 0.000198 ***
## dowTuesday -0.4589532 0.0746317 -6.150 7.77e-10 ***
## dowWednesday -0.2014242 0.0773529 -2.604 0.009215 **
## novel_locations:lockdown -0.0005953 0.0010494 -0.567 0.570502
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.1576 0.1194 -34.823
## (10,20]|(20,30] -3.2809 0.1101 -29.787
## (20,30]|(30,40] -2.4281 0.1056 -22.999
## (30,40]|(40,50] -1.5439 0.1032 -14.963
## (40,50]|(50,60] -0.5850 0.1020 -5.736
## (50,60]|(60,70] 0.6749 0.1021 6.609
## (60,70]|(70,80] 1.6721 0.1040 16.083
## (70,80]|(80,90] 2.7603 0.1088 25.374
## (80,90]|(90,100] 4.1121 0.1243 33.084
## (20375 observations deleted due to missingness)
PA ~ RE * Pre_COVID_RE Post-Lockdown
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ roaming_entropy * pre_covid_re + dow + (1 | subject)
## data: lockdown_df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 3313 -6340.07 12718.14 3530(31738) 1.09e-02 2.2e+04
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.375 1.541
## Number of groups: subject 232
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.02973 0.38135 0.078 0.9379
## pre_covid_re 0.16097 0.31949 0.504 0.6144
## dowMonday -0.14634 0.11580 -1.264 0.2063
## dowSaturday 0.25678 0.12024 2.135 0.0327 *
## dowSunday 0.32322 0.12532 2.579 0.0099 **
## dowThursday -0.05268 0.11736 -0.449 0.6536
## dowTuesday 0.06819 0.12336 0.553 0.5804
## dowWednesday 0.07790 0.12242 0.636 0.5245
## roaming_entropy:pre_covid_re 0.09863 0.15656 0.630 0.5287
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.2122 0.7626 -4.212
## (10,20]|(20,30] -2.2470 0.7601 -2.956
## (20,30]|(30,40] -1.3769 0.7591 -1.814
## (30,40]|(40,50] -0.4269 0.7587 -0.563
## (40,50]|(50,60] 0.6645 0.7585 0.876
## (50,60]|(60,70] 2.0562 0.7584 2.711
## (60,70]|(70,80] 3.1662 0.7588 4.173
## (70,80]|(80,90] 4.3624 0.7612 5.731
## (80,90]|(90,100] 5.7666 0.7710 7.479
## (7082 observations deleted due to missingness)
PA ~ NL * Pre_COVID_RE Post-Lockdown
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ novel_locations * pre_covid_re + dow + (1 | subject)
## data: lockdown_df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 3313 -6344.21 12726.41 2268(20065) 1.42e-01 3.0e+07
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.359 1.536
## Number of groups: subject 232
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.009453 0.009437 1.002 0.3165
## pre_covid_re 0.295803 0.285377 1.037 0.3000
## dowMonday -0.156613 0.115768 -1.353 0.1761
## dowSaturday 0.266242 0.120413 2.211 0.0270 *
## dowSunday 0.321993 0.125369 2.568 0.0102 *
## dowThursday -0.052323 0.117529 -0.445 0.6562
## dowTuesday 0.074890 0.123265 0.608 0.5435
## dowWednesday 0.069199 0.122544 0.565 0.5723
## novel_locations:pre_covid_re -0.001960 0.003785 -0.518 0.6045
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.0390 0.6870 -4.424
## (10,20]|(20,30] -2.0746 0.6843 -3.032
## (20,30]|(30,40] -1.2058 0.6833 -1.765
## (30,40]|(40,50] -0.2571 0.6830 -0.376
## (40,50]|(50,60] 0.8323 0.6832 1.218
## (50,60]|(60,70] 2.2224 0.6840 3.249
## (60,70]|(70,80] 3.3325 0.6854 4.862
## (70,80]|(80,90] 4.5263 0.6889 6.570
## (80,90]|(90,100] 5.9264 0.6994 8.474
## (7082 observations deleted due to missingness)
Notch Analysis
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ bs(novel_locations, knots = c(1, 26), degree = 1) +
## dow + (1 | subject)
## data: prelockdown_df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 4722 -9252.14 18542.27 3216(19529) 5.40e-03 1.0e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.268 1.126
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -0.24789 0.16712 -1.483
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.01486 0.14847 0.100
## bs(novel_locations, knots = c(1, 26), degree = 1)3 3.12261 0.49083 6.362
## dowMonday -0.22045 0.10144 -2.173
## dowSaturday 0.42421 0.10596 4.003
## dowSunday -0.04585 0.10139 -0.452
## dowThursday -0.44556 0.09887 -4.507
## dowTuesday -0.74466 0.09540 -7.805
## dowWednesday -0.39703 0.10112 -3.926
## Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.1380
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.9203
## bs(novel_locations, knots = c(1, 26), degree = 1)3 1.99e-10 ***
## dowMonday 0.0298 *
## dowSaturday 6.24e-05 ***
## dowSunday 0.6511
## dowThursday 6.59e-06 ***
## dowTuesday 5.93e-15 ***
## dowWednesday 8.63e-05 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.4476 0.2002 -22.219
## (10,20]|(20,30] -3.5722 0.1874 -19.059
## (20,30]|(30,40] -2.6463 0.1808 -14.634
## (30,40]|(40,50] -1.7425 0.1779 -9.797
## (40,50]|(50,60] -0.8162 0.1765 -4.624
## (50,60]|(60,70] 0.4314 0.1762 2.448
## (60,70]|(70,80] 1.4136 0.1775 7.965
## (70,80]|(80,90] 2.4908 0.1812 13.748
## (80,90]|(90,100] 3.8521 0.1943 19.823
## (6960 observations deleted due to missingness)
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ bs(novel_locations, knots = c(1, 26), degree = 1) +
## dow + (1 | subject)
## data: lockdown_df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 3313 -6337.51 12713.02 3182(29041) 1.12e-02 3.2e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.386 1.545
## Number of groups: subject 232
##
## Coefficients:
## Estimate Std. Error z value
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.21012 0.09741 2.157
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.45981 0.09727 4.727
## bs(novel_locations, knots = c(1, 26), degree = 1)3 2.53107 1.33026 1.903
## dowMonday -0.13214 0.11593 -1.140
## dowSaturday 0.26788 0.12041 2.225
## dowSunday 0.33390 0.12538 2.663
## dowThursday -0.03889 0.11742 -0.331
## dowTuesday 0.07744 0.12330 0.628
## dowWednesday 0.07814 0.12264 0.637
## Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.03100 *
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.00000228 ***
## bs(novel_locations, knots = c(1, 26), degree = 1)3 0.05708 .
## dowMonday 0.25437
## dowSaturday 0.02610 *
## dowSunday 0.00775 **
## dowThursday 0.74051
## dowTuesday 0.52993
## dowWednesday 0.52403
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.5677 0.1691 -21.093
## (10,20]|(20,30] -2.6032 0.1561 -16.677
## (20,30]|(30,40] -1.7325 0.1504 -11.521
## (30,40]|(40,50] -0.7809 0.1477 -5.289
## (40,50]|(50,60] 0.3123 0.1472 2.122
## (50,60]|(60,70] 1.7071 0.1502 11.369
## (60,70]|(70,80] 2.8196 0.1563 18.036
## (70,80]|(80,90] 4.0156 0.1701 23.602
## (80,90]|(90,100] 5.4195 0.2087 25.963
## (7082 observations deleted due to missingness)
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ bs(novel_locations, knots = c(1, 26), degree = 1) *
## lockdown + dow + (1 | subject)
## data: df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 8035 -15633.47 31312.93 5178(31070) 1.18e-02 9.5e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.389 1.178
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -0.28045 0.16189
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.06716 0.14291
## bs(novel_locations, knots = c(1, 26), degree = 1)3 3.41389 0.50312
## lockdown -0.82908 0.14446
## dowMonday -0.19971 0.07553
## dowSaturday 0.32878 0.07868
## dowSunday 0.09739 0.07834
## dowThursday -0.27311 0.07498
## dowTuesday -0.45339 0.07471
## dowWednesday -0.19954 0.07740
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.43665 0.18323
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.39924 0.16278
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -2.35273 1.37021
## z value Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -1.732 0.08320
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.470 0.63837
## bs(novel_locations, knots = c(1, 26), degree = 1)3 6.785 1.16e-11
## lockdown -5.739 9.51e-09
## dowMonday -2.644 0.00819
## dowSaturday 4.179 2.93e-05
## dowSunday 1.243 0.21379
## dowThursday -3.643 0.00027
## dowTuesday -6.069 1.29e-09
## dowWednesday -2.578 0.00994
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 2.383 0.01717
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 2.453 0.01418
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -1.717 0.08597
##
## bs(novel_locations, knots = c(1, 26), degree = 1)1 .
## bs(novel_locations, knots = c(1, 26), degree = 1)2
## bs(novel_locations, knots = c(1, 26), degree = 1)3 ***
## lockdown ***
## dowMonday **
## dowSaturday ***
## dowSunday
## dowThursday ***
## dowTuesday ***
## dowWednesday **
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown *
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown *
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.3346 0.1784 -24.294
## (10,20]|(20,30] -3.4570 0.1721 -20.085
## (20,30]|(30,40] -2.6033 0.1691 -15.397
## (30,40]|(40,50] -1.7181 0.1675 -10.255
## (40,50]|(50,60] -0.7580 0.1668 -4.545
## (50,60]|(60,70] 0.5031 0.1667 3.017
## (60,70]|(70,80] 1.5012 0.1677 8.950
## (70,80]|(80,90] 2.5903 0.1706 15.180
## (80,90]|(90,100] 3.9426 0.1810 21.788
## (20375 observations deleted due to missingness)
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ bs(novel_locations, knots = c(1, 26), degree = 1) *
## lockdown + dow + (1 | subject)
## data: distance_less_than_100_km
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 7608 -14726.91 29499.82 5453(32935) 3.57e-03 2.2e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.454 1.206
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -0.21976 0.16316
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.22193 0.14554
## bs(novel_locations, knots = c(1, 26), degree = 1)3 2.09923 0.28533
## lockdown -0.84067 0.14520
## dowMonday -0.20175 0.07760
## dowSaturday 0.33090 0.08189
## dowSunday 0.13422 0.08157
## dowThursday -0.23576 0.07703
## dowTuesday -0.39519 0.07683
## dowWednesday -0.22528 0.07970
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.37026 0.18465
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.60401 0.16784
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -1.86581 0.64544
## z value Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -1.347 0.17802
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -1.525 0.12728
## bs(novel_locations, knots = c(1, 26), degree = 1)3 7.357 1.88e-13
## lockdown -5.790 7.05e-09
## dowMonday -2.600 0.00932
## dowSaturday 4.041 5.32e-05
## dowSunday 1.645 0.09989
## dowThursday -3.061 0.00221
## dowTuesday -5.144 2.69e-07
## dowWednesday -2.826 0.00471
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 2.005 0.04495
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 3.599 0.00032
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -2.891 0.00384
##
## bs(novel_locations, knots = c(1, 26), degree = 1)1
## bs(novel_locations, knots = c(1, 26), degree = 1)2
## bs(novel_locations, knots = c(1, 26), degree = 1)3 ***
## lockdown ***
## dowMonday **
## dowSaturday ***
## dowSunday .
## dowThursday **
## dowTuesday ***
## dowWednesday **
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown *
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown ***
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.3964 0.1813 -24.249
## (10,20]|(20,30] -3.4833 0.1745 -19.967
## (20,30]|(30,40] -2.6108 0.1713 -15.244
## (30,40]|(40,50] -1.7157 0.1697 -10.111
## (40,50]|(50,60] -0.7482 0.1689 -4.430
## (50,60]|(60,70] 0.5374 0.1689 3.182
## (60,70]|(70,80] 1.5592 0.1700 9.173
## (70,80]|(80,90] 2.6656 0.1733 15.382
## (80,90]|(90,100] 4.0537 0.1853 21.880
## (20802 observations deleted due to missingness)
## df BIC
## PA_NL_PreLD_notch 19 18665.01
## PA_NL_PostLD_notch 19 12829.03
## PA_NL_LD_notch 23 31473.74
## PA_NL_LD_notch_controldistance_100KM 23 29659.37
Z Notch Analysis
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula:
## PA_avg ~ bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1) *
## lockdown + dow + (1 | subject)
## data: distance_less_than_100_km
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 7608 -14733.13 29512.26 4951(29978) 1.27e-02 9.5e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.459 1.208
## Number of groups: subject 236
##
## Coefficients:
## Estimate
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.30426
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.19963
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 2.35244
## lockdown -0.51970
## dowMonday -0.20781
## dowSaturday 0.32220
## dowSunday 0.12877
## dowThursday -0.24111
## dowTuesday -0.40560
## dowWednesday -0.22993
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown -0.06417
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.27430
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -2.00498
## Std. Error
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.38561
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.34071
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.42261
## lockdown 0.38979
## dowMonday 0.07755
## dowSaturday 0.08194
## dowSunday 0.08153
## dowThursday 0.07706
## dowTuesday 0.07676
## dowWednesday 0.07962
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.44588
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.39304
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.73957
## z value
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.789
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.586
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 5.567
## lockdown -1.333
## dowMonday -2.680
## dowSaturday 3.932
## dowSunday 1.580
## dowThursday -3.129
## dowTuesday -5.284
## dowWednesday -2.888
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown -0.144
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.698
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -2.711
## Pr(>|z|)
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.43008
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.55792
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.000000026
## lockdown 0.18244
## dowMonday 0.00737
## dowSaturday 0.000084161
## dowSunday 0.11422
## dowThursday 0.00176
## dowTuesday 0.000000127
## dowWednesday 0.00388
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.88557
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.48525
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.00671
##
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 ***
## lockdown
## dowMonday **
## dowSaturday ***
## dowSunday
## dowThursday **
## dowTuesday ***
## dowWednesday **
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.9623 0.3604 -10.994
## (10,20]|(20,30] -3.0502 0.3572 -8.539
## (20,30]|(30,40] -2.1780 0.3558 -6.121
## (30,40]|(40,50] -1.2836 0.3552 -3.613
## (40,50]|(50,60] -0.3171 0.3551 -0.893
## (50,60]|(60,70] 0.9677 0.3551 2.725
## (60,70]|(70,80] 1.9888 0.3557 5.592
## (70,80]|(80,90] 3.0943 0.3574 8.658
## (80,90]|(90,100] 4.4811 0.3635 12.327
## (20802 observations deleted due to missingness)
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ bs(within_subjects_lockdown_z_NL, knots = c(-0.725,
## -0.3), degree = 1) * lockdown + dow + (1 | subject)
## data: distance_less_than_100_km
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 7602 -14716.69 29479.38 5175(31543) 2.28e-02 9.3e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.459 1.208
## Number of groups: subject 236
##
## Coefficients:
## Estimate
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 -0.08973
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 -0.10431
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 3.32821
## lockdown -1.48194
## dowMonday -0.21346
## dowSaturday 0.31912
## dowSunday 0.12493
## dowThursday -0.24157
## dowTuesday -0.41367
## dowWednesday -0.23486
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.90163
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 1.03381
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -1.54221
## Std. Error
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.33945
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.30588
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.51410
## lockdown 0.55156
## dowMonday 0.07760
## dowSaturday 0.08196
## dowSunday 0.08161
## dowThursday 0.07703
## dowTuesday 0.07675
## dowWednesday 0.07965
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.60146
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.54885
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.81483
## z value
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 -0.264
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 -0.341
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 6.474
## lockdown -2.687
## dowMonday -2.751
## dowSaturday 3.894
## dowSunday 1.531
## dowThursday -3.136
## dowTuesday -5.390
## dowWednesday -2.949
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 1.499
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 1.884
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -1.893
## Pr(>|z|)
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.79152
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.73309
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 9.55e-11
## lockdown 0.00721
## dowMonday 0.00594
## dowSaturday 9.87e-05
## dowSunday 0.12584
## dowThursday 0.00171
## dowTuesday 7.05e-08
## dowWednesday 0.00319
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.13385
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.05962
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.05840
##
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 ***
## lockdown **
## dowMonday **
## dowSaturday ***
## dowSunday
## dowThursday **
## dowTuesday ***
## dowWednesday **
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown .
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.3402 0.3278 -13.242
## (10,20]|(20,30] -3.4253 0.3242 -10.566
## (20,30]|(30,40] -2.5520 0.3226 -7.912
## (30,40]|(40,50] -1.6565 0.3218 -5.148
## (40,50]|(50,60] -0.6897 0.3215 -2.146
## (50,60]|(60,70] 0.5958 0.3215 1.853
## (60,70]|(70,80] 1.6191 0.3220 5.028
## (70,80]|(80,90] 2.7264 0.3237 8.423
## (80,90]|(90,100] 4.1100 0.3302 12.447
## (20808 observations deleted due to missingness)
Symptom Analysis
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ depression_mean * lockdown + (1 | subject)
## data: df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 8257 -16193.03 32412.07 2020(15823) 4.85e-02 3.3e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.129 1.063
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## depression_mean -0.090700 0.015632 -5.802 0.00000000654 ***
## lockdown -0.522186 0.059158 -8.827 < 2e-16 ***
## depression_mean:lockdown -0.018987 0.008702 -2.182 0.0291 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.6121 0.1269 -36.330
## (10,20]|(20,30] -3.7365 0.1185 -31.538
## (20,30]|(30,40] -2.8937 0.1142 -25.342
## (30,40]|(40,50] -2.0216 0.1117 -18.096
## (40,50]|(50,60] -1.0776 0.1102 -9.776
## (50,60]|(60,70] 0.1533 0.1097 1.397
## (60,70]|(70,80] 1.1166 0.1108 10.081
## (70,80]|(80,90] 2.1900 0.1146 19.105
## (80,90]|(90,100] 3.5202 0.1283 27.432
## (20153 observations deleted due to missingness)
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ novel_locations * depression_mean * lockdown + (1 |
## subject)
## data: df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 8035 -15817.81 31669.61 645(3220) 1.07e+03 3.7e+07
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1 1
## Number of groups: subject 236
##
## Coefficients:
## Estimate Std. Error z value
## novel_locations 0.00250722 0.00062425 4.016
## depression_mean -0.08815745 0.01565996 -5.629
## lockdown -0.40739418 0.07242211 -5.625
## novel_locations:depression_mean 0.00021639 0.00009255 2.338
## novel_locations:lockdown 0.00329196 0.00142246 2.314
## depression_mean:lockdown 0.00376538 0.01083984 0.347
## novel_locations:depression_mean:lockdown -0.00071825 0.00021192 -3.389
## Pr(>|z|)
## novel_locations 0.0000591000 ***
## depression_mean 0.0000000181 ***
## lockdown 0.0000000185 ***
## novel_locations:depression_mean 0.019382 *
## novel_locations:lockdown 0.020652 *
## depression_mean:lockdown 0.728317
## novel_locations:depression_mean:lockdown 0.000701 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.9119 0.1240 -31.553
## (10,20]|(20,30] -3.0984 0.1162 -26.653
## (20,30]|(30,40] -2.3257 0.1125 -20.672
## (30,40]|(40,50] -1.5476 0.1105 -14.004
## (40,50]|(50,60] -0.7306 0.1095 -6.673
## (50,60]|(60,70] 0.3040 0.1093 2.780
## (60,70]|(70,80] 1.1027 0.1103 9.997
## (70,80]|(80,90] 1.9939 0.1132 17.612
## (80,90]|(90,100] 3.1308 0.1235 25.357
## (20375 observations deleted due to missingness)
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ novel_locations * anx_change + (1 | subject)
## data: lockdown_df
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 2970 -5709.03 11444.07 1183(11059) 1.72e+00 1.7e+06
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.194 1.481
## Number of groups: subject 205
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.0055372 0.0011091 4.992 0.000000596 ***
## anx_change -0.0273779 0.0274564 -0.997 0.3187
## novel_locations:anx_change -0.0005719 0.0002773 -2.063 0.0391 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.76984 0.14736 -25.583
## (10,20]|(20,30] -2.77980 0.12922 -21.512
## (20,30]|(30,40] -1.91722 0.12107 -15.835
## (30,40]|(40,50] -0.98637 0.11681 -8.444
## (40,50]|(50,60] 0.06696 0.11552 0.580
## (50,60]|(60,70] 1.45902 0.11867 12.294
## (60,70]|(70,80] 2.59237 0.12712 20.393
## (70,80]|(80,90] 3.74911 0.14582 25.710
## (80,90]|(90,100] 5.15776 0.19616 26.294
## (7425 observations deleted due to missingness)
## [1] "BIC: DEP_on_PA_by_lockdown, NL_on_PA_by_DEP_by_lockdown, NL_on_PA_by_anxiety_change"
## df BIC
## d1 13 32503.31
## d2 17 31788.47
## d3 13 11522.02
BIC Analysis
## df BIC
## complex_PA_RE 21 19788.91
## simple_PA_RE 17 19594.43
## complex_PA_NL 21 18848.48
## simple_PA_NL 17 18651.44
## RE_LD 10 56146.58
## NL_LD 10 222993.24
## PA_LD 17 32403.85
## NA_LD 17 32070.65
## PA_RE_LD 19 32400.29
## PA_NL_LD 19 31451.70
## PA_RE_PreLDRE_PostLD 19 12834.15
## PA_RE_PreLDNL_PostLD 19 12842.42
## PA_NL_PreLD_notch 19 18665.01
## PA_NL_PostLD_notch 19 12829.03
## PA_NL_LD_notch 23 31473.74
## PA_NL_LD_notch_controldistance_100KM 23 29659.37
## Z_within_sub 23 29671.81
## Z_within_sub_LD 23 29638.91
## d1 13 32503.31
## d2 17 31788.47
## d3 13 11522.02