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 3463 -6917.43 13876.86 1001(6008) 1.52e+02 7.4e+07
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1 1
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.0926975 0.0443223 2.091 0.036488 *
## time_of_day 0.0036969 0.0017215 2.148 0.031749 *
## precipi 0.1402720 0.2181033 0.643 0.520130
## mean_temp -0.0122209 0.0044783 -2.729 0.006354 **
## distance 0.0005219 0.0001460 3.574 0.000352 ***
## dowMonday -0.2186004 0.1295690 -1.687 0.091577 .
## dowSaturday 0.5007026 0.1235249 4.053 0.0000504671 ***
## dowSunday -0.0793837 0.1378905 -0.576 0.564817
## dowThursday -0.3166364 0.1209959 -2.617 0.008873 **
## dowTuesday -0.6378163 0.1169935 -5.452 0.0000000499 ***
## dowWednesday -0.6673580 0.1264936 -5.276 0.0000001322 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.4545 0.3998 -11.142
## (10,20]|(20,30] -3.5563 0.3918 -9.077
## (20,30]|(30,40] -2.6860 0.3884 -6.916
## (30,40]|(40,50] -1.9066 0.3871 -4.925
## (40,50]|(50,60] -1.1095 0.3864 -2.871
## (50,60]|(60,70] -0.2445 0.3861 -0.633
## (60,70]|(70,80] 0.5715 0.3864 1.479
## (70,80]|(80,90] 1.5239 0.3875 3.933
## (80,90]|(90,100] 2.7257 0.3934 6.929
## (5439 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 3484 -6897.04 13828.08 3151(21075) 5.20e-03 6.7e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.359 1.166
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.17821 0.04475 3.983 6.82e-05 ***
## dowMonday -0.15062 0.11985 -1.257 0.20882
## dowSaturday 0.62457 0.12483 5.003 5.63e-07 ***
## dowSunday 0.07824 0.13324 0.587 0.55710
## dowThursday -0.37160 0.12068 -3.079 0.00208 **
## dowTuesday -0.66879 0.11429 -5.852 4.86e-09 ***
## dowWednesday -0.85400 0.12545 -6.807 9.95e-12 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.1497 0.2087 -19.884
## (10,20]|(20,30] -3.1664 0.1900 -16.663
## (20,30]|(30,40] -2.1776 0.1811 -12.026
## (30,40]|(40,50] -1.2631 0.1774 -7.121
## (40,50]|(50,60] -0.2876 0.1758 -1.636
## (50,60]|(60,70] 0.7870 0.1762 4.468
## (60,70]|(70,80] 1.8079 0.1788 10.110
## (70,80]|(80,90] 2.9630 0.1848 16.031
## (80,90]|(90,100] 4.3287 0.2029 21.334
## (5418 observations deleted due to missingness)
## df BIC
## complex_PA_RE 21 14006.01
## simple_PA_RE 17 13932.74
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 3250 -6381.74 12805.49 3267(16946) 6.30e+00 7.0e+07
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.401 1.184
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.0056565 0.0007727 7.320 2.47e-13 ***
## time_of_day 0.0062189 0.0018704 3.325 0.000884 ***
## precipi -0.0253801 0.2291122 -0.111 0.911794
## mean_temp -0.0078137 0.0048367 -1.616 0.106200
## distance 0.0003744 0.0001596 2.346 0.018979 *
## dowMonday -0.1785746 0.1449645 -1.232 0.218005
## dowSaturday 0.4041384 0.1326614 3.046 0.002316 **
## dowSunday -0.1303932 0.1462393 -0.892 0.372584
## dowThursday -0.3579095 0.1349437 -2.652 0.007995 **
## dowTuesday -0.6728670 0.1311988 -5.129 2.92e-07 ***
## dowWednesday -0.8008243 0.1447149 -5.534 3.13e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.8114 0.4190 -11.482
## (10,20]|(20,30] -3.8141 0.4085 -9.336
## (20,30]|(30,40] -2.7927 0.4037 -6.917
## (30,40]|(40,50] -1.8527 0.4019 -4.610
## (40,50]|(50,60] -0.8626 0.4011 -2.151
## (50,60]|(60,70] 0.2335 0.4010 0.582
## (60,70]|(70,80] 1.2936 0.4017 3.220
## (70,80]|(80,90] 2.4492 0.4038 6.065
## (80,90]|(90,100] 3.8295 0.4122 9.291
## (5652 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 3267 -6423.77 12881.54 1894(9546) 8.96e-02 5.3e+05
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.394 1.181
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.006573 0.000715 9.193 < 2e-16 ***
## dowMonday -0.095298 0.127132 -0.750 0.453499
## dowSaturday 0.426480 0.127817 3.337 0.000848 ***
## dowSunday -0.016011 0.134135 -0.119 0.904986
## dowThursday -0.332058 0.127749 -2.599 0.009341 **
## dowTuesday -0.596107 0.121750 -4.896 0.0000009774 ***
## dowWednesday -0.776281 0.136046 -5.706 0.0000000116 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.3408 0.1817 -23.894
## (10,20]|(20,30] -3.3466 0.1566 -21.371
## (20,30]|(30,40] -2.3313 0.1442 -16.166
## (30,40]|(40,50] -1.3968 0.1389 -10.055
## (40,50]|(50,60] -0.4090 0.1366 -2.995
## (50,60]|(60,70] 0.6801 0.1368 4.971
## (60,70]|(70,80] 1.7356 0.1403 12.371
## (70,80]|(80,90] 2.8946 0.1483 19.514
## (80,90]|(90,100] 4.2733 0.1715 24.919
## (5635 observations deleted due to missingness)
## df BIC
## complex_PA_NL 21 12933.3
## simple_PA_NL 17 12985.1
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: 53792.8
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.5757 -0.6437 -0.0393 0.5547 8.1706
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 0.1139 0.3375
## Residual 0.7119 0.8437
## Number of obs: 21265, groups: subject, 230
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.68473 0.02795 498.74456 96.055 < 2e-16 ***
## lockdown -1.31063 0.01198 21110.06033 -109.390 < 2e-16 ***
## dowMonday -0.20937 0.02168 21033.73020 -9.656 < 2e-16 ***
## dowSaturday -0.19673 0.02195 21033.98418 -8.961 < 2e-16 ***
## dowSunday -0.37510 0.02188 21034.50705 -17.144 < 2e-16 ***
## dowThursday -0.03920 0.02139 21033.99174 -1.832 0.066925 .
## dowTuesday -0.19120 0.02147 21034.24854 -8.905 < 2e-16 ***
## dowWednesday -0.07390 0.02147 21033.60065 -3.443 0.000577 ***
## ---
## 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.258
## dowMonday -0.387 0.013
## dowSaturday -0.381 0.007 0.489
## dowSunday -0.385 0.017 0.491 0.485
## dowThursday -0.386 -0.012 0.502 0.496 0.497
## dowTuesday -0.395 0.029 0.501 0.494 0.496 0.507
## dowWednesdy -0.386 -0.008 0.500 0.494 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: 216094
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.0030 -0.5144 -0.2301 0.1980 16.7386
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 307.2 17.53
## Residual 3331.8 57.72
## Number of obs: 19693, groups: subject, 230
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 64.5273 1.6858 807.0832 38.276 < 2e-16 ***
## lockdown -20.8291 0.8736 19585.6003 -23.842 < 2e-16 ***
## dowMonday -15.2610 1.5420 19466.6824 -9.897 < 2e-16 ***
## dowSaturday 3.6121 1.5621 19467.1129 2.312 0.020768 *
## dowSunday -6.0837 1.5568 19468.2680 -3.908 9.35e-05 ***
## dowThursday -5.5457 1.5195 19467.1574 -3.650 0.000263 ***
## dowTuesday -13.0101 1.5260 19467.6760 -8.526 < 2e-16 ***
## dowWednesday -10.2652 1.5251 19466.4204 -6.731 1.74e-11 ***
## ---
## 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.337
## dowMonday -0.457 0.014
## dowSaturday -0.449 0.007 0.489
## dowSunday -0.454 0.018 0.491 0.484
## dowThursday -0.455 -0.012 0.502 0.496 0.497
## dowTuesday -0.468 0.032 0.501 0.494 0.496 0.507
## dowWednesdy -0.455 -0.008 0.500 0.494 0.496 0.508 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 9141 -18082.57 36199.15 3611(22311) 2.94e-02 5.7e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.424 1.193
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## lockdown -0.48320 0.03855 -12.534 < 2e-16 ***
## dowMonday -0.31123 0.07008 -4.441 8.97e-06 ***
## dowSaturday 0.22214 0.07238 3.069 0.00215 **
## dowSunday 0.13866 0.07400 1.874 0.06098 .
## dowThursday -0.28215 0.07079 -3.986 6.72e-05 ***
## dowTuesday -0.46187 0.07059 -6.543 6.02e-11 ***
## dowWednesday -0.21268 0.07248 -2.934 0.00334 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.28275 0.11430 -37.469
## (10,20]|(20,30] -3.37236 0.10623 -31.746
## (20,30]|(30,40] -2.51554 0.10239 -24.569
## (30,40]|(40,50] -1.61767 0.10034 -16.121
## (40,50]|(50,60] -0.64164 0.09927 -6.464
## (50,60]|(60,70] 0.41331 0.09918 4.167
## (60,70]|(70,80] 1.41850 0.10050 14.114
## (70,80]|(80,90] 2.49759 0.10447 23.907
## (80,90]|(90,100] 3.83174 0.11809 32.448
## (17955 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 9141 -17809.70 35653.40 3387(33663) 9.42e-03 7.7e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.341 1.53
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## lockdown 0.31440 0.03847 8.172 3.02e-16 ***
## dowMonday 0.19761 0.07035 2.809 0.00497 **
## dowSaturday -0.15430 0.07287 -2.117 0.03423 *
## dowSunday -0.12521 0.07431 -1.685 0.09198 .
## dowThursday 0.23112 0.07147 3.234 0.00122 **
## dowTuesday 0.44665 0.07100 6.291 3.15e-10 ***
## dowWednesday 0.15549 0.07235 2.149 0.03163 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -2.9099 0.1235 -23.564
## (10,20]|(20,30] -1.7013 0.1196 -14.230
## (20,30]|(30,40] -0.7293 0.1182 -6.169
## (30,40]|(40,50] 0.1566 0.1179 1.328
## (40,50]|(50,60] 1.1278 0.1184 9.526
## (50,60]|(60,70] 2.0678 0.1197 17.270
## (60,70]|(70,80] 2.9379 0.1221 24.067
## (70,80]|(80,90] 3.9875 0.1277 31.228
## (80,90]|(90,100] 4.9555 0.1385 35.784
## (17955 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 7651 -15154.66 30347.33 4173(25019) 9.21e-03 1.0e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.377 1.173
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.12884 0.04043 3.187 0.001440 **
## lockdown -0.44045 0.11403 -3.863 0.000112 ***
## dowMonday -0.24380 0.07683 -3.173 0.001507 **
## dowSaturday 0.31771 0.08155 3.896 9.78e-05 ***
## dowSunday 0.15779 0.08171 1.931 0.053463 .
## dowThursday -0.29273 0.07680 -3.811 0.000138 ***
## dowTuesday -0.48027 0.07640 -6.287 3.25e-10 ***
## dowWednesday -0.23404 0.07832 -2.988 0.002805 **
## roaming_entropy:lockdown 0.18717 0.05013 3.734 0.000189 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.8989 0.1596 -24.434
## (10,20]|(20,30] -3.0291 0.1528 -19.823
## (20,30]|(30,40] -2.1814 0.1495 -14.593
## (30,40]|(40,50] -1.2848 0.1478 -8.692
## (40,50]|(50,60] -0.3103 0.1471 -2.109
## (50,60]|(60,70] 0.7325 0.1473 4.972
## (60,70]|(70,80] 1.7411 0.1488 11.705
## (70,80]|(80,90] 2.8483 0.1522 18.718
## (80,90]|(90,100] 4.1904 0.1633 25.654
## (19445 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 7434 -14696.62 29431.24 2672(15817) 5.11e-01 1.0e+06
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.38 1.175
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.0062848 0.0006501 9.667 < 2e-16 ***
## lockdown -0.2195308 0.0568865 -3.859 0.000114 ***
## dowMonday -0.2486938 0.0785489 -3.166 0.001545 **
## dowSaturday 0.2204203 0.0821296 2.684 0.007279 **
## dowSunday 0.1291163 0.0818943 1.577 0.114883
## dowThursday -0.2842230 0.0785596 -3.618 0.000297 ***
## dowTuesday -0.4369405 0.0783031 -5.580 0.000000024 ***
## dowWednesday -0.1723225 0.0807807 -2.133 0.032907 *
## novel_locations:lockdown -0.0023953 0.0008218 -2.915 0.003561 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.9201 0.1255 -31.231
## (10,20]|(20,30] -3.0531 0.1166 -26.192
## (20,30]|(30,40] -2.2037 0.1121 -19.658
## (30,40]|(40,50] -1.3018 0.1098 -11.859
## (40,50]|(50,60] -0.3246 0.1087 -2.986
## (50,60]|(60,70] 0.7217 0.1090 6.620
## (60,70]|(70,80] 1.7452 0.1110 15.725
## (70,80]|(80,90] 2.8514 0.1157 24.635
## (80,90]|(90,100] 4.1980 0.1305 32.161
## (19662 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 4167 -8166.63 16371.27 3852(33251) 5.73e-02 3.6e+04
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.021 1.422
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## roaming_entropy 0.420812 0.237103 1.775 0.075929 .
## pre_covid_re 0.133628 0.267607 0.499 0.617536
## dowMonday -0.337879 0.102395 -3.300 0.000968 ***
## dowSaturday -0.003673 0.110292 -0.033 0.973431
## dowSunday 0.241586 0.104983 2.301 0.021381 *
## dowThursday -0.261812 0.101233 -2.586 0.009703 **
## dowTuesday -0.289497 0.106418 -2.720 0.006521 **
## dowWednesday 0.127037 0.101584 1.251 0.211094
## roaming_entropy:pre_covid_re -0.030658 0.091620 -0.335 0.737913
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.1827 0.6881 -4.626
## (10,20]|(20,30] -2.3054 0.6860 -3.360
## (20,30]|(30,40] -1.4770 0.6852 -2.156
## (30,40]|(40,50] -0.5191 0.6849 -0.758
## (40,50]|(50,60] 0.5346 0.6850 0.780
## (50,60]|(60,70] 1.6385 0.6854 2.390
## (60,70]|(70,80] 2.7107 0.6863 3.950
## (70,80]|(80,90] 3.8327 0.6882 5.569
## (80,90]|(90,100] 5.2210 0.6938 7.525
## (14027 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 4167 -8189.90 16417.81 2170(16590) 8.13e+00 1.1e+08
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.999 1.414
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.006912 0.004420 1.564 0.117850
## pre_covid_re 0.186453 0.254909 0.731 0.464506
## dowMonday -0.389343 0.102038 -3.816 0.000136 ***
## dowSaturday -0.019145 0.110121 -0.174 0.861977
## dowSunday 0.219924 0.104898 2.097 0.036032 *
## dowThursday -0.287855 0.101300 -2.842 0.004489 **
## dowTuesday -0.298304 0.106194 -2.809 0.004969 **
## dowWednesday 0.116675 0.101723 1.147 0.251386
## novel_locations:pre_covid_re -0.001091 0.001677 -0.651 0.515178
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.3050 0.6568 -5.032
## (10,20]|(20,30] -2.4300 0.6547 -3.712
## (20,30]|(30,40] -1.6051 0.6537 -2.455
## (30,40]|(40,50] -0.6522 0.6533 -0.998
## (40,50]|(50,60] 0.3953 0.6534 0.605
## (50,60]|(60,70] 1.4923 0.6537 2.283
## (60,70]|(70,80] 2.5570 0.6545 3.907
## (70,80]|(80,90] 3.6695 0.6563 5.591
## (80,90]|(90,100] 5.0487 0.6620 7.626
## (14027 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 3267 -6423.65 12885.30 3087(15817) 1.84e-03 2.1e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.394 1.181
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.10691 0.28491 0.375
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.20827 0.26872 0.775
## bs(novel_locations, knots = c(1, 26), degree = 1)3 6.68181 0.79461 8.409
## dowMonday -0.09450 0.12711 -0.743
## dowSaturday 0.42374 0.12838 3.301
## dowSunday -0.02093 0.13603 -0.154
## dowThursday -0.33018 0.12780 -2.584
## dowTuesday -0.59699 0.12175 -4.903
## dowWednesday -0.77356 0.13619 -5.680
## Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.707490
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.438317
## bs(novel_locations, knots = c(1, 26), degree = 1)3 < 2e-16 ***
## dowMonday 0.457195
## dowSaturday 0.000964 ***
## dowSunday 0.877712
## dowThursday 0.009780 **
## dowTuesday 0.0000009418 ***
## dowWednesday 0.0000000135 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.2918 0.3164 -13.563
## (10,20]|(20,30] -3.2975 0.3026 -10.897
## (20,30]|(30,40] -2.2820 0.2962 -7.703
## (30,40]|(40,50] -1.3473 0.2937 -4.587
## (40,50]|(50,60] -0.3592 0.2928 -1.227
## (50,60]|(60,70] 0.7300 0.2929 2.493
## (60,70]|(70,80] 1.7855 0.2944 6.065
## (70,80]|(80,90] 2.9444 0.2983 9.870
## (80,90]|(90,100] 4.3232 0.3105 13.925
## (5635 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 4167 -8163.85 16365.70 3285(26396) 6.81e-03 8.4e+02
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 2.023 1.422
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.16665 0.08930 1.866
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.69152 0.08240 8.392
## bs(novel_locations, knots = c(1, 26), degree = 1)3 2.36076 0.58456 4.038
## dowMonday -0.33374 0.10241 -3.259
## dowSaturday -0.01364 0.11034 -0.124
## dowSunday 0.24314 0.10501 2.315
## dowThursday -0.25287 0.10134 -2.495
## dowTuesday -0.28839 0.10647 -2.709
## dowWednesday 0.13816 0.10183 1.357
## Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.06203 .
## bs(novel_locations, knots = c(1, 26), degree = 1)2 < 2e-16 ***
## bs(novel_locations, knots = c(1, 26), degree = 1)3 0.0000538 ***
## dowMonday 0.00112 **
## dowSaturday 0.90163
## dowSunday 0.02059 *
## dowThursday 0.01259 *
## dowTuesday 0.00676 **
## dowWednesday 0.17485
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.4923 0.1522 -22.941
## (10,20]|(20,30] -2.6165 0.1416 -18.484
## (20,30]|(30,40] -1.7882 0.1364 -13.107
## (30,40]|(40,50] -0.8285 0.1339 -6.188
## (40,50]|(50,60] 0.2279 0.1334 1.708
## (50,60]|(60,70] 1.3344 0.1350 9.881
## (60,70]|(70,80] 2.4066 0.1389 17.320
## (70,80]|(80,90] 3.5259 0.1472 23.952
## (80,90]|(90,100] 4.9110 0.1716 28.622
## (14027 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 7434 -14674.41 29394.82 5016(30295) 1.27e-02 7.0e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.385 1.177
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -0.03836 0.26474
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.08386 0.24597
## bs(novel_locations, knots = c(1, 26), degree = 1)3 6.85294 0.75173
## lockdown -0.64335 0.24385
## dowMonday -0.22101 0.07869
## dowSaturday 0.21674 0.08244
## dowSunday 0.12593 0.08241
## dowThursday -0.26259 0.07864
## dowTuesday -0.42931 0.07843
## dowWednesday -0.15647 0.08089
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.20111 0.27647
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.66340 0.25410
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -4.41136 0.93406
## z value Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 -0.145 0.88478
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.341 0.73314
## bs(novel_locations, knots = c(1, 26), degree = 1)3 9.116 < 2e-16
## lockdown -2.638 0.00833
## dowMonday -2.808 0.00498
## dowSaturday 2.629 0.00856
## dowSunday 1.528 0.12647
## dowThursday -3.339 0.00084
## dowTuesday -5.474 0.000000044
## dowWednesday -1.934 0.05307
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.727 0.46697
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 2.611 0.00903
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -4.723 0.000002326
##
## 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.1120 0.2685 -15.313
## (10,20]|(20,30] -3.2442 0.2641 -12.283
## (20,30]|(30,40] -2.3923 0.2619 -9.134
## (30,40]|(40,50] -1.4869 0.2608 -5.700
## (40,50]|(50,60] -0.5055 0.2605 -1.941
## (50,60]|(60,70] 0.5449 0.2606 2.091
## (60,70]|(70,80] 1.5718 0.2613 6.015
## (70,80]|(80,90] 2.6807 0.2634 10.179
## (80,90]|(90,100] 4.0291 0.2702 14.912
## (19662 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 7038 -13836.79 27719.58 5208(30536) 2.06e-02 4.6e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.467 1.211
## Number of groups: subject 230
##
## Coefficients:
## Estimate
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.007711
## bs(novel_locations, knots = c(1, 26), degree = 1)2 -0.123586
## bs(novel_locations, knots = c(1, 26), degree = 1)3 2.443094
## lockdown -0.624915
## dowMonday -0.214249
## dowSaturday 0.257875
## dowSunday 0.150284
## dowThursday -0.238754
## dowTuesday -0.386549
## dowWednesday -0.213381
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.188974
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.589059
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown -0.047602
## Std. Error z value
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.265877 0.029
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.248037 -0.498
## bs(novel_locations, knots = c(1, 26), degree = 1)3 0.388523 6.288
## lockdown 0.244701 -2.554
## dowMonday 0.080815 -2.651
## dowSaturday 0.085778 3.006
## dowSunday 0.085675 1.754
## dowThursday 0.080744 -2.957
## dowTuesday 0.080454 -4.805
## dowWednesday 0.083018 -2.570
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.277748 0.680
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.258186 2.282
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown 0.544378 -0.087
## Pr(>|z|)
## bs(novel_locations, knots = c(1, 26), degree = 1)1 0.97686
## bs(novel_locations, knots = c(1, 26), degree = 1)2 0.61831
## bs(novel_locations, knots = c(1, 26), degree = 1)3 3.21e-10 ***
## lockdown 0.01066 *
## dowMonday 0.00802 **
## dowSaturday 0.00264 **
## dowSunday 0.07941 .
## dowThursday 0.00311 **
## dowTuesday 1.55e-06 ***
## dowWednesday 0.01016 *
## bs(novel_locations, knots = c(1, 26), degree = 1)1:lockdown 0.49626
## bs(novel_locations, knots = c(1, 26), degree = 1)2:lockdown 0.02252 *
## bs(novel_locations, knots = c(1, 26), degree = 1)3:lockdown 0.93032
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.1365 0.2711 -15.260
## (10,20]|(20,30] -3.2435 0.2664 -12.177
## (20,30]|(30,40] -2.3789 0.2641 -9.008
## (30,40]|(40,50] -1.4621 0.2630 -5.559
## (40,50]|(50,60] -0.4730 0.2626 -1.801
## (50,60]|(60,70] 0.6005 0.2627 2.286
## (60,70]|(70,80] 1.6492 0.2636 6.257
## (70,80]|(80,90] 2.7712 0.2659 10.420
## (80,90]|(90,100] 4.1513 0.2739 15.155
## (20058 observations deleted due to missingness)
## df BIC
## PA_NL_PreLD_notch 19 13001.04
## PA_NL_PostLD_notch 19 16486.07
## PA_NL_LD_notch 23 29553.83
## PA_NL_LD_notch_controldistance_100KM 23 27877.34
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 7038 -13830.90 27707.80 5396(31509) 4.01e-03 1.1e+04
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.473 1.214
## Number of groups: subject 230
##
## Coefficients:
## Estimate
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.76564
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.58226
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 2.83159
## lockdown -0.13327
## dowMonday -0.21143
## dowSaturday 0.25914
## dowSunday 0.16563
## dowThursday -0.23949
## dowTuesday -0.38525
## dowWednesday -0.21750
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown -0.35639
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.03120
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -0.11457
## Std. Error
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.42622
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.37930
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.46806
## lockdown 0.42087
## dowMonday 0.08082
## dowSaturday 0.08586
## dowSunday 0.08556
## dowThursday 0.08081
## dowTuesday 0.08050
## dowWednesday 0.08298
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.48179
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.42404
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.59836
## z value
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 1.796
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 1.535
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 6.050
## lockdown -0.317
## dowMonday -2.616
## dowSaturday 3.018
## dowSunday 1.936
## dowThursday -2.963
## dowTuesday -4.786
## dowWednesday -2.621
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown -0.740
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.074
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -0.191
## Pr(>|z|)
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.07244
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.12477
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.00000000145
## lockdown 0.75151
## dowMonday 0.00890
## dowSaturday 0.00254
## dowSunday 0.05287
## dowThursday 0.00304
## dowTuesday 0.00000170513
## dowWednesday 0.00877
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.45948
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.94134
## bs(within_subjects_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.84816
##
## 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.4214 0.3968 -8.622
## (10,20]|(20,30] -2.5283 0.3939 -6.418
## (20,30]|(30,40] -1.6636 0.3929 -4.235
## (30,40]|(40,50] -0.7464 0.3925 -1.902
## (40,50]|(50,60] 0.2429 0.3923 0.619
## (50,60]|(60,70] 1.3172 0.3924 3.357
## (60,70]|(70,80] 2.3675 0.3930 6.024
## (70,80]|(80,90] 3.4914 0.3948 8.843
## (80,90]|(90,100] 4.8729 0.4004 12.170
## (20058 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 7031 -13853.22 27752.43 5135(26294) 1.32e-02 1.9e+04
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.455 1.206
## Number of groups: subject 230
##
## Coefficients:
## Estimate
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 -0.36647
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 -0.39673
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 3.39363
## lockdown -0.45623
## dowMonday -0.26589
## dowSaturday 0.20890
## dowSunday 0.12246
## dowThursday -0.26874
## dowTuesday -0.41713
## dowWednesday -0.22940
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.01849
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.23714
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -2.00161
## Std. Error
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.52879
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.48559
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.70114
## lockdown 0.57737
## dowMonday 0.08067
## dowSaturday 0.08586
## dowSunday 0.08559
## dowThursday 0.08071
## dowTuesday 0.08035
## dowWednesday 0.08320
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.63038
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.57515
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.86126
## z value
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 -0.693
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 -0.817
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 4.840
## lockdown -0.790
## dowMonday -3.296
## dowSaturday 2.433
## dowSunday 1.431
## dowThursday -3.329
## dowTuesday -5.192
## dowWednesday -2.757
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.029
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.412
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown -2.324
## Pr(>|z|)
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1 0.488290
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2 0.413923
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3 0.000001297
## lockdown 0.429419
## dowMonday 0.000981
## dowSaturday 0.014976
## dowSunday 0.152495
## dowThursday 0.000870
## dowTuesday 0.000000209
## dowWednesday 0.005831
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)1:lockdown 0.976602
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)2:lockdown 0.680108
## bs(within_subjects_lockdown_z_NL, knots = c(-0.725, -0.3), degree = 1)3:lockdown 0.020123
##
## 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.4856 0.5024 -8.928
## (10,20]|(20,30] -3.5915 0.5000 -7.183
## (20,30]|(30,40] -2.7288 0.4989 -5.470
## (30,40]|(40,50] -1.8157 0.4984 -3.643
## (40,50]|(50,60] -0.8325 0.4980 -1.672
## (50,60]|(60,70] 0.2350 0.4979 0.472
## (60,70]|(70,80] 1.2786 0.4981 2.567
## (70,80]|(80,90] 2.3967 0.4991 4.802
## (80,90]|(90,100] 3.7674 0.5033 7.485
## (20065 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 9141 -18138.40 36302.79 2126(13066) 1.97e-02 4.5e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.167 1.08
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## depression_mean -0.085925 0.016292 -5.274 0.0000001335 ***
## lockdown -0.317932 0.056762 -5.601 0.0000000213 ***
## depression_mean:lockdown -0.026334 0.008252 -3.191 0.00142 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -4.5429 0.1295 -35.074
## (10,20]|(20,30] -3.6333 0.1223 -29.696
## (20,30]|(30,40] -2.7778 0.1189 -23.371
## (30,40]|(40,50] -1.8836 0.1169 -16.119
## (40,50]|(50,60] -0.9163 0.1157 -7.920
## (50,60]|(60,70] 0.1251 0.1154 1.084
## (60,70]|(70,80] 1.1163 0.1163 9.599
## (70,80]|(80,90] 2.1826 0.1195 18.260
## (80,90]|(90,100] 3.5056 0.1313 26.705
## (17955 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 7434 -14855.55 29745.10 644(3217) 2.70e+03 6.0e+07
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1 1
## Number of groups: subject 230
##
## Coefficients:
## Estimate Std. Error z value
## novel_locations 0.0052478 0.0008990 5.837
## depression_mean -0.0780872 0.0166997 -4.676
## lockdown 0.0027775 0.0797833 0.035
## novel_locations:depression_mean 0.0001800 0.0001272 1.415
## novel_locations:lockdown -0.0022914 0.0011493 -1.994
## depression_mean:lockdown -0.0257101 0.0116002 -2.216
## novel_locations:depression_mean:lockdown -0.0001016 0.0001655 -0.614
## Pr(>|z|)
## novel_locations 0.00000000531 ***
## depression_mean 0.00000292574 ***
## lockdown 0.9722
## novel_locations:depression_mean 0.1571
## novel_locations:lockdown 0.0462 *
## depression_mean:lockdown 0.0267 *
## novel_locations:depression_mean:lockdown 0.5391
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.6612 0.1323 -27.670
## (10,20]|(20,30] -2.8611 0.1250 -22.881
## (20,30]|(30,40] -2.0950 0.1215 -17.248
## (30,40]|(40,50] -1.3038 0.1195 -10.909
## (40,50]|(50,60] -0.4763 0.1186 -4.016
## (50,60]|(60,70] 0.3773 0.1187 3.180
## (60,70]|(70,80] 1.1959 0.1197 9.993
## (70,80]|(80,90] 2.1027 0.1225 17.165
## (80,90]|(90,100] 3.2418 0.1321 24.548
## (19662 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 3721 -7336.83 14699.67 1118(8847) 1.12e+01 4.6e+06
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 1.882 1.372
## Number of groups: subject 203
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## novel_locations 0.0046479 0.0005848 7.948 1.9e-15 ***
## anx_change -0.0390810 0.0248931 -1.570 0.116
## novel_locations:anx_change -0.0001814 0.0001503 -1.207 0.227
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## [0,10]|(10,20] -3.66870 0.13285 -27.614
## (10,20]|(20,30] -2.77721 0.11808 -23.520
## (20,30]|(30,40] -1.96050 0.11088 -17.681
## (30,40]|(40,50] -1.02058 0.10692 -9.545
## (40,50]|(50,60] 0.01294 0.10571 0.122
## (50,60]|(60,70] 1.10051 0.10725 10.261
## (60,70]|(70,80] 2.18821 0.11248 19.454
## (70,80]|(80,90] 3.27198 0.12361 26.471
## (80,90]|(90,100] 4.64658 0.15554 29.874
## (14473 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 36395.36
## d2 17 29862.64
## d3 13 14780.55
BIC Analysis
## df BIC
## complex_PA_RE 21 14006.01
## simple_PA_RE 17 13932.74
## complex_PA_NL 21 12933.30
## simple_PA_NL 17 12985.10
## RE_LD 10 53892.41
## NL_LD 10 216192.91
## PA_LD 17 36320.20
## NA_LD 17 35774.45
## PA_RE_LD 19 30479.24
## PA_NL_LD 19 29562.60
## PA_RE_PreLDRE_PostLD 19 16491.63
## PA_RE_PreLDNL_PostLD 19 16538.17
## PA_NL_PreLD_notch 19 13001.04
## PA_NL_PostLD_notch 19 16486.07
## PA_NL_LD_notch 23 29553.83
## PA_NL_LD_notch_controldistance_100KM 23 27877.34
## Z_within_sub 23 27865.56
## Z_within_sub_LD 23 27910.17
## d1 13 36395.36
## d2 17 29862.64
## d3 13 14780.55