R |mixed effect model

---
title: "mixed effect model"
author: "gavin"
date: "2019年8月16日"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

```{r}
library(lme4)
politeness=read.csv("/Users/Gavin/Downloads/politeness_data.csv")
# head(politeness), tail(politeness),summary(politeness), str(politeness), colnames(politeness)
which(is.na(politeness)==T)
boxplot(frequency ~ attitude*gender,
        col=c("white","lightgray"),politeness)
```

```{r}
NEWdata <- data.frame(t(politeness))
NEWdata <- NEWdata[,-39]
NEW_MODEL = lmer(as.numeric(NEWdata[5,]) ~ as.numeric(NEWdata[6,]) + as.numeric((1|NEWdata[7,])))
NEW_MODEL
```


```{r}
politeness.model = lmer(frequency ~ attitude +
                          (1|subject) + (1|scenario), data=politeness)
politeness.model

pol