Basic Models
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: moving_depression ~ PA_avg + (1 | subject)
## Data: df
##
## REML criterion at convergence: 31954.7
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.8659 -0.4592 -0.0273 0.4247 5.6507
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 23.742 4.873
## Residual 7.343 2.710
## Number of obs: 6386, groups: subject, 251
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 5.519697 0.325739 304.306515 16.95 < 2e-16 ***
## PA_avg -0.005220 0.001977 6198.256809 -2.64 0.00831 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr)
## PA_avg -0.308
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: interpolated_depression ~ PA_avg + (1 | subject)
## Data: df
##
## REML criterion at convergence: 27345.8
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.7911 -0.4215 -0.0186 0.4072 6.0763
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 24.992 4.999
## Residual 3.458 1.859
## Number of obs: 6386, groups: subject, 251
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 5.683644 0.324060 273.991242 17.539 < 2e-16 ***
## PA_avg -0.003803 0.001359 6164.188301 -2.799 0.00515 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr)
## PA_avg -0.213
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: moving_depression ~ novel_locations + (1 | subject)
## Data: df
##
## REML criterion at convergence: 63147.7
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.9939 -0.4512 -0.0297 0.3775 7.1246
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 24.308 4.930
## Residual 7.136 2.671
## Number of obs: 12878, groups: subject, 252
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 5.3019093 0.3122608 252.7151119 16.979 < 2e-16 ***
## novel_locations -0.0014034 0.0004075 12641.1851967 -3.444 0.000576 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr)
## novel_lctns -0.060
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: interpolated_depression ~ novel_locations + (1 | subject)
## Data: df
##
## REML criterion at convergence: 55360.2
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.7750 -0.3802 -0.0126 0.3670 6.6998
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 25.728 5.072
## Residual 3.847 1.961
## Number of obs: 12878, groups: subject, 252
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 5.5042747 0.3204037 251.8607517 17.18 <2e-16 ***
## novel_locations -0.0004100 0.0002993 12633.3426533 -1.37 0.171
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr)
## novel_lctns -0.043
Interaction Models
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: moving_depression ~ PA_avg * lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 25010.5
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.0933 -0.4119 -0.0513 0.3687 6.8894
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 24.178 4.917
## Residual 7.092 2.663
## Number of obs: 4995, groups: subject, 250
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 4.883881 0.342632 350.860445 14.254 < 2e-16 ***
## PA_avg 0.001126 0.002549 4803.553742 0.442 0.6586
## lockdown1 1.322040 0.229832 4766.634377 5.752 0.00000000936 ***
## PA_avg:lockdown1 -0.008587 0.004251 4761.897278 -2.020 0.0434 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) PA_avg lckdw1
## PA_avg -0.394
## lockdown1 -0.224 0.507
## PA_vg:lckd1 0.193 -0.491 -0.920
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: interpolated_depression ~ PA_avg * lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 21602
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.8551 -0.3652 -0.0129 0.3732 6.0937
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 25.24 5.024
## Residual 3.45 1.858
## Number of obs: 4995, groups: subject, 250
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 5.671312 0.333276 295.630089 17.017 <2e-16 ***
## PA_avg -0.004136 0.001781 4772.535541 -2.322 0.0203 *
## lockdown1 -0.084715 0.160421 4754.004540 -0.528 0.5975
## PA_avg:lockdown1 0.004799 0.002967 4751.703310 1.618 0.1058
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) PA_avg lckdw1
## PA_avg -0.283
## lockdown1 -0.161 0.507
## PA_vg:lckd1 0.139 -0.491 -0.920
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: moving_depression ~ PA_avg * who_lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 31748
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.1598 -0.4478 -0.0694 0.4127 5.9774
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 24.109 4.910
## Residual 7.082 2.661
## Number of obs: 6386, groups: subject, 251
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 4.571002 0.343284 363.774280 13.316 < 2e-16 ***
## PA_avg 0.003836 0.002643 6183.723944 1.452 0.146638
## who_lockdown1 1.605597 0.191215 6158.216736 8.397 < 2e-16 ***
## PA_avg:who_lockdown1 -0.011416 0.003447 6154.255466 -3.312 0.000931 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) PA_avg wh_lc1
## PA_avg -0.405
## who_lckdwn1 -0.300 0.657
## PA_vg:wh_l1 0.273 -0.673 -0.931
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: interpolated_depression ~ PA_avg * who_lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 27331.2
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.8672 -0.4126 -0.0241 0.4114 6.1101
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 25.055 5.005
## Residual 3.442 1.855
## Number of obs: 6386, groups: subject, 251
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 5.663809 0.332319 301.389246 17.043 < 2e-16 ***
## PA_avg -0.005567 0.001845 6157.414856 -3.018 0.00256 **
## who_lockdown1 -0.008724 0.133381 6144.722236 -0.065 0.94785
## PA_avg:who_lockdown1 0.004869 0.002404 6142.751175 2.026 0.04284 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) PA_avg wh_lc1
## PA_avg -0.292
## who_lckdwn1 -0.217 0.656
## PA_vg:wh_l1 0.197 -0.673 -0.931
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: moving_depression ~ novel_locations * lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 57755.3
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -5.1270 -0.3950 -0.0537 0.3830 7.2374
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 24.21 4.921
## Residual 7.08 2.661
## Number of obs: 11779, groups: subject, 252
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 4.9736511 0.3134909 257.9542314 15.865
## novel_locations 0.0004269 0.0005449 11542.0868789 0.784
## lockdown1 0.6732695 0.0663731 11553.3791763 10.144
## novel_locations:lockdown1 -0.0028375 0.0009345 11537.4069927 -3.036
## Pr(>|t|)
## (Intercept) <2e-16 ***
## novel_locations 0.4333
## lockdown1 <2e-16 ***
## novel_locations:lockdown1 0.0024 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) nvl_lc lckdw1
## novel_lctns -0.096
## lockdown1 -0.093 0.461
## nvl_lctns:1 0.055 -0.570 -0.554
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: interpolated_depression ~ novel_locations * lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 50652.6
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.8112 -0.3657 -0.0142 0.3474 6.6863
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 25.559 5.056
## Residual 3.819 1.954
## Number of obs: 11779, groups: subject, 252
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.5090914 0.3203398 254.6421937 17.198
## novel_locations 0.0001592 0.0004003 11533.6077780 0.398
## lockdown1 0.0226232 0.0487594 11539.4010493 0.464
## novel_locations:lockdown1 -0.0022632 0.0006864 11531.0300874 -3.297
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## novel_locations 0.690884
## lockdown1 0.642674
## novel_locations:lockdown1 0.000979 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) nvl_lc lckdw1
## novel_lctns -0.069
## lockdown1 -0.067 0.461
## nvl_lctns:1 0.040 -0.570 -0.554
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: moving_depression ~ novel_locations * who_lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 62876
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -5.2059 -0.4220 -0.0740 0.4156 7.3300
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 24.536 4.953
## Residual 6.975 2.641
## Number of obs: 12878, groups: subject, 252
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 4.7394831 0.3164944 261.8836403 14.975
## novel_locations 0.0009501 0.0006819 12643.0689446 1.393
## who_lockdown1 0.9492635 0.0657723 12649.8861152 14.433
## novel_locations:who_lockdown1 -0.0020609 0.0008397 12637.6027873 -2.454
## Pr(>|t|)
## (Intercept) <2e-16 ***
## novel_locations 0.1635
## who_lockdown1 <2e-16 ***
## novel_locations:who_lockdown1 0.0141 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) nvl_lc wh_lc1
## novel_lctns -0.121
## who_lckdwn1 -0.133 0.588
## nvl_lctn:_1 0.098 -0.802 -0.650
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: interpolated_depression ~ novel_locations * who_lockdown + (1 |
## subject)
## Data: df
##
## REML criterion at convergence: 55373.8
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.7952 -0.3802 -0.0174 0.3651 6.6902
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 25.748 5.074
## Residual 3.846 1.961
## Number of obs: 12878, groups: subject, 252
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.4422234 0.3220608 256.6980999 16.898
## novel_locations 0.0001441 0.0005064 12633.7277858 0.285
## who_lockdown1 0.0981531 0.0488546 12637.4062803 2.009
## novel_locations:who_lockdown1 -0.0007396 0.0006236 12630.7748211 -1.186
## Pr(>|t|)
## (Intercept) <2e-16 ***
## novel_locations 0.7760
## who_lockdown1 0.0445 *
## novel_locations:who_lockdown1 0.2357
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) nvl_lc wh_lc1
## novel_lctns -0.088
## who_lckdwn1 -0.097 0.588
## nvl_lctn:_1 0.071 -0.802 -0.650
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: interpolated_depression ~ novel_locations * lockdown + (1 | subject)
## Data: df[which(df$novel_locations < 200), ]
##
## REML criterion at convergence: 49595.9
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.8494 -0.3700 -0.0174 0.3405 6.8212
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 25.571 5.057
## Residual 3.815 1.953
## Number of obs: 11530, groups: subject, 252
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.5437122 0.3210675 256.6605243 17.266
## novel_locations -0.0006192 0.0006086 11291.1277397 -1.017
## lockdown1 0.0559301 0.0545678 11291.4780056 1.025
## novel_locations:lockdown1 -0.0042722 0.0010087 11285.4588236 -4.236
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## novel_locations 0.309
## lockdown1 0.305
## novel_locations:lockdown1 0.000023 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) nvl_lc lckdw1
## novel_lctns -0.093
## lockdown1 -0.083 0.542
## nvl_lctns:1 0.054 -0.580 -0.648
Double Interaction Models
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: moving_depression ~ novel_locations * lockdown + PA_avg * lockdown +
## (1 | subject)
## Data: df
##
## REML criterion at convergence: 24017.5
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.0634 -0.4074 -0.0485 0.4017 6.8045
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 24.329 4.932
## Residual 7.249 2.692
## Number of obs: 4767, groups: subject, 248
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 4.93069784 0.34885461 362.48332415 14.134
## novel_locations -0.00005831 0.00083158 4534.27380061 -0.070
## lockdown1 1.30545816 0.24115856 4539.53377848 5.413
## PA_avg 0.00066153 0.00268271 4575.50184822 0.247
## novel_locations:lockdown1 -0.00323547 0.00241962 4531.82458813 -1.337
## lockdown1:PA_avg -0.00756361 0.00435684 4534.75725840 -1.736
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## novel_locations 0.9441
## lockdown1 0.0000000651 ***
## PA_avg 0.8052
## novel_locations:lockdown1 0.1812
## lockdown1:PA_avg 0.0826 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) nvl_lc lckdw1 PA_avg nvl_:1
## novel_lctns -0.084
## lockdown1 -0.243 0.131
## PA_avg -0.394 -0.109 0.493
## nvl_lctns:1 0.032 -0.332 -0.200 0.028
## lckdwn1:PA_ 0.199 0.056 -0.884 -0.506 -0.050
## [[1]]
##
## [[2]]
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: interpolated_depression ~ novel_locations * lockdown + PA_avg *
## lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 20601.8
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.8237 -0.3691 -0.0158 0.3766 6.1060
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 25.747 5.074
## Residual 3.393 1.842
## Number of obs: 4767, groups: subject, 248
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.7771997 0.3390016 296.6595734 17.042
## novel_locations -0.0002556 0.0005693 4523.3469282 -0.449
## lockdown1 -0.1443090 0.1651219 4525.7214264 -0.874
## PA_avg -0.0049054 0.0018391 4542.9735920 -2.667
## novel_locations:lockdown1 -0.0036818 0.0016563 4522.0231496 -2.223
## lockdown1:PA_avg 0.0063331 0.0029827 4523.4911305 2.123
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## novel_locations 0.65350
## lockdown1 0.38219
## PA_avg 0.00767 **
## novel_locations:lockdown1 0.02627 *
## lockdown1:PA_avg 0.03378 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) nvl_lc lckdw1 PA_avg nvl_:1
## novel_lctns -0.059
## lockdown1 -0.171 0.131
## PA_avg -0.278 -0.108 0.493
## nvl_lctns:1 0.022 -0.332 -0.200 0.028
## lckdwn1:PA_ 0.140 0.056 -0.884 -0.506 -0.051
## [[1]]
##
## [[2]]
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: moving_depression ~ novel_locations * who_lockdown + PA_avg *
## who_lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 25756.2
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.0968 -0.4225 -0.0481 0.4197 7.0443
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 24.508 4.951
## Residual 6.989 2.644
## Number of obs: 5158, groups: subject, 249
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 4.4888946 0.3534391 383.0710088 12.701
## novel_locations 0.0001941 0.0011059 4930.7280405 0.175
## who_lockdown1 1.9443943 0.2161076 4928.9422548 8.997
## PA_avg 0.0051515 0.0028916 4956.7585544 1.782
## novel_locations:who_lockdown1 -0.0010226 0.0014016 4921.7086727 -0.730
## who_lockdown1:PA_avg -0.0151254 0.0038319 4925.5967148 -3.947
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## novel_locations 0.8607
## who_lockdown1 < 2e-16 ***
## PA_avg 0.0749 .
## novel_locations:who_lockdown1 0.4657
## who_lockdown1:PA_avg 0.0000802 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) nvl_lc wh_lc1 PA_avg nv_:_1
## novel_lctns -0.090
## who_lckdwn1 -0.314 0.156
## PA_avg -0.402 -0.178 0.587
## nvl_lctn:_1 0.074 -0.778 -0.177 0.132
## wh_lck1:PA_ 0.263 0.125 -0.881 -0.657 -0.138
## [[1]]
##
## [[2]]
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: interpolated_depression ~ novel_locations * who_lockdown + PA_avg *
## who_lockdown + (1 | subject)
## Data: df
##
## REML criterion at convergence: 22220.3
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.9448 -0.3835 -0.0273 0.3905 6.1598
##
## Random effects:
## Groups Name Variance Std.Dev.
## subject (Intercept) 25.739 5.073
## Residual 3.392 1.842
## Number of obs: 5158, groups: subject, 249
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.7564345 0.3409770 307.9754663 16.882
## novel_locations -0.0001059 0.0007710 4916.6748695 -0.137
## who_lockdown1 -0.0342530 0.1506519 4915.8841119 -0.227
## PA_avg -0.0057388 0.0020175 4929.6670357 -2.845
## novel_locations:who_lockdown1 -0.0002739 0.0009769 4912.3205951 -0.280
## who_lockdown1:PA_avg 0.0052888 0.0026710 4914.2635642 1.980
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## novel_locations 0.89071
## who_lockdown1 0.82015
## PA_avg 0.00447 **
## novel_locations:who_lockdown1 0.77917
## who_lockdown1:PA_avg 0.04775 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) nvl_lc wh_lc1 PA_avg nv_:_1
## novel_lctns -0.065
## who_lckdwn1 -0.227 0.156
## PA_avg -0.291 -0.178 0.587
## nvl_lctn:_1 0.053 -0.778 -0.177 0.132
## wh_lck1:PA_ 0.190 0.125 -0.881 -0.657 -0.138
## [[1]]
##
## [[2]]
