#Import libraries and dataset:

Tided data:

#Create durations:
tided.CF <- CF %>%
  mutate(
    WordDuration  = WordEnd - WordStart,
    PhoneDuration = PhoneEnd - PhoneStart
  )

# Keep only rows where annotators made decisions
TP <- tided.CF %>% filter(!is.na(Prominence))
TB <- tided.CF %>% filter(!is.na(Boundary))

#Scale predictors
TP <- TP %>%
  mutate(across(c(Pitch_Prominence, Intensity_Prominence, WordDuration), scale))

TB <- TB %>%
  mutate(across(c(Pitch_Boundary, Intensity_Boundary, Silence_Boundary, PhoneDuration), scale))

Mixed Model for Prominence:

prom.model <- lmer(
  Prominence ~ Pitch_Prominence + Intensity_Prominence + WordDuration +
    (1 | Speaker) + (1 | Annotator),
  data = TP
)

summary(prom.model)
Linear mixed model fit by REML ['lmerMod']
Formula: Prominence ~ Pitch_Prominence + Intensity_Prominence + WordDuration +  
    (1 | Speaker) + (1 | Annotator)
   Data: TP

REML criterion at convergence: -490593.8

Scaled residuals: 
     Min       1Q   Median       3Q      Max 
-1.78378 -0.86425 -0.00702  0.86096  1.81577 

Random effects:
 Groups    Name        Variance  Std.Dev. 
 Speaker   (Intercept) 1.664e-08 1.290e-04
 Annotator (Intercept) 9.897e-10 3.146e-05
 Residual              8.344e-06 2.889e-03
Number of obs: 55409, groups:  Speaker, 32; Annotator, 5

Fixed effects:
                       Estimate Std. Error  t value
(Intercept)           5.642e-01  3.672e-05 15363.91
Pitch_Prominence      1.741e-01  2.444e-05  7123.95
Intensity_Prominence  1.555e-01  2.514e-05  6184.03
WordDuration         -2.055e-05  1.995e-05    -1.03

Correlation of Fixed Effects:
            (Intr) Ptch_P Intn_P
Ptch_Prmnnc -0.191              
Intnsty_Prm  0.070 -0.639       
WordDuratin  0.056 -0.175 -0.465

Mixed Model for Boundary:

bound.model <- lmer(
  Prominence ~ Pitch_Prominence + Intensity_Prominence + WordDuration +
    (1 | Speaker) + (1 | Annotator),
  data = TB
)

summary(bound.model)
Linear mixed model fit by REML ['lmerMod']
Formula: Prominence ~ Pitch_Prominence + Intensity_Prominence + WordDuration +  
    (1 | Speaker) + (1 | Annotator)
   Data: TB

REML criterion at convergence: -490601.8

Scaled residuals: 
     Min       1Q   Median       3Q      Max 
-1.78378 -0.86425 -0.00702  0.86096  1.81577 

Random effects:
 Groups    Name        Variance  Std.Dev. 
 Speaker   (Intercept) 1.664e-08 1.290e-04
 Annotator (Intercept) 9.897e-10 3.146e-05
 Residual              8.344e-06 2.889e-03
Number of obs: 55409, groups:  Speaker, 32; Annotator, 5

Fixed effects:
                       Estimate Std. Error  t value
(Intercept)          -1.323e-05  4.745e-05   -0.279
Pitch_Prominence      5.001e-01  7.020e-05 7123.949
Intensity_Prominence  4.999e-01  8.085e-05 6184.027
WordDuration         -1.231e-04  1.195e-04   -1.030

Correlation of Fixed Effects:
            (Intr) Ptch_P Intn_P
