---
title: "english"
output: html_document
date: "2023-05-23"
---

```{r}
table <- read.delim ("english.txt", stringsAsFactors=TRUE)
table
```

```{r}
contrasts <- cbind (c(+0.5, -0.5))
colnames (contrasts) <- c("-N+L")
contrasts (table$PROFICIENCY) <- contrasts
contrasts (table$PROFICIENCY)
```

```{r}
model <- lmerTest::lmer (F1..HZ. ~ PROFICIENCY + (PROFICIENCY | VOWEL) + (1 | SPEAKER), data=table)
summary (model)
confint (model)
```

```{r}
model2 <- lmerTest::lmer (F2..HZ. ~ PROFICIENCY + (PROFICIENCY | VOWEL) + (1 | SPEAKER), data=table)
summary (model2)
confint (model2)
```