Introduction
In all variables, I did the following: For each muscle, I did a two-way “analysis of variance”, one direction, the side of the operation (OP.side) which is a between-subjects factor with two levels, “OP”, “Nonop”, and the second direction, also two levels within subject factor, when the measurment was taken (Type), “before” or “after” to operation.
digastric.coronal.plain.long.axis
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 88.456 88.456 1 115.57 39.865 5.221e-09 ***
Time 106.888 106.888 1 115.57 48.172 2.416e-10 ***
OP.side:Time 72.121 72.121 1 115.57 32.503 9.270e-08 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: digastric.coronal.plain.long.axis.in.mm ~ OP.side * Time + (1 |
ID)
Data: temp
REML criterion at convergence: 645.7
Scaled residuals:
Min 1Q Median 3Q Max
-1.8842 -0.5269 -0.1670 0.5498 3.5231
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 4.007 2.002
Residual 2.219 1.490
Number of obs: 156, groups: ID, 38
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 12.8400 0.4031 66.1294 31.851 < 2e-16 ***
OP.sideOP -2.8659 0.3373 115.5733 -8.496 7.83e-14 ***
TimeBefore 0.2956 0.3373 115.5733 0.876 0.383
OP.sideOP:TimeBefore 2.7197 0.4771 115.5733 5.701 9.27e-08 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 12.84 0.403 65.4 12.03 13.6
OP 9.97 0.403 65.4 9.17 10.8
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 13.14 0.403 65.4 12.33 13.9
OP 12.99 0.403 65.4 12.18 13.8
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP 2.866 0.337 115 8.496 <.0001
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP 0.146 0.337 115 0.433 0.6656
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before -0.296 0.337 115 -0.876 0.3826
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before -3.015 0.337 115 -8.939 <.0001
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop -0.296 0.337 115 -0.876 0.8171
After Nonop - After OP 2.866 0.337 115 8.496 <.0001
After Nonop - Before OP -0.149 0.337 115 -0.443 0.9708
Before Nonop - After OP 3.162 0.337 115 9.372 <.0001
Before Nonop - Before OP 0.146 0.337 115 0.433 0.9726
After OP - Before OP -3.015 0.337 115 -8.939 <.0001
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
digastric.coronal.plain.short.axis
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 66.535 66.535 1 117.21 51.509 6.999e-11 ***
Time 108.533 108.533 1 117.21 84.021 1.971e-15 ***
OP.side:Time 59.422 59.422 1 117.21 46.002 5.082e-10 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: digastric.coronal.plain.short.axis.in.mm ~ OP.side * Time + (1 |
ID)
Data: temp
REML criterion at convergence: 531.7
Scaled residuals:
Min 1Q Median 3Q Max
-2.8883 -0.5451 -0.1279 0.5299 2.3255
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 0.8049 0.8972
Residual 1.2917 1.1365
Number of obs: 156, groups: ID, 38
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 6.8952 0.2332 104.7107 29.565 < 2e-16 ***
OP.sideOP -2.5405 0.2574 117.2113 -9.871 < 2e-16 ***
TimeBefore 0.4338 0.2574 117.2113 1.686 0.0945 .
OP.sideOP:TimeBefore 2.4687 0.3640 117.2113 6.782 5.08e-10 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 6.90 0.233 102 6.43 7.36
OP 4.35 0.233 102 3.89 4.82
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 7.33 0.233 102 6.87 7.79
OP 7.26 0.233 102 6.79 7.72
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP 2.5405 0.257 115 9.871 <.0001
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP 0.0718 0.257 115 0.279 0.7808
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before -0.434 0.257 115 -1.686 0.0946
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before -2.903 0.257 115 -11.277 <.0001
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop -0.4338 0.257 115 -1.686 0.3359
After Nonop - After OP 2.5405 0.257 115 9.871 <.0001
After Nonop - Before OP -0.3621 0.257 115 -1.407 0.4978
Before Nonop - After OP 2.9744 0.257 115 11.556 <.0001
Before Nonop - Before OP 0.0718 0.257 115 0.279 0.9924
After OP - Before OP -2.9026 0.257 115 -11.277 <.0001
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
digastric.area.in.coronal.plain
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 12535 12535 1 115.84 53.799 3.276e-11 ***
Time 17015 17015 1 115.84 73.024 5.943e-14 ***
OP.side:Time 11548 11548 1 115.84 49.562 1.457e-10 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: digastric.area.in.coronal.plain ~ OP.side * Time + (1 | ID)
Data: temp
REML criterion at convergence: 1355.5
Scaled residuals:
Min 1Q Median 3Q Max
-2.62339 -0.58435 -0.03063 0.59457 2.64715
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 452.7 21.28
Residual 233.0 15.26
Number of obs: 156, groups: ID, 38
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 78.815 4.231 64.652 18.626 < 2e-16 ***
OP.sideOP -35.136 3.457 115.838 -10.165 < 2e-16 ***
TimeBefore 3.679 3.457 115.838 1.064 0.289
OP.sideOP:TimeBefore 34.415 4.889 115.838 7.040 1.46e-10 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 78.8 4.23 63.5 70.4 87.3
OP 43.7 4.23 63.5 35.2 52.1
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 82.5 4.23 63.5 74.0 90.9
OP 81.8 4.23 63.5 73.3 90.2
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP 35.136 3.46 115 10.165 <.0001
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP 0.721 3.46 115 0.208 0.8353
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before -3.68 3.46 115 -1.064 0.2894
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before -38.09 3.46 115 -11.021 <.0001
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop -3.679 3.46 115 -1.064 0.7118
After Nonop - After OP 35.136 3.46 115 10.165 <.0001
After Nonop - Before OP -2.959 3.46 115 -0.856 0.8274
Before Nonop - After OP 38.815 3.46 115 11.229 <.0001
Before Nonop - Before OP 0.721 3.46 115 0.208 0.9968
After OP - Before OP -38.095 3.46 115 -11.021 <.0001
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
mylohyoid.coronal.plain.width.in.mm
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 4.5105 4.5105 1 101.41 33.933 6.768e-08 ***
Time 11.8053 11.8053 1 101.07 88.811 1.657e-15 ***
OP.side:Time 5.8908 5.8908 1 101.07 44.316 1.472e-09 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: mylohyoid.coronal.plain.width.in.mm ~ OP.side * Time + (1 | ID)
Data: temp
REML criterion at convergence: 214.1
Scaled residuals:
Min 1Q Median 3Q Max
-2.67264 -0.57979 -0.04548 0.56869 2.32702
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 0.4763 0.6902
Residual 0.1329 0.3646
Number of obs: 138, groups: ID, 34
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 3.29727 0.13411 46.87997 24.587 < 2e-16 ***
OP.sideOP -0.77717 0.08807 101.24049 -8.824 3.40e-14 ***
TimeBefore 0.17176 0.08843 101.07353 1.942 0.0549 .
OP.sideOP:TimeBefore 0.82652 0.12416 101.07353 6.657 1.47e-09 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 3.30 0.134 46.8 3.03 3.57
OP 2.52 0.133 46.0 2.25 2.79
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 3.47 0.134 46.8 3.20 3.74
OP 3.52 0.133 46.0 3.25 3.79
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP 0.7772 0.0881 101 8.824 <.0001
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP -0.0493 0.0881 101 -0.560 0.5765
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before -0.172 0.0884 101 -1.942 0.0549
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before -0.998 0.0872 101 -11.454 <.0001
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop -0.1718 0.0884 101 -1.942 0.2171
After Nonop - After OP 0.7772 0.0881 101 8.824 <.0001
After Nonop - Before OP -0.2211 0.0881 101 -2.510 0.0643
Before Nonop - After OP 0.9489 0.0881 101 10.774 <.0001
Before Nonop - Before OP -0.0493 0.0881 101 -0.560 0.9435
After OP - Before OP -0.9983 0.0872 101 -11.454 <.0001
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
infrahyoid.axial.plain.long.axis.in.mm
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 268.37 268.37 1 112.62 52.528 5.678e-11 ***
Time 351.76 351.76 1 112.62 68.850 2.587e-13 ***
OP.side:Time 170.47 170.47 1 112.62 33.366 6.883e-08 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: infrahyoid.axial.plain.long.axis.in.mm ~ OP.side * Time + (1 |
ID)
Data: temp
REML criterion at convergence: 759.3
Scaled residuals:
Min 1Q Median 3Q Max
-2.10193 -0.52977 -0.08592 0.54426 2.55214
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 11.428 3.381
Residual 5.109 2.260
Number of obs: 152, groups: ID, 37
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 22.1585 0.6661 59.4556 33.265 < 2e-16 ***
OP.sideOP -4.7755 0.5186 112.6159 -9.209 2.12e-15 ***
TimeBefore 0.9245 0.5186 112.6159 1.783 0.0773 .
OP.sideOP:TimeBefore 4.2361 0.7333 112.6159 5.776 6.88e-08 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 22.2 0.666 58.6 20.8 23.5
OP 17.4 0.666 58.6 16.0 18.7
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 23.1 0.666 58.6 21.7 24.4
OP 22.5 0.666 58.6 21.2 23.9
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP 4.776 0.519 112 9.209 <.0001
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP 0.539 0.519 112 1.040 0.3004
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before -0.924 0.519 112 -1.783 0.0773
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before -5.161 0.519 112 -9.952 <.0001
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop -0.924 0.519 112 -1.783 0.2870
After Nonop - After OP 4.776 0.519 112 9.209 <.0001
After Nonop - Before OP -0.385 0.519 112 -0.742 0.8797
Before Nonop - After OP 5.700 0.519 112 10.992 <.0001
Before Nonop - Before OP 0.539 0.519 112 1.040 0.7261
After OP - Before OP -5.161 0.519 112 -9.952 <.0001
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
infrahyoid.axial.plain.short.axis.in.mm
In this variable, you have an incorrect figure in the second measurement of patient number 1. The value is 11.3, much higher than the first measurement, 7.4, and in the opposite direction compared to the other patients.
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 9.5016 9.5016 1 114.84 19.961 1.862e-05 ***
Time 8.4467 8.4467 1 114.84 17.745 5.049e-05 ***
OP.side:Time 13.3878 13.3878 1 114.84 28.125 5.588e-07 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: infrahyoid.axial.plain.short.axis.in.mm ~ OP.side * Time + (1 |
ID)
Data: temp
REML criterion at convergence: 404.8
Scaled residuals:
Min 1Q Median 3Q Max
-2.6470 -0.4775 -0.1013 0.4093 3.8946
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 0.6925 0.8321
Residual 0.4760 0.6899
Number of obs: 156, groups: ID, 38
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 4.9889 0.1745 71.1495 28.583 < 2e-16 ***
OP.sideOP -1.0795 0.1562 114.8448 -6.909 2.88e-10 ***
TimeBefore -0.1205 0.1562 114.8448 -0.771 0.442
OP.sideOP:TimeBefore 1.1718 0.2210 114.8448 5.303 5.59e-07 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 4.99 0.175 71.4 4.64 5.34
OP 3.91 0.175 71.4 3.56 4.26
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 4.87 0.175 71.4 4.52 5.22
OP 4.96 0.175 71.4 4.61 5.31
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP 1.0795 0.156 115 6.909 <.0001
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP -0.0923 0.156 115 -0.591 0.5558
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before 0.121 0.156 115 0.771 0.4421
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before -1.051 0.156 115 -6.729 <.0001
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop 0.1205 0.156 115 0.771 0.8672
After Nonop - After OP 1.0795 0.156 115 6.909 <.0001
After Nonop - Before OP 0.0282 0.156 115 0.181 0.9979
Before Nonop - After OP 0.9590 0.156 115 6.138 <.0001
Before Nonop - Before OP -0.0923 0.156 115 -0.591 0.9347
After OP - Before OP -1.0513 0.156 115 -6.729 <.0001
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
infrahyoid.axial.plain.area
In this variable, you have an incorrect figure in the second measurement of patient number 1. The value is 240, much higher than the first measurement, 199, and in the opposite direction compared to the other patients.
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 6281.1 6281.1 1 112.04 41.891 2.637e-09 ***
Time 7544.9 7544.9 1 112.04 50.320 1.252e-10 ***
OP.side:Time 7918.6 7918.6 1 112.04 52.812 5.249e-11 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: infrahyoid.axial.plain.area ~ OP.side * Time + (1 | ID)
Data: temp
REML criterion at convergence: 1284.9
Scaled residuals:
Min 1Q Median 3Q Max
-2.3476 -0.4795 -0.1041 0.5317 4.2059
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 719.7 26.83
Residual 149.9 12.24
Number of obs: 152, groups: ID, 37
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 90.0861 4.8383 46.8986 18.620 < 2e-16 ***
OP.sideOP -27.2921 2.8092 112.0358 -9.715 < 2e-16 ***
TimeBefore -0.3447 2.8092 112.0358 -0.123 0.903
OP.sideOP:TimeBefore 28.8711 3.9728 112.0358 7.267 5.25e-11 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 90.1 4.84 46.9 80.4 99.8
OP 62.8 4.84 46.9 53.1 72.5
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 89.7 4.84 46.9 80.0 99.5
OP 91.3 4.84 46.9 81.6 101.1
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP 27.29 2.81 112 9.715 <.0001
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP -1.58 2.81 112 -0.562 0.5752
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before 0.345 2.81 112 0.123 0.9026
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before -28.526 2.81 112 -10.155 <.0001
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop 0.345 2.81 112 0.123 0.9993
After Nonop - After OP 27.292 2.81 112 9.715 <.0001
After Nonop - Before OP -1.234 2.81 112 -0.439 0.9715
Before Nonop - After OP 26.947 2.81 112 9.593 <.0001
Before Nonop - Before OP -1.579 2.81 112 -0.562 0.9431
After OP - Before OP -28.526 2.81 112 -10.155 <.0001
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
Rectus.capitis.posterior.major…long.axis
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 0.01525 0.01525 1 99.965 0.0148 0.9033
Time 0.04095 0.04095 1 99.965 0.0398 0.8422
OP.side:Time 1.44942 1.44942 1 99.965 1.4103 0.2378
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: Rectus.capitis.posterior.major...long.axis ~ OP.side * Time +
(1 | ID)
Data: temp
REML criterion at convergence: 495.3
Scaled residuals:
Min 1Q Median 3Q Max
-2.01762 -0.60176 -0.03084 0.59731 2.37767
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 6.029 2.455
Residual 1.028 1.014
Number of obs: 136, groups: ID, 33
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 13.2924 0.4616 39.8531 28.799 <2e-16 ***
OP.sideOP 0.2276 0.2459 99.9648 0.926 0.357
TimeBefore 0.1718 0.2459 99.9648 0.699 0.486
OP.sideOP:TimeBefore -0.4129 0.3477 99.9648 -1.188 0.238
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 13.3 0.462 39.9 12.4 14.2
OP 13.5 0.462 39.9 12.6 14.5
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 13.5 0.462 39.9 12.5 14.4
OP 13.3 0.462 39.9 12.3 14.2
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP -0.228 0.246 100 -0.926 0.3568
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP 0.185 0.246 100 0.754 0.4529
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before -0.172 0.246 100 -0.699 0.4864
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before 0.241 0.246 100 0.981 0.3290
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop -0.1718 0.246 100 -0.699 0.8974
After Nonop - After OP -0.2276 0.246 100 -0.926 0.7911
After Nonop - Before OP 0.0135 0.246 100 0.055 0.9999
Before Nonop - After OP -0.0559 0.246 100 -0.227 0.9958
Before Nonop - Before OP 0.1853 0.246 100 0.754 0.8749
After OP - Before OP 0.2412 0.246 100 0.981 0.7606
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
Rectus.capitis.posterior.major…short.axis
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 0.19125 0.19125 1 99.929 0.2883 0.5925
Time 1.50360 1.50360 1 99.929 2.2665 0.1354
OP.side:Time 0.16243 0.16243 1 99.929 0.2448 0.6218
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: Rectus.capitis.posterior.major...short.axis ~ OP.side * Time +
(1 | ID)
Data: temp
REML criterion at convergence: 418.3
Scaled residuals:
Min 1Q Median 3Q Max
-2.05263 -0.58282 -0.01263 0.49617 2.42003
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 2.0600 1.4353
Residual 0.6634 0.8145
Number of obs: 136, groups: ID, 33
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 8.180608 0.286362 46.614103 28.567 <2e-16 ***
OP.sideOP 0.005882 0.197543 99.928834 0.030 0.976
TimeBefore 0.141176 0.197543 99.928834 0.715 0.476
OP.sideOP:TimeBefore 0.138235 0.279368 99.928834 0.495 0.622
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 8.18 0.286 46.7 7.60 8.76
OP 8.19 0.286 46.7 7.61 8.76
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 8.32 0.286 46.7 7.75 8.90
OP 8.47 0.286 46.7 7.89 9.04
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP -0.00588 0.198 100 -0.030 0.9763
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP -0.14412 0.198 100 -0.730 0.4674
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before -0.141 0.198 100 -0.715 0.4765
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before -0.279 0.198 100 -1.414 0.1603
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop -0.14118 0.198 100 -0.715 0.8911
After Nonop - After OP -0.00588 0.198 100 -0.030 1.0000
After Nonop - Before OP -0.28529 0.198 100 -1.444 0.4751
Before Nonop - After OP 0.13529 0.198 100 0.685 0.9026
Before Nonop - Before OP -0.14412 0.198 100 -0.730 0.8850
After OP - Before OP -0.27941 0.198 100 -1.414 0.4934
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
Rectus.capitis.posterior.major…area
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
OP.side 12.421 12.421 1 99.98 0.2052 0.6515
Time 27.271 27.271 1 99.98 0.4506 0.5036
OP.side:Time 10.451 10.451 1 99.98 0.1727 0.6786
Linear mixed model fit by REML. t-tests use Satterthwaite's method [
lmerModLmerTest]
Formula: Rectus.capitis.posterior.major...area ~ OP.side * Time + (1 | ID)
Data: temp
REML criterion at convergence: 1055.1
Scaled residuals:
Min 1Q Median 3Q Max
-2.53223 -0.58265 -0.01463 0.55591 2.37070
Random effects:
Groups Name Variance Std.Dev.
ID (Intercept) 716.29 26.764
Residual 60.52 7.779
Number of obs: 136, groups: ID, 33
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 91.079 4.847 35.908 18.791 <2e-16 ***
OP.sideOP -0.050 1.887 99.980 -0.027 0.979
TimeBefore 1.450 1.887 99.980 0.769 0.444
OP.sideOP:TimeBefore -1.109 2.668 99.980 -0.416 0.679
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Time = After:
OP.side emmean SE df lower.CL upper.CL
Nonop 91.1 4.85 35.9 81.2 101
OP 91.0 4.85 35.9 81.2 101
Time = Before:
OP.side emmean SE df lower.CL upper.CL
Nonop 92.5 4.85 35.9 82.7 102
OP 91.4 4.85 35.9 81.5 101
Degrees-of-freedom method: kenward-roger
Confidence level used: 0.95
Time = After:
contrast estimate SE df t.ratio p.value
Nonop - OP 0.05 1.89 100 0.027 0.9789
Time = Before:
contrast estimate SE df t.ratio p.value
Nonop - OP 1.16 1.89 100 0.614 0.5405
Degrees-of-freedom method: kenward-roger
OP.side = Nonop:
contrast estimate SE df t.ratio p.value
After - Before -1.450 1.89 100 -0.769 0.4440
OP.side = OP:
contrast estimate SE df t.ratio p.value
After - Before -0.341 1.89 100 -0.181 0.8569
Degrees-of-freedom method: kenward-roger
contrast estimate SE df t.ratio p.value
After Nonop - Before Nonop -1.450 1.89 100 -0.769 0.8684
After Nonop - After OP 0.050 1.89 100 0.027 1.0000
After Nonop - Before OP -0.291 1.89 100 -0.154 0.9987
Before Nonop - After OP 1.500 1.89 100 0.795 0.8565
Before Nonop - Before OP 1.159 1.89 100 0.614 0.9273
After OP - Before OP -0.341 1.89 100 -0.181 0.9979
Degrees-of-freedom method: kenward-roger
P value adjustment: tukey method for comparing a family of 4 estimates