Ptch_Prmnnc -0.389              
Intnsty_Prm  0.116 -0.639       
WordDuratin -0.118 -0.175 -0.465
LS0tCnRpdGxlOiAiTWl4ZWQgbGluZWFyIHN0YXRzIC0gQXV0b1JQVCIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKI0ltcG9ydCBsaWJyYXJpZXMgYW5kIGRhdGFzZXQ6CgpgYGB7cn0KbGlicmFyeSh0aWR5dmVyc2UpCmxpYnJhcnkobG1lNCkKCkNGIDwtIHJlYWQuY3N2KCIvVXNlcnMvZmVybmFuZGFsb3VyZW5jby9Eb3dubG9hZHMvZmluYWxfY29tYmluZWRfZmlsZSAtIGZpbmFsX2NvbWJpbmVkX2ZpbGUuY3N2IiwKICAgICAgICAgICAgICAgaGVhZGVyID0gVFJVRSwKICAgICAgICAgICAgICAgcm93Lm5hbWVzID0gTlVMTCwKICAgICAgICAgICAgICAgY2hlY2submFtZXMgPSBGQUxTRSkKQ0YKI0NGOiBDb21iaW5lZCBGaWxlCmBgYAoKIyBUaWRlZCBkYXRhOgoKYGBge3J9CiNDcmVhdGUgZHVyYXRpb25zOgp0aWRlZC5DRiA8LSBDRiAlPiUKICBtdXRhdGUoCiAgICBXb3JkRHVyYXRpb24gID0gV29yZEVuZCAtIFdvcmRTdGFydCwKICAgIFBob25lRHVyYXRpb24gPSBQaG9uZUVuZCAtIFBob25lU3RhcnQKICApCgojIEtlZXAgb25seSByb3dzIHdoZXJlIGFubm90YXRvcnMgbWFkZSBkZWNpc2lvbnMKVFAgPC0gdGlkZWQuQ0YgJT4lIGZpbHRlcighaXMubmEoUHJvbWluZW5jZSkpClRCIDwtIHRpZGVkLkNGICU+JSBmaWx0ZXIoIWlzLm5hKEJvdW5kYXJ5KSkKCiNTY2FsZSBwcmVkaWN0b3JzClRQIDwtIFRQICU+JQogIG11dGF0ZShhY3Jvc3MoYyhQaXRjaF9Qcm9taW5lbmNlLCBJbnRlbnNpdHlfUHJvbWluZW5jZSwgV29yZER1cmF0aW9uKSwgc2NhbGUpKQoKVEIgPC0gVEIgJT4lCiAgbXV0YXRlKGFjcm9zcyhjKFBpdGNoX0JvdW5kYXJ5LCBJbnRlbnNpdHlfQm91bmRhcnksIFNpbGVuY2VfQm91bmRhcnksIFBob25lRHVyYXRpb24pLCBzY2FsZSkpCgojVFA6IFRpZGVkIGZvciBQcm9taW5lbmNlLyBUQjogVGlkZWQgZm9yIEJvdW5kYXJ5CmBgYAoKIyBNaXhlZCBNb2RlbCBmb3IgUHJvbWluZW5jZToKCmBgYHtyfQpwcm9tLm1vZGVsIDwtIGxtZXIoCiAgUHJvbWluZW5jZSB+IFBpdGNoX1Byb21pbmVuY2UgKyBJbnRlbnNpdHlfUHJvbWluZW5jZSArIFdvcmREdXJhdGlvbiArCiAgICAoMSB8IFNwZWFrZXIpICsgKDEgfCBBbm5vdGF0b3IpLAogIGRhdGEgPSBUUAopCgpzdW1tYXJ5KHByb20ubW9kZWwpCmBgYAoKIyBNaXhlZCBNb2RlbCBmb3IgQm91bmRhcnk6CgpgYGB7cn0KYm91bmQubW9kZWwgPC0gbG1lcigKICBQcm9taW5lbmNlIH4gUGl0Y2hfUHJvbWluZW5jZSArIEludGVuc2l0eV9Qcm9taW5lbmNlICsgV29yZER1cmF0aW9uICsKICAgICgxIHwgU3BlYWtlcikgKyAoMSB8IEFubm90YXRvciksCiAgZGF0YSA9IFRCCikKCnN1bW1hcnkoYm91bmQubW9kZWwpCmBgYAoK