1 Overview
This supplementary document complements the paper Distance vs. time. Articulatory and acoustic consequences of reduced vowel duration in Polish.
Only relevant R code is shown in this document.
The full code can be inspected on the OSF repository at https://osf.io/hy7nt/?view_only=d06775a2cb1d4786a0442b0f15d73296, in the file code/analysis.Rmd
.
2 Atriculatory data processing
The figure below illustrates the PC1-PC2 rotation for all individual speakers.
This is Figure 1 from the paper, which illustrates the articulatory data processing, as well as vowel formants. See paper for details.
3 Descriptive statistics plots
3.1 Articulatory undershoot
Tongue contours from ultrasound tongue imaging of speaker PL01. Tongue tip on the right.
3.2 Duration
Figure 4 from the paper with boxplots and means of vowel duration depending on speech rate and stress.
3.3 Vowel spaces
4 Inferential statistics
4.1 VIF (multicollinearity)
VIF values below 3 indicate absence of multicollinearity (Zuur, Ieno, and Elphick 2010). For all predictors, VIF < 3.
Variables | VIF |
---|---|
duration | 2.169934 |
f1.z | 1.379401 |
f2.z | 1.208097 |
f0.z | 1.122029 |
rate | 1.639475 |
stress | 1.451903 |
4.2 Duration
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
dur.lmer3 | 31 | -291.7616 | -128.3173 | 176.8808 | -353.7616 | NA | NA | NA |
dur.lmer2 | 33 | -341.1601 | -167.1710 | 203.5801 | -407.1601 | 53.39852 | 2 | 0 |
dur.lmer.full | 48 | -375.6100 | -122.5349 | 235.8050 | -471.6100 | 64.44983 | 15 | 0 |
4.3 F1
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
f1.lmer.base | 15 | 1760.472 | 1839.558 | -865.2361 | 1730.472 | NA | NA | NA |
f1.lmer.stress1 | 16 | 1697.351 | 1781.709 | -832.6755 | 1665.351 | 65.12115 | 1 | 0e+00 |
f1.lmer.stress2 | 21 | 1671.028 | 1781.748 | -814.5139 | 1629.028 | 36.32313 | 5 | 8e-07 |
## # R2 for Mixed Models
##
## Conditional R2: 0.830
## Marginal R2: 0.772
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.828
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.832
## V1 = a:
## contrast estimate SE df t.ratio p.value
## unstressed - stressed -0.1947 0.0803 70.9 -2.426 0.0178
##
## V1 = e:
## contrast estimate SE df t.ratio p.value
## unstressed - stressed -0.3938 0.0721 46.9 -5.462 <.0001
##
## V1 = i:
## contrast estimate SE df t.ratio p.value
## unstressed - stressed -0.0232 0.0708 43.8 -0.328 0.7444
##
## V1 = o:
## contrast estimate SE df t.ratio p.value
## unstressed - stressed -0.4242 0.0750 54.8 -5.655 <.0001
##
## V1 = u:
## contrast estimate SE df t.ratio p.value
## unstressed - stressed -0.0926 0.0759 57.4 -1.219 0.2277
##
## V1 = ɨ:
## contrast estimate SE df t.ratio p.value
## unstressed - stressed -0.1340 0.0739 51.7 -1.813 0.0756
##
## Degrees-of-freedom method: kenward-roger
V1 | duration.trend | SE | df | t.ratio | p.value |
---|---|---|---|---|---|
a | 0.0188220 | 0.0012463 | 1387.754 | 15.1027779 | 0.0000000 |
e | 0.0089742 | 0.0013664 | 1410.827 | 6.5679103 | 0.0000000 |
i | -0.0034263 | 0.0013199 | 1385.225 | -2.5959280 | 0.0095333 |
o | 0.0027981 | 0.0011339 | 1408.425 | 2.4675926 | 0.0137207 |
u | 0.0012137 | 0.0013415 | 1398.657 | 0.9047526 | 0.3657523 |
ɨ | 0.0009668 | 0.0010876 | 1406.486 | 0.8889355 | 0.3741898 |
4.4 F2
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
f2.lmer.base | 15 | 686.0166 | 765.1026 | -328.0083 | 656.0166 | NA | NA | NA |
f2.lmer.stress1 | 16 | 685.7771 | 770.1354 | -326.8885 | 653.7771 | 2.239533 | 1 | 0.1345216 |
f2.lmer.stress2 | 21 | 660.6134 | 771.3337 | -309.3067 | 618.6134 | 35.163703 | 5 | 0.0000014 |
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.906
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.907
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.909
V1 | duration.trend | SE | df | t.ratio | p.value |
---|---|---|---|---|---|
a | 0.0010623 | 0.0007896 | 71.72053 | 1.345283 | 0.1827721 |
e | 0.0036083 | 0.0009683 | 153.87450 | 3.726346 | 0.0002723 |
i | 0.0064001 | 0.0009630 | 149.92347 | 6.645943 | 0.0000000 |
o | -0.0046407 | 0.0007924 | 72.27819 | -5.856616 | 0.0000001 |
u | -0.0061957 | 0.0008945 | 114.65049 | -6.926301 | 0.0000000 |
ɨ | 0.0035699 | 0.0007826 | 67.64431 | 4.561831 | 0.0000219 |
4.5 Z1
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
z1.lmer.base | 15 | 1079.825 | 1158.911 | -524.9127 | 1049.825 | NA | NA | NA |
z1.lmer.stress1 | 16 | 1078.523 | 1162.882 | -523.2617 | 1046.523 | 3.301855 | 1 | 0.0692017 |
z1.lmer.stress2 | 21 | 1080.686 | 1191.407 | -519.3431 | 1038.686 | 7.837281 | 5 | 0.1654343 |
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.731
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.732
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.732
V1 | duration.trend | SE | df | t.ratio | p.value |
---|---|---|---|---|---|
a | 0.0088917 | 0.0011077 | 27.49535 | 8.0269993 | 0.0000000 |
e | 0.0041344 | 0.0012787 | 48.22526 | 3.2334337 | 0.0022094 |
i | -0.0009243 | 0.0012737 | 47.31917 | -0.7256695 | 0.4716154 |
o | 0.0094759 | 0.0011100 | 27.75403 | 8.5371785 | 0.0000000 |
u | 0.0000687 | 0.0012072 | 38.57507 | 0.0569243 | 0.9548993 |
ɨ | -0.0004451 | 0.0011012 | 26.79854 | -0.4042152 | 0.6892624 |
4.6 Z2
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
z2.lmer.base | 15 | 753.9717 | 833.0577 | -361.9858 | 723.9717 | NA | NA | NA |
z2.lmer.stress1 | 16 | 755.9562 | 840.3146 | -361.9781 | 723.9562 | 0.0154439 | 1 | 0.9010988 |
z2.lmer.stress2 | 21 | 760.4627 | 871.1831 | -359.2314 | 718.4627 | 5.4935092 | 5 | 0.3586581 |
## # R2 for Mixed Models
##
## Conditional R2: 0.896
## Marginal R2: 0.721
## # R2 for Mixed Models
##
## Conditional R2: 0.895
## Marginal R2: 0.721
## # R2 for Mixed Models
##
## Conditional R2: 0.896
## Marginal R2: 0.721
V1 | duration.trend | SE | df | t.ratio | p.value |
---|---|---|---|---|---|
a | -0.0031410 | 0.0009973 | 25.64689 | -3.1493657 | 0.0041252 |
e | -0.0006506 | 0.0011446 | 44.07179 | -0.5684129 | 0.5726397 |
i | 0.0040653 | 0.0011405 | 43.36942 | 3.5644529 | 0.0009023 |
o | -0.0066949 | 0.0009992 | 25.86902 | -6.7001979 | 0.0000004 |
u | -0.0035357 | 0.0010826 | 35.41626 | -3.2660975 | 0.0024243 |
ɨ | 0.0035870 | 0.0009917 | 25.02169 | 3.6172118 | 0.0013131 |
5 Model summaries
5.1 Duration
duration.log | ||||
---|---|---|---|---|
Predictors | Estimates | CI | p | df |
(Intercept) | 4.3187 | 4.1263 – 4.5111 | <0.001 | 14.8100 |
rate [normal] | -0.2209 | -0.3372 – -0.1046 | 0.00041 | 43.2044 |
rate [fast] | -0.4046 | -0.5810 – -0.2281 | 0.00016 | 16.4521 |
V1 [e] | -0.1048 | -0.3544 – 0.1449 | 0.37010 | 9.6185 |
V1 [i] | -0.1178 | -0.3944 – 0.1588 | 0.35959 | 8.8555 |
V1 [o] | -0.0733 | -0.1675 – 0.0208 | 0.12679 | 1367.8525 |
V1 [u] | -0.3184 | -0.5950 – -0.0418 | 0.02863 | 8.8555 |
V1 [ɨ] | 0.0229 | -0.2537 – 0.2995 | 0.85529 | 8.8555 |
stress [stressed] | 0.4710 | 0.3768 – 0.5652 | <0.001 | 238.9499 |
rate [normal] * V1 [e] | -0.0473 | -0.1707 – 0.0761 | 0.45222 | 1361.0048 |
rate [fast] * V1 [e] | -0.0470 | -0.1705 – 0.0764 | 0.45501 | 1361.0048 |
rate [normal] * V1 [i] | 0.0429 | -0.0805 – 0.1663 | 0.49544 | 1361.0048 |
rate [fast] * V1 [i] | -0.0223 | -0.1457 – 0.1011 | 0.72313 | 1361.0048 |
rate [normal] * V1 [o] | -0.0607 | -0.1842 – 0.0627 | 0.33475 | 1361.0048 |
rate [fast] * V1 [o] | -0.0667 | -0.1901 – 0.0567 | 0.28941 | 1361.0048 |
rate [normal] * V1 [u] | -0.0806 | -0.2040 – 0.0429 | 0.20060 | 1361.0048 |
rate [fast] * V1 [u] | -0.0502 | -0.1736 – 0.0733 | 0.42553 | 1361.0048 |
rate [normal] * V1 [ɨ] | -0.0581 | -0.1815 – 0.0653 | 0.35608 | 1361.0048 |
rate [fast] * V1 [ɨ] | -0.0838 | -0.2073 – 0.0396 | 0.18305 | 1361.0048 |
rate [normal] * stress [stressed] |
-0.0019 | -0.1253 – 0.1215 | 0.97590 | 1361.0048 |
rate [fast] * stress [stressed] |
-0.1241 | -0.2475 – -0.0007 | 0.04880 | 1361.0048 |
V1 [e] * stress [stressed] |
-0.1580 | -0.2815 – -0.0346 | 0.01213 | 1361.0048 |
V1 [i] * stress [stressed] |
-0.2180 | -0.3415 – -0.0946 | 0.00055 | 1361.0048 |
V1 [o] * stress [stressed] |
-0.0017 | -0.1252 – 0.1217 | 0.97795 | 1361.0048 |
V1 [u] * stress [stressed] |
0.0419 | -0.0816 – 0.1653 | 0.50578 | 1361.0048 |
V1 [ɨ] * stress [stressed] |
-0.0527 | -0.1761 – 0.0707 | 0.40247 | 1361.0048 |
(rate [normal] * V1 [e]) * stress [stressed] |
-0.0292 | -0.2037 – 0.1454 | 0.74299 | 1361.0048 |
(rate [fast] * V1 [e]) * stress [stressed] |
-0.0129 | -0.1875 – 0.1617 | 0.88478 | 1361.0048 |
(rate [normal] * V1 [i]) * stress [stressed] |
-0.0241 | -0.1987 – 0.1505 | 0.78666 | 1361.0048 |
(rate [fast] * V1 [i]) * stress [stressed] |
0.0225 | -0.1521 – 0.1970 | 0.80075 | 1361.0048 |
(rate [normal] * V1 [o]) * stress [stressed] |
-0.0624 | -0.2369 – 0.1122 | 0.48338 | 1361.0048 |
(rate [fast] * V1 [o]) * stress [stressed] |
-0.0897 | -0.2643 – 0.0848 | 0.31347 | 1361.0048 |
(rate [normal] * V1 [u]) * stress [stressed] |
-0.0273 | -0.2018 – 0.1473 | 0.75924 | 1361.0048 |
(rate [fast] * V1 [u]) * stress [stressed] |
-0.1707 | -0.3452 – 0.0039 | 0.05532 | 1361.0048 |
(rate [normal] * V1 [ɨ]) * stress [stressed] |
-0.0381 | -0.2127 – 0.1365 | 0.66853 | 1361.0048 |
(rate [fast] * V1 [ɨ]) * stress [stressed] |
-0.0843 | -0.2589 – 0.0903 | 0.34366 | 1361.0048 |
Random Effects | ||||
σ2 | 0.0396 | |||
τ00 frame | 0.0154 | |||
τ00 speaker | 0.0193 | |||
τ11 speaker.ratenormal | 0.0135 | |||
τ11 speaker.ratefast | 0.0498 | |||
τ11 speaker.stressstressed | 0.0031 | |||
ρ01 speaker.ratenormal | -0.6410 | |||
ρ01 speaker.ratefast | -0.5484 | |||
ρ01 speaker.stressstressed | -0.9139 | |||
N speaker | 10 | |||
N frame | 12 | |||
Observations | 1440 | |||
Marginal R2 / Conditional R2 | 0.706 / NA |
5.2 F1
f1.z | ||||
---|---|---|---|---|
Predictors | Estimates | CI | p | df |
(Intercept) | 0.1655 | -0.2195 – 0.5504 | 0.36633 | 11.6159 |
V1 [e] | -0.8117 | -1.3554 – -0.2681 | 0.00709 | 11.4664 |
V1 [i] | -0.9752 | -1.5739 – -0.3766 | 0.00452 | 10.2636 |
V1 [o] | -0.3559 | -0.5738 – -0.1381 | 0.00138 | 1405.5179 |
V1 [u] | -0.9014 | -1.4958 – -0.3071 | 0.00717 | 9.6528 |
V1 [ɨ] | -0.7690 | -1.3653 – -0.1726 | 0.01660 | 9.9580 |
duration | 0.0188 | 0.0164 – 0.0213 | <0.001 | 1387.7540 |
stress [stressed] | 0.1947 | 0.0346 – 0.3547 | 0.01783 | 70.9286 |
V1 [e] * duration | -0.0098 | -0.0135 – -0.0062 | <0.001 | 1412.7627 |
V1 [i] * duration | -0.0222 | -0.0258 – -0.0187 | <0.001 | 1413.7521 |
V1 [o] * duration | -0.0160 | -0.0193 – -0.0127 | <0.001 | 1410.9211 |
V1 [u] * duration | -0.0176 | -0.0212 – -0.0140 | <0.001 | 1411.8737 |
V1 [ɨ] * duration | -0.0179 | -0.0211 – -0.0146 | <0.001 | 1409.6337 |
V1 [e] * stress [stressed] |
0.1991 | 0.0228 – 0.3753 | 0.02687 | 1404.6665 |
V1 [i] * stress [stressed] |
-0.1714 | -0.3457 – 0.0029 | 0.05387 | 1404.0928 |
V1 [o] * stress [stressed] |
0.2295 | 0.0487 – 0.4104 | 0.01288 | 1403.3462 |
V1 [u] * stress [stressed] |
-0.1021 | -0.2845 – 0.0803 | 0.27231 | 1404.7296 |
V1 [ɨ] * stress [stressed] |
-0.0607 | -0.2398 – 0.1183 | 0.50609 | 1402.6657 |
Random Effects | ||||
σ2 | 0.1740 | |||
τ00 frame | 0.0697 | |||
τ00 speaker | 0.0029 | |||
τ11 speaker.stressstressed | 0.0178 | |||
ρ01 speaker | -1.0000 | |||
N speaker | 10 | |||
N frame | 12 | |||
Observations | 1440 | |||
Marginal R2 / Conditional R2 | 0.836 / NA |
5.3 F2
f2.z | ||||
---|---|---|---|---|
Predictors | Estimates | CI | p | df |
(Intercept) | -0.3897 | -0.7957 – 0.0163 | 0.05794 | 8.8492 |
V1 [e] | 0.6104 | 0.0373 – 1.1836 | 0.03937 | 8.7240 |
V1 [i] | 1.6015 | 0.9638 – 2.2391 | 0.00038 | 8.2679 |
V1 [o] | -0.2358 | -0.3885 – -0.0831 | 0.00249 | 1411.1814 |
V1 [u] | -0.3495 | -0.9855 – 0.2866 | 0.24094 | 8.0210 |
V1 [ɨ] | 0.4526 | -0.1841 – 1.0894 | 0.14020 | 8.1408 |
duration | 0.0011 | -0.0005 – 0.0026 | 0.18277 | 71.7205 |
V1 [e] * duration | 0.0025 | 0.0003 – 0.0048 | 0.02587 | 1419.9932 |
V1 [i] * duration | 0.0053 | 0.0031 – 0.0076 | <0.001 | 1419.3853 |
V1 [o] * duration | -0.0057 | -0.0076 – -0.0038 | <0.001 | 1412.3655 |
V1 [u] * duration | -0.0073 | -0.0094 – -0.0051 | <0.001 | 1417.6616 |
V1 [ɨ] * duration | 0.0025 | 0.0006 – 0.0044 | 0.01082 | 1412.6052 |
Random Effects | ||||
σ2 | 0.0889 | |||
τ00 frame | 0.0842 | |||
τ00 speaker | 0.0043 | |||
τ11 speaker.duration | 0.0000 | |||
ρ01 speaker | -1.0000 | |||
N speaker | 10 | |||
N frame | 12 | |||
Observations | 1440 | |||
Marginal R2 / Conditional R2 | 0.907 / NA |
5.4 Z1
z1 | ||||
---|---|---|---|---|
Predictors | Estimates | CI | p | df |
(Intercept) | -0.0013 | -0.2825 – 0.2799 | 0.99230 | 16.5061 |
V1 [e] | -0.4182 | -0.7976 – -0.0389 | 0.03314 | 13.3278 |
V1 [i] | -0.8547 | -1.2696 – -0.4398 | 0.00078 | 11.5918 |
V1 [o] | -0.0092 | -0.1832 – 0.1648 | 0.91747 | 1411.4250 |
V1 [u] | -0.0972 | -0.5074 – 0.3130 | 0.61118 | 10.6587 |
V1 [ɨ] | -0.0935 | -0.5058 – 0.3188 | 0.62769 | 11.1077 |
duration | 0.0089 | 0.0066 – 0.0112 | <0.001 | 27.4954 |
V1 [e] * duration | -0.0048 | -0.0073 – -0.0022 | 0.00026 | 1418.6184 |
V1 [i] * duration | -0.0098 | -0.0124 – -0.0073 | <0.001 | 1419.1450 |
V1 [o] * duration | 0.0006 | -0.0016 – 0.0028 | 0.60526 | 1413.4468 |
V1 [u] * duration | -0.0088 | -0.0112 – -0.0064 | <0.001 | 1415.6416 |
V1 [ɨ] * duration | -0.0093 | -0.0115 – -0.0071 | <0.001 | 1412.0221 |
Random Effects | ||||
σ2 | 0.1154 | |||
τ00 frame | 0.0320 | |||
τ00 speaker | 0.0239 | |||
τ11 speaker.duration | 0.0000 | |||
ρ01 speaker | -0.9948 | |||
ICC | 0.2402 | |||
N speaker | 10 | |||
N frame | 12 | |||
Observations | 1440 | |||
Marginal R2 / Conditional R2 | 0.680 / 0.757 |
5.5 Z2
z2 | ||||
---|---|---|---|---|
Predictors | Estimates | CI | p | df |
(Intercept) | -0.1205 | -0.6668 – 0.4257 | 0.62680 | 8.4109 |
V1 [e] | 0.5192 | -0.2493 – 1.2877 | 0.15759 | 7.9250 |
V1 [i] | 1.2666 | 0.4094 – 2.1238 | 0.00950 | 7.6823 |
V1 [o] | -0.2783 | -0.4311 – -0.1255 | 0.00036 | 1408.2080 |
V1 [u] | -0.5101 | -1.3663 – 0.3462 | 0.20469 | 7.5513 |
V1 [ɨ] | 0.1881 | -0.6686 – 1.0448 | 0.62395 | 7.6137 |
duration | -0.0031 | -0.0052 – -0.0011 | 0.00413 | 25.6469 |
V1 [e] * duration | 0.0025 | 0.0002 – 0.0047 | 0.02945 | 1412.6676 |
V1 [i] * duration | 0.0072 | 0.0050 – 0.0094 | <0.001 | 1412.9083 |
V1 [o] * duration | -0.0036 | -0.0055 – -0.0016 | 0.00035 | 1408.6364 |
V1 [u] * duration | -0.0004 | -0.0025 – 0.0017 | 0.71441 | 1412.3145 |
V1 [ɨ] * duration | 0.0067 | 0.0048 – 0.0087 | <0.001 | 1408.1702 |
Random Effects | ||||
σ2 | 0.0889 | |||
τ00 frame | 0.1548 | |||
τ00 speaker | 0.0190 | |||
τ11 speaker.duration | 0.0000 | |||
ρ01 speaker | -0.9900 | |||
ICC | 0.6412 | |||
N speaker | 10 | |||
N frame | 12 | |||
Observations | 1440 | |||
Marginal R2 / Conditional R2 | 0.722 / 0.900 |
6 Further exploratory analyses
6.1 F0
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
f0.lmer1 | 11 | 3832.230 | 3890.196 | -1905.115 | 3810.230 | NA | NA | NA |
f0.lmer2 | 16 | 3816.583 | 3900.897 | -1892.292 | 3784.583 | 25.6465446 | 5 | 0.0001045 |
f0.lmer3 | 21 | 3825.897 | 3936.559 | -1891.948 | 3783.897 | 0.6864549 | 5 | 0.9837038 |
f0.z | |||
---|---|---|---|
Predictors | Estimates | CI | p |
(Intercept) | -0.65 | -0.82 – -0.49 | <0.001 |
V11 | 0.11 | -0.23 – 0.46 | 0.520 |
V12 | 0.09 | -0.28 – 0.46 | 0.629 |
V13 | -0.26 | -0.64 – 0.11 | 0.171 |
V14 | 0.04 | -0.28 – 0.36 | 0.820 |
V15 | 0.39 | 0.05 – 0.74 | 0.026 |
duration | 0.01 | 0.01 – 0.01 | <0.001 |
stress [stressed] | 0.29 | 0.19 – 0.40 | <0.001 |
V11 * duration | -0.01 | -0.01 – -0.00 | 0.001 |
V12 * duration | -0.00 | -0.01 – 0.00 | 0.055 |
V13 * duration | 0.00 | -0.00 – 0.01 | 0.053 |
V14 * duration | -0.00 | -0.01 – 0.00 | 0.187 |
V15 * duration | 0.01 | 0.00 – 0.01 | 0.018 |
Random Effects | |||
σ2 | 0.81 | ||
τ00 frame | 0.03 | ||
τ00 speaker | 0.00 | ||
N speaker | 10 | ||
N frame | 12 | ||
Observations | 1436 | ||
Marginal R2 / Conditional R2 | 0.186 / NA |
6.2 Modelling after excluding slow speech rate data
6.2.1 F1
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
f1.lmer.base | 15 | 1239.295 | 1312.299 | -604.6477 | 1209.295 | NA | NA | NA |
f1.lmer.stress1 | 16 | 1166.619 | 1244.490 | -567.3097 | 1134.619 | 74.67616 | 1 | 0 |
f1.lmer.stress2 | 21 | 1126.087 | 1228.293 | -542.0436 | 1084.087 | 50.53219 | 5 | 0 |
## # R2 for Mixed Models
##
## Conditional R2: 0.781
## Marginal R2: 0.718
## # R2 for Mixed Models
##
## Conditional R2: 0.798
## Marginal R2: 0.732
## # R2 for Mixed Models
##
## Conditional R2: 0.808
## Marginal R2: 0.743
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
f2.lmer.base | 15 | 473.2132 | 546.2172 | -221.6066 | 443.2132 | NA | NA | NA |
f2.lmer.stress1 | 16 | 472.4649 | 550.3358 | -220.2324 | 440.4649 | 2.748322 | 1 | 0.0973565 |
f2.lmer.stress2 | 21 | 430.0225 | 532.2281 | -194.0113 | 388.0225 | 52.442374 | 5 | 0.0000000 |
## # R2 for Mixed Models
##
## Conditional R2: 0.904
## Marginal R2: 0.818
## # R2 for Mixed Models
##
## Conditional R2: 0.904
## Marginal R2: 0.819
## # R2 for Mixed Models
##
## Conditional R2: 0.909
## Marginal R2: 0.826
V1 | duration.trend | SE | df | t.ratio | p.value |
---|---|---|---|---|---|
a | 0.0020044 | 0.0013460 | 38.26200 | 1.489133 | 0.1446494 |
e | 0.0025912 | 0.0016847 | 90.40466 | 1.538116 | 0.1275131 |
i | 0.0076520 | 0.0015817 | 69.63964 | 4.837880 | 0.0000076 |
o | -0.0041391 | 0.0014093 | 46.03407 | -2.936982 | 0.0051610 |
u | -0.0063933 | 0.0016050 | 74.53133 | -3.983481 | 0.0001566 |
ɨ | 0.0030065 | 0.0014465 | 49.48094 | 2.078507 | 0.0428693 |
6.2.2 Z1
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
z1.lmer.base | 15 | 496.0472 | 569.0512 | -233.0236 | 466.0472 | NA | NA | NA |
z1.lmer.stress1 | 16 | 497.7308 | 575.6017 | -232.8654 | 465.7308 | 0.3164545 | 1 | 0.5737462 |
z1.lmer.stress2 | 21 | 499.9211 | 602.1267 | -228.9605 | 457.9211 | 7.8096519 | 5 | 0.1670428 |
## # R2 for Mixed Models
##
## Conditional R2: 0.749
## Marginal R2: 0.642
## # R2 for Mixed Models
##
## Conditional R2: 0.749
## Marginal R2: 0.642
## # R2 for Mixed Models
##
## Conditional R2: 0.751
## Marginal R2: 0.642
V1 | duration.trend | SE | df | t.ratio | p.value |
---|---|---|---|---|---|
a | 0.0062979 | 0.0013590 | 38.41946 | 4.6341868 | 0.0000405 |
e | 0.0024480 | 0.0017051 | 92.59074 | 1.4356780 | 0.1544638 |
i | 0.0008437 | 0.0016004 | 70.19398 | 0.5271550 | 0.5997481 |
o | 0.0057680 | 0.0014241 | 46.54966 | 4.0502032 | 0.0001922 |
u | -0.0008461 | 0.0016260 | 76.80410 | -0.5203748 | 0.6042992 |
ɨ | -0.0005280 | 0.0014609 | 49.66675 | -0.3614091 | 0.7193266 |
6.3 Z2
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
z2.lmer.base | 15 | 202.6340 | 275.6380 | -86.31698 | 172.6340 | NA | NA | NA |
z2.lmer.stress1 | 16 | 204.3809 | 282.2518 | -86.19046 | 172.3809 | 0.2530395 | 1 | 0.6149429 |
z2.lmer.stress2 | 21 | 198.0688 | 300.2744 | -78.03441 | 156.0688 | 16.3120910 | 5 | 0.0060073 |
## # R2 for Mixed Models
##
## Conditional R2: 0.914
## Marginal R2: 0.720
## # R2 for Mixed Models
##
## Conditional R2: 0.914
## Marginal R2: 0.720
## # R2 for Mixed Models
##
## Conditional R2: 0.915
## Marginal R2: 0.720
`
V1 | duration.trend | SE | df | t.ratio | p.value |
---|---|---|---|---|---|
a | -0.0028013 | 0.0011663 | 35.80871 | -2.4019719 | 0.0216210 |
e | -0.0006214 | 0.0014538 | 84.24792 | -0.4274275 | 0.6701585 |
i | 0.0024638 | 0.0013669 | 64.69699 | 1.8025080 | 0.0761235 |
o | -0.0058029 | 0.0012198 | 42.96368 | -4.7572397 | 0.0000224 |
u | -0.0057657 | 0.0013876 | 69.89782 | -4.1552199 | 0.0000906 |
ɨ | 0.0010473 | 0.0012515 | 46.08212 | 0.8368970 | 0.4069697 |
6.4 Include more PCs
This function takes a wide spline dataset as input, where spline data are in Cartesian coordinates named X_1:Y_42. This function drops fanlines with missing data, and performs a PCA on the rest of the data
This function applies linear regression performing formants, based on PC1 and PC2. Predicted formant values are added to the data (comp1 and comp2).
Run PCA for each speaker in dataset
Join with acoustic data and rotate by-speaker.
mean.comb |
---|
0.979386 |
Plot space.
6.4.1 Z1
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
z1.lmer.base | 15 | 1450.505 | 1529.591 | -710.2527 | 1420.505 | NA | NA | NA |
z1.lmer.stress1 | 16 | 1449.534 | 1533.892 | -708.7670 | 1417.534 | 2.971243 | 1 | 0.0847567 |
z1.lmer.stress2 | 21 | 1452.380 | 1563.100 | -705.1899 | 1410.380 | 7.154220 | 5 | 0.2094211 |
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.719
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.719
## Random effect variances not available. Returned R2 does not account for random effects.
## # R2 for Mixed Models
##
## Conditional R2: NA
## Marginal R2: 0.720
No significant effect of stress emerges.
## Linear mixed model fit by REML ['lmerModLmerTest']
## Formula: z1 ~ (1 + duration | speaker) + (1 | frame) + V1 * duration
## Data: midpoint_rotated_4PCs
## REML criterion at convergence: 1507.675
## Random effects:
## Groups Name Std.Dev. Corr
## frame (Intercept) 0.1850794
## speaker (Intercept) 0.0725332
## duration 0.0009817 -1.00
## Residual 0.3918978
## Number of obs: 1440, groups: frame, 12; speaker, 10
## Fixed Effects:
## (Intercept) V1e V1i V1o V1u
## -0.0071179 -0.3859335 -0.8904109 -0.1072723 -0.2021698
## V1ɨ duration V1e:duration V1i:duration V1o:duration
## -0.1597102 0.0113978 -0.0066102 -0.0117384 -0.0006565
## V1u:duration V1ɨ:duration
## -0.0127163 -0.0115170
6.4.2 Z2
npar | AIC | BIC | logLik | deviance | Chisq | Df | Pr(>Chisq) | |
---|---|---|---|---|---|---|---|---|
z2.lmer.base | 15 | 667.4576 | 746.5436 | -318.7288 | 637.4576 | NA | NA | NA |
z2.lmer.stress1 | 16 | 669.3933 | 753.7517 | -318.6967 | 637.3933 | 0.064281 | 1 | 0.7998533 |
z2.lmer.stress2 | 21 | 671.7338 | 782.4542 | -314.8669 | 629.7338 | 7.659469 | 5 | 0.1760285 |
## # R2 for Mixed Models
##
## Conditional R2: 0.905
## Marginal R2: 0.753
## # R2 for Mixed Models
##
## Conditional R2: 0.905
## Marginal R2: 0.753
## # R2 for Mixed Models
##
## Conditional R2: 0.905
## Marginal R2: 0.752
## Linear mixed model fit by REML ['lmerModLmerTest']
## Formula: z2 ~ (1 + duration | speaker) + (1 | frame) + V1 * duration
## Data: midpoint_rotated_4PCs
## REML criterion at convergence: 629.3296
## Random effects:
## Groups Name Std.Dev. Corr
## frame (Intercept) 0.374359
## speaker (Intercept) 0.174768
## duration 0.002869 -1.00
## Residual 0.284728
## Number of obs: 1440, groups: frame, 12; speaker, 10
## Fixed Effects:
## (Intercept) V1e V1i V1o V1u
## -0.1968455 0.6291168 1.3587379 -0.3125603 -0.5222861
## V1ɨ duration V1e:duration V1i:duration V1o:duration
## 0.2731720 -0.0015114 0.0008838 0.0060753 -0.0040087
## V1u:duration V1ɨ:duration
## -0.0018357 0.0047848
7 Session info
## ─ Session info ───────────────────────────────────────────────────────────────
## setting value
## version R version 4.1.0 (2021-05-18)
## os macOS Big Sur 10.16
## system x86_64, darwin17.0
## ui X11
## language (EN)
## collate en_US.UTF-8
## ctype en_GB.UTF-8
## tz Europe/Berlin
## date 2021-06-11
##
## ─ Packages ───────────────────────────────────────────────────────────────────
## ! package * version date lib
## P abind 1.4-5 2016-07-21 [?]
## P assertthat 0.2.1 2019-03-21 [?]
## P backports 1.2.1 2020-12-09 [?]
## P bayestestR 0.8.2 2021-01-26 [?]
## P bookdown 0.21 2020-10-13 [?]
## P boot 1.3-27 2021-02-12 [?]
## P broom 0.7.5 2021-02-19 [?]
## P cachem 1.0.4 2021-02-13 [?]
## P callr 3.5.1 2020-10-13 [?]
## P car 3.0-10 2020-09-29 [?]
## P carData * 3.0-4 2020-05-22 [?]
## P cellranger 1.1.0 2016-07-27 [?]
## P cli 2.3.1 2021-02-23 [?]
## P codetools 0.2-18 2020-11-04 [4]
## P colorspace 2.0-0 2020-11-11 [?]
## P crayon 1.4.1 2021-02-08 [?]
## P curl 4.3 2019-12-02 [?]
## P data.table 1.14.0 2021-02-21 [?]
## P DBI 1.1.1 2021-01-15 [?]
## P dbplyr 2.1.0 2021-02-03 [?]
## P desc 1.3.0 2021-03-05 [?]
## P devtools 2.3.2 2020-09-18 [?]
## P digest 0.6.27 2020-10-24 [?]
## P dplyr * 1.0.5 2021-03-05 [?]
## P effects * 4.2-0 2020-08-11 [?]
## P effectsize 0.4.4 2021-03-14 [?]
## P ellipsis 0.3.1 2020-05-15 [?]
## P emmeans * 1.5.5-1 2021-03-21 [?]
## P estimability 1.3 2018-02-11 [?]
## P evaluate 0.14 2019-05-28 [?]
## P fansi 0.4.2 2021-01-15 [?]
## P farver 2.1.0 2021-02-28 [?]
## P fastmap 1.1.0 2021-01-25 [?]
## P forcats * 0.5.1 2021-01-27 [?]
## P foreign 0.8-81 2020-12-22 [4]
## P fs 1.5.0 2020-07-31 [?]
## P generics 0.1.0 2020-10-31 [?]
## P ggeffects 1.0.2 2021-03-17 [?]
## P ggplot2 * 3.3.3 2020-12-30 [?]
## P ggpubr * 0.4.0 2020-06-27 [?]
## P ggrepel * 0.9.1 2021-01-15 [?]
## P ggsignif 0.6.1 2021-02-23 [?]
## P glue 1.4.2 2020-08-27 [?]
## P gtable 0.3.0 2019-03-25 [?]
## P haven 2.3.1 2020-06-01 [?]
## P here 1.0.1 2020-12-13 [?]
## P highr 0.8 2019-03-20 [?]
## P hms 1.0.0 2021-01-13 [?]
## P htmltools 0.5.1.1 2021-01-22 [?]
## P httr 1.4.2 2020-07-20 [?]
## P insight 0.13.1 2021-02-22 [?]
## P jsonlite 1.7.2 2020-12-09 [?]
## P knitr 1.31 2021-01-27 [?]
## P labeling 0.4.2 2020-10-20 [?]
## P lattice 0.20-41 2020-04-02 [?]
## P lifecycle 1.0.0 2021-02-15 [?]
## P lme4 * 1.1-26 2020-12-01 [?]
## P lmerTest * 3.1-3 2020-10-23 [?]
## P lubridate 1.7.10 2021-02-26 [?]
## P magrittr 2.0.1 2020-11-17 [?]
## P MASS 7.3-53.1 2021-02-12 [?]
## P Matrix * 1.3-2 2021-01-06 [?]
## P memoise 2.0.0 2021-01-26 [?]
## P minqa 1.2.4 2014-10-09 [?]
## P mitools 2.4 2019-04-26 [?]
## P modelr 0.1.8 2020-05-19 [?]
## P munsell 0.5.0 2018-06-12 [?]
## P mvtnorm 1.1-1 2020-06-09 [?]
## P nlme 3.1-152 2021-02-04 [4]
## P nloptr 1.2.2.2 2020-07-02 [?]
## P nnet 7.3-15 2021-01-24 [?]
## P numDeriv 2016.8-1.1 2019-06-06 [?]
## P openxlsx 4.2.3 2020-10-27 [?]
## P optimx * 2020-4.2 2020-04-08 [?]
## P parameters 0.12.0 2021-02-21 [?]
## P patchwork * 1.1.1 2020-12-17 [?]
## P pbkrtest 0.5.1 2021-03-09 [?]
## P performance * 0.7.0 2021-02-03 [?]
## P pillar 1.5.1 2021-03-05 [?]
## P pkgbuild 1.2.0 2020-12-15 [?]
## P pkgconfig 2.0.3 2019-09-22 [?]
## P pkgload 1.2.0 2021-02-23 [?]
## P plyr 1.8.6 2020-03-03 [?]
## P prettyunits 1.1.1 2020-01-24 [?]
## P printr * 0.1.1 2021-01-27 [?]
## P processx 3.5.0 2021-03-23 [?]
## P ps 1.6.0 2021-02-28 [?]
## P purrr * 0.3.4 2020-04-17 [?]
## P R6 2.5.0 2020-10-28 [?]
## P raster * 3.4-5 2020-11-14 [?]
## P Rcpp 1.0.6 2021-01-15 [?]
## P readr * 1.4.0 2020-10-05 [?]
## P readxl 1.3.1 2019-03-13 [?]
## P remotes 2.2.0 2020-07-21 [?]
## P reprex 1.0.0 2021-01-27 [?]
## P rio 0.5.26 2021-03-01 [?]
## P rlang 0.4.10 2020-12-30 [?]
## P rmarkdown 2.7 2021-02-19 [?]
## P rmdformats 1.0.1 2021-01-13 [?]
## P rprojroot 2.0.2 2020-11-15 [?]
## P rstatix 0.7.0 2021-02-13 [?]
## P rstudioapi 0.13 2020-11-12 [?]
## P rticulate * 1.6.0 2021-06-04 [?]
## P rvest 1.0.0 2021-03-09 [?]
## P scales 1.1.1 2020-05-11 [?]
## P sessioninfo 1.1.1 2018-11-05 [?]
## P sjlabelled 1.1.7 2020-09-24 [?]
## P sjmisc 2.8.6 2021-01-07 [?]
## P sjPlot * 2.8.7 2021-01-10 [?]
## P sjstats 0.18.1 2021-01-09 [?]
## P sp * 1.4-5 2021-01-10 [?]
## P statmod 1.4.35 2020-10-19 [?]
## P stringi 1.5.3 2020-09-09 [?]
## P stringr * 1.4.0 2019-02-10 [?]
## P survey 4.0 2020-04-03 [?]
## P survival 3.2-10 2021-03-16 [?]
## P testthat 3.0.2 2021-02-14 [?]
## P tibble * 3.1.0 2021-02-25 [?]
## P tidyr * 1.1.3 2021-03-03 [?]
## P tidyselect 1.1.0 2020-05-11 [?]
## P tidyverse * 1.3.0 2019-11-21 [?]
## P usdm * 1.1-18 2017-06-25 [?]
## P usethis 2.0.1 2021-02-10 [?]
## P utf8 1.2.1 2021-03-12 [?]
## P vctrs 0.3.6 2020-12-17 [?]
## P withr 2.4.1 2021-01-26 [?]
## P xfun 0.22 2021-03-11 [?]
## P xml2 1.3.2 2020-04-23 [?]
## P xtable 1.8-4 2019-04-21 [?]
## P yaml 2.2.1 2020-02-01 [?]
## P zip 2.1.1 2020-08-27 [?]
## source
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## Github (stefanocoretta/rticulate@7e8ec81)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
## CRAN (R 4.1.0)
##
## [1] /Users/ste/repos/research/polish-reduction/renv/library/R-4.1/x86_64-apple-darwin17.0
## [2] /private/var/folders/yj/b5zvw9lj2zg6tp920wyg14y80000gn/T/RtmpmPGmjh/renv-system-library
## [3] /Library/Frameworks/R.framework/Versions/4.1/Resources/site-library
## [4] /Library/Frameworks/R.framework/Versions/4.1/Resources/library
##
## P ── Loaded and on-disk path mismatch.