Multilevel Analysis
Carson Dover
2026-04-20
32-day Graphs
All Time Points Multilevel Models
Feeling Down
Intercept Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID) + (1 | FamID)
## Data: core1_data
##
## REML criterion at convergence: 2105
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.5978 -0.4830 -0.0851 0.1528 4.0634
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 2.432e-01 0.4931956
## FamID (Intercept) 6.254e-08 0.0002501
## Residual 3.486e-01 0.5904398
## Number of obs: 1106, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 1.47415 0.07574 43.89434 19.46 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Intercept Model (No Family Nesting)
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID)
## Data: core1_data
##
## REML criterion at convergence: 2105
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.5978 -0.4830 -0.0851 0.1528 4.0634
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.2432 0.4932
## Residual 0.3486 0.5904
## Number of obs: 1106, groups: ID, 45
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 1.47415 0.07574 43.92405 19.46 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1ICC Estimates
## [1] "ICC for ID is 0.411"## [1] "ICC for residuals is 0.589"Effects Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ Role * Day + Role * I(Day^2) + (1 + Day | ID)
## Data: core1_data
##
## REML criterion at convergence: 2124.1
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.3312 -0.4460 -0.0907 0.1663 4.0095
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.1911649 0.43722
## Day 0.0002019 0.01421 -0.09
## Residual 0.3256072 0.57062
## Number of obs: 1106, groups: ID, 45
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 1.522e+00 1.528e-01 6.472e+01 9.965 1.07e-14 ***
## RoleMother -1.048e-01 2.021e-01 6.331e+01 -0.518 0.6060
## RoleFather -2.808e-01 2.081e-01 6.389e+01 -1.349 0.1820
## Day 5.353e-03 1.477e-02 8.286e+02 0.362 0.7172
## I(Day^2) 2.108e-04 4.506e-04 1.059e+03 0.468 0.6399
## RoleMother:Day 2.874e-02 1.923e-02 8.033e+02 1.495 0.1354
## RoleFather:Day -2.619e-03 1.990e-02 8.152e+02 -0.132 0.8954
## RoleMother:I(Day^2) -1.397e-03 5.822e-04 1.057e+03 -2.399 0.0166 *
## RoleFather:I(Day^2) -4.220e-04 6.057e-04 1.058e+03 -0.697 0.4862
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) RlMthr RlFthr Day I(D^2) RlMt:D RlFt:D RM:I(D
## RoleMother -0.756
## RoleFather -0.734 0.555
## Day -0.519 0.392 0.381
## I(Day^2) 0.441 -0.333 -0.324 -0.930
## RoleMthr:Dy 0.399 -0.514 -0.293 -0.768 0.714
## RoleFthr:Dy 0.385 -0.291 -0.516 -0.742 0.690 0.570
## RlMt:I(D^2) -0.341 0.436 0.250 0.719 -0.774 -0.929 -0.534
## RlFt:I(D^2) -0.328 0.248 0.437 0.691 -0.744 -0.531 -0.928 0.576
## optimizer (nloptwrap) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.00871004 (tol = 0.002, component 1)
Feeling Lonely
Intercept Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID) + (1 | FamID)
## Data: core2_data
##
## REML criterion at convergence: 2181.5
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.5523 -0.4624 -0.1717 0.1154 4.2105
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.18056 0.4249
## FamID (Intercept) 0.03997 0.1999
## Residual 0.37651 0.6136
## Number of obs: 1106, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 1.46867 0.08222 14.11493 17.86 4.34e-11 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1ICC Estimates
## [1] "ICC for ID is 0.3024"## [1] "ICC for famID is 0.0669"## [1] "ICC for residuals is 0.6306"Effects Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ Role * Day + Role * I(Day^2) + (1 + Day | ID) + (1 |
## FamID)
## Data: core2_data
##
## REML criterion at convergence: 2192.3
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.7197 -0.5145 -0.1804 0.1840 4.3539
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.1298918 0.36041
## Day 0.0001453 0.01205 0.24
## FamID (Intercept) 0.0157877 0.12565
## Residual 0.3518598 0.59318
## Number of obs: 1106, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 1.287e+00 1.430e-01 6.504e+01 9.000 4.93e-13 ***
## RoleMother -1.501e-01 1.840e-01 4.193e+01 -0.816 0.4191
## RoleFather -5.409e-02 1.896e-01 4.269e+01 -0.285 0.7769
## Day 3.836e-02 1.513e-02 8.876e+02 2.536 0.0114 *
## I(Day^2) -7.964e-04 4.658e-04 1.057e+03 -1.710 0.0876 .
## RoleMother:Day 1.994e-02 1.967e-02 8.668e+02 1.014 0.3110
## RoleFather:Day -2.915e-02 2.037e-02 8.762e+02 -1.431 0.1527
## RoleMother:I(Day^2) -8.280e-04 6.023e-04 1.058e+03 -1.375 0.1695
## RoleFather:I(Day^2) 5.254e-04 6.263e-04 1.058e+03 0.839 0.4017
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) RlMthr RlFthr Day I(D^2) RlMt:D RlFt:D RM:I(D
## RoleMother -0.742
## RoleFather -0.719 0.557
## Day -0.523 0.406 0.394
## I(Day^2) 0.489 -0.380 -0.369 -0.941
## RoleMthr:Dy 0.402 -0.530 -0.303 -0.769 0.724
## RoleFthr:Dy 0.388 -0.302 -0.532 -0.743 0.699 0.571
## RlMt:I(D^2) -0.379 0.497 0.285 0.728 -0.773 -0.940 -0.540
## RlFt:I(D^2) -0.364 0.283 0.499 0.700 -0.744 -0.538 -0.940 0.575
## optimizer (nloptwrap) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.00302669 (tol = 0.002, component 1)
Physical Discomfort
Intercept Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID) + (1 | FamID)
## Data: core3_data
##
## REML criterion at convergence: 2386.3
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.3890 -0.4549 -0.1594 0.2975 4.1816
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.1249 0.3534
## FamID (Intercept) 0.1708 0.4132
## Residual 0.4566 0.6757
## Number of obs: 1105, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 1.6519 0.1143 16.4385 14.45 9.1e-11 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1ICC Estimates
## [1] "ICC for ID is 0.166"## [1] "ICC for famID is 0.227"## [1] "ICC for residuals is 0.607"Effects Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ Role * Day + Role * I(Day^2) + (1 + Day | ID) + (1 |
## FamID)
## Data: core3_data
##
## REML criterion at convergence: 2367.2
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.8014 -0.4710 -0.1721 0.2256 4.4382
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.0711992 0.26683
## Day 0.0005927 0.02434 -0.50
## FamID (Intercept) 0.1934300 0.43981
## Residual 0.4051311 0.63650
## Number of obs: 1105, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 1.776e+00 1.680e-01 5.707e+01 10.567 4.72e-15 ***
## RoleMother -1.453e-01 1.713e-01 5.706e+01 -0.848 0.400
## RoleFather -9.947e-02 1.763e-01 5.870e+01 -0.564 0.575
## Day 3.290e-04 1.738e-02 5.997e+02 0.019 0.985
## I(Day^2) 2.440e-04 5.076e-04 1.058e+03 0.481 0.631
## RoleMother:Day -4.777e-03 2.259e-02 5.712e+02 -0.211 0.833
## RoleFather:Day -1.584e-02 2.337e-02 5.801e+02 -0.678 0.498
## RoleMother:I(Day^2) -1.620e-06 6.534e-04 1.054e+03 -0.002 0.998
## RoleFather:I(Day^2) 8.607e-05 6.805e-04 1.058e+03 0.126 0.899
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) RlMthr RlFthr Day I(D^2) RlMt:D RlFt:D RM:I(D
## RoleMother -0.601
## RoleFather -0.578 0.558
## Day -0.578 0.566 0.551
## I(Day^2) 0.453 -0.443 -0.431 -0.887
## RoleMthr:Dy 0.447 -0.743 -0.423 -0.768 0.682
## RoleFthr:Dy 0.431 -0.421 -0.745 -0.743 0.659 0.571
## RlMt:I(D^2) -0.354 0.581 0.335 0.688 -0.776 -0.885 -0.511
## RlFt:I(D^2) -0.339 0.330 0.582 0.661 -0.745 -0.508 -0.884 0.578
Feeling Anxious
Intercept Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID) + (1 | FamID)
## Data: core7_data
##
## REML criterion at convergence: 2405.9
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.6365 -0.4986 -0.0609 0.5002 3.9398
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.2230 0.4722
## FamID (Intercept) 0.2084 0.4565
## Residual 0.4568 0.6759
## Number of obs: 1105, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 1.7705 0.1323 18.0404 13.38 8.27e-11 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1ICC Estimates
## [1] "ICC for ID is 0.251"## [1] "ICC for famID is 0.2346"## [1] "ICC for residuals is 0.5143"Effects Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ Role * Day + Role * I(Day^2) + (1 + Day | ID) + (1 |
## FamID)
## Data: core7_data
##
## REML criterion at convergence: 2350.3
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.3416 -0.5219 -0.0741 0.4317 4.3994
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.4668469 0.68326
## Day 0.0005766 0.02401 -0.79
## FamID (Intercept) 0.2551135 0.50509
## Residual 0.3877941 0.62273
## Number of obs: 1105, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 1.833e+00 2.503e-01 5.174e+01 7.320 1.57e-09 ***
## RoleMother 5.061e-01 2.887e-01 4.579e+01 1.753 0.08627 .
## RoleFather -1.923e-01 2.967e-01 4.574e+01 -0.648 0.52010
## Day -2.150e-02 1.708e-02 5.856e+02 -1.259 0.20850
## I(Day^2) 8.036e-04 4.973e-04 1.051e+03 1.616 0.10639
## RoleMother:Day -1.322e-02 2.223e-02 5.544e+02 -0.595 0.55224
## RoleFather:Day 4.124e-02 2.297e-02 5.661e+02 1.795 0.07316 .
## RoleMother:I(Day^2) -4.690e-04 6.411e-04 1.048e+03 -0.732 0.46454
## RoleFather:I(Day^2) -1.881e-03 6.668e-04 1.050e+03 -2.821 0.00488 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) RlMthr RlFthr Day I(D^2) RlMt:D RlFt:D RM:I(D
## RoleMother -0.666
## RoleFather -0.644 0.553
## Day -0.556 0.482 0.468
## I(Day^2) 0.298 -0.259 -0.250 -0.888
## RoleMthr:Dy 0.428 -0.641 -0.360 -0.768 0.682
## RoleFthr:Dy 0.414 -0.359 -0.638 -0.743 0.660 0.571
## RlMt:I(D^2) -0.232 0.340 0.195 0.688 -0.775 -0.885 -0.512
## RlFt:I(D^2) -0.222 0.194 0.337 0.662 -0.745 -0.508 -0.885 0.578
## optimizer (nloptwrap) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.0248291 (tol = 0.002, component 1)
Feeling Thankful
Intercept Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID) + (1 | FamID)
## Data: core8_data
##
## REML criterion at convergence: 2553.9
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.2249 -0.5813 -0.0882 0.7063 2.7713
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.26855 0.5182
## FamID (Intercept) 0.07861 0.2804
## Residual 0.52534 0.7248
## Number of obs: 1106, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.7931 0.1056 14.4930 26.46 1.12e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1ICC Estimates
## [1] "ICC for ID is 0.3078"## [1] "ICC for famID is 0.0901"## [1] "ICC for residuals is 0.6021"Effects Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ Role * Day + Role * I(Day^2) + (1 + Day | ID) + (1 |
## FamID)
## Data: core8_data
##
## REML criterion at convergence: 2551.2
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.5521 -0.5848 -0.0434 0.6826 3.1250
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.336410 0.58001
## Day 0.000445 0.02109 -0.49
## FamID (Intercept) 0.086960 0.29489
## Residual 0.474474 0.68882
## Number of obs: 1106, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.975e+00 2.119e-01 5.959e+01 14.040 <2e-16 ***
## RoleMother 1.967e-01 2.626e-01 3.949e+01 0.749 0.458
## RoleFather 3.841e-02 2.701e-01 3.961e+01 0.142 0.888
## Day -2.182e-02 1.829e-02 7.336e+02 -1.193 0.233
## I(Day^2) 6.235e-04 5.479e-04 1.055e+03 1.138 0.255
## RoleMother:Day -9.442e-03 2.378e-02 7.023e+02 -0.397 0.691
## RoleFather:Day -2.075e-02 2.459e-02 7.154e+02 -0.844 0.399
## RoleMother:I(Day^2) 4.297e-05 7.068e-04 1.053e+03 0.061 0.952
## RoleFather:I(Day^2) 1.457e-04 7.355e-04 1.055e+03 0.198 0.843
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) RlMthr RlFthr Day I(D^2) RlMt:D RlFt:D RM:I(D
## RoleMother -0.716
## RoleFather -0.693 0.556
## Day -0.553 0.446 0.433
## I(Day^2) 0.389 -0.314 -0.305 -0.913
## RoleMthr:Dy 0.426 -0.588 -0.333 -0.769 0.702
## RoleFthr:Dy 0.411 -0.332 -0.587 -0.743 0.679 0.572
## RlMt:I(D^2) -0.303 0.411 0.237 0.708 -0.775 -0.912 -0.526
## RlFt:I(D^2) -0.290 0.234 0.411 0.680 -0.745 -0.523 -0.911 0.577
## optimizer (nloptwrap) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.00742022 (tol = 0.002, component 1)
Thinking About Freedom
Intercept Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID) + (1 | FamID)
## Data: core9_data
##
## REML criterion at convergence: 2488.2
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.7398 -0.6976 0.0454 0.6432 3.5700
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.63295 0.7956
## FamID (Intercept) 0.04406 0.2099
## Residual 0.48119 0.6937
## Number of obs: 1105, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.3693 0.1311 15.9300 18.07 4.89e-12 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1ICC Estimates
## [1] "ICC for ID is 0.5465"## [1] "ICC for famID is 0.038"## [1] "ICC for residuals is 0.4155"Effects Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ Role * Day + Role * I(Day^2) + (1 + Day | ID) + (1 |
## FamID)
## Data: core9_data
##
## REML criterion at convergence: 2407.5
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.4480 -0.6517 0.0480 0.5548 3.9005
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.3363644 0.57997
## Day 0.0005468 0.02338 -0.66
## FamID (Intercept) 0.1995140 0.44667
## Residual 0.4115694 0.64154
## Number of obs: 1105, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.672e+00 2.237e-01 5.523e+01 11.946 < 2e-16 ***
## RoleMother -9.494e-01 2.581e-01 4.755e+01 -3.679 0.000596 ***
## RoleFather -9.073e-01 2.652e-01 4.759e+01 -3.421 0.001291 **
## Day 9.303e-02 1.744e-02 5.765e+02 5.333 1.39e-07 ***
## I(Day^2) -2.856e-03 5.125e-04 1.050e+03 -5.573 3.18e-08 ***
## RoleMother:Day -6.122e-02 2.270e-02 5.447e+02 -2.697 0.007212 **
## RoleFather:Day -7.780e-02 2.346e-02 5.566e+02 -3.316 0.000972 ***
## RoleMother:I(Day^2) 2.580e-03 6.605e-04 1.047e+03 3.906 9.98e-05 ***
## RoleFather:I(Day^2) 2.413e-03 6.873e-04 1.050e+03 3.510 0.000467 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) RlMthr RlFthr Day I(D^2) RlMt:D RlFt:D RM:I(D
## RoleMother -0.669
## RoleFather -0.646 0.554
## Day -0.552 0.479 0.465
## I(Day^2) 0.344 -0.298 -0.289 -0.895
## RoleMthr:Dy 0.426 -0.634 -0.358 -0.768 0.687
## RoleFthr:Dy 0.411 -0.356 -0.632 -0.744 0.665 0.571
## RlMt:I(D^2) -0.268 0.391 0.225 0.694 -0.776 -0.893 -0.516
## RlFt:I(D^2) -0.257 0.223 0.389 0.667 -0.745 -0.512 -0.892 0.578
## optimizer (nloptwrap) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.0218905 (tol = 0.002, component 1)
Meaningful Contact
Intercept Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID) + (1 | FamID)
## Data: core11_data
##
## REML criterion at convergence: 2866.5
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.4691 -0.6586 -0.1229 0.7180 2.4953
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.11128 0.3336
## FamID (Intercept) 0.05435 0.2331
## Residual 0.72665 0.8524
## Number of obs: 1105, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.31483 0.07993 16.83735 28.96 8.4e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1ICC Estimates
## [1] "ICC for ID is 0.1247"## [1] "ICC for famID is 0.0609"## [1] "ICC for residuals is 0.8144"Effects Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ Role * Day + Role * I(Day^2) + (1 + Day | ID) + (1 |
## FamID)
## Data: core11_data
##
## REML criterion at convergence: 2811.3
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.5861 -0.6544 -0.1205 0.6378 2.9070
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.340475 0.58350
## Day 0.000748 0.02735 -0.86
## FamID (Intercept) 0.060814 0.24660
## Residual 0.623036 0.78933
## Number of obs: 1105, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.645e+00 2.180e-01 6.156e+01 12.131 <2e-16 ***
## RoleMother 2.588e-01 2.763e-01 5.521e+01 0.937 0.3531
## RoleFather 2.327e-01 2.840e-01 5.556e+01 0.819 0.4160
## Day -3.805e-02 2.129e-02 6.425e+02 -1.787 0.0743 .
## I(Day^2) 8.742e-04 6.275e-04 1.057e+03 1.393 0.1639
## RoleMother:Day -1.408e-02 2.769e-02 6.107e+02 -0.508 0.6113
## RoleFather:Day -3.916e-02 2.863e-02 6.225e+02 -1.368 0.1719
## RoleMother:I(Day^2) 5.332e-05 8.094e-04 1.055e+03 0.066 0.9475
## RoleFather:I(Day^2) 6.838e-04 8.419e-04 1.055e+03 0.812 0.4169
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) RlMthr RlFthr Day I(D^2) RlMt:D RlFt:D RM:I(D
## RoleMother -0.731
## RoleFather -0.709 0.557
## Day -0.708 0.559 0.543
## I(Day^2) 0.436 -0.344 -0.334 -0.901
## RoleMthr:Dy 0.545 -0.738 -0.418 -0.769 0.693
## RoleFthr:Dy 0.527 -0.416 -0.736 -0.744 0.670 0.572
## RlMt:I(D^2) -0.339 0.450 0.259 0.698 -0.775 -0.899 -0.519
## RlFt:I(D^2) -0.325 0.257 0.449 0.671 -0.745 -0.516 -0.899 0.577
Being Kind to Self
Intercept Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID) + (1 | FamID)
## Data: core17_data
##
## REML criterion at convergence: 2498.8
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -2.8549 -0.7657 -0.0290 0.4691 3.4109
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.258200 0.50813
## FamID (Intercept) 0.001423 0.03772
## Residual 0.505598 0.71105
## Number of obs: 1104, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.27514 0.07941 15.00841 28.65 1.62e-14 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1ICC Estimates
## [1] "ICC for ID is 0.3374"## [1] "ICC for famID is 0.0019"## [1] "ICC for residuals is 0.6607"Effects Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ Role * Day + Role * I(Day^2) + (1 + Day | ID) + (1 |
## FamID)
## Data: core17_data
##
## REML criterion at convergence: 2504.3
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.1599 -0.6223 -0.0339 0.5259 3.5472
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.2381623 0.48802
## Day 0.0004065 0.02016 -0.44
## FamID (Intercept) 0.0363623 0.19069
## Residual 0.4626990 0.68022
## Number of obs: 1104, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.651e+00 1.832e-01 6.380e+01 14.470 <2e-16 ***
## RoleMother -9.062e-02 2.334e-01 4.812e+01 -0.388 0.6995
## RoleFather -4.383e-01 2.402e-01 4.851e+01 -1.825 0.0741 .
## Day -1.773e-02 1.804e-02 7.432e+02 -0.983 0.3260
## I(Day^2) 1.030e-04 5.424e-04 1.055e+03 0.190 0.8494
## RoleMother:Day 6.214e-03 2.343e-02 7.116e+02 0.265 0.7909
## RoleFather:Day -8.606e-03 2.423e-02 7.239e+02 -0.355 0.7226
## RoleMother:I(Day^2) -2.788e-05 6.987e-04 1.053e+03 -0.040 0.9682
## RoleFather:I(Day^2) 4.022e-04 7.270e-04 1.054e+03 0.553 0.5802
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) RlMthr RlFthr Day I(D^2) RlMt:D RlFt:D RM:I(D
## RoleMother -0.735
## RoleFather -0.713 0.557
## Day -0.598 0.470 0.456
## I(Day^2) 0.445 -0.349 -0.339 -0.917
## RoleMthr:Dy 0.461 -0.617 -0.352 -0.770 0.706
## RoleFthr:Dy 0.446 -0.350 -0.617 -0.745 0.683 0.573
## RlMt:I(D^2) -0.346 0.456 0.263 0.712 -0.776 -0.915 -0.530
## RlFt:I(D^2) -0.332 0.261 0.456 0.684 -0.746 -0.527 -0.915 0.579
## optimizer (nloptwrap) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.00265103 (tol = 0.002, component 1)
Taking Care of Self
Intercept Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ 1 + (1 | ID) + (1 | FamID)
## Data: core18_data
##
## REML criterion at convergence: 2716
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.2379 -0.6074 -0.0543 0.6333 2.9563
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.2440 0.4940
## FamID (Intercept) 0.1288 0.3588
## Residual 0.6184 0.7864
## Number of obs: 1101, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.3757 0.1166 17.4412 20.37 1.33e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1ICC Estimates
## [1] "ICC for ID is 0.2462"## [1] "ICC for famID is 0.1299"## [1] "ICC for residuals is 0.6239"Effects Model
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: value ~ Role * Day + Role * I(Day^2) + (1 + Day | ID) + (1 |
## FamID)
## Data: core18_data
##
## REML criterion at convergence: 2763.7
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.2625 -0.5995 -0.0613 0.6279 3.0173
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.2469683 0.49696
## Day 0.0002661 0.01631 -0.27
## FamID (Intercept) 0.1333922 0.36523
## Residual 0.5957962 0.77188
## Number of obs: 1101, groups: ID, 45; FamID, 18
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.822e+00 2.118e-01 6.837e+01 13.324 <2e-16 ***
## RoleMother -4.676e-01 2.522e-01 4.870e+01 -1.854 0.0698 .
## RoleFather -4.631e-01 2.590e-01 4.852e+01 -1.788 0.0800 .
## Day -4.867e-02 2.009e-02 8.953e+02 -2.423 0.0156 *
## I(Day^2) 1.344e-03 6.149e-04 1.053e+03 2.185 0.0291 *
## RoleMother:Day 3.490e-02 2.610e-02 8.669e+02 1.337 0.1815
## RoleFather:Day 4.816e-02 2.691e-02 8.748e+02 1.790 0.0739 .
## RoleMother:I(Day^2) -8.479e-04 7.943e-04 1.051e+03 -1.067 0.2860
## RoleFather:I(Day^2) -1.377e-03 8.234e-04 1.052e+03 -1.672 0.0947 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) RlMthr RlFthr Day I(D^2) RlMt:D RlFt:D RM:I(D
## RoleMother -0.696
## RoleFather -0.672 0.559
## Day -0.551 0.462 0.450
## I(Day^2) 0.448 -0.376 -0.366 -0.941
## RoleMthr:Dy 0.425 -0.605 -0.347 -0.770 0.724
## RoleFthr:Dy 0.411 -0.345 -0.604 -0.746 0.702 0.575
## RlMt:I(D^2) -0.348 0.491 0.284 0.728 -0.774 -0.940 -0.543
## RlFt:I(D^2) -0.335 0.281 0.489 0.702 -0.746 -0.540 -0.939 0.578
## optimizer (nloptwrap) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.00672121 (tol = 0.002, component 1)