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 4704 -9350.94 18743.88 957(5773) 3.38e+02 NaN
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1 1
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.0463139 NA NA NA
## time_of_day 0.0034675 NA NA NA
## precipi 0.1410273 NA NA NA
## mean_temp -0.0136288 NA NA NA
## distance 0.0001739 NA NA NA
## dowMonday -0.3313003 NA NA NA
## dowSaturday 0.4500323 NA NA NA
## dowSunday -0.1091738 NA NA NA
## dowThursday -0.3357451 NA NA NA
## dowTuesday -0.7124834 NA NA NA
## dowWednesday -0.4044147 NA NA NA
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.6343 NA NA
## (10,20]|(20,30] -3.7971 NA NA
## (20,30]|(30,40] -2.9789 NA NA
## (30,40]|(40,50] -2.2027 NA NA
## (40,50]|(50,60] -1.4357 NA NA
## (50,60]|(60,70] -0.4421 NA NA
## (60,70]|(70,80] 0.3265 NA NA
## (70,80]|(80,90] 1.2347 NA NA
## (80,90]|(90,100] 2.4177 NA NA
## (6400 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 4704 -9274.42 18582.84 3232(19515) 5.25e-03 5.1e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.257 1.121
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.079193 0.030762 2.574 0.01004 *
## dowMonday -0.261272 0.100388 -2.603 0.00925 **
## dowSaturday 0.520515 0.107365 4.848 1.25e-06 ***
## dowSunday 0.002201 0.102991 0.021 0.98295
## dowThursday -0.394568 0.098389 -4.010 6.06e-05 ***
## dowTuesday -0.746825 0.094102 -7.936 2.08e-15 ***
## dowWednesday -0.444648 0.099599 -4.464 8.03e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.2872 0.1613 -26.578
## (10,20]|(20,30] -3.3824 0.1454 -23.269
## (20,30]|(30,40] -2.4696 0.1373 -17.982
## (30,40]|(40,50] -1.5765 0.1337 -11.789
## (40,50]|(50,60] -0.6605 0.1321 -5.000
## (50,60]|(60,70] 0.5608 0.1320 4.249
## (60,70]|(70,80] 1.5166 0.1339 11.327
## (70,80]|(80,90] 2.6106 0.1392 18.758
## (80,90]|(90,100] 3.9695 0.1563 25.390
## (6400 observations deleted due to missingness)
## df BIC
## complex_PA_RE 21 18879.46
## simple_PA_RE 17 18692.60
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 4486 -8890.54 17823.08 957(5737) 7.75e+02 1.2e+08
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1 1
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.00246667 0.00048941 5.040 4.65e-07 ***
## time_of_day 0.00443046 0.00149934 2.955 0.003127 **
## precipi 0.12721492 0.17416526 0.730 0.465130
## mean_temp -0.01191610 0.00305673 -3.898 9.69e-05 ***
## distance 0.00006105 0.00007375 0.828 0.407742
## dowMonday -0.24975818 0.10719956 -2.330 0.019814 *
## dowSaturday 0.41888978 0.10613504 3.947 7.92e-05 ***
## dowSunday -0.10031150 0.10359235 -0.968 0.332880
## dowThursday -0.31718113 0.10051644 -3.156 0.001602 **
## dowTuesday -0.66771765 0.09682937 -6.896 5.36e-12 ***
## dowWednesday -0.36610888 0.10308538 -3.552 0.000383 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.4980 0.2690 -16.720
## (10,20]|(20,30] -3.6603 0.2605 -14.053
## (20,30]|(30,40] -2.8203 0.2565 -10.997
## (30,40]|(40,50] -2.0324 0.2548 -7.975
## (40,50]|(50,60] -1.2600 0.2541 -4.960
## (50,60]|(60,70] -0.2589 0.2535 -1.021
## (60,70]|(70,80] 0.5229 0.2536 2.062
## (70,80]|(80,90] 1.4212 0.2547 5.579
## (80,90]|(90,100] 2.6060 0.2617 9.960
## (6618 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 4486 -8813.27 17660.55 1786(10853) 3.79e-01 6.4e+05
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.252 1.119
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.0035268 0.0004926 7.159 8.12e-13 ***
## dowMonday -0.1794433 0.1049249 -1.710 0.087228 .
## dowSaturday 0.4651903 0.1090316 4.267 1.99e-05 ***
## dowSunday -0.0239454 0.1040002 -0.230 0.817903
## dowThursday -0.3726155 0.1017977 -3.660 0.000252 ***
## dowTuesday -0.6966520 0.0978759 -7.118 1.10e-12 ***
## dowWednesday -0.3877805 0.1043200 -3.717 0.000201 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.3218 0.1489 -29.031
## (10,20]|(20,30] -3.4176 0.1298 -26.338
## (20,30]|(30,40] -2.4815 0.1198 -20.714
## (30,40]|(40,50] -1.5745 0.1152 -13.669
## (40,50]|(50,60] -0.6532 0.1131 -5.777
## (50,60]|(60,70] 0.5773 0.1130 5.111
## (60,70]|(70,80] 1.5517 0.1154 13.450
## (70,80]|(80,90] 2.6357 0.1216 21.673
## (80,90]|(90,100] 4.0018 0.1416 28.252
## (6618 observations deleted due to missingness)
## df BIC
## complex_PA_NL 21 17957.67
## simple_PA_NL 17 17769.49
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: 50657.6
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.8257 -0.6328 0.0263 0.6198 8.3853
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 0.1144 0.3382
## Residual 0.6973 0.8351
## Number of obs: 20179, groups: subject, 236
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.50502 0.02751 484.85536 91.046 < 2e-16 ***
## lockdown -1.27637 0.01223 20140.08462 -104.364 < 2e-16 ***
## dowMonday -0.18130 0.02205 19944.67274 -8.222 < 2e-16 ***
## dowSaturday -0.18098 0.02236 19944.74709 -8.093 6.14e-16 ***
## dowSunday -0.41082 0.02220 19946.22733 -18.502 < 2e-16 ***
## dowThursday -0.06820 0.02172 19943.51170 -3.140 0.00169 **
## dowTuesday -0.15661 0.02176 19944.87451 -7.198 6.35e-13 ***
## dowWednesday -0.12033 0.02180 19943.46055 -5.520 3.44e-08 ***
## ---
## 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.181
## dowMonday -0.394 -0.007
## dowSaturday -0.386 -0.024 0.487
## dowSunday -0.395 0.006 0.491 0.484
## dowThursday -0.399 -0.011 0.501 0.495 0.498
## dowTuesday -0.401 0.000 0.501 0.493 0.497 0.508
## dowWednesdy -0.401 0.005 0.500 0.492 0.496 0.507 0.506
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: 201341.8
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.0560 -0.5151 -0.2065 0.2194 18.1264
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 287.2 16.95
## Residual 2878.1 53.65
## Number of obs: 18594, groups: subject, 236
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 65.4615 1.5660 702.9386 41.802 < 2e-16 ***
## lockdown -29.9202 0.8134 18586.0000 -36.783 < 2e-16 ***
## dowMonday -15.9583 1.4763 18360.9852 -10.810 < 2e-16 ***
## dowSaturday 3.0628 1.4973 18362.4381 2.046 0.0408 *
## dowSunday -8.3229 1.4870 18363.3419 -5.597 0.0000000221 ***
## dowThursday -7.5872 1.4519 18359.8449 -5.226 0.0000001753 ***
## dowTuesday -13.5934 1.4553 18361.3992 -9.341 < 2e-16 ***
## dowWednesday -12.2980 1.4581 18359.8400 -8.434 < 2e-16 ***
## ---
## 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.231
## dowMonday -0.462 -0.008
## dowSaturday -0.451 -0.027 0.486
## dowSunday -0.463 0.006 0.489 0.483
## dowThursday -0.469 -0.012 0.501 0.494 0.497
## dowTuesday -0.471 0.000 0.500 0.493 0.497 0.509
## dowWednesdy -0.471 0.005 0.499 0.492 0.495 0.507 0.506
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 7600 -14867.69 29769.38 3447(20982) 4.33e-03 4.5e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.39 1.179
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## lockdown -0.63442 0.04392 -14.446 < 2e-16 ***
## dowMonday -0.26516 0.07720 -3.435 0.000593 ***
## dowSaturday 0.34884 0.08113 4.300 1.71e-05 ***
## dowSunday 0.06021 0.08006 0.752 0.452048
## dowThursday -0.27182 0.07667 -3.545 0.000392 ***
## dowTuesday -0.51421 0.07593 -6.772 1.27e-11 ***
## dowWednesday -0.30702 0.07877 -3.898 9.71e-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.31902 0.11769 -36.697
## (10,20]|(20,30] -3.45427 0.10810 -31.954
## (20,30]|(30,40] -2.59811 0.10314 -25.189
## (30,40]|(40,50] -1.73118 0.10052 -17.222
## (40,50]|(50,60] -0.77450 0.09907 -7.817
## (50,60]|(60,70] 0.46241 0.09888 4.676
## (60,70]|(70,80] 1.43093 0.10048 14.241
## (70,80]|(80,90] 2.52825 0.10530 24.010
## (80,90]|(90,100] 3.86978 0.12133 31.894
## (20810 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 7600 -14728.70 29491.40 3147(24948) 3.24e-02 6.3e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.259 1.503
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## lockdown 0.36413 0.04369 8.334 < 2e-16 ***
## dowMonday 0.13688 0.07739 1.769 0.076956 .
## dowSaturday -0.31120 0.08152 -3.817 0.000135 ***
## dowSunday -0.09212 0.08015 -1.149 0.250416
## dowThursday 0.21003 0.07716 2.722 0.006489 **
## dowTuesday 0.41859 0.07590 5.515 0.0000000349 ***
## dowWednesday 0.18908 0.07858 2.406 0.016120 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -2.951914 0.123081 -23.983
## (10,20]|(20,30] -1.765735 0.118439 -14.908
## (20,30]|(30,40] -0.835310 0.116820 -7.150
## (30,40]|(40,50] 0.007914 0.116338 0.068
## (40,50]|(50,60] 0.961772 0.116755 8.238
## (50,60]|(60,70] 2.165868 0.118838 18.225
## (60,70]|(70,80] 2.998266 0.121813 24.614
## (70,80]|(80,90] 4.061267 0.129392 31.387
## (80,90]|(90,100] 5.031805 0.143870 34.975
## (20810 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 7600 -14862.26 29762.52 4220(25554) 8.90e-03 6.9e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.39 1.179
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.05323 0.02966 1.795 0.072724 .
## lockdown -0.63438 0.09120 -6.956 3.50e-12 ***
## dowMonday -0.25323 0.07731 -3.276 0.001054 **
## dowSaturday 0.34511 0.08156 4.231 2.32e-05 ***
## dowSunday 0.08592 0.08117 1.058 0.289849
## dowThursday -0.26915 0.07667 -3.511 0.000447 ***
## dowTuesday -0.51081 0.07596 -6.725 1.76e-11 ***
## dowWednesday -0.29746 0.07884 -3.773 0.000161 ***
## roaming_entropy:lockdown 0.08498 0.05580 1.523 0.127799
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.1888 0.1387 -30.204
## (10,20]|(20,30] -3.3243 0.1306 -25.455
## (20,30]|(30,40] -2.4679 0.1265 -19.510
## (30,40]|(40,50] -1.6002 0.1245 -12.858
## (40,50]|(50,60] -0.6427 0.1234 -5.209
## (50,60]|(60,70] 0.5949 0.1233 4.823
## (60,70]|(70,80] 1.5644 0.1248 12.537
## (70,80]|(80,90] 2.6631 0.1289 20.663
## (80,90]|(90,100] 4.0057 0.1425 28.108
## (20810 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 7382 -14399.39 28836.78 2613(15958) 4.24e-01 6.5e+05
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.394 1.18
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.0035726 0.0004885 7.313 2.60e-13 ***
## lockdown -0.4888271 0.0554096 -8.822 < 2e-16 ***
## dowMonday -0.2028731 0.0793864 -2.556 0.01060 *
## dowSaturday 0.3277272 0.0825058 3.972 7.12e-05 ***
## dowSunday 0.0870023 0.0816766 1.065 0.28678
## dowThursday -0.2503839 0.0783401 -3.196 0.00139 **
## dowTuesday -0.4681178 0.0778627 -6.012 1.83e-09 ***
## dowWednesday -0.2512366 0.0811878 -3.095 0.00197 **
## novel_locations:lockdown -0.0014190 0.0011371 -1.248 0.21207
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.1354 0.1229 -33.646
## (10,20]|(20,30] -3.2735 0.1134 -28.867
## (20,30]|(30,40] -2.4087 0.1085 -22.201
## (30,40]|(40,50] -1.5348 0.1060 -14.480
## (40,50]|(50,60] -0.5731 0.1047 -5.473
## (50,60]|(60,70] 0.6720 0.1048 6.410
## (60,70]|(70,80] 1.6561 0.1067 15.517
## (70,80]|(80,90] 2.7501 0.1117 24.614
## (80,90]|(90,100] 4.0999 0.1277 32.095
## (21028 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 2896 -5555.81 11149.61 2839(27959) 5.72e+00 5.8e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.486 1.577
## Number of groups: subject 208
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.069111 0.003050 22.661 < 2e-16 ***
## pre_covid_re 0.169707 0.003037 55.871 < 2e-16 ***
## dowMonday -0.207506 0.003089 -67.175 < 2e-16 ***
## dowSaturday 0.181929 0.003087 58.929 < 2e-16 ***
## dowSunday 0.303363 0.003085 98.328 < 2e-16 ***
## dowThursday -0.061897 0.003087 -20.052 < 2e-16 ***
## dowTuesday 0.020698 0.003086 6.708 1.98e-11 ***
## dowWednesday -0.025345 0.003087 -8.212 < 2e-16 ***
## roaming_entropy:pre_covid_re 0.073653 0.002970 24.795 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.206273 0.003068 -1045.16
## (10,20]|(20,30] -2.295589 0.003079 -745.49
## (20,30]|(30,40] -1.422244 0.003076 -462.31
## (30,40]|(40,50] -0.514941 0.003081 -167.13
## (40,50]|(50,60] 0.590722 0.042265 13.98
## (50,60]|(60,70] 1.969253 0.060396 32.61
## (60,70]|(70,80] 3.055603 0.077136 39.61
## (70,80]|(80,90] 4.253363 0.104544 40.69
## (80,90]|(90,100] 5.630289 0.160115 35.16
## (6179 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 2896 -5559.35 11156.71 2196(21200) 4.21e-01 3.0e+07
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.471 1.572
## Number of groups: subject 208
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.011635 0.009999 1.164 0.2446
## pre_covid_re 0.278096 0.313261 0.888 0.3747
## dowMonday -0.215775 0.124486 -1.733 0.0830 .
## dowSaturday 0.191241 0.129188 1.480 0.1388
## dowSunday 0.301928 0.134201 2.250 0.0245 *
## dowThursday -0.063954 0.125219 -0.511 0.6095
## dowTuesday 0.024394 0.132021 0.185 0.8534
## dowWednesday -0.037227 0.131669 -0.283 0.7774
## novel_locations:pre_covid_re -0.003046 0.003997 -0.762 0.4459
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.0939 0.7529 -4.109
## (10,20]|(20,30] -2.1832 0.7503 -2.910
## (20,30]|(30,40] -1.3104 0.7492 -1.749
## (30,40]|(40,50] -0.4024 0.7489 -0.537
## (40,50]|(50,60] 0.7007 0.7491 0.935
## (50,60]|(60,70] 2.0801 0.7497 2.775
## (60,70]|(70,80] 3.1670 0.7510 4.217
## (70,80]|(80,90] 4.3647 0.7545 5.785
## (80,90]|(90,100] 5.7404 0.7645 7.509
## (6179 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_distance_less_than_100_km
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 2850 -5439.55 10917.11 3159(28307) 4.40e-03 6.8e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.601 1.613
## Number of groups: subject 208
##
## Coefficients:
## Estimate Std. Error z value
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.19757 0.10553 1.872
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.52514 0.11011 4.769
## bs(novel_locations, knots = c(1, 26), degree = 1)3 0.47886 0.66594 0.719
## dowMonday -0.17769 0.12541 -1.417
## dowSaturday 0.24295 0.13114 1.853
## dowSunday 0.31499 0.13518 2.330
## dowThursday -0.03433 0.12615 -0.272
## dowTuesday 0.04480 0.13279 0.337
## dowWednesday -0.02407 0.13270 -0.181
## Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.0612 .
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.00000185 ***
## bs(novel_locations, knots = c(1, 26), degree = 1)3 0.4721
## dowMonday 0.1565
## dowSaturday 0.0639 .
## dowSunday 0.0198 *
## dowThursday 0.7855
## dowTuesday 0.7358
## dowWednesday 0.8560
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.6158 0.1843 -19.614
## (10,20]|(20,30] -2.6825 0.1711 -15.680
## (20,30]|(30,40] -1.7857 0.1647 -10.841
## (30,40]|(40,50] -0.8659 0.1619 -5.349
## (40,50]|(50,60] 0.2566 0.1612 1.592
## (50,60]|(60,70] 1.6600 0.1641 10.113
## (60,70]|(70,80] 2.7598 0.1704 16.197
## (70,80]|(80,90] 3.9836 0.1850 21.533
## (80,90]|(90,100] 5.3199 0.2219 23.974
## (5889 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_distance_less_than_100_km
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 4157 -8114.15 16266.30 3443(27077) 8.16e-03 9.4e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.327 1.152
## Number of groups: subject 234
##
## Coefficients:
## Estimate Std. Error z value
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -0.2022 0.1745 -1.158
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.1294 0.1577 -0.821
## bs(novel_locations, knots = c(1, 26), degree = 1)3 1.8106 0.2873 6.302
## dowMonday -0.1772 0.1095 -1.618
## dowSaturday 0.4511 0.1163 3.877
## dowSunday 0.0474 0.1116 0.425
## dowThursday -0.3167 0.1062 -2.982
## dowTuesday -0.6185 0.1025 -6.032
## dowWednesday -0.4207 0.1093 -3.850
## Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.246712
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.411863
## bs(novel_locations, knots = c(1, 26), degree = 1)3 2.95e-10 ***
## dowMonday 0.105581
## dowSaturday 0.000106 ***
## dowSunday 0.671046
## dowThursday 0.002864 **
## dowTuesday 1.62e-09 ***
## dowWednesday 0.000118 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.5084 0.2130 -21.165
## (10,20]|(20,30] -3.5497 0.1979 -17.934
## (20,30]|(30,40] -2.5880 0.1908 -13.563
## (30,40]|(40,50] -1.6690 0.1878 -8.889
## (40,50]|(50,60] -0.7442 0.1864 -3.993
## (50,60]|(60,70] 0.5182 0.1862 2.783
## (60,70]|(70,80] 1.5266 0.1879 8.126
## (70,80]|(80,90] 2.6285 0.1925 13.656
## (80,90]|(90,100] 4.0732 0.2095 19.440
## (6146 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 7382 -14391.94 28829.87 5006(30186) 5.18e-03 1.0e+04
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.397 1.182
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -0.28684 0.16737
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.09069 0.14829
## bs(novel_locations, knots = c(1, 26), degree = 1)3 3.37253 0.54039
## lockdown -0.81185 0.15026
## dowMonday -0.19465 0.07945
## dowSaturday 0.33461 0.08259
## dowSunday 0.09388 0.08210
## dowThursday -0.24270 0.07837
## dowTuesday -0.46085 0.07796
## dowWednesday -0.24776 0.08126
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.43648 0.19095
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.43881 0.17029
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -3.69524 1.47291
## z value Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -1.714 0.08657
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.612 0.54082
## bs(novel_locations, knots = c(1, 26), degree = 1)3 6.241 4.35e-10
## lockdown -5.403 6.56e-08
## dowMonday -2.450 0.01429
## dowSaturday 4.051 5.09e-05
## dowSunday 1.144 0.25283
## dowThursday -3.097 0.00196
## dowTuesday -5.912 3.39e-09
## dowWednesday -3.049 0.00230
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 2.286 0.02226
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 2.577 0.00997
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -2.509 0.01211
##
## 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.3292 0.1847 -23.442
## (10,20]|(20,30] -3.4662 0.1782 -19.450
## (20,30]|(30,40] -2.6004 0.1750 -14.863
## (30,40]|(40,50] -1.7255 0.1734 -9.953
## (40,50]|(50,60] -0.7625 0.1725 -4.419
## (50,60]|(60,70] 0.4840 0.1725 2.807
## (60,70]|(70,80] 1.4692 0.1735 8.470
## (70,80]|(80,90] 2.5643 0.1765 14.532
## (80,90]|(90,100] 3.9148 0.1870 20.930
## (21028 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 + distance + (1 | subject)
## data: distance_less_than_100_km
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 7007 -13577.13 27202.26 3995(21252) 9.12e-02 4.3e+05
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.482 1.217
## Number of groups: subject 234
##
## Coefficients:
## Estimate
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -0.2255064
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.2412445
## bs(novel_locations, knots = c(1, 26), degree = 1)3 2.0718822
## lockdown -0.8248059
## dowMonday -0.1913897
## dowSaturday 0.3348037
## dowSunday 0.1286502
## dowThursday -0.2075564
## dowTuesday -0.3981668
## dowWednesday -0.2686236
## distance -0.0001218
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.3753480
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.6375293
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -2.1984363
## Std. Error z value
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.1688531 -1.336
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.1521086 -1.586
## bs(novel_locations, knots = c(1, 26), degree = 1)3 0.4880736 4.245
## lockdown 0.1511686 -5.456
## dowMonday 0.0816188 -2.345
## dowSaturday 0.0860328 3.892
## dowSunday 0.0854470 1.506
## dowThursday 0.0805565 -2.577
## dowTuesday 0.0802596 -4.961
## dowWednesday 0.0836655 -3.211
## distance 0.0023445 -0.052
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.1926185 1.949
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.1762922 3.616
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown 0.6889188 -3.191
## Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.181707
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.112739
## bs(novel_locations, knots = c(1, 26), degree = 1)3 0.0000218574 ***
## lockdown 0.0000000486 ***
## dowMonday 0.019031 *
## dowSaturday 0.0000995931 ***
## dowSunday 0.132166
## dowThursday 0.009980 **
## dowTuesday 0.0000007014 ***
## dowWednesday 0.001324 **
## distance 0.958556
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.051336 .
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.000299 ***
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown 0.001417 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.3931 0.1882 -23.345
## (10,20]|(20,30] -3.4965 0.1812 -19.297
## (20,30]|(30,40] -2.6086 0.1778 -14.674
## (30,40]|(40,50] -1.7233 0.1761 -9.784
## (40,50]|(50,60] -0.7506 0.1753 -4.282
## (50,60]|(60,70] 0.5230 0.1752 2.985
## (60,70]|(70,80] 1.5355 0.1764 8.707
## (70,80]|(80,90] 2.6532 0.1798 14.758
## (80,90]|(90,100] 4.0472 0.1922 21.060
## (13899 observations deleted due to missingness)
## df BIC
## PA_NL_PreLD_notch 19 11030.26
## PA_NL_PostLD_notch 19 16386.62
## PA_NL_LD_notch 23 28988.73
## PA_NL_LD_notch_controldistance_100KM 24 27366.77
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 + distance + dow + (1 | subject)
## data: distance_less_than_100_km
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 7007 -13586.72 27221.45 4036(21172) 1.72e-02 2.4e+06
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.477 1.215
## Number of groups: subject 234
##
## Coefficients:
## Estimate
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 -0.033710
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 -0.113235
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 2.079031
## lockdown -0.690859
## distance 0.004018
## dowMonday -0.205631
## dowSaturday 0.323723
## dowSunday 0.110579
## dowThursday -0.218383
## dowTuesday -0.416785
## dowWednesday -0.281256
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.125099
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.389118
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -3.504748
## Std. Error
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.444823
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.409181
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.721311
## lockdown 0.472898
## distance 0.001980
## dowMonday 0.081523
## dowSaturday 0.086105
## dowSunday 0.085327
## dowThursday 0.080546
## dowTuesday 0.080172
## dowWednesday 0.083498
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.522230
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.474230
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 1.367077
## z value
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 -0.076
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 -0.277
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 2.882
## lockdown -1.461
## distance 2.030
## dowMonday -2.522
## dowSaturday 3.760
## dowSunday 1.296
## dowThursday -2.711
## dowTuesday -5.199
## dowWednesday -3.368
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.240
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.821
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -2.564
## Pr(>|z|)
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.939591
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.781983
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.003948
## lockdown 0.144041
## distance 0.042399
## dowMonday 0.011657
## dowSaturday 0.000170
## dowSunday 0.194998
## dowThursday 0.006702
## dowTuesday 0.000000201
## dowWednesday 0.000756
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.810680
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.411917
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.010357
##
## 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
## distance *
## 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] -4.2668 0.4254 -10.030
## (10,20]|(20,30] -3.3723 0.4225 -7.982
## (20,30]|(30,40] -2.4866 0.4210 -5.907
## (30,40]|(40,50] -1.6030 0.4203 -3.814
## (40,50]|(50,60] -0.6314 0.4200 -1.503
## (50,60]|(60,70] 0.6407 0.4200 1.526
## (60,70]|(70,80] 1.6518 0.4204 3.929
## (70,80]|(80,90] 2.7685 0.4218 6.563
## (80,90]|(90,100] 4.1619 0.4272 9.742
## (13899 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 + distance + (1 | subject)
## data: distance_less_than_100_km
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 7007 -13576.59 27201.18 4254(22364) 4.69e-02 1.7e+06
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.482 1.218
## Number of groups: subject 234
##
## Coefficients:
## Estimate
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 -0.268014
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 -0.256007
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 2.900872
## lockdown -1.615995
## dowMonday -0.199366
## dowSaturday 0.315171
## dowSunday 0.119004
## dowThursday -0.216217
## dowTuesday -0.411796
## dowWednesday -0.277158
## distance 0.001183
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 1.087866
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 1.204510
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -1.417898
## Std. Error
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.345981
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.314009
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.632054
## lockdown 0.568925
## dowMonday 0.081596
## dowSaturday 0.086129
## dowSunday 0.085465
## dowThursday 0.080533
## dowTuesday 0.080203
## dowWednesday 0.083531
## distance 0.001788
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.621001
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.566071
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.852220
## z value
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 -0.775
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 -0.815
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 4.590
## lockdown -2.840
## dowMonday -2.443
## dowSaturday 3.659
## dowSunday 1.392
## dowThursday -2.685
## dowTuesday -5.134
## dowWednesday -3.318
## distance 0.661
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 1.752
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 2.128
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -1.664
## Pr(>|z|)
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.438547
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.414909
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.000004441
## lockdown 0.004505
## dowMonday 0.014553
## dowSaturday 0.000253
## dowSunday 0.163790
## dowThursday 0.007257
## dowTuesday 0.000000283
## dowWednesday 0.000907
## distance 0.508336
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.079809
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.033350
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.096158
##
## 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 ***
## distance
## 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.4601 0.3346 -13.329
## (10,20]|(20,30] -3.5644 0.3310 -10.769
## (20,30]|(30,40] -2.6769 0.3292 -8.132
## (30,40]|(40,50] -1.7916 0.3283 -5.457
## (40,50]|(50,60] -0.8190 0.3280 -2.497
## (50,60]|(60,70] 0.4548 0.3279 1.387
## (60,70]|(70,80] 1.4679 0.3284 4.470
## (70,80]|(80,90] 2.5866 0.3300 7.838
## (80,90]|(90,100] 3.9825 0.3368 11.826
## (13899 observations deleted due to missingness)
## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula: PA_avg ~ bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1) *
## lockdown + distance + dow + (1 | subject)
## data: distance_less_than_100_km
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 7007 -13579.00 27201.99 3376(18194) 1.20e-01 2.0e+06
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.479 1.216
## Number of groups: subject 234
##
## Coefficients:
## Estimate
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1 -2.2744396
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2 -2.3718410
## lockdown -3.1557158
## distance -0.0002255
## dowMonday -0.1914640
## dowSaturday 0.3326636
## dowSunday 0.1261414
## dowThursday -0.2081687
## dowTuesday -0.4015485
## dowWednesday -0.2660535
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1:lockdown 2.5268934
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2:lockdown 3.0033808
## Std. Error
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1 0.4664470
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2 0.4709670
## lockdown 0.6702249
## distance 0.0023458
## dowMonday 0.0815850
## dowSaturday 0.0860273
## dowSunday 0.0854357
## dowThursday 0.0805504
## dowTuesday 0.0802038
## dowWednesday 0.0836347
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1:lockdown 0.6685571
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2:lockdown 0.7205404
## z value
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1 -4.876
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2 -5.036
## lockdown -4.708
## distance -0.096
## dowMonday -2.347
## dowSaturday 3.867
## dowSunday 1.476
## dowThursday -2.584
## dowTuesday -5.007
## dowWednesday -3.181
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1:lockdown 3.780
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2:lockdown 4.168
## Pr(>|z|)
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1 0.000001082
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2 0.000000475
## lockdown 0.000002496
## distance 0.923413
## dowMonday 0.018935
## dowSaturday 0.000110
## dowSunday 0.139823
## dowThursday 0.009757
## dowTuesday 0.000000554
## dowWednesday 0.001467
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1:lockdown 0.000157
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2:lockdown 0.000030697
##
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1 ***
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2 ***
## lockdown ***
## distance
## dowMonday *
## dowSaturday ***
## dowSunday
## dowThursday **
## dowTuesday ***
## dowWednesday **
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)1:lockdown ***
## bs(znovel_locations, knots = c(-0.73, -0.3), degree = 1)2:lockdown ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -6.5130 0.4790 -13.596
## (10,20]|(20,30] -5.6172 0.4763 -11.792
## (20,30]|(30,40] -4.7298 0.4748 -9.962
## (30,40]|(40,50] -3.8447 0.4739 -8.113
## (40,50]|(50,60] -2.8721 0.4732 -6.070
## (50,60]|(60,70] -1.5989 0.4724 -3.385
## (60,70]|(70,80] -0.5868 0.4722 -1.243
## (70,80]|(80,90] 0.5305 0.4729 1.122
## (80,90]|(90,100] 1.9241 0.4769 4.034
## (13899 observations deleted due to missingness)
## df BIC
## Z_within_sub 24 27385.96
## Z_within_sub_LD 24 27365.69
## Z_NL 22 27352.80
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 7600 -14928.17 29882.35 2090(16156) 1.20e-02 3.0e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.129 1.062
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## depression_mean -0.09050 0.01572 -5.756 8.61e-09 ***
## lockdown -0.46866 0.06265 -7.480 7.41e-14 ***
## depression_mean:lockdown -0.02420 0.00920 -2.630 0.00854 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.5745 0.1288 -35.508
## (10,20]|(20,30] -3.7118 0.1200 -30.940
## (20,30]|(30,40] -2.8582 0.1153 -24.792
## (30,40]|(40,50] -1.9958 0.1126 -17.717
## (40,50]|(50,60] -1.0482 0.1110 -9.441
## (50,60]|(60,70] 0.1697 0.1105 1.536
## (60,70]|(70,80] 1.1206 0.1116 10.041
## (70,80]|(80,90] 2.2011 0.1157 19.023
## (80,90]|(90,100] 3.5280 0.1300 27.129
## (20810 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 7382 -14562.88 29159.76 642(3210) 3.80e+03 3.5e+07
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1 1
## Number of groups: subject 235
##
## Coefficients:
## Estimate Std. Error z value
## novel_locations 0.00280109 0.00067576 4.145
## depression_mean -0.08519778 0.01577509 -5.401
## lockdown -0.27595852 0.07643967 -3.610
## novel_locations:depression_mean 0.00013694 0.00009661 1.417
## novel_locations:lockdown 0.00225339 0.00151777 1.485
## depression_mean:lockdown -0.01494193 0.01138438 -1.312
## novel_locations:depression_mean:lockdown -0.00063002 0.00022090 -2.852
## Pr(>|z|)
## novel_locations 0.0000339666 ***
## depression_mean 0.0000000664 ***
## lockdown 0.000306 ***
## novel_locations:depression_mean 0.156365
## novel_locations:lockdown 0.137629
## depression_mean:lockdown 0.189353
## novel_locations:depression_mean:lockdown 0.004344 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.8728 0.1261 -30.710
## (10,20]|(20,30] -3.0738 0.1180 -26.039
## (20,30]|(30,40] -2.2918 0.1140 -20.110
## (30,40]|(40,50] -1.5235 0.1118 -13.624
## (40,50]|(50,60] -0.7038 0.1107 -6.360
## (50,60]|(60,70] 0.3168 0.1105 2.868
## (60,70]|(70,80] 1.1038 0.1115 9.903
## (70,80]|(80,90] 1.9998 0.1145 17.465
## (80,90]|(90,100] 3.1320 0.1252 25.023
## (21028 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 2557 -4931.45 9888.90 1132(9950) 2.04e+00 1.9e+06
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.317 1.522
## Number of groups: subject 182
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.0049387 0.0012196 4.050 0.0000513 ***
## anx_change -0.0223949 0.0295411 -0.758 0.4484
## novel_locations:anx_change -0.0005875 0.0002843 -2.066 0.0388 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.72124 0.15813 -23.533
## (10,20]|(20,30] -2.78403 0.14060 -19.801
## (20,30]|(30,40] -1.91157 0.13195 -14.487
## (30,40]|(40,50] -1.02699 0.12764 -8.046
## (40,50]|(50,60] 0.03444 0.12619 0.273
## (50,60]|(60,70] 1.41614 0.12935 10.948
## (60,70]|(70,80] 2.52694 0.13780 18.338
## (70,80]|(80,90] 3.68205 0.15692 23.464
## (80,90]|(90,100] 5.06263 0.20665 24.499
## (6518 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 29972.512
## d2 17 29277.171
## d3 13 9964.901
BIC Analysis
## df BIC
## complex_PA_RE 21 18879.458
## simple_PA_RE 17 18692.598
## complex_PA_NL 21 17957.667
## simple_PA_NL 17 17769.494
## RE_LD 10 50756.718
## NL_LD 10 201440.118
## PA_LD 17 29887.290
## NA_LD 17 29609.310
## PA_RE_LD 19 29894.301
## PA_NL_LD 19 28968.009
## PA_RE_PreLDRE_PostLD 19 11263.065
## PA_RE_PreLDNL_PostLD 19 11270.157
## PA_NL_PreLD_notch 19 11030.255
## PA_NL_PostLD_notch 19 16386.619
## PA_NL_LD_notch 23 28988.729
## PA_NL_LD_notch_controldistance_100KM 24 27366.767
## Z_within_sub 24 27385.958
## Z_within_sub_LD 24 27365.694
## Z_NL 22 27352.797
## d1 13 29972.512
## d2 17 29277.171
## d3 13 9964.901